Merge branch 'fix/scaling-parameter' into 'master'

Divide the regularization parameter through the scaling parameter in the PN...

See merge request !133

dune/gfe/riemannianpnsolver.cc

@@ -281,7 +281,7 @@ void RiemannianProximalNewtonSolver<Basis,TargetSpace,Assembler>::solve()
             // Add the regularization - Identity Matrix for now
             for (std::size_t i=0; i<stiffnessMatrix.N(); i++)
               for(int j=0; j<blocksize; j++)
-                stiffnessMatrix[i][i][j][j] += regularization*scaling_[j];
+                stiffnessMatrix[i][i][j][j] += regularization/scaling_[j];
 
             innerSolver_->setProblem(stiffnessMatrix,corr_global,rhs_global);
             innerSolver_->preprocess();