hydro_iact.h File Reference

#include "adaptive_softening_iact.h"
#include "adiabatic_index.h"
#include "hydro_parameters.h"
#include "signal_velocity.h"

Include dependency graph for hydro_iact.h:

This graph shows which files directly or indirectly include this file:

Go to the source code of this file.

Functions
	__attribute__ ((always_inline)) INLINE static void runner_iact_density(const float r2
	Density interaction between two particles.
	kernel_deval (ui, &wi, &wi_dx)
	adaptive_softening_add_correction_term (pi, ui, hi_inv, mj)
	kernel_deval (uj, &wj, &wj_dx)
	adaptive_softening_add_correction_term (pj, uj, hj_inv, mi)
	adaptive_softening_add_correction_term (pi, ui, h_inv, mj)

Variables
const float	dx [3]
const float const float	hi
const float const float const float	hj
const float const float const float struct part *restrict	pi
const float const float const float struct part *restrict struct part *restrict	pj
const float const float const float struct part *restrict struct part *restrict const float	a
const float const float const float struct part *restrict struct part *restrict const float const float	H
const float	r = sqrtf(r2)
const float	r_inv = r ? 1.0f / r : 0.0f
const float	mi = pi->mass
const float	mj = pj->mass
const float	hi_inv = 1.f / hi
const float	ui = r * hi_inv
pi	rho = mj * wi
pi density	rho_dh = mj * (hydro_dimension * wi + ui * wi_dx)
pi density	wcount = wi
pi density	wcount_dh = (hydro_dimension * wi + ui * wi_dx)
const float	hj_inv = 1.f / hj
const float	uj = r * hj_inv
float	dv [3] = pi->v[0] - pj->v[0]
float	curlvr [3] = dv[1] * dx[2] - dv[2] * dx[1]
const float	faci = mj * wi_dx * r_inv
const float	facj = mi * wj_dx * r_inv
const float	dvdr = dv[0] * dx[0] + dv[1] * dx[1] + dv[2] * dx[2]
pi density	div_v = faci * dvdr
pi density	rot_v [0] = faci * curlvr[0]
const float	h_inv = 1.f / hi

Function Documentation

__attribute__()

__attribute__

(

(always_inline)

)

const

Density interaction between two particles.

Calculate the gradient interaction between particle i and particle j.

Density interaction between two particles (non-symmetric).

Parameters

r2	Comoving square distance between the two particles.
dx	Comoving vector separating both particles (pi - pj).
hi	Comoving smoothing-length of particle i.
hj	Comoving smoothing-length of particle j.
pi	First particle.
pj	Second particle.
a	Current scale factor.
H	Current Hubble parameter.
r2	Comoving square distance between the two particles.
dx	Comoving vector separating both particles (pi - pj).
hi	Comoving smoothing-length of particle i.
hj	Comoving smoothing-length of particle j.
pi	First particle.
pj	Second particle (not updated).
a	Current scale factor.
H	Current Hubble parameter.

Nothing to do here in this scheme.

Parameters

r2	Comoving squared distance between particle i and particle j.
dx	Comoving distance vector between the particles (dx = pi->x - pj->x).
hi	Comoving smoothing-length of particle i.
hj	Comoving smoothing-length of particle j.
pi	Particle i.
pj	Particle j.
a	Current scale factor.
H	Current Hubble parameter.

References INLINE.

adaptive_softening_add_correction_term() [1/3]

adaptive_softening_add_correction_term	(	pi	,
		ui	,
		h_inv	,
		mj	)

References curlvr, dv, dvdr, dx, faci, h_inv, mj, pi, pj, r_inv, and ui.

adaptive_softening_add_correction_term() [2/3]

adaptive_softening_add_correction_term	(	pi	,
		ui	,
		hi_inv	,
		mj	)

References hi_inv, mj, pi, and ui.

adaptive_softening_add_correction_term() [3/3]

adaptive_softening_add_correction_term	(	pj	,
		uj	,
		hj_inv	,
		mi	)

References hj_inv, mi, pj, and uj.

kernel_deval() [1/2]

kernel_deval	(	ui	,
		&	wi,
		&	wi_dx )

References ui.

Referenced by swift2::kernels::legacy::density_kernel(), and hydro_kernel_dump().

Here is the caller graph for this function:

kernel_deval() [2/2]

kernel_deval	(	uj	,
		&	wj,
		&	wj_dx )

References mi, pj, and uj.

Variable Documentation

a

const float const float const float struct part *restrict struct part *restrict const float a

Definition at line 55 of file hydro_iact.h.

Referenced by applications::exahype2::euler::sphericalaccretion::addGravitationalSource_AlphaCDM(), tl::compose_fn< F, G >::call(), TP::Utilities::chebev(), checkUserDefinedType(), tl::compose_fn< F, G >::compose_fn(), do_density(), do_density(), TP::Utilities::fourft(), applications::exahype2::euler::sphericalaccretion::SSInfall::getForceDensityNorm(), TP::TwoPunctures::interpol(), toolbox::blockstructured::interpolateCellDataAssociatedToVolumesIntoOverlappingCell_secondOrder(), applications::exahype2::euler::sphericalaccretion::SSInfall::mass_interpolate(), tarch::la::max(), applications::exahype2::CompressibleNavierStokes::NavierStokesSolver::nonconservativeProduct(), TP::Utilities::nrerror(), applications::exahype2::ccz4::Psi4Calc(), TP::TwoPunctures::PunctTaylorExpandAtArbitPosition(), toolbox::finiteelements::stencilProduct(), toolbox::finiteelements::stencilProduct(), toolbox::finiteelements::stencilProduct(), toolbox::finiteelements::stencilProduct(), toolbox::finiteelements::stencilProduct(), toolbox::finiteelements::stencilProduct(), toolbox::finiteelements::stencilProduct(), toolbox::finiteelements::stencilProduct(), toolbox::finiteelements::stencilProduct(), toolbox::finiteelements::stencilProduct(), toolbox::curvi::Face::subface(), tarch::la::tests::ScalarTest::testAbs(), tarch::la::tests::ScalarTest::testComparison(), toolbox::particles::assignmentchecks::tests::TestHelpers::testTruthTableSearchAndIDKeys(), TP::TwoPunctures::TestVector_w(), and tarch::la::tests::VectorTest::testVectorVectorOperations().

curlvr

float curlvr = dv[1] * dx[2] - dv[2] * dx[1]

Definition at line 98 of file hydro_iact.h.

Referenced by adaptive_softening_add_correction_term(), and swift2::kernels::legacy::density_kernel().

div_v

pi density div_v = faci * dvdr

Definition at line 108 of file hydro_iact.h.

dv

float dv = pi->v[0] - pj->v[0]

Definition at line 98 of file hydro_iact.h.

Referenced by adaptive_softening_add_correction_term(), TP::TwoPunctures::bicgstab(), swift2::kernels::legacy::density_kernel(), toolbox::particles::assignmentchecks::tests::internal::dropParticle(), TP::TwoPunctures::J_times_dv(), TP::TwoPunctures::JFD_times_dv(), toolbox::particles::assignmentchecks::tests::internal::liftParticle(), TP::TwoPunctures::LineRelax_al(), TP::TwoPunctures::LineRelax_be(), TP::TwoPunctures::Newton(), TP::TwoPunctures::relax(), TP::TwoPunctures::resid(), TP::TwoPunctures::SetMatrix_JFD(), TP::TwoPunctures::TestRelax(), and TP::TwoPunctures::TestVector_w().

dvdr

const float dvdr = dv[0] * dx[0] + dv[1] * dx[1] + dv[2] * dx[2]

Definition at line 106 of file hydro_iact.h.

Referenced by adaptive_softening_add_correction_term(), and swift2::kernels::legacy::density_kernel().

dx

const float dx

Definition at line 54 of file hydro_iact.h.

Referenced by adaptive_softening_add_correction_term(), ExaSeis::Derivatives< Shortcuts, num_nodes >::computeDerivatives_x_3D(), toolbox::curvi::Transformation::computeDerivatives_x_3D(), swift2::kernels::legacy::density_kernel(), swift2::kernels::legacy::densityKernelEvaluatePairInteractionPredicate(), do_density(), do_density(), Refinement::CoarseBoundaryLayer< Shortcuts >::eval(), Refinement::RefineBetweenPositions< Shortcuts >::eval(), Refinement::RefineCubeAroundPosition< Shortcuts >::eval(), Refinement::RefineDownToPosition< Shortcuts >::eval(), Refinement::RefineDownToPositionCustomCoordinates< Shortcuts >::eval(), Refinement::RefineFilterCube< Shortcuts >::eval(), Refinement::RefinementCriterion< Shortcuts >::eval(), Refinement::RefinePosition< Shortcuts >::eval(), Refinement::RefinePositionCustomCoordinates< Shortcuts >::eval(), Refinement::StaticAMR< Shortcuts >::eval(), Refinement::trackVelocity< Shortcuts >::eval(), swift2::kernels::forceKernelDistanceCheck(), ContextCurvilinear< Shortcuts, basisSize >::getElementSize(), toolbox::curvi::Block::getMetricDerivatives(), toolbox::curvi::Index::getRange(), Context< Shortcuts, basisSize >::initUnknownsPatch(), ContextCartesian< Shortcuts, basisSize >::initUnknownsPatch(), ContextCurvilinear< Shortcuts, basisSize >::initUnknownsPatch(), ContextDiffuse< Shortcuts, basisSize >::initUnknownsPatch(), toolbox::curvi::Interpolation::Lagrange::interpolate3D(), main(), ExaSeis::Derivatives< Shortcuts, num_nodes >::metricDerivatives(), toolbox::curvi::Transformation::metricDerivativesAndJacobian3D(), applications::exahype2::ccz4::internal::recomputeAuxiliaryVariablesFD4_centralDifferences_LoopBody(), applications::exahype2::ccz4::internal::recomputeAuxiliaryVariablesFD4_leftDifferences_LoopBody(), applications::exahype2::ccz4::internal::recomputeAuxiliaryVariablesFD4_rightDifferences_LoopBody(), toolbox::particles::assignmentchecks::tests::TestHelpers::testAddingParticleMovingEvents(), toolbox::particles::assignmentchecks::tests::TestHelpers::testParticleWalk(), toolbox::particles::assignmentchecks::tests::TestHelpers::testParticleWalkSameTreeLevel(), toolbox::particles::assignmentchecks::tests::TestHelpers::testPeriodicBoundaryConditions(), toolbox::particles::assignmentchecks::tests::TestHelpers::testPeriodicBoundaryConditionsCorner(), and toolbox::particles::assignmentchecks::tests::TestHelpers::testSieveSet().

faci

const float faci = mj * wi_dx * r_inv

Definition at line 100 of file hydro_iact.h.

Referenced by adaptive_softening_add_correction_term().

facj

const float facj = mi * wj_dx * r_inv

Definition at line 101 of file hydro_iact.h.

H

const float const float const float struct part *restrict const struct part *restrict const float const float H

Initial value:

{
 
  float wi, wj, wi_dx, wj_dx

Definition at line 56 of file hydro_iact.h.

Referenced by applications::exahype2::ccz4::diagonal_gaugeWave(), do_density(), and do_density().

h_inv

const float h_inv = 1.f / hi

Definition at line 158 of file hydro_iact.h.

Referenced by adaptive_softening_add_correction_term(), swift2::kernels::legacy::hydro_end_density(), swift2::kernels::legacy::hydro_predict_extra(), and swift2::kernels::legacy::hydro_prepare_force().

hi

const float const float hi

Definition at line 54 of file hydro_iact.h.

Referenced by swift2::kernels::legacy::density_kernel(), swift2::kernels::legacy::densityKernelEvaluatePairInteractionPredicate(), do_density(), and do_density().

hi_inv

const float hi_inv = 1.f / hi

Definition at line 76 of file hydro_iact.h.

hj

const float const float const float hj

Definition at line 54 of file hydro_iact.h.

Referenced by do_density(), and do_density().

hj_inv

const float hj_inv = 1.f / hj

Definition at line 87 of file hydro_iact.h.

Referenced by adaptive_softening_add_correction_term().

mi

const float mi = pi->mass

Definition at line 72 of file hydro_iact.h.

Referenced by adaptive_softening_add_correction_term(), and kernel_deval().

mj

const float mj = pj->mass

Definition at line 73 of file hydro_iact.h.

Referenced by adaptive_softening_add_correction_term(), adaptive_softening_add_correction_term(), and swift2::kernels::legacy::density_kernel().

pi

const float const float const float struct part* restrict pi

Definition at line 55 of file hydro_iact.h.

pj

const float const float const float struct part *restrict struct part *restrict pj

Definition at line 55 of file hydro_iact.h.

Referenced by adaptive_softening_add_correction_term(), adaptive_softening_add_correction_term(), do_density(), do_density(), and kernel_deval().

r

const float r = sqrtf(r2)

Definition at line 68 of file hydro_iact.h.

Referenced by TP::TwoPunctures::bicgstab(), tbb::detail::d1::blocked_nd_range_impl(), tbb::detail::d1::blocked_nd_range_impl(), TP::TwoPunctures::C_To_c(), DiffuseInterface< Shortcuts, basisSize >::computeAlphaPatch(), swift2::kernels::gravity::computeGravitationalForce(), Region< dim, fp, topBoundaryInclusive >::contains(), swift2::kernels::legacy::density_kernel(), TP::TwoPunctures::Derivatives_AB3(), tbb::detail::d1::do_split(), TP::TwoPunctures::F_of_v(), DiffuseInterface< Shortcuts, basisSize >::getAlpha(), toolbox::curvi::Interpolation::Transfinite::interpolate2D(), Region< dim, fp, topBoundaryInclusive >::intersection(), TP::TwoPunctures::J_times_dv(), TP::TwoPunctures::JFD_times_dv(), TP::TwoPunctures::LinEquations(), tarch::la::matVec(), tarch::la::multiply(), TP::TwoPunctures::NonLinEquations(), Region< dim, fp, topBoundaryInclusive >::operator!=(), tl::views::detail::stride_fn_base::operator()(), Region< dim, fp, topBoundaryInclusive >::operator==(), Region< dim, fp, topBoundaryInclusive >::overlaps(), Numerics::riemannSolver(), Numerics::riemannSolver_BC0(), Numerics::riemannSolver_BCn(), Numerics::riemannSolver_boundary(), Numerics::riemannSolver_boundary(), TP::TwoPunctures::rx3_To_xyz(), TP::TwoPunctures::set_initial_guess(), Region< dim, fp, topBoundaryInclusive >::sharesBoundary(), and TP::TwoPunctures::TestVector_w().

r_inv

const float r_inv = r ? 1.0f / r : 0.0f

Definition at line 69 of file hydro_iact.h.

Referenced by adaptive_softening_add_correction_term(), swift2::kernels::gravity::computeGravitationalForce(), swift2::kernels::legacy::density_kernel(), and TP::TwoPunctures::rx3_To_xyz().

rho

pi rho = mj * wi

Definition at line 80 of file hydro_iact.h.

Referenced by applications::exahype2::euler::sphericalaccretion::SSInfall::add_mass(), TP::TwoPunctures::bicgstab(), Numerics::compute_parameters(), Numerics::compute_parameters(), EasiScen< Shortcuts, basisSize >::initUnknownsPointwise(), GaussianWave< Variables, basisSize >::initUnknownsPointwise(), WholeSpaceProblem< VariableShortcuts, basisSize >::initUnknownsPointwise(), Zugspitze< Shortcuts, basisSize >::initUnknownsPointwise(), Zugspitze_noTopo< Shortcuts, basisSize >::initUnknownsPointwise(), main(), Numerics::right_eigenvectors(), Numerics::right_eigenvectors_inverse(), and applications::exahype2::euler::selfgravitation::Poisson::sourceTerm().

rho_dh

pi density rho_dh = mj * (hydro_dimension * wi + ui * wi_dx)

Definition at line 81 of file hydro_iact.h.

rot_v

pi density rot_v = faci * curlvr[0]

Definition at line 116 of file hydro_iact.h.

ui

const float ui = r * hi_inv

Definition at line 77 of file hydro_iact.h.

uj

const float uj = r * hj_inv

Definition at line 88 of file hydro_iact.h.

Referenced by adaptive_softening_add_correction_term(), and kernel_deval().

wcount

pi density wcount = wi

Definition at line 82 of file hydro_iact.h.

Referenced by swift2::kernels::legacy::hydro_update_smoothing_length_and_rerun_if_required().

wcount_dh

pi density wcount_dh = (hydro_dimension * wi + ui * wi_dx)

Definition at line 83 of file hydro_iact.h.