From 6cc923eafe0f879e8ae56abf475d8171ce60c9e9 Mon Sep 17 00:00:00 2001
From: Oliver Sander <sander@igpm.rwth-aachen.de>
Date: Thu, 24 Aug 2006 09:53:41 +0000
Subject: [PATCH] minor cleanup
[[Imported from SVN: r910]]
---
rod3d.cc | 18 ++++++++----------
1 file changed, 8 insertions(+), 10 deletions(-)
diff --git a/rod3d.cc b/rod3d.cc
index 85e51b62..289f1395 100644
--- a/rod3d.cc
+++ b/rod3d.cc
@@ -15,7 +15,6 @@
#include "../solver/iterativesolver.hh"
#include "../common/geomestimator.hh"
-#include "../common/refinegrid.hh"
#include "../common/energynorm.hh"
#include <dune/common/configparser.hh>
@@ -88,7 +87,7 @@ int main (int argc, char *argv[]) try
Array<SimpleVector<BoxConstraint<blocksize> > > trustRegionObstacles(1);
Array<BitField> hasObstacle(1);
- Array<BitField> dirichletNodes(1);
+ std::vector<BitField> dirichletNodes(1);
// ////////////////////////////////
// Create a multigrid solver
@@ -148,14 +147,13 @@ int main (int argc, char *argv[]) try
x[i].q = Quaternion<double>::identity();
}
- x[x.size()-1].r[0] = 0;
- x[x.size()-1].r[1] = 0.1;
- x[x.size()-1].r[2] = 0;
-#if 0
- x[x.size()-1].q[0] = 0;
- x[x.size()-1].q[1] = 0;
- x[x.size()-1].q[2] = 1/sqrt(2);
- x[x.size()-1].q[3] = 1/sqrt(2);
+// x[x.size()-1].r[0] = 0;
+// x[x.size()-1].r[1] = 0;
+// x[x.size()-1].r[2] = 0;
+#if 1
+ FieldVector<double,3> zAxis(0);
+ zAxis[2] = 1;
+ x[x.size()-1].q = Quaternion<double>(zAxis, M_PI);
#endif
std::cout << "Left boundary orientation:" << std::endl;
--
GitLab