source: src/comm/comm.cpp@ dfed1c

/**
 * @file   comm.hpp
 * @author Julian Iseringhausen <isering@ins.uni-bonn.de>
 * @date   Wed Jun 16 13:21:06 2010
 *
 * @brief  Base class for communication.
 *
 */

#ifdef HAVE_CONFIG_H
#include <config.h>
#endif

#ifdef HAVE_BOOST_FILESYSTEM
#include <boost/filesystem.hpp>
namespace fs = boost::filesystem;
#endif

#include "base/discretization.hpp"
#include "base/helper.hpp"
#include "base/stencil.hpp"
#include "comm/comm.hpp"
#include "comm/domain_decomposition.hpp"
#include "grid/grid.hpp"
#include "mg.hpp"

using namespace VMG;

vmg_float Comm::ComputeResidualNorm(Multigrid& sol, Multigrid& rhs)
{
#ifdef DEBUG_MATRIX_CHECKS
  sol().IsCompatible(rhs());
  sol().IsConsistent();
  rhs().IsConsistent();
#endif

  Grid::iterator grid_iter;
  Stencil::iterator stencil_iter;
  vmg_float val;
  vmg_float norm = 0.0;

  const vmg_float prefactor = MG::GetDiscretization()->OperatorPrefactor(sol());
  const Stencil& A = MG::GetDiscretization()->GetStencil();

  this->CommToGhosts(sol());

  if (sol().Global().BoundaryType() == LocallyRefined)
    MG::GetDiscretization()->SetInnerBoundary(sol(), rhs(), sol(sol.Level()-1));

  for (grid_iter=rhs().Iterators().Local().Begin(); grid_iter!=rhs().Iterators().Local().End(); ++grid_iter) {

    val = rhs().GetVal(*grid_iter) - prefactor * A.GetDiag() * sol().GetVal(*grid_iter);

    for (stencil_iter=A.begin(); stencil_iter!=A.end(); ++stencil_iter)
      val -= prefactor * stencil_iter->Val() * sol().GetVal(*grid_iter + stencil_iter->Disp());

    norm += val*val;

  }

  norm = sqrt( Helper::pow_3(sol().MeshWidth()) * GlobalSum(norm) );

  return norm;
}

void Comm::InitComm()
{
}

Comm::~Comm()
{
  delete dd;
}

const std::string& Comm::OutputPath()
{
  if (!output_directory_is_created) {
    output_path_str = CreateOutputDirectory();
    output_directory_is_created = true;
  }

  return output_path_str;
}