public static double OneNorm(double[] store, int offset, int rowStride, int colStride, int nRows, int nColumns)
{
double norm = 0.0;
for (int c = 0; c < nColumns; c++)
{
double csum = Blas1.dNrm1(store, offset + colStride * c, rowStride, nRows);
if (csum > norm)
{
norm = csum;
}
}
return(norm);
}
public static double InfinityNorm(double[] store, int offset, int rowStride, int colStride, int nRows, int nColumns)
{
double norm = 0.0;
for (int r = 0; r < nRows; r++)
{
double rsum = Blas1.dNrm1(store, offset + rowStride * r, colStride, nColumns);
if (rsum > norm)
{
norm = rsum;
}
}
return(norm);
}
