ufunction.hpp

FreeFem++可以生成高质量的有限元网格。可以用于流体力学

// ORIG-DATE:     Dec 2007
// -*- Mode : c++ -*-
//
// SUMMARY  :  Model of usefull function generic function   
// USAGE    : LGPL      
// ORG      : LJLL Universite Pierre et Marie Curi, Paris,  FRANCE 
// AUTHOR   : Frederic Hecht
// E-MAIL   : [e-mail]
//

/*
 
 This file is part of Freefem++
 
 Freefem++ is free software; you can redistribute it and/or modify
 it under the terms of the GNU Lesser General Public License as published by
 the Free Software Foundation; either version 2.1 of the License, or
 (at your option) any later version.
 
 Freefem++  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 Lesser General Public License for more details.
 
 You should have received a copy of the GNU Lesser General Public License
 along with Freefem++; if not, write to the Free Software
 Foundation, Inc., 51 Franklin St, Fifth Floor, Boston, MA  02110-1301  USA

 Thank to the ARN ()  FF2A3 grant
 ref:ANR-07-CIS7-002-01 
 */

#ifndef UFUNCTION_HPP_
#define UFUNCTION_HPP_
// some usefull function
template<class T> inline T Min (const T &a,const T &b){return a < b ? a : b;}
template<class T> inline T Max (const T &a,const T & b){return a > b ? a : b;}
template<class T> inline T Abs (const T &a){return a <0 ? -a : a;}
template<class T> inline double Norme (const T &a){return sqrt(a*a);}
template<class T> inline void Exchange (T& a,T& b) {T c=a;a=b;b=c;}
template<class T> inline T Max (const T &a,const T & b,const T & c){return Max(Max(a,b),c);}
template<class T> inline T Min (const T &a,const T & b,const T & c){return Min(Min(a,b),c);}


template<int d,typename A,typename B>
struct Copy {
static A * copy(A a[d],B * const b)
{
  for (int i=0;i<d;++i)
    a[i]=b[i];
  return a;
}};

template<typename A,typename B>
struct Copy<1,A,B> {
static A * copy(A a[1],B * const b  )
{
    a[0]=b[0];
  return a;
}};

template<typename A,typename B>
struct Copy<2,A,B> {
static A * copy(A a[2],B * const b)
{
    a[0]=b[0];
    a[1]=b[1];

  return a;
}};

template<typename A,typename B>
struct Copy<3,A,B> {
static A * copy(A a[3],B * const b)
{
    a[0]=b[0];
    a[1]=b[1];
    a[2]=b[2];
  return a;
}};

template<typename A,typename B>
struct Copy<4,A,B> {
static A * copy(A a[4],B * const b)
{
    a[0]=b[0];
    a[1]=b[1];
    a[2]=b[2];
    a[3]=b[3];
  return a;
}};

struct UniqueffId {
private:
  static int count ;
  int id;
public:
  UniqueffId() : id(++count) {}
  bool operator==(UniqueffId u) const {return id==u.id;}
  bool operator!=(UniqueffId u) const {return id!=u.id;}
  void init(){id=++count;}
};

#endif