<div dir="ltr"><div class="gmail_default" style="font-family:tahoma,sans-serif;font-size:small">Thanks René.</div><div class="gmail_default" style="font-family:tahoma,sans-serif">The pdes that I'm working on are highly nonlinear and I'm trying to reduce the assembly time.</div><div class="gmail_default" style="font-family:tahoma,sans-serif">I'm wondering whether it is possible to implement jacobian_volume,i.e. analytic jacobian and for the skeleton and boundary parts just use the numerical jacobian(default)?  </div><div class="gmail_default" style="font-family:tahoma,sans-serif">Best</div><div class="gmail_default" style="font-family:tahoma,sans-serif;font-size:small">Amir</div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Tue, Jul 7, 2020 at 5:12 PM René Heß <<a href="mailto:rene.hk-edv@gmx.de">rene.hk-edv@gmx.de</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">Hi Amir,<br>
<br>
<br>
the jacobian methods will assemble the Jacobian matrix J. This matrix<br>
will be used by the linear solvers to solve your system.<br>
<br>
Some solvers only need J to calculate the application of J to a vector<br>
x. This means the solvers only needs to be able to compute Jx for a<br>
given vector x. With the jacobian apply methods it is possible to<br>
calculate this application of the Jacobian without ever assembling J.<br>
<br>
This can be used to implement matrix-free solvers that could be faster<br>
for certain cases. In your case you can ignore this issue and you won't<br>
need the jacobian_apply methods (or if the compiler complains you can<br>
just use the numerical methods, they won't be used anyway).<br>
<br>
<br>
Best regards,<br>
René<br>
<br>
<br>
Amir Peiraviminaei <<a href="mailto:apminaei@gmail.com" target="_blank">apminaei@gmail.com</a>> writes:<br>
<br>
> Hi everyone,<br>
> I'm using pdelab for a FV scheme and tutorial 02 as an example. So far I<br>
> have used only numerical jacobian, Now I want to implement analytical<br>
> jacobian myself, my question is<br>
> 1. What is the difference between jacobian_volume and jacobian_apply_volume?<br>
> (tutorial 02 uses both)<br>
><br>
> Best<br>
> Amir<br>
> _______________________________________________<br>
> dune-pdelab mailing list<br>
> <a href="mailto:dune-pdelab@lists.dune-project.org" target="_blank">dune-pdelab@lists.dune-project.org</a><br>
> <a href="https://lists.dune-project.org/mailman/listinfo/dune-pdelab" rel="noreferrer" target="_blank">https://lists.dune-project.org/mailman/listinfo/dune-pdelab</a><br>
</blockquote></div>