Re: Jacobi Method

Searching for "jacobi method" finds passing mention in PDLA::MatrixOps, but if you perhaps asked a more specific question with sample code showing what you're attempting to do you might get a more germane (or possibly even useful) reply.

As for your second question, yes; see perlxs or Inline::C for example. In general you'd want to search CPAN first because if you're wanting to use an open source C library there's a good chance someone's written a module wrapping it already.

The cake is a lie.
The cake is a lie.
The cake is a lie.

Comment on Re: Jacobi Method

Replies are listed 'Best First'.
Re^2: Jacobi Method by Lucero (Novice) on Sep 30, 2021 at 00:35 UTC
Thank you. What I have to do is find eigenvalues and eigenvectors using that method but I have been asked to use certain routine from the book Numerical Recipes in C.	[reply]
Re^3: Jacobi Method by etj (Priest) on Oct 31, 2021 at 19:34 UTC
I think https://metacpan.org/pod/PDL::MatrixOps#eigens looks like what you want. Please don't use PDLA, it's effectively deprecated (I should update it accordingly).	[reply]
Re^3: Jacobi Method by etj (Priest) on Apr 18, 2022 at 19:27 UTC
An additional PDL module that uses Jacobians within a Levenberg-Marquardt approach is PDL::Fit::Levmar.	[reply]
Re^3: Jacobi Method by perlfan (Parson) on Oct 12, 2021 at 15:58 UTC
FWIR, jacobi (and SOR, Gauss-Seidel) are basically just averaging values over time steps. How many dimensions are you dealing with and what is the PDE you're solving? PDL almost surely can help you out; but even if not it'd be pretty easy to implement this in pure Perl^tm.	[reply]