On 3/1/08, Pauli Virtanen <[EMAIL PROTECTED]> wrote: > A silly question: did you check directly that the pure-numpy code and > the F90 code give the same results for the Jacobian-vector product > J(z0) z for some randomly chosen vectors z0, z?
No, I did not do that. However, I've checked the output of of the finite diferencing routines for random X input of 32*32 and alpha=6.8, and the maximum difference is always 4.4408920985e-16. At first, this seems good. -- Lisandro Dalcín --------------- Centro Internacional de Métodos Computacionales en Ingeniería (CIMEC) Instituto de Desarrollo Tecnológico para la Industria Química (INTEC) Consejo Nacional de Investigaciones Científicas y Técnicas (CONICET) PTLC - Güemes 3450, (3000) Santa Fe, Argentina Tel/Fax: +54-(0)342-451.1594 _______________________________________________ Numpy-discussion mailing list Numpy-discussion@scipy.org http://projects.scipy.org/mailman/listinfo/numpy-discussion