amroc/amr/StatCPTLevelSet.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@cacr.caltech.edu

#ifndef AMROC_STATCPTLEVELSET_H
#define AMROC_STATCPTLEVELSET_H

#include "CPTLevelSet.h"

template <class DataType, int dim>
class StatCPTLevelSet : public CPTLevelSet<DataType,dim> {
  typedef CPTLevelSet<DataType,dim> base;

public:
  typedef typename base::grid_fct_type grid_fct_type;   
  typedef typename base::grid_data_type grid_data_type;   
  typedef typename cpt::State<dim,DataType> cpt_type;
  typedef ads::FixedArray<dim,DataType> point_type;
  typedef ads::FixedArray<dim,int> multi_index_type;

  StatCPTLevelSet() : base(), brep_filetype(0), vertices(0), connections(0) {
    base::_Stationary = 1;
  }

  virtual ~StatCPTLevelSet() {    
    if (vertices!=0) delete[] vertices;
    if (connections!=0) delete[] connections;
  }

  virtual void register_at(ControlDevice& Ctrl, const std::string& prefix) {
    base::LocCtrl = Ctrl.getSubDevice(prefix+"StatCPT");
    RegisterAt(base::LocCtrl,"Brep",brep_filename);
    RegisterAt(base::LocCtrl,"FileType",brep_filetype);
    RegisterAt(base::LocCtrl,"PlotPhi",base::_PlotPhi);
    RegisterAt(base::LocCtrl,"Inverse",base::inverse);
    RegisterAt(base::LocCtrl,"Unsigned",base::unsign);
    RegisterAt(base::LocCtrl,"FillWidth",base::_FillWidth);
  } 
  virtual void register_at(ControlDevice& Ctrl) {
    register_at(Ctrl, "");
  }

  virtual void SetupData(GridHierarchy* gh, const int& ghosts) {
    base::SetupData(gh, ghosts);

    std::ifstream brep_file( brep_filename.c_str() );
    if ( ! brep_file ) {
      std::cerr << "Could not read the b-rep file.  Exiting...\n";
      std::exit( 1 );
    }

    int num_vertices, num_connections;

    // ISS file
    if (brep_filetype==1) {
      int fdim, mdim;
      brep_file >> fdim >> mdim;
      if ( fdim!=dim || mdim!=dim-1 ) {
        std::cerr << "Mesh in b-rep file has wrong dimension.  Exiting...\n";
        std::exit( 1 );
      }
      brep_file >> num_vertices;
      vertices = new point_type[ num_vertices ];
      const point_type* point_iter_end = vertices + num_vertices;
      for ( point_type* piter = vertices; piter != point_iter_end; ++piter )
        brep_file >> *piter;
      brep_file >> num_connections;
      connections = new multi_index_type[ num_connections ];
    }
    else {
      brep_file >> num_vertices >> num_connections;
      vertices = new point_type[ num_vertices ];
      connections = new multi_index_type[ num_connections ];
      const point_type* point_iter_end = vertices + num_vertices;
      for ( point_type* piter = vertices; piter != point_iter_end; ++piter )
        brep_file >> *piter;
    }
      
    const multi_index_type* connections_iter_end = connections + num_connections;
    for ( multi_index_type* eiter = connections; eiter != connections_iter_end; ++eiter ) 
      brep_file >> *eiter;

    assert( brep_file.good() );
    brep_file.close();
    
    base::SetBrep(num_vertices, reinterpret_cast<const DataType*>(vertices), 
                  num_connections, reinterpret_cast<const int*>(connections));

  }

  inline void SetFilename(const std::string filename) { brep_filename = filename; }
  inline std::string Filename() { return brep_filename; }

protected:
  std::string brep_filename;
  int brep_filetype;
  point_type* vertices;
  multi_index_type* connections;
};


#endif

---