I’m solving a 3D heat problem. This is a minimal working version of the code: load "msh3"; func T0 = x < .5; int p = 10; mesh3 dominio = cube(p, p, p); fespace Vh(dominio, P2, periodic = [ [1, x, z], [3, x, z], [2, y, z], [4, y, z], [5, x, y], [6, x, y] ] ); Vh T = T0, v, Tant; real dt = 0.…

Indeed, using UMFPACK64 solves the problem. Thanks!

[image] prj: UMFPACK64. solver UMFPACK64，May it be used to solve PDE in complex fields?

Umpfack errors with periodic boundary and slightly large mesh in 3D problem

prj March 8, 2020, 11:11am 2

Use another solver, like MUMPS, or even UMFPACK64.

Topic		Replies	Views
Simple problem gives "Error umfpack_di_solve status -3" General Discussion	1	1204	May 5, 2021
3D error and warning General Discussion	1	419	March 13, 2021
Error occurs in UMFPACK solver 3D problem General Discussion	2	179	February 8, 2024
Umpfack error in heat conduction problem General Discussion	1	735	June 2, 2021
Problem with periodic bc General Discussion	4	461	January 31, 2021