From: Wolfgang Bangerth Date: Wed, 18 Feb 2015 12:46:37 +0000 (-0600) Subject: Further fixes to accommodate new Point vs Tensor distinction. X-Git-Tag: v8.3.0-rc1~449^2 X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=refs%2Fpull%2F562%2Fhead;p=dealii.git Further fixes to accommodate new Point vs Tensor distinction. --- diff --git a/source/opencascade/boundary_lib.cc b/source/opencascade/boundary_lib.cc index c5b3a3a80b..ce4873000a 100644 --- a/source/opencascade/boundary_lib.cc +++ b/source/opencascade/boundary_lib.cc @@ -171,19 +171,23 @@ namespace OpenCASCADE { Tensor<1,3> u = surrounding_points[1]-surrounding_points[0]; Tensor<1,3> v = surrounding_points[2]-surrounding_points[0]; - Tensor<1,3> n1(u(1)*v(2)-u(2)*v(1),u(2)*v(0)-u(0)*v(2),u(0)*v(1)-u(1)*v(0)); + const double n1_coords[3] = {u(1)*v(2)-u(2)*v(1),u(2)*v(0)-u(0)*v(2),u(0)*v(1)-u(1)*v(0)}; + Tensor<1,3> n1(n1_coords); n1 = n1/n1.norm(); u = surrounding_points[2]-surrounding_points[3]; v = surrounding_points[1]-surrounding_points[3]; - Tensor<1,3> n2(u(1)*v(2)-u(2)*v(1),u(2)*v(0)-u(0)*v(2),u(0)*v(1)-u(1)*v(0)); + const double n2_coords[3] = {u(1)*v(2)-u(2)*v(1),u(2)*v(0)-u(0)*v(2),u(0)*v(1)-u(1)*v(0)}; + Tensor<1,3> n2(n2_coords); n2 = n2/n2.norm(); u = surrounding_points[4]-surrounding_points[7]; v = surrounding_points[6]-surrounding_points[7]; - Point<3> n3(u(1)*v(2)-u(2)*v(1),u(2)*v(0)-u(0)*v(2),u(0)*v(1)-u(1)*v(0)); + const double n3_coords[3] = {u(1)*v(2)-u(2)*v(1),u(2)*v(0)-u(0)*v(2),u(0)*v(1)-u(1)*v(0)}; + Tensor<1,3> n3(n3_coords); n3 = n3/n3.norm(); u = surrounding_points[6]-surrounding_points[7]; v = surrounding_points[5]-surrounding_points[7]; - Point<3> n4(u(1)*v(2)-u(2)*v(1),u(2)*v(0)-u(0)*v(2),u(0)*v(1)-u(1)*v(0)); + const double n4_coords[3] = {u(1)*v(2)-u(2)*v(1),u(2)*v(0)-u(0)*v(2),u(0)*v(1)-u(1)*v(0)}; + Tensor<1,3> n4(n4_coords); n4 = n4/n4.norm(); //for (unsigned int i=0; i