diff --git a/src/averagedistanceassembler.hh b/src/averagedistanceassembler.hh index 18de1a9e1f7adf4ca076640fffcee2068ffc8c20..82b04e83ea52212128b1e65c64aca9ac309f52d3 100644 --- a/src/averagedistanceassembler.hh +++ b/src/averagedistanceassembler.hh @@ -57,8 +57,8 @@ public: Dune::FieldMatrix<double,size,size>& matrix) const { matrix = 0; - for (int i=0; i<dim+1; i++) - matrix.axpy(w[i], TargetSpace::secondDerivativeOfDistanceSquaredWRTSecondArgument(coefficients_[i], q)); + for (int i=0; i<coefficients_.size(); i++) + matrix.axpy(weights_[i], TargetSpace::secondDerivativeOfDistanceSquaredWRTSecondArgument(coefficients_[i], x)); } const std::vector<TargetSpace> coefficients_;