513 amrex::Real b1 = ( sm3 * sm3 * 6649./2880.0
514 - sm3 * sm2 * 2623./160.0
515 + sm3 * sm1 * 9449./480.0
516 - sm3 * s * 11389./1440.0
517 + sm2 * sm2 * 28547./960.0
518 - sm2 * sm1 * 35047./480.0
519 + sm2 * s * 14369./480.0
520 + sm1 * sm1 * 44747./960.0
521 - sm1 * s * 6383./160.0
522 + s * s * 25729./2880.0 );
523 amrex::Real b2 = ( sm2 * sm2 * 3169/2880.0
524 - sm2 * sm1 * 3229/480.0
525 + sm2 * s * 3169/480.0
526 - sm2 * sp1 * 2989/1440.0
527 + sm1 * sm1 * 11147/960.0
528 - sm1 * s * 11767/480.0
529 + sm1 * sp1 * 1283/160.0
530 + s * s * 13667/960.0
531 - s * sp1 * 5069/480.0
532 + sp1 * sp1 * 6649/2880.0 );
533 amrex::Real b3 = ( sm1 * sm1 * 6649./2880.0
534 - sm1 * s * 5069./480.0
535 + sm1 * sp1 * 1283./160.0
536 - sm1 * sp2 * 2989./1440.0
537 + s * s * 13667./960.0
538 - s * sp1 * 11767./480.0
539 + s * sp2 * 3169./480.0
540 + sp1 * sp1 * 11147./960.0
541 - sp1 * sp2 * 3229./480.0
542 + sp2 * sp2 * 3169./2880.0 );
543 amrex::Real b4 = ( s * s * 25729./2880.
544 - s * sp1 * 6383./160.
545 + s * sp2 * 14369./480.
546 - s * sp3 * 11389./1440.
547 + sp1 * sp1 * 44747./960.
548 - sp1 * sp2 * 35047./480.
549 + sp1 * sp3 * 9449./480.
550 + sp2 * sp2 * 28547./960.
551 - sp2 * sp3 * 2623./160.
552 + sp3 * sp3 * 6649./2880. );
555 amrex::Real t5 = std::abs(b1 - b2 - b3 + b4);
556 amrex::Real w1 =
g1 * ( 1.0 + (t5*t5) / ((
eps + b1) * (
eps + b1)) );
557 amrex::Real w2 =
g2 * ( 1.0 + (t5*t5) / ((
eps + b2) * (
eps + b2)) );
558 amrex::Real w3 =
g3 * ( 1.0 + (t5*t5) / ((
eps + b3) * (
eps + b3)) );
559 amrex::Real w4 =
g4 * ( 1.0 + (t5*t5) / ((
eps + b4) * (
eps + b4)) );
562 amrex::Real wsum = w1 + w2 + w3 + w4;
565 amrex::Real v1 = (-0.3125)*sm3 + ( 1.3125)*sm2 + (-2.1875)*sm1 + ( 2.1875)*s;
566 amrex::Real v2 = ( 0.0625)*sm2 + (-0.3125)*sm1 + ( 0.9375)*s + ( 0.3125)*sp1;
567 amrex::Real v3 = (-0.0625)*sm1 + ( 0.5625)*s + ( 0.5625)*sp1 + (-0.0625)*sp2;
568 amrex::Real v4 = ( 0.3125)*s + ( 0.9375)*sp1 + (-0.3125)*sp2 + ( 0.0625)*sp3;
571 return ( (w1 * v1 + w2 * v2 + w3 * v3 + w4 * v4) / (wsum) );
static constexpr amrex::Real g1
Definition: ERF_Interpolation_WENO_Z.H:577
const amrex::Real eps
Definition: ERF_Interpolation_WENO_Z.H:576
static constexpr amrex::Real g4
Definition: ERF_Interpolation_WENO_Z.H:580
static constexpr amrex::Real g2
Definition: ERF_Interpolation_WENO_Z.H:578
static constexpr amrex::Real g3
Definition: ERF_Interpolation_WENO_Z.H:579