> First problem I face is that I cannot explicitly define the size of
> the matrix, so I cannot write
> template <int dim>
> struct QuadratureHistory
> {
> dealii::SymmetricTensor<2,dim> effective_eigenstrain;
> dealii::SymmetricTensor<2,dim> strain;
> dealii::FullMatrix<double> effective_stiffness(6,6);
> };
> it cause a compiler error:
Yes, what you have is the *declaration* of a structure and you can't
initialize members there like you want to do. You'll have to write a
constructor in which you can initialize the member variable in a way
similar to the above.
> quadrature_history[i].effective_stiffness(6,6);
All indices in deal.II are zero-based, so if your matrix has size 6x6, then
valid indices are 0...5. That said, from this error...
> Index 6 is not in [0,0[
...it seems that you are accessing a matrix that hasn't even be set to size
6x6 yet.
W.
-------------------------------------------------------------------------
Wolfgang Bangerth email: [email protected]
www: http://www.math.tamu.edu/~bangerth/
_______________________________________________
dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
