Skip to content
Snippets Groups Projects
Commit ecaab53c authored by graeser's avatar graeser
Browse files

Change entries from Field*<double,1,1> to plain double

parent 6aba9157
No related branches found
No related tags found
No related merge requests found
...@@ -87,8 +87,8 @@ int main(int argc, char** argv) ...@@ -87,8 +87,8 @@ int main(int argc, char** argv)
using GridView = decltype(gridView); using GridView = decltype(gridView);
using Matrix = Dune::BCRSMatrix<Dune::FieldMatrix<double,1,1>>; using Matrix = Dune::BCRSMatrix<double>;
using Vector = Dune::BlockVector<Dune::FieldVector<double,1>>; using Vector = Dune::BlockVector<double>;
Matrix A; Matrix A;
Vector rhs; Vector rhs;
......
...@@ -132,17 +132,17 @@ void assembleLocalRHS( ...@@ -132,17 +132,17 @@ void assembleLocalRHS(
template<class GridView, class RHSFunction> template<class GridView, class RHSFunction>
void assemblePoissonProblemPQ1( void assemblePoissonProblemPQ1(
const GridView& gridView, const GridView& gridView,
Dune::BCRSMatrix<Dune::FieldMatrix<double,1,1>>& matrix, Dune::BCRSMatrix<double>& matrix,
Dune::BlockVector<Dune::FieldVector<double,1>>& rhs, Dune::BlockVector<double>& rhs,
const RHSFunction& rhsFunction) const RHSFunction& rhsFunction)
{ {
static const int dim = GridView::dimension; static const int dim = GridView::dimension;
const auto& indexSet = gridView.indexSet(); const auto& indexSet = gridView.indexSet();
std::size_t size = indexSet.size(dim); std::size_t size = indexSet.size(dim);
using Matrix = Dune::BCRSMatrix<Dune::FieldMatrix<double,1,1>>; using Matrix = Dune::BCRSMatrix<double>;
using ElementMatrix = Dune::Matrix<Dune::FieldMatrix<double,1,1>>; using ElementMatrix = Dune::Matrix<double>;
using ElementRhs = Dune::BlockVector<Dune::FieldVector<double,1>>; using ElementRhs = Dune::BlockVector<double>;
rhs.resize(size); rhs.resize(size);
rhs = 0; rhs = 0;
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment