error.cpp

OFELI is an object oriented library of C++ classes for development of finite element codes. Its main

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

void error(const Mesh &ms, User &ud, const Vect<complex<double> > &u)
{
   double ee=0, e;
   const Node *nd;
   for (ms.TopNode(); (nd=ms.GetNode());) {
      complex<double> U = ud.ExactSolution(nd->Coord());
      e = abs(U - u(nd->Label()));
      ee += e*e;
   }
   cout << "Discrete L2 error = " << sqrt(ee/ms.NbNodes()) << endl;
}