examples.diffusion.nthOrder.input4thOrder_line

>>> eq.solve(var,
...          boundaryConditions = BCs,
...          solver = solver)

Using the Pysparse solvers, the answer is totally inaccurate. This is due to the 4th order term having a high matrix condition number. In this particular example, multigrid preconditioners such as those provided by Trilinos allow a more accurate solution.

>>> print(var.allclose(mesh.cellCenters[0], atol = 10))
1

Last updated on Jun 26, 2024. Created using Sphinx 7.1.2.