[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
changeset: fsolve
From: |
Jaroslav Hajek |
Subject: |
changeset: fsolve |
Date: |
Fri, 31 Oct 2008 12:24:30 +0100 |
hello,
please consider the attached changeset, reimplementing fsolve as
m-file function and removing the old fsolve code, including MINPACK.
The advantages we get using an m-file include:
1. more understandable, modifiable and maintainable code (hopefully :)
2. reentrancy
3. will compute in single precision if given single precision inputs.
4. possibly faster for large problems, because it uses blocked QR from
LAPACK instead of MINPACK's own implementation.
The method is essentially the same as in MINPACK, i.e. Powell's hybrid
method (trust-region approach, updating jacobian by broyden updates if
possible).
regards
--
RNDr. Jaroslav Hajek
computing expert
Aeronautical Research and Test Institute (VZLU)
Prague, Czech Republic
url: www.highegg.matfyz.cz
fsolve.diff
Description: Text Data
- changeset: fsolve,
Jaroslav Hajek <=