src/expde/formulas/subdiv.h

Go to the documentation of this file.

//    expde: expression templates for partial differential equations.
//    Copyright (C) 2001  Christoph Pflaum
//    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; either version 2 of the License, or
//    (at your option) any later version.
//
//    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.
//
//                 SEE  Notice1.doc made by 
//                 LAWRENCE LIVERMORE NATIONAL LABORATORY
//

// ------------------------------------------------------------
//
// subdiv.h
//
//
// ------------------------------------------------------------

#ifndef SUBDIV_H_
#define SUBDIV_H_

inline void Tell_about_round_of_problem(double a, double b) {
  if(tell_round_off_problem) {
    if(ABS(a - b) < eps_round_off_problem)
      cout << " there may be a round off problem"
             << " in the subdivision of the boundary cells! " << endl;
  }
}


// Klasse zur Berechnung von FE auf Randzellen


class Calculator_FE {
public:
  void Start(Bo_description* desc, BoCeData* boc) {
    bo_desc = desc;
    bo_cell = boc;
  } 
  void Tetraeder(Edge_Corner_point p0, Edge_Corner_point p1,
                 Edge_Corner_point p2, Edge_Corner_point p3);
  // local_coord_normi ist lokale Koordinate in [-0.5,0.5]^3 Zelle
  void Add_Bo_freedom(D3vector local_coord_normi) {
    bo_cell->Add_Bo_freedom(local_coord_normi*bo_cell->Meshsize());
  };
private:
  Bo_description* bo_desc;
  BoCeData* bo_cell;
};


void Act_on_subdivision(Bo_description* Description, Calculator_FE* Act);


#endif
   

---