From 161930f9a5db246d7ce25ddcb59ec937ae95c462 Mon Sep 17 00:00:00 2001
From: Oliver Sander <sander@igpm.rwth-aachen.de>
Date: Fri, 6 Jun 2008 15:10:36 +0000
Subject: [PATCH] remove lapackcpp code
[[Imported from SVN: r2377]]
---
dirneucoupling.cc | 24 ------------------------
1 file changed, 24 deletions(-)
diff --git a/dirneucoupling.cc b/dirneucoupling.cc
index 62e61ac0..2630adad 100644
--- a/dirneucoupling.cc
+++ b/dirneucoupling.cc
@@ -399,7 +399,6 @@ int main (int argc, char *argv[]) try
std::cout << "resultant force: " << resultantForce << std::endl;
std::cout << "resultant torque: " << resultantTorque << std::endl;
-#if 1
VectorType neumannValues(rhs3d.size());
// Using that index 0 is always the left boundary for a uniformly refined OneDGrid
@@ -413,29 +412,6 @@ int main (int argc, char *argv[]) try
neumannValues,
rhs3d);
-#else
-#ifndef HAVE_LAPACKPP
-#error You need LaPack++ for this!
-#endif
- // For the time being the Neumann data coming from the rod is a dg function (== not continuous)
- // Maybe that is not necessary
- DGIndexSet<GridType> dgIndexSet(grid,grid.maxLevel());
- dgIndexSet.setup(grid,grid.maxLevel());
-
- VectorType neumannValues(dgIndexSet.size());
-
- // Using that index 0 is always the left boundary for a uniformly refined OneDGrid
- computeAveragePressure<GridType>(resultantForce, resultantTorque,
- interfaceBoundary[grid.maxLevel()],
- rodX[0],
- neumannValues);
-
- rhs3d = 0;
- assembleAndAddNeumannTerm<GridType, VectorType>(interfaceBoundary[grid.maxLevel()],
- dgIndexSet,
- neumannValues,
- rhs3d);
-#endif
// ///////////////////////////////////////////////////////////
// Solve the Neumann problem for the 3d body
// ///////////////////////////////////////////////////////////
--
GitLab