WdgG2Gradient.h

Go to the documentation of this file.

#define size_tri ((order + 1) * (order + 2) / 2)
#define size_lin (order + 1)
 
WORKSPACE(RealT, term_tri, size_tri);
WORKSPACE(RealT, dterm_dr_tri, size_tri);
WORKSPACE(RealT, dterm_ds_tri, size_tri);
 
WORKSPACE(RealT, term_lin, size_lin);
WORKSPACE(RealT, dterm_dt_lin, size_lin);
 
// clang-format off
term_tri[0] = 7.6156408160581188226956328435335 * ss * rr - 5.4599331074840975475126469973475 * rr -
  5.459933107484095771155807597097 * ss + 3.7248334213839786599464787286706 * ss * ss +
  3.7248334213839773276788491784828 * rr * rr + 1.8736592735119068642291040305281;
term_tri[1] = 0.16597397329018215295093341410393 * ss - 1.989733735284602289539179764688 * rr -
  0.16597397329016608247265196496301 * ss * rr - 0.16597397329016555511671526801365 * ss * ss +
  3.7248334213839764395004294783575 * rr * rr + 0.13855958741198567962804588660219;
term_tri[2] = 0.1659739732901820974397821828461 * rr - 1.9897337352846020674945748396567 * ss -
  0.16597397329016638778398373688105 * ss * rr + 3.7248334213839768835896393284202 * ss * ss -
  0.16597397329016555511671526801365 * rr * rr + 0.13855958741198570738362150223111;
term_tri[3] = 7.9963036868787469302333192899823 * ss + 0.40859490520043550976225787962903 * rr -
  7.9963036868779449051203300768975 * ss * rr - 7.9963036868779484578340088773984 * ss * ss +
  0.3563054087000795400186348160787 * rr * rr - 0.63855958741216845009347480299766;
term_tri[4] = 0.40859490520043528771765295459772 * ss + 7.9963036868787469302333192899823 * rr -
  7.9963036868779449051203300768975 * ss * rr + 0.35630540870007948450748358482087 * ss * ss -
  7.9963036868779475696555891772732 * rr * rr - 0.63855958741216845009347480299766;
term_tri[5] = 8.7089145042781037631129947840236 * ss * rr - 1.1212057226006655330508010592894 * rr -
  1.1212057226006653110061961342581 * ss + 0.35630540870007937348518112230522 * ss * ss +
  0.35630540870007937348518112230522 * rr * rr + 0.12634072648845873221290503352066;
 
dterm_dr_tri[0] = 7.6156408160581188226956328435335 * ss - 5.4599331074840975475126469973475 +
  2. * 3.7248334213839773276788491784828 * rr;
dterm_dr_tri[1] = - 1.989733735284602289539179764688 -
  0.16597397329016608247265196496301 * ss +
  2. * 3.7248334213839764395004294783575 * rr;
dterm_dr_tri[2] = 0.1659739732901820974397821828461 -
  0.16597397329016638778398373688105 * ss -
  2. * 0.16597397329016555511671526801365 * rr;
dterm_dr_tri[3] = 0.40859490520043550976225787962903 -
  7.9963036868779449051203300768975 * ss +
  2. * 0.3563054087000795400186348160787 * rr;
dterm_dr_tri[4] = 7.9963036868787469302333192899823 -
  7.9963036868779449051203300768975 * ss -
  2. * 7.9963036868779475696555891772732 * rr;
dterm_dr_tri[5] = 8.7089145042781037631129947840236 * ss - 1.1212057226006655330508010592894 +
  2. * 0.35630540870007937348518112230522 * rr;
 
dterm_ds_tri[0] = 7.6156408160581188226956328435335 * rr -
  5.459933107484095771155807597097 + 2. * 3.7248334213839786599464787286706 * ss;
dterm_ds_tri[1] = 0.16597397329018215295093341410393 -
  0.16597397329016608247265196496301 * rr - 2. * 0.16597397329016555511671526801365 * ss;
dterm_ds_tri[2] = - 1.9897337352846020674945748396567 -
  0.16597397329016638778398373688105 * rr + 2. * 3.7248334213839768835896393284202 * ss;
dterm_ds_tri[3] = 7.9963036868787469302333192899823 -
  7.9963036868779449051203300768975 * rr - 2. * 7.9963036868779484578340088773984 * ss;
dterm_ds_tri[4] = 0.40859490520043528771765295459772 -
  7.9963036868779449051203300768975 * rr + 2. * 0.35630540870007948450748358482087 * ss;
dterm_ds_tri[5] = 8.7089145042781037631129947840236 * rr -
  1.1212057226006653110061961342581 + 2. * 0.35630540870007937348518112230522 * ss;
 
term_lin[0] = 0.83333333333333333333333333333333 * tt * (tt - 0.77459666924148337703585307995648);
term_lin[1] = 1.0 - 1.6666666666666666666666666666667 * tt * tt;
term_lin[2] = 0.83333333333333333333333333333333 * tt * (tt + 0.77459666924148337703585307995648);
 
dterm_dt_lin[0] =
  2. * 0.83333333333333333333333333333333 * tt - 0.77459666924148337703585307995648 * 0.83333333333333333333333333333333;
dterm_dt_lin[1] =
 -2. * 1.6666666666666666666666666666667 * tt;
dterm_dt_lin[2] =
  2. * 0.83333333333333333333333333333333 * tt + 0.77459666924148337703585307995648 * 0.83333333333333333333333333333333;
 
for (int ii = 0; ii < size_lin; ++ii)
{
  for (int jj = 0; jj < size_tri; ++jj)
  {
    int i_sp = ii * size_tri + jj;
    // ∂/∂r, ∂/∂s, ∂/∂t:
    basisGradient[3 * i_sp    ] = dterm_dr_tri[jj] * term_lin[ii]                  /* zero */;
    basisGradient[3 * i_sp + 1] = dterm_ds_tri[jj] * term_lin[ii]                  /* zero */;
    basisGradient[3 * i_sp + 2] = /* zero */                        term_tri[jj] * dterm_dt_lin[ii];
  }
}
// clang-format on