Mailing Lists: Apple Mailing Lists

Image of Mac OS face in stamp
 
[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

[apple scitech] Determinant refinement



Daniel Farrell:

If you want to try Martin Costabel's variant for 4x4 det, as:

======================================
  det = G(0,1,2,3) + G(0,2,1,3) + G(0,3,1,2) +
         G(1,2,0,3) + G(1,3,0,2) + G(2,3,0,1);


double G(int j, int k, int m, int n) { return( ( a[j][0] * a[k][1] - a[k][0] * a[j][1] ) * ( a[m][2] * a[n][3] - a[n][2] * a[m][3] ) ); }

======================================

then observe that function overhead could be reduced by unrolling the whole thing, instead of calling G six times...or in-lining G, something such. It would be good to know if the F- or G- scheme gives you the speed you want, for there are yet more speedups possible...it all becomes, as is typical, a tradeoff between research time and execution time, i.e. quality-of-life vs. quality-of-machine-life!

Richard Crandall
Advanced Computation Group



_______________________________________________
Do not post admin requests to the list. They will be ignored.
Scitech mailing list      (email@hidden)
Help/Unsubscribe/Update your Subscription:
http://lists.apple.com/mailman/options/scitech/email@hidden

This email sent to email@hidden


Visit the Apple Store online or at retail locations.
1-800-MY-APPLE

Contact Apple | Terms of Use | Privacy Policy

Copyright © 2007 Apple Inc. All rights reserved.