I just saw a paper on ArXiv about making the complex division routine
more robust:
http://arxiv.org/abs/1210.4539
Second author is actually the original inventor of the current algorithm
in Lapack. I have attached my modified DLADIV routine, which passes all
the tests in the build process.
--
Victor Liu
Stanford University
Ph.D. Candidate, Dept. of Electrical Engineering
Ginzton Laboratory, Nano Building
Mail Code 4088
348 Via Pueblo Mall, Stanford, CA 94305-4088
Phone: (408) 482-6910
Web: http://cylith.dyndns.org/
-------------- next part --------------
A non-text attachment was scrubbed...
Name: dladiv_robust.f
Type: text/x-fortran
Size: 6116 bytes
Desc: not available
Url :
http://lists.eecs.utk.edu/mailman/private/lapack/attachments/20121017/2f8eab2c/dladiv_robust.bin
|