From 59b577fd9c35fdc5c2858bb564995b58319c2de8 Mon Sep 17 00:00:00 2001
From: Oliver Sander <sander@igpm.rwth-aachen.de>
Date: Thu, 11 Feb 2010 21:37:53 +0000
Subject: [PATCH] implement matrix assembly by simply returning the identity
matrix. This should a sort of projected gradient descent. It is too late at
night for me to compute the actual matrix
[[Imported from SVN: r5545]]
---
src/averagedistanceassembler.hh | 4 +++-
1 file changed, 3 insertions(+), 1 deletion(-)
diff --git a/src/averagedistanceassembler.hh b/src/averagedistanceassembler.hh
index 5c54f707..17f9c3ee 100644
--- a/src/averagedistanceassembler.hh
+++ b/src/averagedistanceassembler.hh
@@ -41,7 +41,9 @@ public:
void assembleMatrix(const TargetSpace& x,
Dune::FieldMatrix<double,size,size>& matrix) const
{
- DUNE_THROW(Dune::NotImplemented, "assembleMatrix");
+ for (int i=0; i<size; i++)
+ for (int j=0; j<size; j++)
+ matrix[i][j] = (i==j);
}
const std::vector<TargetSpace> coefficients_;
--
GitLab