Livermore Loop loop09.sac sparse matrix multiply sensitive to argument ordering
| Bugzilla Link | 1051 |
| Created on | Mar 08, 2013 18:44 |
| Version | svn |
| OS | Linux |
| Architecture | PC |
Extended Description
This sparse matrix kernel in loop09.sac exhibits an undesirable sensitivity to
argument ordering:
inline double[.] vecmatmul( double[.] wts, double[.,.] PX) {
/* Sparse wts matrix multiply: wts f.g PX */
colsx = shape(wts)[0];
colsy = shape(PX)[1];
z = genarray([colsy], 0.0d);
for (colx=0; colx