e1[0] = 1.;
e2[1] = 1.;
e3[2] = 1.;
- Tensor<1,3> result = cross_product(e1, e2);
+ Tensor<1,3> result = cross_product_3d(e1, e2);
deallog << '\t' << result[0]
<< '\t' << result[1]
<< '\t' << result[2] << std::endl;
- result = cross_product(e2, e3);
+ result = cross_product_3d(e2, e3);
deallog << '\t' << result[0]
<< '\t' << result[1] << '\t'
<< result[2] << std::endl;
- result = cross_product(e3, e1);
+ result = cross_product_3d(e3, e1);
deallog << '\t' << result[0]
<< '\t' << result[1]
<< '\t' << result[2] << std::endl;
e1[0] = 1.;
e2[1] = 1.;
e3[2] = 1.;
- Tensor<1,3,std::complex<double> > result = cross_product(e1, e2);
+ Tensor<1,3,std::complex<double> > result = cross_product_3d(e1, e2);
deallog << '\t' << result[0]
<< '\t' << result[1]
<< '\t' << result[2] << std::endl;
- result = cross_product(e2, e3);
+ result = cross_product_3d(e2, e3);
deallog << '\t' << result[0]
<< '\t' << result[1] << '\t'
<< result[2] << std::endl;
- result = cross_product(e3, e1);
+ result = cross_product_3d(e3, e1);
deallog << '\t' << result[0]
<< '\t' << result[1]
<< '\t' << result[2] << std::endl;
e1[0] = 1.;
e2[1] = 1.;
e3[2] = 1.;
- Tensor<1,3,float> result = cross_product(e1, e2);
+ Tensor<1,3,float> result = cross_product_3d(e1, e2);
deallog << '\t' << static_cast<double>(result[0])
<< '\t' << static_cast<double>(result[1])
<< '\t' << static_cast<double>(result[2]) << std::endl;
- result = cross_product(e2, e3);
+ result = cross_product_3d(e2, e3);
deallog << '\t' << static_cast<double>(result[0])
<< '\t' << static_cast<double>(result[1]) << '\t'
<< static_cast<double>(result[2]) << std::endl;
- result = cross_product(e3, e1);
+ result = cross_product_3d(e3, e1);
deallog << '\t' << static_cast<double>(result[0])
<< '\t' << static_cast<double>(result[1])
<< '\t' << static_cast<double>(result[2]) << std::endl;
const Tensor<1,3> &n,
const double &rn_c)
{
- Tensor<1, 3> ra1 = cross_product(r, a1);
- Tensor<1, 3> ra2 = cross_product(r, a2);
- Tensor<1, 3> a12 = cross_product(a1, a2);
+ Tensor<1, 3> ra1 = cross_product_3d(r, a1);
+ Tensor<1, 3> ra2 = cross_product_3d(r, a2);
+ Tensor<1, 3> a12 = cross_product_3d(a1, a2);
double integral =
-1./2./numbers::PI
const Tensor<1,3> &a1,
const Tensor<1,3> &a2)
{
- Tensor<1, 3> ra1 = cross_product(r, a1);
- Tensor<1, 3> ra2 = cross_product(r, a2);
- Tensor<1, 3> a12 = cross_product(a1, a2);
+ Tensor<1, 3> ra1 = cross_product_3d(r, a1);
+ Tensor<1, 3> ra2 = cross_product_3d(r, a2);
+ Tensor<1, 3> a12 = cross_product_3d(a1, a2);
double integral = 1./2./numbers::PI
*atan2( ra1*ra2, (r.norm()* (r*a12)));