Function N_VWrmsNorm

externalpackages/PVODE/source/nvector.cpp:473–493

real N_VWrmsNorm(N_Vector x, N_Vector w)
{
  integer N, N_global;
  real sum = ZERO, *xd, *wd, gsum;
  machEnvType machenv;

  N = x->length;
  N_global = x->global_length;
  xd = x->data;
  wd = w->data;
  machenv = x->machEnv;

  #pragma omp parallel for reduction(+: sum)
  for (integer i=0; i < N; i++) {
    real prodi = xd[i] * wd[i];
    sum += prodi * prodi;
  }

  gsum = PVecAllReduce(sum, 1, machenv);
  return(RSqrt(gsum / N_global));
}

real N_VMin(N_Vector x)
{