
<html dir="ltr">

<head>

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

<style id="owaParaStyle" type="text/css">P {margin-top:0;margin-bottom:0;}</style>

</head>

<body ocsi="0" fpstyle="1">

<div style="direction: ltr;font-family: Tahoma;color: #000000;font-size: 10pt;">Hi all,<br>

<br>

I have a quick question about the preconditioning.<br>

<br>

So far I have written a linear operator for a system I need to solve.<br>

Suppose I want to use a BiCGSTAB or a GMRES to solve the system; one option is to create an instance of Dune::Fem::OEMBICGSTABOp.<br>

Since I want to pass to the solver my own preconditioner (derived from one paper in literature), how can I provide it correctly?<br>

In addition, since it's a block preconditioner, I suspect I'll have to provide a solve () method, i.e. how to "invert" it numerically.<br>

<br>

How can I do all these things in Dune-Fem?<br>

Is there any special restriction with the type of grid functions (Adaptive vs. BlockVector type)?<br>

<br>

Thanks in advance!<br>

<br>

Andrea<br>

<div>

<div style="font-family:Tahoma; font-size:13px">__________________________________________________________<br>

<br>

Andrea Sacconi<br>

PhD student, Applied Mathematics<br>

AMMP Section, Department of Mathematics, Imperial College London,<br>

London SW7 2AZ, UK<br>

a.sacconi11@imperial.ac.uk<br>

</div>

</div>

</div>

</body>

</html>