[Dune] C++ question

[Mike Rosing]

I'm trying to create a BCRSMatrix of runtime sized Matrix. I tried

   typedef Matrix<double> occupy(np, np);
   BCRSMatrix<occupy> stiffns(4*numelement, 4*numelement, 9, .12,
                              BCRSMatrix<occupy>::implicit);


and several other iterations - but I obviously don't understand what I'm 
doing.  The "occupy" matrix is mostly full, the "stiffns" matrix is mostly
empty.

In case it matters, the list of includes is
#include <iostream>
#include <fstream>
#include <dune/common/parallel/mpihelper.hh> // An initializer of MPI
// #include <dune/common/exceptions.hh> // We use exceptions
#include <dune/grid/uggrid.hh>
#include <dune/grid/io/file/gmshreader.hh>
#include <dune/geometry/affinegeometry.hh>
#include <dune/istl/matrix.hh>

// { include_matrix_vector_begin }
#include <dune/istl/bcrsmatrix.hh>
#include <dune/istl/bvector.hh>
// { include_matrix_vector_end }

#include <dune/istl/matrixindexset.hh>
#include <dune/istl/preconditioners.hh>
#include <dune/istl/solvers.hh>
#include <dune/istl/matrixmarket.hh>

What's the correct way to set up this kind of structure?

Mike Rosing
Engineering Research
Madison WI USA
Elliptic Curve Crypto:  http://mng.bz/D9NA