error.cpp

来自「OFELI is an object oriented library of C」· C++ 代码 · 共 17 行

#include "Mesh.h"
#include "Vect.h"
using namespace OFELI;

void error(double time, const Mesh &ms, const Vect<double> &u)
{
   double ee, e=0;
   double a = 10, U;
   const Node *nd;
   for (ms.TopNode(); (nd=ms.GetNode());) {
      U = tanh(a*nd->Coord(2))*(exp(time)-1);
      ee = U - u(nd->Label());
      e += ee*ee;
   }
   cout << "Discrete L2 error = " << sqrt(e/ms.NbNodes()) << endl;
}