Re: Matrix Library Message #2 Posted by Paul Dale on 19 Oct 2011, 4:30 a.m., in response to message #1 by John Abbott (S. Africa)
Quote: I have found an error (finger trouble probably) in the matrix library - namely in wp34s-1.dat.
Won't be finger trouble -- these are generated by the assembler from the sources in the library. Still, the source only has one 'M-' label so something has gone awry.
Quote: For Pauli - in the 'M-y' routine you execute a [cmplx]DROP, which does not leave a descriptor on the stack. Should this not be just a DROP?
The matrix descriptor is saved in L. The value on the stack isn't the descriptor anymore -- it is a loop index for ISG that has finished. Thus, complex DROP is correct here.
Quote: BTW- Thanks for these routines. The LU routine (if included) would be great :)
The LU routine will be back if we can find enough space for the wrapper. We're working on it.
- Pauli
|