1 #ifndef DUNE_GRIDGLUE_COMMON_CROSSPRODUCT_HH 2 #define DUNE_GRIDGLUE_COMMON_CROSSPRODUCT_HH 1 12 template <
class T,
int dim>
13 static Dune::FieldVector<T,dim>
crossProduct(
const Dune::FieldVector<T,dim>& a,
14 const Dune::FieldVector<T,dim>& b)
17 DUNE_THROW(Dune::NotImplemented,
"crossProduct does not work for dimension " << dim);
19 Dune::FieldVector<T,dim> c;
20 c[0] = a[1]*b[2] - a[2]*b[1];
21 c[1] = a[2]*b[0] - a[0]*b[2];
22 c[2] = a[0]*b[1] - a[1]*b[0];
Definition: gridglue.hh:33
static Dune::FieldVector< T, dim > crossProduct(const Dune::FieldVector< T, dim > &a, const Dune::FieldVector< T, dim > &b)
compute cross product
Definition: crossproduct.hh:13