For a mostly unused pure PDL approach, see PDL::Fit::Householder<http://search.cpan.org/%7Edcmertens/PDL-Fit-Householder-v0.0.2/lib/PDL/Fit/Householder.pd>. You'll need to install Module::Build::PDL; at this point, if I really wanted to polish off this module, I'd rewrite it to only need Module::Build.
The documentation discusses how to solve a linear combination of functions on a single dataset. You need to diagonalize a matrix with a different constructions, in which case you'll probably use the _Householder method instead of the Householder method. I'm not in a position to discuss how you would go about applying this method to solve you're problem, but at least you'll be aware that such a tool exists to solve the problem. David On Sun, Jul 15, 2012 at 5:49 PM, Joel Berger <[email protected]>wrote: > Mark, I would think you might want to look into a linear algebra style > solution to your system. If you wanted to do a root-finding solution, > there is a GSL library for multidimensional root-finding > > http://www.gnu.org/software/gsl/manual/html_node/Multidimensional-Root_002dFinding.html > > I haven't evaluated it for inclusion into the PerlGSL namespace yet, > but you are welcome to look at the examples and give it a try. > > Cheers, > Joel > > On Sun, Jul 15, 2012 at 4:35 PM, MARK BAKER <[email protected]> > wrote: > > I made a error in the math in the prior email > > below is the multiplication of the unknowns > > the prior email was of addition of the x,x z variables ... > > > > Hey all, > > > > does any one in the development have any idea on if > > it would be possible to add functionality like this to PDL > > > > > > p $k = pdl([[ 1x+2y+3z=47], [3x+2z+1y=41],[2x+1z+3x=38]]); > > > > > > x = 5 > > y = 6 > > z = 10 > > > > does any one think this might be possible; > > it would make simultaneous equations a snap ... > > > > or can someone point me to the source module so I can play with it.. > > > > Thanks > > > > -Mark > > > > > > _______________________________________________ > > Perldl mailing list > > [email protected] > > http://mailman.jach.hawaii.edu/mailman/listinfo/perldl > > > > > > > > > > > > _______________________________________________ > > Perldl mailing list > > [email protected] > > http://mailman.jach.hawaii.edu/mailman/listinfo/perldl > > > > _______________________________________________ > Perldl mailing list > [email protected] > http://mailman.jach.hawaii.edu/mailman/listinfo/perldl > -- "Debugging is twice as hard as writing the code in the first place. Therefore, if you write the code as cleverly as possible, you are, by definition, not smart enough to debug it." -- Brian Kernighan
_______________________________________________ Perldl mailing list [email protected] http://mailman.jach.hawaii.edu/mailman/listinfo/perldl
