Homec4science
Diffusion LibMultiScale (features/pybind) Browse

libmultiscale/src/filter/filter.cc
b05d5f2aba78features/pybind

Last Change
Raw File

filter.cc
View Options

/**
* @file filter.cc
*
* @author Guillaume Anciaux <guillaume.anciaux@epfl.ch>
*
* @date Tue Jul 29 17:03:26 2014
*
* @brief This is the mother class of all filters
*
* @section LICENSE
*
* Copyright (©) 2010-2011 EPFL (Ecole Polytechnique Fédérale de Lausanne)
* Laboratory (LSMS - Laboratoire de Simulation en Mécanique des Solides)
*
* LibMultiScale is free software: you can redistribute it and/or modify it under the
* terms of the GNU Lesser General Public License as published by the Free
* Software Foundation, either version 3 of the License, or (at your option) any
* later version.
*
* LibMultiScale 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 Lesser General Public License for more
* details.
*
* You should have received a copy of the GNU Lesser General Public License
* along with LibMultiScale. If not, see <http://www.gnu.org/licenses/>.
*
*/
#include "lm_common.hh"
#include "lib_dumper.hh"
#include "lib_filter.hh"
#include "lib_stimulation.hh"
#include "filter_manager.hh"
#include "action_manager.hh"
#include "factory_multiscale_inline_impl.hh"
#include "lib_md.hh"
#include "lib_continuum.hh"
#include "lib_dd.hh"
/* -------------------------------------------------------------------------- */
__BEGIN_LIBMULTISCALE__
template <>
void Filter<ContainerArray<Real> >::connect(FilterManager & f){
LM_FATAL("internal error");
}
/* -------------------------------------------------------------------------- */
template <>
void Filter<ContainerArray<Real> >::connect(ActionManager & f){
LM_FATAL("internal error");
}
/* -------------------------------------------------------------------------- */
template <>
void Filter<ContainerArray<RefPointData> >::connect(ActionManager & f){
LM_FATAL("internal error");
}
/* -------------------------------------------------------------------------- */
template <>
void Filter<ContainerArray<RefPointData> >::connect(FilterManager & f){
LM_FATAL("internal error");
}
/* -------------------------------------------------------------------------- */
template <typename _Input>
void Filter<_Input>::connect(FilterManager & f){
f.template create<Filter<_Input> >(*this);
}
/* -------------------------------------------------------------------------- */
template <typename _Input>
void Filter<_Input>::connect(ActionManager & f){
f.template create<Filter<_Input> >(*this);
}
/* -------------------------------------------------------------------------- */
template <typename _Input>
void Filter<_Input>::visit(ComponentLMOut<_Input> & obj){
this->copyContainerInfo(obj.getOutput());
if (obj.getOutput().getRelease() != UINT_MAX &&
obj.getOutput().getRelease() <= this->getRelease()) return;
this->copyReleaseInfo(obj);
Communicator & comm = Communicator::getCommunicator();
CommGroup group = obj.getOutput().getCommGroup();
if (comm.amIinGroup(group))
this->build(obj);
}
/* -------------------------------------------------------------------------- */
template <typename _Input>
void Filter<_Input>::build(){
if (!this->input) {
DUMP("this filter has no input",DBG_WARNING);
}
else this->input->accept(*this);
}
/* -------------------------------------------------------------------------- */
template <typename _Input>
void Filter<_Input>::manualBuild(_Input & obj){
this->copyContainerInfo(obj);
if (obj.getRelease() != UINT_MAX &&
obj.getRelease() <= this->getRelease()) return;
this->copyReleaseInfo(obj);
Communicator & comm = Communicator::getCommunicator();
CommGroup group = obj.getCommGroup();
if (comm.amIinGroup(group))
this->build(obj);
}
/* -------------------------------------------------------------------------- */
/* LMDESC FILTER
This is a generic internal filterf
*/
template <typename _Input>
void Filter<_Input>::declareParams(){
/* LMKEYWORD DOF_TYPE
Allow to specify the dof type to iterate on: useful only for parallelism.
The possible choices are:
\begin{itemize}
\item dt\_local\_master : only the local or master dofs
\item dt\_local\_slave : only the slave dofs
\item dt\_local : all local dofs
\item dt\_ghost : only the (pure) ghost dofs
\item dt\_all : all dofs
\end{itemize}
*/
this->parseKeyword("DOF_TYPE", dt, dt_local);
}
/* -------------------------------------------------------------------------- */
#define BOOST_INSTANCIATE_FILTER_CLASS_TEMPLATE(r,data,x) \
template class Filter<BOOST_PP_TUPLE_ELEM(3,0,x)::ContainerSubset>; \
template class Filter<BOOST_PP_TUPLE_ELEM(3,0,x)::ContainerPoints>;
#define DECLARE_FILTER_CLASS(list) \
BOOST_PP_SEQ_FOR_EACH(BOOST_INSTANCIATE_FILTER_CLASS_TEMPLATE,"",list)
/* -------------------------------------------------------------------------- */
DECLARE_FILTER_CLASS(LIST_ATOM_MODEL);
DECLARE_FILTER_CLASS(LIST_CONTINUUM_MODEL);
DECLARE_FILTER_CLASS(LIST_DD_MODEL);
template class Filter<ContainerArray<RefPointData> >;
template class Filter<ContainerGenericMesh>;
__END_LIBMULTISCALE__
c4science · Help