Hi,

What is the best way to solve a system of positive-definite linear equationsusing smath( as used in structural analysis of space frames)?

Elsewhere I have used software that makes use of public domain solver libraries for this. Can these be adapted for use in Smath?

Thanks

farlanw

Iterative methods are faster than direct for large systems,
for positive definite linear matrices I've read here that PCG ( "preconditioned conjugate gradient method" ) is quite efficient in general.
Wikipedia
Algorithm pseudocode

Wrote

Iterative methods are faster than direct for large systems,
for positive definite linear matrices I've read here that PCG ( "preconditioned conjugate gradient method" ) is quite efficient in general.
Wikipedia
Algorithm pseudocode

The sm code for linear systems.


regards,

w3b5urf3r

---

ConjugateGradient.sm (7 KiB) downloaded 222 time(s).