src/expde/domain/sca_dom.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
//

// ------------------------------------------------------------
//
// sca_dom.h
//
// ------------------------------------------------------------

#ifndef SCADOM_H_
#define SCADOM_H_

#ifdef COMP_GNUOLD

#else
  using std::ifstream;
#endif 

// 1. domains described by scattered data which
//    are points at the boundary and normal vectors
// 2. domains described by voxels
class Scat_domain : public All_Domains  {
 protected:
 protected:
  // Konstruktor
  Scat_domain(int number);

  // storage for scattered data points (samples)
  int punkte;           // number of samples
  D3vector * sample, * normal;

  // interior bounding box
  D3vector bordermax;
  D3vector bordermin;

 public:
  // Konstruktor 
  Scat_domain(ifstream& Datei);

  // Groesse eines uniformen Grundgitters, Angabe in Levels:
  int  Give_n_uniform() { return n_uniform; } ;
  // Kann berechnet werden ob Kante im Gebiet? 
  bool calc_edge(D3vector V1, D3vector V2) { return true; }
  // ist Kante im Gebiet?       
  bool edge(D3vector V1, D3vector V2);
  // Kann der Abstand zum Rand in Richting d berechnet werden?
  calc_dis calc_distance(D3vector V, dir_3D d) {return yes; }
  // Abstand zum Rand in Richting d  
  double distance(D3vector V, dir_3D d);
  // ist Punkt im Gebiet? Diese Funktion sollte nur amAnfang verwendet werden.
  bool point_in_domain(D3vector V);

  // Beschreibung des Umgebenden Quadrates 
  D3vector GiveA() const { return A; }
  double   GiveH() const { return H; }
  D3vector GiveVecH() const { return VecH; }

  //
  int Give_number_of_scattered_data() { return punkte; };
};

class Six_corner_pyramid : public Scat_domain {
 public:
  Six_corner_pyramid(D3vector top1, D3vector top2, D3vector top3, 
                     D3vector top4, D3vector top5, D3vector top6,
                     D3vector down1, D3vector down2, D3vector down3, 
                     D3vector down4, D3vector down5, D3vector down6);
};



// some elementary member functions
// --------------------------------------

inline bool Scat_domain::edge(D3vector V1, D3vector V2) {
  return (point_in_domain(V1) && point_in_domain(V2));
}

#endif

        

---