aboutsummaryrefslogtreecommitdiff
path: root/CarpetDev/CarpetIOF5/src/topology.cc
diff options
context:
space:
mode:
authorErik Schnetter <schnetter@cct.lsu.edu>2010-08-24 12:04:51 -0400
committerBarry Wardell <barry.wardell@gmail.com>2011-12-14 18:21:12 +0000
commita62e8ee3f1e86100cda633f51c22ee063e28df76 (patch)
tree96c4522bba984cfb4a1e53ce75733d6cdbf08c01 /CarpetDev/CarpetIOF5/src/topology.cc
parenteb813afcef1be0aba7154fe41d21e4a7657371ad (diff)
CarpetIOF5: Rename old thorn
Diffstat (limited to 'CarpetDev/CarpetIOF5/src/topology.cc')
-rw-r--r--CarpetDev/CarpetIOF5/src/topology.cc200
1 files changed, 0 insertions, 200 deletions
diff --git a/CarpetDev/CarpetIOF5/src/topology.cc b/CarpetDev/CarpetIOF5/src/topology.cc
deleted file mode 100644
index bc05c73a7..000000000
--- a/CarpetDev/CarpetIOF5/src/topology.cc
+++ /dev/null
@@ -1,200 +0,0 @@
-#include <cassert>
-#include <sstream>
-
-#include <hdf5.h>
-
-#include "cctk.h"
-
-#include "defs.hh"
-#include "vect.hh"
-
-#include "topology.hh"
-
-
-
-namespace CarpetIOF5 {
-
- namespace F5 {
-
- topology_t::
- topology_t (simulation_t & simulation)
- : m_simulation (simulation)
- {
- }
-
-
-
- topology_t::
- ~ topology_t ()
- {
- }
-
-
-
- simulation_t & topology_t::
- get_simulation ()
- const
- {
- return m_simulation;
- }
-
-
-
- hid_t topology_t::
- get_hdf5_topology()
- const
- {
- return m_hdf5_topology;
- }
-
-
-
- bool topology_t::
- invariant()
- const
- {
- return m_hdf5_topology >= 0;
- }
-
-
-
- unigrid_topology_t::
- unigrid_topology_t (simulation_t & simulation)
- : topology_t (simulation)
- {
- char const * const name = "uniform";
- m_name = string (name);
- m_hdf5_topology
- = open_or_create_group (m_simulation.get_hdf5_simulation(), name);
- assert (m_hdf5_topology >= 0);
-
- assert (invariant());
- }
-
-
-
- unigrid_topology_t::
- ~ unigrid_topology_t ()
- {
- herr_t const herr = H5Gclose (m_hdf5_topology);
- assert (not herr);
- }
-
-
-
- bool unigrid_topology_t::
- invariant()
- const
- {
- return topology_t::invariant();
- }
-
-
-
- mesh_refinement_topology_t::
- mesh_refinement_topology_t (simulation_t & simulation,
- int const map,
- int const maps,
- int const refinement_level,
- int const max_refinement_levels,
- vect<int, dim> const & level_refinement_factor,
- vect<int, dim> const & max_refinement_factor)
- : topology_t (simulation),
- m_refinement_level (refinement_level),
- m_max_refinement_levels (max_refinement_levels),
- m_level_refinement_factor (level_refinement_factor),
- m_max_refinement_factor (max_refinement_factor)
- {
- assert (refinement_level >= 0);
- assert (refinement_level < max_refinement_levels);
- assert (all (level_refinement_factor > 0));
- assert (all (level_refinement_factor <= max_refinement_factor));
-
- ostringstream namebuf;
- namebuf << "Vertices";
- if (maps > 1)
- {
- namebuf << "-map" << map;
- }
- if (max_refinement_levels > 1)
- {
- namebuf << "-level" << refinement_level;
- }
- string const namestr = namebuf.str();
- m_name = namestr;
- char const * const name = namestr.c_str();
-
- m_hdf5_topology
- = open_or_create_group (m_simulation.get_hdf5_simulation(), name);
- assert (m_hdf5_topology >= 0);
-
- assert (invariant());
- }
-
-
-
- void mesh_refinement_topology_t::
- calculate_level_origin_delta (vect<CCTK_REAL, dim> const & coarse_origin,
- vect<CCTK_REAL, dim> const & coarse_delta,
- vect<int, dim> const & level_offset,
- vect<int, dim> const & level_offset_denominator,
- vect<CCTK_REAL, dim> & level_origin,
- vect<CCTK_REAL, dim> & level_delta)
- const
- {
- assert (all (coarse_delta > (CCTK_REAL) 0));
- assert (all (level_offset_denominator > 0));
-
- level_delta
- = coarse_delta / vect<CCTK_REAL, dim> (m_level_refinement_factor);
- level_origin
- = (coarse_origin
- + (vect<CCTK_REAL, dim> (level_offset)
- / vect<CCTK_REAL, dim> (level_offset_denominator)
- / vect<CCTK_REAL, dim> (m_max_refinement_factor)));
- }
-
-
-
- mesh_refinement_topology_t::
- ~ mesh_refinement_topology_t ()
- {
- herr_t const herr = H5Gclose (m_hdf5_topology);
- assert (not herr);
- }
-
-
-
- void topology_t::
- get_link_destination (int const proc,
- string & filename,
- string & objectname)
- const
- {
- get_simulation().get_link_destination (proc, filename, objectname);
- if (objectname.empty())
- {
- objectname = m_name;
- }
- else
- {
- objectname += string ("/") + m_name;
- }
- }
-
-
-
- bool mesh_refinement_topology_t::
- invariant()
- const
- {
- return (topology_t::invariant()
- and m_refinement_level >= 0
- and m_refinement_level < m_max_refinement_levels
- and all (m_level_refinement_factor > 0)
- and all (m_level_refinement_factor <= m_max_refinement_factor));
- }
-
- } // namespace F5
-
-} // namespace CarpetIOF5