<html>
<head>
<meta content="text/html; charset=ISO-8859-1"
http-equiv="Content-Type">
</head>
<body text="#000000" bgcolor="#FFFFFF">
<div class="moz-cite-prefix">Hey Steffen, <br>
<br>
thank you for the fast and helpful answer. Indeed, I overlooked
that I reimplement the jacobian_volume method where I explicitly
use a LocalVector to call the residual methods. That caused the
warnings. <br>
<br>
Kind regards <br>
Bernd <br>
<br>
On 08/20/2012 02:18 PM, Steffen Müthing wrote:<br>
</div>
<blockquote
cite="mid:31D1C2A6-3374-4DDA-945C-B58D42113258@ipvs.uni-stuttgart.de"
type="cite">
<pre wrap="">Hi Bernd,
Am 20.08.2012 um 14:01 schrieb Bernd Flemisch:
</pre>
<blockquote type="cite">
<pre wrap="">Dear PDELab,
I am using the 1.0 release. I get warnings like
../../../dumux/boxmodels/linearelasticity/2p/el2plocaloperator.hh:453:33: warning: ‘void Dune::PDELab::LocalVector<T, LFSFlavorTag, W>::accumulate(const LFS&, Dune::PDELab::LocalVector<T, LFSFlavorTag, W>::size_type, Dune::PDELab::LocalVector<T, LFSFlavorTag, W>::value_type) [with LFS = Dune::PDELab::LeafLocalFunctionSpaceNode<Dune::PDELab::GridFunctionSpace<Dune::GridView<Dune::DefaultLeafGridViewTraits<const Dune::ALUCubeGrid<3, 3>, (Dune::PartitionIteratorType)4u> >, Dune::PDELab::Q1LocalFiniteElementMap<double, double, 3>, Dune::PDELab::NoConstraints, Dune::PDELab::ISTLVectorBackend<1>, Dune::PDELab::GridFunctionGeneralMapper>, Dune::PDELab::MultiIndex<long unsigned int, 3u> >, T = double, LFSFlavorTag = Dune::PDELab::AnySpaceTag, W = double, Dune::PDELab::LocalVector<T, LFSFlavorTag, W>::size_type = long unsigned int, Dune::PDELab::LocalVector<T, LFSFlavorTag, W>::value_type = double]
’ i
s deprecated (declared at /temp/bernd/DUMUX22/dune-pdelab/dune/pdelab/gridfunctionspace/localvector.hh:349) [-Wdeprecated-declarations]
>From taking a look at gridfunctionspace/localvector.hh, I would think that I should use the accumulate method of the WeightedVectorAccumulationView rather than of the LocalVector. However, I do not know how to enforce this, also not by comparing with non-warning producing tests. Can you please help me? I am using the types given in the message above and also listed below.
</pre>
</blockquote>
<pre wrap="">
hmm, that's strange. Did you create the LocalVector yourself, or is it an object that the GridOperator passed to your LocalOperator?
Because in that case you shouldn't get a warning… ;-)
If you created the LocalVector yourself, you can just extract an AccumulationView using the nested type LocalVector::WeightedAccumulationView
and the member method localVector.weightedAccumulationView(double). The double parameter is the weight that will be applied to the accumulated
entries (just set it to 1 if you don't need it). For an example of this mechanism in action take a look at dune/pdelab/localoperator/defaultimp.hh. The
numerical differentiation code in there creates its own residual vectors from the jacobian passed in by the GridOperator, wraps them in
WeightedAccumulationViews and passes those on to the alpha_*() methods of the LocalOperator.
By the way, you also don't get a LocalMatrix for the jacobian from the GridOperator, but a LocalMatrix::WeightedAccumulationView. ;-)
</pre>
<blockquote type="cite">
<pre wrap="">
Thank you. Kind regards
Bernd
</pre>
</blockquote>
<pre wrap="">
I hope that clears things up a little! ;-)
Greetings,
Steffen
</pre>
<blockquote type="cite">
<pre wrap="">
enum{dim = GridView::dimension};
typedef double Scalar;
typedef Dune::PDELab::NoConstraints Constraints;
typedef Dune::PDELab::Q1LocalFiniteElementMap<Scalar,Scalar,dim> DisplacementFEM;
typedef Dune::PDELab::GridFunctionSpace<GridView,
DisplacementFEM,
Constraints,
Dune::PDELab::ISTLVectorBackend<1>
</pre>
<blockquote type="cite">
<pre wrap="">ScalarDisplacementGFS;
</pre>
</blockquote>
<pre wrap="">typedef Dune::PDELab::PowerGridFunctionSpace<ScalarDisplacementGFS,
dim,
Dune::PDELab::GridFunctionSpaceBlockwiseMapper
</pre>
<blockquote type="cite">
<pre wrap="">DisplacementGFS;
</pre>
</blockquote>
<pre wrap="">
typedef Dune::PDELab::Q1LocalFiniteElementMap<Scalar,Scalar,dim> PressureFEM;
typedef Dune::PDELab::GridFunctionSpace<GridView,
PressureFEM,
Constraints,
Dune::PDELab::ISTLVectorBackend<1>
</pre>
<blockquote type="cite">
<pre wrap="">ScalarPressureGFS;
</pre>
</blockquote>
<pre wrap="">typedef Dune::PDELab::PowerGridFunctionSpace<ScalarPressureGFS,
2,
Dune::PDELab::GridFunctionSpaceBlockwiseMapper
</pre>
<blockquote type="cite">
<pre wrap="">PressureGFS;
</pre>
</blockquote>
<pre wrap="">
typedef Dune::PDELab::GridFunctionSpaceLexicographicMapper GFSMapper;
typedef Dune::PDELab::CompositeGridFunctionSpace<GFSMapper,
PressureGFS,
DisplacementGFS> GridFunctionSpace;
typedef typename GridFunctionSpace::template ConstraintsContainer<Scalar>::Type ConstraintsTrafo;
typedef Dumux::PDELab::El2PLocalOperator<TypeTag> LocalOperator;
typedef Dune::PDELab::GridOperator<GridFunctionSpace,
GridFunctionSpace,
LocalOperator,
Dune::PDELab::ISTLBCRSMatrixBackend<1, 1>,
Scalar, Scalar, Scalar,
ConstraintsTrafo,
ConstraintsTrafo,
true
</pre>
<blockquote type="cite">
<pre wrap="">GridOperator;
</pre>
</blockquote>
<pre wrap="">
The local operator looks like
template<class TypeTag>
class El2PLocalOperator
:
public Dune::PDELab::FullVolumePattern,
public Dune::PDELab::LocalOperatorDefaultFlags
{...};
--
_____________________________________________________________________
Bernd Flemisch phone: +49 711 685 69162
IWS, Universität Stuttgart fax: +49 711 685 60430
Pfaffenwaldring 61 email: <a class="moz-txt-link-abbreviated" href="mailto:bernd@iws.uni-stuttgart.de">bernd@iws.uni-stuttgart.de</a>
D-70569 Stuttgart url: <a class="moz-txt-link-abbreviated" href="http://www.hydrosys.uni-stuttgart.de">www.hydrosys.uni-stuttgart.de</a>
_____________________________________________________________________
_______________________________________________
dune-pdelab mailing list
<a class="moz-txt-link-abbreviated" href="mailto:dune-pdelab@dune-project.org">dune-pdelab@dune-project.org</a>
<a class="moz-txt-link-freetext" href="http://lists.dune-project.org/mailman/listinfo/dune-pdelab">http://lists.dune-project.org/mailman/listinfo/dune-pdelab</a>
</pre>
</blockquote>
<pre wrap="">
Steffen Müthing
Universität Stuttgart
Institut für Parallele und Verteilte Systeme
Universitätsstr. 38
70569 Stuttgart
Tel: +49 711 685 88429
Fax: +49 711 685 88340
Email: <a class="moz-txt-link-abbreviated" href="mailto:steffen.muething@ipvs.uni-stuttgart.de">steffen.muething@ipvs.uni-stuttgart.de</a>
</pre>
<br>
<fieldset class="mimeAttachmentHeader"></fieldset>
<br>
<pre wrap="">_______________________________________________
dune-pdelab mailing list
<a class="moz-txt-link-abbreviated" href="mailto:dune-pdelab@dune-project.org">dune-pdelab@dune-project.org</a>
<a class="moz-txt-link-freetext" href="http://lists.dune-project.org/mailman/listinfo/dune-pdelab">http://lists.dune-project.org/mailman/listinfo/dune-pdelab</a></pre>
</blockquote>
<br>
<br>
<pre class="moz-signature" cols="72">--
_____________________________________________________________________
Bernd Flemisch phone: +49 711 685 69162
IWS, Universität Stuttgart fax: +49 711 685 60430
Pfaffenwaldring 61 email: <a class="moz-txt-link-abbreviated" href="mailto:bernd@iws.uni-stuttgart.de">bernd@iws.uni-stuttgart.de</a>
D-70569 Stuttgart url: <a class="moz-txt-link-abbreviated" href="http://www.hydrosys.uni-stuttgart.de">www.hydrosys.uni-stuttgart.de</a>
_____________________________________________________________________
</pre>
</body>
</html>