From 9f4dfb67b572de5154a14c4d7a011cf688a379dd Mon Sep 17 00:00:00 2001
From: Oliver Sander <sander@igpm.rwth-aachen.de>
Date: Sat, 22 Mar 2014 12:56:14 +0000
Subject: [PATCH] Set the initial iterate by interpolating an analytical
function
[[Imported from SVN: r9657]]
---
cosserat-continuum.cc | 25 +++++++++++++++++++++----
1 file changed, 21 insertions(+), 4 deletions(-)
diff --git a/cosserat-continuum.cc b/cosserat-continuum.cc
index 42510b2e..3ce7e430 100644
--- a/cosserat-continuum.cc
+++ b/cosserat-continuum.cc
@@ -23,6 +23,7 @@
#include <dune/fufem/boundarypatch.hh>
#include <dune/fufem/functiontools/boundarydofs.hh>
+#include <dune/fufem/functiontools/basisinterpolator.hh>
#include <dune/fufem/functionspacebases/p1nodalbasis.hh>
#include <dune/solvers/solvers/iterativesolver.hh>
@@ -46,6 +47,19 @@ const int blocksize = TargetSpace::TangentVector::dimension;
using namespace Dune;
+class Identity
+: public Dune::VirtualFunction<FieldVector<double,dim>, FieldVector<double,3>>
+{
+public:
+ void evaluate(const FieldVector<double,dim>& x, FieldVector<double,3>& y) const
+ {
+ y = 0;
+ for (int i=0; i<dim; i++)
+ y[i] = x[i];
+ }
+};
+
+
#if 1
// Dirichlet boundary data for the shear/wrinkling example
void dirichletValues(const FieldVector<double,dim>& in, FieldVector<double,3>& out,
@@ -228,15 +242,18 @@ int main (int argc, char *argv[]) try
SolutionType x(feBasis.size());
+ Identity identity;
+ std::vector<FieldVector<double,3> > v;
+ Functions::interpolate(feBasis, v, identity);
+
+ for (size_t i=0; i<x.size(); i++)
+ x[i].r = v[i];
+
vIt = gridView.begin<dim>();
for (; vIt!=vEndIt; ++vIt) {
int idx = indexSet.index(*vIt);
- x[idx].r = 0;
- for (int i=0; i<dim; i++)
- x[idx].r[i] = vIt->geometry().corner(0)[i];
-
x[idx].r[2] = 0.002*std::cos(1e4*vIt->geometry().corner(0)[0]);
// x[idx].q is the identity, set by the default constructor
--
GitLab