[Top][All Lists]

[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

Re: [Help-gsl] Using dgbtrf / dgbtrs with odeintv2

From: Patrick Alken
Subject: Re: [Help-gsl] Using dgbtrf / dgbtrs with odeintv2
Date: Mon, 29 Aug 2016 09:05:19 -0600
User-agent: Mozilla/5.0 (X11; Linux x86_64; rv:45.0) Gecko/20100101 Thunderbird/45.2

There are no plans to do this as far as I know, but you're welcome to code it as a library extension and we can put it on the webpage for others to try out

On 08/27/2016 03:34 AM, Björn Dahlgren wrote:
Dear all,

I'd like to use the ODE solvers in GSL for large stiff problems where the
Jacobian matrix is banded (doing dense LU solve is intractable). Currently
this is not possible (unless I've missed something) since the factorization
and solving of the linear system is hardcoded to use `gsl_linalg_LU_decomp`
and `gsl_linalg_LU_solve` in functions with internal linkage (grep:ing for
these reveals use in bsimp.c, modnewton1.c and msbdf.c). So currently, my
only option (it would seem), would be to recompile a modified version of

Are there any plans to allow for banded solve? (it could be an "advanced"
feature without any fancy interface, just somewhere to pass in function
pointers for factorization and solving routines of the linear system, the
responsibility of handling padding and matrix dimensions could be left to
the user).

Best regards,
Björn Dahlgren, KTH - Royal Institute of Technology

reply via email to

[Prev in Thread] Current Thread [Next in Thread]