Skip to content
Snippets Groups Projects
Commit 96008d7c authored by Dr. Felix Tobias Schindler's avatar Dr. Felix Tobias Schindler
Browse files

[local.integrand.abs] enable vector-valued case

parent 83225404
No related branches found
No related tags found
No related merge requests found
......@@ -21,6 +21,7 @@ namespace GDT {
template <class E, size_t r = 1, size_t rC = 1, class R = double, class F = double>
class LocalElementAbsIntegrand : public LocalUnaryElementIntegrandInterface<E, r, rC, R, F>
{
static_assert(rC == 1, "");
using ThisType = LocalElementAbsIntegrand<E, r, rC, R, F>;
using BaseType = LocalUnaryElementIntegrandInterface<E, r, rC, R, F>;
......@@ -52,9 +53,8 @@ public:
// evaluate
basis.evaluate(point_in_reference_element, basis_values_, param);
// compute integrand
using std::abs;
for (size_t ii = 0; ii < size; ++ii)
result[ii] = abs(basis_values_[ii]);
result[ii] = basis_values_[ii].two_norm();
} // ... evaluate(...)
private:
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment