I added some documentation to the header file:
https://github.com/NanoComp/meep/blob/master/src/bicgstab.hpp
<https://github.com/NanoComp/meep/blob/master/src/bicgstab.hpp>
There are probably more sophisticated implementations of this algorithm
elsewhere. It doesn't accept a preconditioner, for example.
(Of course, you could use a left preconditioner M by transforming Ax=b to MAx =
Mb, but there are other ways to include preconditioning in biCGSTAB that may be
better, if I recall correctly.)
> On Sep 14, 2019, at 3:44 PM, Mert TAS <mert....@gmail.com> wrote:
>
> Hi everyone,
>
> I don't know if this is asked before since I cannot see the archives.
> I need to use the bicstabL() method, but I could not figure out how to call
> it:
>
> ptrdiff_t bicgstabL(const int L, const size_t n, realnum *x, bicgstab_op A,
> void *Adata,
> const realnum *b, const double tol,
> int *iters, // input *iters = max iters, output = actual iters
> realnum *work, // if you pass work=NULL, bicgstab
> returns nwork const bool quiet);
>
> Could you please show me an example of its usage and explain the required
> parameters.
> Also is it possible to provide any preconditioners to this method like ILU or
> Jacobi?
>
> Thanks.
>
>
>
> _______________________________________________
> meep-discuss mailing list
> meep-discuss@ab-initio.mit.edu
> http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
_______________________________________________
meep-discuss mailing list
meep-discuss@ab-initio.mit.edu
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
