src/expde/grid/input.h

Go to the documentation of this file.

//    expde: expression templates for partial differential equations.
//    Copyright (C) 2000  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.

// ------------------------------------------------------------
//
// input.h
//
// ------------------------------------------------------------

   
#ifndef INPUT_H_
#define INPUT_H_
                   
#define factor_eps_input 0.00000001
#define search_length 3


// for hashtable

class Input_data_object;

class Point_hashtable_input {
 public:
  int nx, ny, nz;   // inter code

  double wert;

  Point_hashtable_input* next;

  Point_hashtable_input(int Nx, int Ny, int Nz, 
                        double value) 
    : nx(Nx), ny(Ny), nz(Nz), wert(value) {
    next = NULL;
  };
  bool Is_point(int Nx, int Ny,int Nz) { return (Nx==nx && Ny==ny && Nz==nz); }
  double Value() { return wert; }
  ~Point_hashtable_input() {
    delete next;
  };
};

inline int hashtable_input_function(int x, int y, int z, int length) {
  int sp;
  sp = (x + 101*y + 10007*z) % length;
  if(sp>0) return sp;
  return -sp;
 //  return (x + 101*y + 10007*z) % length;
}


// Input_data_object

class Input_data_object {
 public:
  Input_data_object(D3vector* coordinates, double* value, int num_points,
                    double mesh_x, double mesh_y, double mesh_z);

  double Interpolate_value(D3vector coordinates);
  
 private:
  D3vector* ort;
  double*   werte;
  
  double org_x;
  double org_y;
  double org_z;

  int anz_punkte;

  double h_x;             // Maschenweite x
  double h_y;             // Maschenweite y
  double h_z;             // Maschenweite z

  double epsx;
  double epsy;
  double epsz;

  // Punkt hinzufuegen usw.
  void Add_point(double x, double y, double z, double value);
  bool   Exists_point(D3vector coordinates);
  double Give_value(D3vector coordinates);
  D3vector Staggered_coordinate(D3vector coordinates);

  // Speicher des hashtables:
  Point_hashtable_input* point;              // fuer Iterationen
  int hashtable_input_leng;                  // Laenge des hashtable
  int hashtable_input_occ;                   // Besetzungszahl des hashtable
  Point_hashtable_input** hashtable_input_start;   // Zeiger auf Tabelle

  // gehe zu hashtable
  Point_hashtable_input* hashtable_input_point(int Nx, int Ny, int Nz) const;
};


inline Point_hashtable_input* 
Input_data_object::hashtable_input_point(int Nx, int Ny, int Nz) const {
  Point_hashtable_input* point;
  point = hashtable_input_start[hashtable_input_function(Nx,Ny,Nz,
                                 hashtable_input_leng)];
  while(point!=NULL) {
    if(point->nx==Nx && point->ny==Ny && point->nz==Nz) return point;
    point = point->next;
  }
  return NULL;
}


#endif

---