TestDensityCalculation.py

Go to the documentation of this file.

class TestDensityCalculation:
    """
    Initialisation snippet for testing the calculation of the SPH density.
    Returns a string with C++ code.
    """
 
    def __init__(self):
        self.initialisation_snippet = """
    // Start of template in src/exaparty/api/scenario/TestDensityCalculation.py
 
    /*
     * Most of the initialised variables are dummy (i.e. not used)
     */
 
    // Internal search radious used by Peano and with max value R = grid_size/2
 
    // Internal search radius used by Peano and with max value R = grid_size/2
    tarch::la::Vector<Dimensions,double> GridSize = marker.h();
    const double GridSizeHalf_x = 0.5 * GridSize(0);
 
    const double h_max = particle->getSmlMax();
    const double H_max = h_max * kernel_gamma;
    double search_radius = 0.999 * std::min(H_max, GridSizeHalf_x);
 
    particle->setSearchRadius( search_radius);
 
    particle->setDensity ( 1.0 );
    particle->setMass( particle->getDensity() / std::pow(HYDRO_PART_NUMBER, HYDRO_DIMENSIONS) );
 
    particle->setU( 1e-1 );
    particle->setUDot( 0 );
    particle->setU_full( 1e-1 );
    particle->setPressure( 1e-1 );
    particle->setSoundSpeed( 1e-1 );
 
    /* Compute initial smoothing length */
    const double eta_fac = particle->getEtaFactor();
    particle->setSmoothingLength( eta_fac * std::pow(particle->getMass()/particle->getDensity(), 1.0 / HYDRO_DIMENSIONS) );
 
    /* Initialize remaining fields  */
 
    // Vectors
    for(int d=0; d<Dimensions; d++){
       particle->setV          (d,0);
       particle->setV_full     (d,0);
       particle->setA          (d,0);
    }
 
    particle->setRot_v(0);
    particle->setHDot  (0.0);
    particle->setWcount(0.0);
    particle->setF     (0.0);
    particle->setRho_dh(0.0);
    particle->setWcount_dh(0.0);
 
    particle->setSmlIterLeftBound(particle->getSmlMin());
    particle->setSmlIterRightBound(particle->getSmlMax());
 
    particle->setHasNoNeighbours(false);
    particle->setCellHasUpdatedParticle(false);
    particle->setSmoothingLengthIterCount(0);
 
    // Artificial viscosity scheme
    particle->setDiv_v   (0.0);
    particle->setV_sig_AV(2.0 * particle->getSoundSpeed());
    particle->setBalsara (0.0);
 
#if PeanoDebug > 0
    particle->setPartid(1);
#endif
 
    // End of inserted block
 
    """