Skip to content
Snippets Groups Projects
Commit e0fdc851 authored by Anna's avatar Anna
Browse files

Remove unused chimera transforms

parent 89cab54a
No related branches found
No related tags found
1 merge request!122Merge changes from develop into open source
......@@ -37,62 +37,6 @@
using namespace vf::lbm::constant;
////////////////////////////////////////////////////////////////////////////////
//! \brief forward chimera transformation \ref forwardInverseChimeraWithK
//! Transformation from distributions to central moments according to Eq. (6)-(14) in \ref
//! <a href="https://doi.org/10.1016/j.jcp.2017.05.040"><b>[ M. Geier et al. (2017), DOI:10.1016/j.jcp.2017.05.040 ]</b></a>
//! Modified for lower round-off errors.
inline __device__ void forwardInverseChimeraWithKincompressible(real &mfa, real &mfb, real &mfc, real vv, real v2, real Kinverse, real K, real oneMinusRho) {
real m2 = mfa + mfc;
real m1 = mfc - mfa;
real m0 = m2 + mfb;
mfa = m0;
m0 *= Kinverse;
m0 += oneMinusRho;
mfb = (m1*Kinverse - m0 * vv) * K;
mfc = ((m2 - c2o1* m1 * vv)*Kinverse + v2 * m0) * K;
}
////////////////////////////////////////////////////////////////////////////////
//! \brief backward chimera transformation \ref backwardInverseChimeraWithK
//! Transformation from central moments to distributions according to Eq. (57)-(65) in \ref
//! <a href="https://doi.org/10.1016/j.jcp.2017.05.040"><b>[ M. Geier et al. (2017), DOI:10.1016/j.jcp.2017.05.040 ]</b></a>
//! Modified for lower round-off errors.
inline __device__ void backwardInverseChimeraWithKincompressible(real &mfa, real &mfb, real &mfc, real vv, real v2, real Kinverse, real K, real oneMinusRho) {
real m0 = (((mfc - mfb) * c1o2 + mfb * vv)*Kinverse + (mfa*Kinverse + oneMinusRho) * (v2 - vv) * c1o2) * K;
real m1 = (((mfa - mfc) - c2o1 * mfb * vv)*Kinverse + (mfa*Kinverse + oneMinusRho) * (-v2)) * K;
mfc = (((mfc + mfb) * c1o2 + mfb * vv)*Kinverse + (mfa*Kinverse + oneMinusRho) * (v2 + vv) * c1o2) * K;
mfa = m0;
mfb = m1;
}
////////////////////////////////////////////////////////////////////////////////
//! \brief forward chimera transformation \ref forwardChimera
//! - Chimera transform from distributions to central moments as defined in Eq. (43)-(45) in \ref
//! <a href="https://doi.org/10.1016/j.camwa.2015.05.001"><b>[ M. Geier et al. (2015), DOI:10.1016/j.camwa.2015.05.001 ]</b></a>
inline __device__ void forwardChimera(real &mfa, real &mfb, real &mfc, real vv, real v2) {
real m1 = (mfa + mfc) + mfb;
real m2 = mfc - mfa;
mfc = (mfc + mfa) + (v2*m1 - c2o1*vv*m2);
mfb = m2 - vv*m1;
mfa = m1;
}
////////////////////////////////////////////////////////////////////////////////
//! \brief backward chimera transformation \ref backwardChimera
//! - Chimera transform from central moments to distributions as defined in Eq. (88)-(96) in \ref
//! <a href="https://doi.org/10.1016/j.camwa.2015.05.001"><b>[ M. Geier et al. (2015), DOI:10.1016/j.camwa.2015.05.001 ]</b></a>
inline __device__ void backwardChimera(real &mfa, real &mfb, real &mfc, real vv, real v2) {
real ma = (mfc + mfa*(v2 - vv))*c1o2 + mfb*(vv - c1o2);
real mb = ((mfa - mfc) - mfa*v2) - c2o1*mfb*vv;
mfc = (mfc + mfa*(v2 + vv))*c1o2 + mfb*(vv + c1o2);
mfb = mb;
mfa = ma;
}
////////////////////////////////////////////////////////////////////////////////
inline __device__ real calcDistributionBC_AD_interpol(real q, real weight, real v, real v_sq, real f, real finf, real omegaDiffusivity, real jTangential, real concentration) {
real feq = weight * concentration * (c1o1 + c3o1 * v + c9o2 * v * v * concentration - v_sq * concentration);
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment