MultiscaleTransitionsTest.h

Go to the documentation of this file.

// This file is part of the Peano project. For conditions of distribution and
// use, please see the copyright notice at www.peano-framework.org
#pragma once
 
#include "tarch/logging/Log.h"
#include "tarch/tests/TestCase.h"
 
namespace toolbox {
  namespace particles {
    namespace tests {
      class MultiscaleTransitionsTest;
    } // namespace tests
  }   // namespace particles
} // namespace toolbox
 
class toolbox::particles::tests::MultiscaleTransitionsTest: public tarch::tests::TestCase {
private:
  static tarch::logging::Log _log;
 
  void testSievePredicate();
 
  void testLiftDropOfParticleAssociatedWithVertex01();
 
 
  /*
   * Another test case
   *
   * We start with this lift operation:
   *
   * ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  77929644137  00:01:17     rank:0       core:1       info         benchmarks::exaparty::noh::observers::StepHydroPartDrift_sweep12peano4_toolbox_particles_api_UpdateParticleGridAssociation_LiftDrop1::touchVertexLastTime(...)                   have to lift particle (debugX=[0.505,0.165],debugH=[0,0],x=[0.504
 975,0.1667],cellH=[0.037037,0.037037],searchRadius=0.0185185,ParallelState=Local,NewParallelState=Local,MoveState=Moved,CellHasUpdatedParticle=1,cfl=0.1,initialTimeStepSize=0.0001,adjustTimeStepSize=0,hydroDimensions=2,etaFactor=1.2348,smlMin=1e-06,smlMax=0.166667,smlTolerance=0.0001,smlMaxIterations=30,a
 lphaAV=0.8,betaAV=3,mass=2.5e-05,v=[-0.0149237,0.999888],a=[7.92017e-06,-0.000467427],density=1.00459,pressure=6.71868e-07,smoothingLength=0.00615988,u=1.0032e-06,uDot=2.00584e-06,f=-1.07003e-07,wcount_dh=-55604.5,rho_dh=-1.39011,wcount=40183.6,hDot=-0.00923729,smoothingLengthIterCount=0,balsara=0.799993,
 rot_v=8.2102e-06,div_v=-2.9864,v_sig_AV=0.0916537,soundSpeed=0.00105578,v_full=[-0.0149237,0.999888],u_full=1.0033e-06,hasNoNeighbours=0,isBoundaryParticle=0,partid=1,dependencyChecksPeanoEventUsedBySwift=touchVertexFirstTime,dependencyChecksAlgorithmStepLastUpdated=FlagBoundaryParticles,dependencyChecksA
 lgorithmStepUpdates=1,...,dependencyChecksAlgorithmStepMaskOuts=0,...,dependencyChecksInitStepLastUpdated=FlagBoundaryParticles,dependencyChecksInitStepUpdates=1,...,dependencyChecksInitStepMaskOuts=0,...) to next coarser vertex 1. Previously assigned to (cell-centre=[0.537037,0.12963],cell-h=[0.037037,0.
 037037],local=1111,hanging=0000,select=2,has-been-refined=0000,will-be-refined=0000,is-parent-vertex-local=1111,is-parent-cell-local=1,rel-pos-within-father=[2,0])
   * ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
   *
   * Then we lift again:
   *
   * ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
 77935688595  00:01:17     rank:0       core:1       info         benchmarks::exaparty::noh::observers::StepHydroPartDrift_sweep12peano4_toolbox_particles_api_UpdateParticleGridAssociation_LiftDrop1::touchVertexLastTime(...)                 have to lift particle (debugX=[0.505,0.165],debugH=[0,0],x=[0.50497
5,0.1667],cellH=[0.111111,0.111111],searchRadius=0.0185185,ParallelState=Local,NewParallelState=Local,MoveState=Moved,CellHasUpdatedParticle=1,cfl=0.1,initialTimeStepSize=0.0001,adjustTimeStepSize=0,hydroDimensions=2,etaFactor=1.2348,smlMin=1e-06,smlMax=0.166667,smlTolerance=0.0001,smlMaxIterations=30,alp
haAV=0.8,betaAV=3,mass=2.5e-05,v=[-0.0149237,0.999888],a=[7.92017e-06,-0.000467427],density=1.00459,pressure=6.71868e-07,smoothingLength=0.00615988,u=1.0032e-06,uDot=2.00584e-06,f=-1.07003e-07,wcount_dh=-55604.5,rho_dh=-1.39011,wcount=40183.6,hDot=-0.00923729,smoothingLengthIterCount=0,balsara=0.799993,ro
t_v=8.2102e-06,div_v=-2.9864,v_sig_AV=0.0916537,soundSpeed=0.00105578,v_full=[-0.0149237,0.999888],u_full=1.0033e-06,hasNoNeighbours=0,isBoundaryParticle=0,partid=1,dependencyChecksPeanoEventUsedBySwift=touchVertexFirstTime,dependencyChecksAlgorithmStepLastUpdated=FlagBoundaryParticles,dependencyChecksAlg
orithmStepUpdates=1,...,dependencyChecksAlgorithmStepMaskOuts=0,...,dependencyChecksInitStepLastUpdated=FlagBoundaryParticles,dependencyChecksInitStepUpdates=1,...,dependencyChecksInitStepMaskOuts=0,...) to next coarser vertex 1. Previously assigned to (cell-centre=[0.611111,0.166667],cell-h=[0.111111,0.1
11111],local=1111,hanging=0000,select=0,has-been-refined=1111,will-be-refined=1111,is-parent-vertex-local=1111,is-parent-cell-local=1,rel-pos-within-father=[2,1])
   * ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
   *
   * This is wrong, as we should have assigned to the vertex 3 in the first
   * step and actually uncovers a potential bug.
   *
   */
  void testLiftDropOfParticleAssociatedWithVertex02();
 
 
  void testLiftDropOfParticleAssociatedWithVertex03();
 
 
  void testLiftDropOfParticleAssociatedWithVertex04();
 
public:
  MultiscaleTransitionsTest();
  virtual void run() override;
};