29 #ifndef __SystemImplicitLinear_H
30 #define __SystemImplicitLinear_H 1
44 template<
typename SpaceType>
63 typedef typename super::value_type value_type;
64 typedef typename super::functionspace_type functionspace_type;
65 typedef typename super::functionspace_type functionspace_ptrtype;
66 typedef typename super::element_type element_type;
69 typedef typename super::backend_ptrtype backend_ptrtype;
71 typedef typename super::sparse_matrix_ptrtype sparse_matrix_ptrtype;
73 typedef typename super::vector_ptrtype vector_ptrtype;
111 sparse_matrix_ptrtype
const&
lhs()
const
117 sparse_matrix_ptrtype&
lhs()
123 vector_ptrtype
const&
rhs()
const
153 vector_ptrtype U( this->
backend->newVector( _test=this->functionSpace() ) );
154 this->
backend->solve( M_lhs, M_lhs, U, M_rhs );
164 sparse_matrix_ptrtype M_lhs;
165 vector_ptrtype M_rhs;
172 template<
typename SpaceType>
173 SystemImplicitLinear<SpaceType>::SystemImplicitLinear( functionspace_ptrtype
const& Xh,
174 po::variables_map
const& vm )
177 M_lhs( this->backend()->newMatrix( Xh, Xh ) ),
178 M_rhs( this->backend()->newVector( Xh ) )
180 template<
typename SpaceType>
181 SystemImplicitLinear<SpaceType>::SystemImplicitLinear( SystemImplicitLinear
const& sil )