}
// Test multiplication of sparse transpose by vector
double[] f = new double[as.nrows];
ranvec(f);
double[] g=as.atx(f); // Sparse transpose multiplied by vector
// Find same result by brute force
double[] h= new double[as.ncols];
for (i=0;i<as.ncols;i++) {
val=0.0;
for (j=0;j<as.nrows;j++) val += a[j][i]*f[j];