After hard coding in some values for a short example script, I see a
warning telling me that the system had no solution. I don't know why that
didn't show up when running the function within my larger script. However,
I do not receive a warning when the inputs are all zeros (save one), which
I assume also has no solution. I've attached the example script in case you
want to see what I'm talking about.
-Bridget
On Tue, Feb 14, 2012 at 9:11 AM, Dennis Shea <shea@ucar.edu> wrote:
> solve_linsys is an interface to the LAPACK subroutines
>
> C LAPACK ROUTINES (USE "AC" AS INPUT)
>
> CALL DGETRF( M, N, AC, LDA, IPIV, INFO )
> CALL DGETRI( N, AC, LDA, IPIV, WORK, LWORK, INFO )
> ===
>
> Not sure why this would return zeros.
>
> ==
>
> You would have to provide a sample data set and a script that
> cleanly calls the solve_linsys
>
> D
>
>
> On 2/13/12 4:03 PM, Bridget Thrasher wrote:
>
>> I am coming across a problem where the solve_linsys function is
>> returning all fill values. I have verified that the input arrays are
>> filled properly, and I am not receiving an error. It appears to be
>> correlated with the A matrix having a lot of values = 0. Is this correct
>> behavior?
>>
>> -Bridget
>>
>> --
>> Bridget Thrasher, PhD
>> Independent Contractor, Research Scientist
>>
>>
>>
>>
>> ______________________________**_________________
>> ncl-talk mailing list
>> List instructions, subscriber options, unsubscribe:
>> http://mailman.ucar.edu/**mailman/listinfo/ncl-talk<http://mailman.ucar.edu/mailman/listinfo/ncl-talk>
>>
>
-- Bridget Thrasher, PhD Independent Contractor, Research Scientist
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
This archive was generated by hypermail 2.1.8 : Fri Feb 17 2012 - 08:50:10 MST