
<html>

<head>

<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">

<style type="text/css" style="display:none;"><!-- P {margin-top:0;margin-bottom:0;} --></style>

</head>

<body dir="ltr">

<div id="divtagdefaultwrapper" style="font-size:12pt;color:#000000;font-family:Calibri,Helvetica,sans-serif;" dir="ltr">

<p>Hello Dune developers and users,</p>

<p><br>

</p>

<p>I'm currently working on a solver for linear elasticity problems and already</p>

<p>implemented a version for "static" loadcases in Dune that works well. Now, I'm in the process</p>

<p>of adding time dependency and thus need to solve a system of ODE's of</p>

<p>the form: <span>y'-Ay = r</span>, where all components are nested datatypes.</p>

<p>Matrix A for example is of type Matrix<BCRSMatrix<FieldMatrix<double, 3, 3>>>,</p>

<p>the other variables are of the corresponding nested block-vector type.</p>

<p>The FIeldMatrix represents the Elementmatrix, the BCRSMatrix the global system</p>

<p>matrix (both for Mass and Stiffness) and the outer Matrix is of size 2x2 that couples the</p>

<p>system matrices and resembles the second order ode transformed into the system</p>

<p>of first order ones.<br>

</p>

<p><br>

</p>

<p>My question is, if there is an interface in Dune that can solve ODE's and handle those</p>

<p>nested datatypes. For example, if i rewrite the system to y'=F(...):</p>

<p>F[0] = y[1];</p>

<p>F[1] = -M^-1*K*y[0]+M^-1*p</p>

<p>that the basic operations can be done with the BCRSMatrix and BlockVectors.<br>

</p>

<p>Or if someone has recommendations on how to solve this sytem in</p>

<p>Dune (in a simple way, that is still performing quite well concerning efficiency).</p>

<p>I've seen that there are implemented time stepping methods in the Dune-FEM and

<br>

</p>

<p>Dune-PDELab modules, but I'm not sure what they are capable of.</p>

<br>

<p>Thanks for your help and best regards,</p>

<p>Max<br>

</p>

<p><br>

</p>

</div>

</body>

</html>