Function LaplaceDL_Triangle

bem/analytic_integrals.cpp:221–232

  double LaplaceDL_Triangle (Vec<3> v0, Vec<3> v1, Vec<3> v2, Vec<3> x,
                             Vec<3> source_normal)
  {
    Vec<3> n = Cross(v1-v0, v2-v0);
    if (InnerProduct(n, source_normal) > 0)
      {
        Swap(v1, v2);
        n *= -1;
      }
    n /= L2Norm(n);
    return LaplaceDL_TriangleOriented(v0, v1, v2, x, n);
  }

}