vorticity.cpp

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

/*==============================================================================

   Copyright (C) 1998 - 2002 Rachid Touzani

   This program is free software; you can redistribute it and/or modify it under
   the terms of the GNU General Public License as published by the Free 
   Software Foundation; Version 2 of the License.

   This program is distributed in the hope that it will be useful, but WITHOUT
   ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS
   FOR A PARTICULAR PURPOSE. See the GNU General Public License for more 
   details.

   You should have received a copy of the GNU General Public License 
   along with this program; if not, write to the :

   Free Software Foundation
   Inc., 59 Temple Place - Suite 330
   Boston, MA  02111-1307, USA

  ==============================================================================*/

#include "Mesh.h"
#include "NodeVect.h"
#include "Vect.h"
#include "LocalVect.h"
#include "CalcElem.h"

using namespace OFELI;

void Vorticity(double q, const Mesh &ms, const Vect<> &u, const Vect<> &M,
               NodeVect<> &w)
//-----------------------------------------------------------------------------
//                            Calculate Vorticity
//-----------------------------------------------------------------------------
{
   const Element *e;
   LocalVect<double,3> dx, dy;
   double omega;
   w = 0;

   for (ms.TopElement(); (e=ms.GetElement());) {
      LocalVect<double,6> ue(e,u);
      CalcElem(*e,dx,dy);
      omega = dx[0]*ue[1] + dx[1]*ue[3] + dx[2]*ue[5] -
              dy[0]*ue[0] - dy[1]*ue[2] - dy[2]*ue[4];
      w(e->NodeLabel(1),1) += OFELI_SIXTH*q*omega;
      w(e->NodeLabel(2),1) += OFELI_SIXTH*q*omega;
      w(e->NodeLabel(3),1) += OFELI_SIXTH*q*omega;
   }

   for (size_t i=1; i<=ms.NbNodes(); i++)
      w(i,1) /= M(i);
}