Function N_VMin

externalpackages/PVODE/source/nvector.cpp:495–528

}

real N_VMin(N_Vector x)
{
  integer N;
  real min, *xd, gmin;
  machEnvType machenv;

  N = x->length;
  xd = x->data;
  min = xd[0];
  machenv = x->machEnv;

#ifndef _OPENMP
  for (integer i=1; i < N; i++, xd++) {
    if ((*xd) < min) min = *xd;
  }
#else
  #pragma omp parallel 
  {
    real min2 = min;
    #pragma omp for nowait
    for (integer i=0; i < N; i++) {
      if (xd[i] < min) min = xd[i];
    }
    #pragma omp critical
    {
      if(min2 < min)
        min = min2;
    }
  }
#endif

  gmin = PVecAllReduce(min, 3, machenv);
  return(gmin);
}


real N_VWL2Norm(N_Vector x, N_Vector w)