method infinityNorm

Source

double infinityNorm() {
  double norm = 0.0;
  {
    double row_norm = 0.0;
    row_norm += _m4storage[0].abs();
    row_norm += _m4storage[1].abs();
    row_norm += _m4storage[2].abs();
    row_norm += _m4storage[3].abs();
    norm = row_norm > norm ? row_norm : norm;
  }
  {
    double row_norm = 0.0;
    row_norm += _m4storage[4].abs();
    row_norm += _m4storage[5].abs();
    row_norm += _m4storage[6].abs();
    row_norm += _m4storage[7].abs();
    norm = row_norm > norm ? row_norm : norm;
  }
  {
    double row_norm = 0.0;
    row_norm += _m4storage[8].abs();
    row_norm += _m4storage[9].abs();
    row_norm += _m4storage[10].abs();
    row_norm += _m4storage[11].abs();
    norm = row_norm > norm ? row_norm : norm;
  }
  {
    double row_norm = 0.0;
    row_norm += _m4storage[12].abs();
    row_norm += _m4storage[13].abs();
    row_norm += _m4storage[14].abs();
    row_norm += _m4storage[15].abs();
    norm = row_norm > norm ? row_norm : norm;
  }
  return norm;
}