On Monday, October 6, 2014 12:31:17 PM UTC-5, Andreas Noack wrote:
>
> There are fwdTriSolve! and bwdTriSolve! in linalg/sparse.jl, but there are
> not A_ldiv_B!(Triangular{SparsesCSC}, StridedVector) versions yet, but they
> should be fairly easy to add.
>
Thanks, Andreas. I thought I had seen them somewhere but I was looking for
such methods in the base/sparse directory instead of base/linalg.
>
> Med venlig hilsen
>
> Andreas Noack
>
> 2014-10-06 13:06 GMT-04:00 Douglas Bates <[email protected] <javascript:>>:
>
>> Is there an existing method for solving sparse triangular systems -
>> taking advantage of the sparseness and the triangularity?
>>
>> It seems that with the templated Triangular type taking the matrix type
>> as well as the element type it should be straightforward to create a
>> Triangular{Float64,SparseMatrixCSC{Float64},:L,false) type and incorporate
>> it in a method signature.
>>
>> I want to avoid reading the code in Csparse or in Tim Davis's book before
>> writing such a method because that code is GPL'd. Can anyone point me to
>> general references on sparse matrix methods other than Tim's book "Direct
>> Methods for Sparse Linear Systems".
>>
>
>
