http://petsc.cs.iit.edu/petsc/petsc-dev/rev/a8a483b98169
This baffles me. I can think of no good reason for this, which gives me the impression that we are optimizing for an implementation quirk. If you have evidence that the performance of BLAS dot() is better than nrm2() across platforms and implementations, then we are witnessing a major implementation failure and people need to be shamed. Aliasing is also *explicitly disallowed* by Fortran, so the result of BLASdot_(&bn,xx,&one,xx,&one); is not defined. -------------- next part -------------- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20120103/25e1da72/attachment.html>
