AmrObject.cpp

Go to the documentation of this file.

//
// Class AmrObject
//   The AMR interface to OPAL. A new AMR library needs
//   to inherit from this class in order to work properly
//   with OPAL. Among other things it specifies the refinement
//   strategies.
//
// Copyright (c) 2016 - 2020, Matthias Frey, Paul Scherrer Institut, Villigen PSI, Switzerland
// All rights reserved
//
// Implemented as part of the PhD thesis
// "Precise Simulations of Multibunches in High Intensity Cyclotrons"
//
// This file is part of OPAL.
//
// OPAL 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.
//
// You should have received a copy of the GNU General Public License
// along with OPAL. If not, see <https://www.gnu.org/licenses/>.
//
#include "AmrObject.h"
 
AmrObject::AmrObject()
    : AmrObject(CHARGE_DENSITY, 0.75, 1.0e-15)
{ }
 
 
AmrObject::AmrObject(TaggingCriteria tagging,
                     double scaling,
                     double chargedensity)
    : tagging_m(tagging)
    , scaling_m(scaling)
    , chargedensity_m(chargedensity)
    , maxNumPart_m(1)
    , minNumPart_m(1)
    , refined_m(false)
    , amrSolveTimer_m(IpplTimings::getTimer("AMR solve"))
    , amrRegridTimer_m(IpplTimings::getTimer("AMR regrid"))
{ }
 
 
AmrObject::~AmrObject()
{ }
 
 
void AmrObject::setTagging(TaggingCriteria tagging) {
    tagging_m = tagging;
}
 
 
void AmrObject::setTagging(const std::string& tagging) {
    if ( tagging == "POTENTIAL" )
        tagging_m = TaggingCriteria::POTENTIAL;
    else if (tagging == "EFIELD" )
        tagging_m = TaggingCriteria::EFIELD;
    else if ( tagging == "MOMENTA" )
        tagging_m = TaggingCriteria::MOMENTA;
    else if ( tagging == "MAX_NUM_PARTICLES" )
        tagging_m = TaggingCriteria::MAX_NUM_PARTICLES;
    else if ( tagging == "MIN_NUM_PARTICLES" )
        tagging_m = TaggingCriteria::MIN_NUM_PARTICLES;
    else if ( tagging == "CHARGE_DENSITY" )
        tagging_m = TaggingCriteria::CHARGE_DENSITY;
    else
        throw OpalException("AmrObject::setTagging(std::string)",
                            "Not supported refinement criteria "
                            "[CHARGE_DENSITY | POTENTIAL | EFIELD | "
                            "MOMENTA | MAX_NUM_PARTICLES | MIN_NUM_PARTICLES].");
}
 
 
void AmrObject::setScalingFactor(double scaling) {
    scaling_m = scaling;
}
 
 
void AmrObject::setChargeDensity(double chargedensity) {
    chargedensity_m = chargedensity;
}
 
 
void AmrObject::setMaxNumParticles(size_t maxNumPart) {
    maxNumPart_m = maxNumPart;
}
 
 
void AmrObject::setMinNumParticles(size_t minNumPart) {
    minNumPart_m = minNumPart;
}
 
 
const bool& AmrObject::isRefined() const {
    return refined_m;
}
 
 
std::string AmrObject::enum2string(int number) {
    std::string tagging;
    switch ( number ) {
        case TaggingCriteria::CHARGE_DENSITY:
            tagging = "CHARGE_DENSITY";
            break;
        case TaggingCriteria::POTENTIAL:
            tagging = "POTENTIAL";
            break;
        case TaggingCriteria::EFIELD:
            tagging = "EFIELD";
            break;
        case TaggingCriteria::MOMENTA:
            tagging = "MOMENTA";
            break;
        case TaggingCriteria::MIN_NUM_PARTICLES:
            tagging = "MIN_NUM_PARTICLES";
            break;
        case TaggingCriteria::MAX_NUM_PARTICLES:
            tagging = "MAX_NUM_PARTICLES";
            break;
        default:
            throw OpalException("AmrObject::enum2string(int)",
                                "Only numbers between 0 and 5 allowed.");
    }
    return tagging;
}