Geom.h

Go to the documentation of this file.

// -*- C++ -*-

// Copyright (C) 2002 Ralf Deiterding
// Brandenburgische Universitaet Cottbus
//
// Copyright (C) 2003-2007 California Institute of Technology
// Ralf Deiterding, ralf@amroc.net

#ifndef AMROC_GEOM_H
#define AMROC_GEOM_H

template <class DataType, int dim>
class Geom {
public:
  typedef GridData<DataType,dim> grid_data_type;   
  typedef Vector<DataType,dim> point_type;

  Geom() {}

  int LocalList(GridHierarchy& GH, const BBox& bb, const int& nc, 
                const point_type* xc, bool* lused) {
    DCoords lbc = GH.worldCoords(bb.lower(), bb.stepsize());
    DCoords ubc = GH.worldCoords(bb.upper()+bb.stepsize(), bb.stepsize());
    int nl = nc;
    for (register int n=0; n<nc; n++) {
      lused[n] = true;
      for (register int d=0; d<dim; d++) 
        if (lbc(d)>xc[n](d) || ubc(d)<xc[n](d)) {
          lused[n] = false;
          nl--;
          break;
        }
    }
    return nl;
  }

  void CellCoords(GridHierarchy& GH, const Coords& ss, const int& nc, 
                  const int* idx, point_type* xc) {
    point_type dxh2;
    GH.worldStep(ss, dxh2.data());
    dxh2 *= static_cast<DataType>(0.5);

    for (register int n=0; n<nc; n++) {
      GH.worldCoords(&(idx[dim*n]), ss(), (DataType*)(xc[n].data()));
      xc[n] += dxh2;
    }
  }

  void CellIndices(GridHierarchy& GH, const Coords& ss, const int& nc, 
                   const point_type* xc, int* idx) {
    for (register int n=0; n<nc; n++) {      
      GH.localCoords(xc[n].data(), &(idx[dim*n]));
      for (register int d=0; d<dim; d++) {
        if (idx[n*dim+d]<0) idx[n*dim+d] -= ss(d)-1;
        idx[n*dim+d] = (idx[n*dim+d]/ss(d))*ss(d);
      }
    }
  }
};

#endif

---