This is a read only copy of the old FEniCS QA forum. Please visit the new QA forum to ask questions

Which solver to use in a case with an exponential?

0 votes

Hello

I am pretty new to fenics, and am looking for advice about solvers. So I have a system with four functions that i want to solve who all depend on each other, for instance

dq=q+p +e^(q-q1)

where q,p, and q1 are all functions. As this problem is not linear I get an error when simply using "solve". So I was wondering what nonlinear solver would be used in this case.

Thanks in advance for any advice.

asked Mar 13, 2017 by haakon555 FEniCS Novice (330 points)
edited Mar 13, 2017 by haakon555

Can you say why this problem is linear?

commented Mar 13, 2017 by nate FEniCS Expert (17,050 points)

Ah, sorry, I of course meant "since this problem is NOT linear".

commented Mar 13, 2017 by haakon555 FEniCS Novice (330 points)

Have you tried this: https://fenicsproject.org/olddocs/dolfin/1.6.0/python/demo/documented/nonlinear-poisson/python/documentation.html ?

commented Apr 16, 2017 by nabarnaf FEniCS User (2,940 points)
...