/* * vmg - a versatile multigrid solver * Copyright (C) 2012 Institute for Numerical Simulation, University of Bonn * * vmg 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 3 of the License, or * (at your option) any later version. * * vmg 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. * * You should have received a copy of the GNU General Public License * along with this program. If not, see . */ /** * @file com_force_discrete_compatibility.cpp * @author Julian Iseringhausen * @date Mon Apr 18 12:37:11 2011 * * @brief May be used to explicitly force the discrete * compatibility condition \sum_{i,j,k} f_{i,j,k} = 0. * */ #ifdef HAVE_CONFIG_H #include #endif #include "base/command.hpp" #include "grid/grid.hpp" #include "grid/multigrid.hpp" #include "mg.hpp" using namespace VMG; class VMGCommandForceDiscreteCompatibility : public Command { public: Request Run(Command::argument_vector arguments) { MPE_EVENT_BEGIN() Grid& rhs = (*MG::GetRhs())(MG::GetRhs()->MaxLevel()); Comm& comm = *MG::GetComm(); if (comm.BoundaryConditions()[0] == Periodic && comm.BoundaryConditions()[1] == Periodic && comm.BoundaryConditions()[2] == Periodic) rhs.ForceDiscreteCompatibilityCondition(); MPE_EVENT_END() return Continue; } static const char* Name() {return "ForceDiscreteCompatibility";} static int Arguments() {return 0;} }; CREATE_INITIALIZER(VMGCommandForceDiscreteCompatibility)