Merge Pull Request #9424 from rppawlo/Trilinos/panzer-no-uvm-core-tests

Automatically Merged using Trilinos Pull Request AutoTester
PR Title: Panzer: remove UVM use from IntegrationValues
PR Author: rppawlo

packages/panzer/disc-fe/src/Panzer_ConvertNormalToRotationMatrix.hpp

@@ -0,0 +1,61 @@
+#ifndef PANZER_CONVERT_NORMAL_TO_ROTATION_MATRIX_HPP
+#define PANZER_CONVERT_NORMAL_TO_ROTATION_MATRIX_HPP
+
+namespace panzer {
+
+template <typename Scalar>
+KOKKOS_INLINE_FUNCTION
+void
+convertNormalToRotationMatrix(const Scalar normal[3], Scalar transverse[3], Scalar binormal[3])
+{
+  using T = Scalar;
+
+  const T n  = sqrt(normal[0]*normal[0]+normal[1]*normal[1]+normal[2]*normal[2]);
+
+  // If this fails then the geometry for this cell is probably undefined
+  if(n > 0.){
+    // Make sure transverse is not parallel to normal within some margin of error
+    transverse[0]=0.;transverse[1]=1.;transverse[2]=0.;
+    if(Kokkos::Experimental::fabs(normal[0]*transverse[0]+normal[1]*transverse[1])>0.9){
+      transverse[0]=1.;transverse[1]=0.;
+    }
+
+    const T nt = normal[0]*transverse[0]+normal[1]*transverse[1]+normal[2]*transverse[2];
+
+    // Note normal has unit length
+    const T mult = nt/(n*n); // = nt
+
+    // Remove normal projection from transverse
+    for(int dim=0;dim<3;++dim){
+      transverse[dim] = transverse[dim] - mult * normal[dim];
+    }
+
+    const T t = sqrt(transverse[0]*transverse[0]+transverse[1]*transverse[1]+transverse[2]*transverse[2]);
+    KOKKOS_ASSERT(t != 0.);
+    for(int dim=0;dim<3;++dim){
+      transverse[dim] /= t;
+    }
+
+    // We assume a right handed system such that b = n \times t
+    binormal[0] = (normal[1] * transverse[2] - normal[2] * transverse[1]);
+    binormal[1] = (normal[2] * transverse[0] - normal[0] * transverse[2]);
+    binormal[2] = (normal[0] * transverse[1] - normal[1] * transverse[0]);
+
+    // Normalize binormal
+    const T b = sqrt(binormal[0]*binormal[0]+binormal[1]*binormal[1]+binormal[2]*binormal[2]);
+    for(int dim=0;dim<3;++dim){
+      binormal[dim] /= b;
+    }
+  } else {
+    transverse[0] = 0.;
+    transverse[1] = 0.;
+    transverse[2] = 0.;
+    binormal[0] = 0.;
+    binormal[1] = 0.;
+    binormal[2] = 0.;
+  }
+}
+
+}
+
+#endif