<html>
<head>
<meta content="text/html; charset=utf-8" http-equiv="Content-Type">
</head>
<body bgcolor="#FFFFFF" text="#000000">
Hey Dominic,<br>
<br>
thank you for your answer. I am solving the Richards Equation
(nonlinear PDE) with the DG method, essentially using:<br>
<br>
typedef
Dune::PDELab::QkDGLocalFiniteElementMap<Real,Real,order,dim>
FEM; as finite element map,<br>
typedef Dune::PDELab::ISTLBackend_BCGS_AMG_SSOR<IGO> LS;
as linear solver,<br>
typedef Dune::PDELab::Newton<IGO,LS,U> PDESOLVER; as
non-linear solver,<br>
Dune::PDELab::Alexander2Parameter<Real> alex2; as
timestepper.<br>
<br>
The grid is set up with<br>
<br>
Dune::FieldVector<double,dim> N;<br>
for (unsigned int i=0; i<dim; i++)<br>
N[i] = elements[i]; <br>
Dune::FieldVector<bool,dim> periodic(false);<br>
int overlap = 2;<br>
const YaspXDirPartition<dim> yp;<br>
Dune::YaspGrid<dim>
grid(helper.getCommunicator(),upperRight,N,periodic,overlap,&yp);<br>
typedef Dune::YaspGrid<dim>::LeafGridView GV;<br>
const GV& gv=grid.leafGridView();<br>
<br>
I have appended the source code responsible for the actual solution,
along with two logfiles - one with barely over 2000 elements in the
grid, one barely below. You can see that the amount of linear
iterations (and thus linear solve time) rises sharply if the
elements exceed 2000.<br>
<br>
I hope this is enough information to address this problem - if not,
please tell me so.<br>
<br>
<br>
Thank you kindly for your time!<br>
<br>
Dion<br>
<br>
<br>
<div class="moz-cite-prefix">Am 21.01.15 um 10:23 schrieb Dominic
Kempf:<br>
</div>
<blockquote
cite="mid:CA+Cg7M5ffcr26RrjHmqxZw_p-NXBnGLz6dn=PKaMiTVxnXpoXQ@mail.gmail.com"
type="cite">
<div dir="ltr">
<div>
<div>Hello Dion,<br>
<br>
</div>
to answer your question we need some information on the exact
problem you are using for this benchmark. Can you provide us
the code? I will refrain from any wild guesses that come to my
mind before having seen the actual problem.<br>
<br>
</div>
Best,<br>
Dominic<br>
</div>
<div class="gmail_extra"><br>
<div class="gmail_quote">On Tue, Jan 20, 2015 at 7:33 PM, Dion
Häfner <span dir="ltr"><<a moz-do-not-send="true"
href="mailto:dhaefner@iup.uni-heidelberg.de"
target="_blank">dhaefner@iup.uni-heidelberg.de</a>></span>
wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0
.8ex;border-left:1px #ccc solid;padding-left:1ex">Dear
PDELab Developers,<br>
<br>
when I was benchmarking my DUNE / PDELab module, I noticed a
sudden jump in runtime (approx. by a factor of 4) as soon as
the number of grid elements exceeded 2000 in a YaspCube
grid. I could not reproduce this behavior with UG Grids, so
I think it is somewhat specific to YaspCube Grids.<br>
<br>
Has anyone encountered a similar behavior before? How can
this be adressed?<br>
<br>
Sorry if my question is not very specific. If you need
further information on my module or the methods used, please
ask and I will provide.<br>
<br>
Thank you in advance,<br>
<br>
Dion Häfner<br>
<br>
_______________________________________________<br>
dune-pdelab mailing list<br>
<a moz-do-not-send="true"
href="mailto:dune-pdelab@dune-project.org" target="_blank">dune-pdelab@dune-project.org</a><br>
<a moz-do-not-send="true"
href="http://lists.dune-project.org/mailman/listinfo/dune-pdelab"
target="_blank">http://lists.dune-project.org/mailman/listinfo/dune-pdelab</a><br>
</blockquote>
</div>
<br>
</div>
</blockquote>
<br>
</body>
</html>