← Back to branch summary

~chaffra/ffc/main-old

« back to all changes in this revision

Viewing changes to test/regression/references/r_auto/MathFunctions.h

  • Committer: Chaffra Affouda
  • Date: 2013-03-17 13:09:45 UTC
  • mfrom: (1675.1.153 work)
  • Revision ID: chaffra@gmail.com-20130317130945-6znd0m92hslnarkb
commit to trunk r1828

Show diffs side-by-side

added added

removed removed

Lines of Context:
test/regression/references/r_auto/MathFunctions.h
30
30
#include <stdexcept>
31
31
#include <fstream>
32
32
#include <ufc.h>
33
 
#include <boost/math/tr1.hpp>
 
33
#include <boost/math/special_functions.hpp>
34
34
 
35
35
/// This class defines the interface for a finite element.
36
36
 
1240
1240
      
1241
1241
      // Number of operations for primary indices: 114
1242
1242
      // Number of operations to compute entry: 114
1243
 
      A[0] += ((boost::math::tr1::cyl_neumann(0.0, F0*F1) + (boost::math::tr1::cyl_neumann(1.0, ((-1.0)*F1 + 3.0*F0)) + ((boost::math::tr1::cyl_bessel_k(0.0, F0/(F1)) + (((boost::math::tr1::cyl_bessel_k(1.0, F0) + ((boost::math::tr1::cyl_bessel_j(0.0, F0*F1) + ((boost::math::tr1::cyl_bessel_j(1.0, F0) + ((((((erf(F0*F1) + ((erf(F0) + (((std::atan(((-1.0)*F1 + 3.0*F0)) + (std::atan(F0) + ((std::asin(F0*F1) + ((std::asin(F0) + (std::acos(F0/(F1)) + ((std::acos(((-1.0)*F1 + 3.0*F0)) + (std::acos(F0) + (std::tan(F0/(F1)) + ((((std::sin(F0/(F1)) + (((std::sin(F0) + ((((((std::log(F0*F1) + (std::log(((-1.0)*F1 + 3.0*F0)) + (std::log(F0) + (((std::exp(((-1.0)*F1 + 3.0*F0)) + ((std::sqrt(F0/(F1)) + ((std::sqrt(((-1.0)*F1 + 3.0*F0)) + std::sqrt(F0)) + std::sqrt(F0*F1))) + std::exp(F0))) + std::exp(F0*F1)) + std::exp(F0/(F1)))))) + std::log(F0/(F1))) + std::cos(F0)) + std::cos(((-1.0)*F1 + 3.0*F0))) + std::cos(F0*F1)) + std::cos(F0/(F1)))) + std::sin(((-1.0)*F1 + 3.0*F0))) + std::sin(F0*F1))) + std::tan(F0)) + std::tan(((-1.0)*F1 + 3.0*F0))) + std::tan(F0*F1))))) + std::acos(F0*F1)))) + std::asin(((-1.0)*F1 + 3.0*F0)))) + std::asin(F0/(F1))))) + std::atan(F0*F1)) + std::atan(F0/(F1)))) + erf(((-1.0)*F1 + 3.0*F0)))) + erf(F0/(F1))) + boost::math::tr1::cyl_bessel_i(1.0, F0)) + boost::math::tr1::cyl_bessel_i(1.0, ((-1.0)*F1 + 3.0*F0))) + boost::math::tr1::cyl_bessel_i(0.0, F0*F1)) + boost::math::tr1::cyl_bessel_i(0.0, F0/(F1)))) + boost::math::tr1::cyl_bessel_j(1.0, ((-1.0)*F1 + 3.0*F0)))) + boost::math::tr1::cyl_bessel_j(0.0, F0/(F1)))) + boost::math::tr1::cyl_bessel_k(1.0, ((-1.0)*F1 + 3.0*F0))) + boost::math::tr1::cyl_bessel_k(0.0, F0*F1))) + boost::math::tr1::cyl_neumann(1.0, F0)))) + boost::math::tr1::cyl_neumann(0.0, F0/(F1)))*W6[ip]*det;
 
1243
      A[0] += ((boost::math::cyl_neumann(0.0, F0*F1) + (boost::math::cyl_neumann(1.0, ((-1.0)*F1 + 3.0*F0)) + (boost::math::cyl_neumann(1.0, F0) + (((boost::math::cyl_bessel_k(1.0, ((-1.0)*F1 + 3.0*F0)) + (boost::math::cyl_bessel_k(1.0, F0) + ((boost::math::cyl_bessel_j(0.0, F0*F1) + (boost::math::cyl_bessel_j(1.0, ((-1.0)*F1 + 3.0*F0)) + (((((((erf(F0*F1) + ((erf(F0) + (((std::atan(((-1.0)*F1 + 3.0*F0)) + (std::atan(F0) + ((std::asin(F0*F1) + ((std::asin(F0) + (std::acos(F0/(F1)) + ((std::acos(((-1.0)*F1 + 3.0*F0)) + (std::acos(F0) + (std::tan(F0/(F1)) + ((((std::sin(F0/(F1)) + (((std::sin(F0) + ((((((std::log(F0*F1) + (std::log(((-1.0)*F1 + 3.0*F0)) + (std::log(F0) + (((std::exp(((-1.0)*F1 + 3.0*F0)) + ((std::sqrt(F0/(F1)) + ((std::sqrt(((-1.0)*F1 + 3.0*F0)) + std::sqrt(F0)) + std::sqrt(F0*F1))) + std::exp(F0))) + std::exp(F0*F1)) + std::exp(F0/(F1)))))) + std::log(F0/(F1))) + std::cos(F0)) + std::cos(((-1.0)*F1 + 3.0*F0))) + std::cos(F0*F1)) + std::cos(F0/(F1)))) + std::sin(((-1.0)*F1 + 3.0*F0))) + std::sin(F0*F1))) + std::tan(F0)) + std::tan(((-1.0)*F1 + 3.0*F0))) + std::tan(F0*F1))))) + std::acos(F0*F1)))) + std::asin(((-1.0)*F1 + 3.0*F0)))) + std::asin(F0/(F1))))) + std::atan(F0*F1)) + std::atan(F0/(F1)))) + erf(((-1.0)*F1 + 3.0*F0)))) + erf(F0/(F1))) + boost::math::cyl_bessel_i(1.0, F0)) + boost::math::cyl_bessel_i(1.0, ((-1.0)*F1 + 3.0*F0))) + boost::math::cyl_bessel_i(0.0, F0*F1)) + boost::math::cyl_bessel_i(0.0, F0/(F1))) + boost::math::cyl_bessel_j(1.0, F0)))) + boost::math::cyl_bessel_j(0.0, F0/(F1))))) + boost::math::cyl_bessel_k(0.0, F0*F1)) + boost::math::cyl_bessel_k(0.0, F0/(F1)))))) + boost::math::cyl_neumann(0.0, F0/(F1)))*W6[ip]*det;
1244
1244
    }// end loop over 'ip'
1245
1245
  }
1246
1246
 
1292
1292
  /// Return a string identifying the form
1293
1293
  virtual const char* signature() const
1294
1294
  {
1295
 
    return "08dbb1086ad455f1da2544c74baf7e535bc15603f3c6251234660e834f2bf4eedaf41987922e99275d699deb744f3665c650cd93431a0536db58f084adde372d";
 
1295
    return "4c61cfefd887576046ad0c94450a06bc13d537048640a066d2532627ced60c1cee1d0922fa36b76d636c4e3783b4ade6c63afe3ce9ce9cb12157c9624c1e2bfe";
1296
1296
  }
1297
1297
 
1298
1298
  /// Return the rank of the global tensor (r)
1416
1416
  /// Create a new point integral on sub domain i
1417
1417
  virtual ufc::point_integral* create_point_integral(std::size_t i) const
1418
1418
  {
1419
 
    std::cerr << "*** FFC warning: " << "create_point_integral not yet implemented." << std::endl;
1420
1419
    return 0;
1421
1420
  }
1422
1421
 
1441
1440
  /// Create a new point integral on everywhere else
1442
1441
  virtual ufc::point_integral* create_default_point_integral() const
1443
1442
  {
1444
 
    std::cerr << "*** FFC warning: " << "create_default_point_integral not yet implemented." << std::endl;
1445
1443
    return 0;
1446
1444
  }
1447
1445