dolfin team mailing list archive

[Siva Shankar Rudraraju <question120034@xxxxxxxxxxxxxxxxxxxxx>]

New question #120034 on DOLFIN:
https://answers.launchpad.net/dolfin/+question/120034


We have a multifield problem set up using assemble(...) and solve(...). 

ffc_options0 = {"quadrature_degree": 2, "optimize": True,  "eliminate_zeros": True, "precompute_ip_const": True}

A = assemble(a, form_compiler_parameters=ffc_options0)
b = assemble(-L, form_compiler_parameters=ffc_options0)
bcb.apply(A, b); bcro.apply(A, b); bcrg.apply(A, b);
solve(A, duc.vector(), b)

Running it in serial gives the desired results, but with mpirun the assemble process is slow, and later its stuck forever in the linear solve stage, as shown below.  Any thoughts on what's going wrong?.

$ mpirun -np 6 python tumor11.py
Process 1: Partitioned mesh, edge cut is 997.
Process 0: Partitioned mesh, edge cut is 997.
Process 2: Partitioned mesh, edge cut is 997.
Process 3: Partitioned mesh, edge cut is 997.
Process 4: Partitioned mesh, edge cut is 997.
Process 5: Partitioned mesh, edge cut is 997.
Process 0: Got expression dimension = 8
Process 5: Got expression dimension = 8
Process 4: Got expression dimension = 8
Process 2: Got expression dimension = 8
Process 1: Got expression dimension = 8
Process 3: Got expression dimension = 8
Time = 3600.0
Time = 3600.0
Time = 3600.0
Time = 3600.0
Time = 3600.0
Time = 3600.0
Reset matrix (B)
Reset matrix (B)
Reset matrix (B)
Reset matrix (B)
Reset matrix (B)
Reset matrix (B)
Process 3: *** Warning: Found no facets matching domain for boundary condition.
Process 3: Applying boundary conditions to linear system.
Process 5: *** Warning: Found no facets matching domain for boundary condition.
Process 5: Applying boundary conditions to linear system.
Process 4: *** Warning: Found no facets matching domain for boundary condition.
Process 4: Applying boundary conditions to linear system.
Process 1: Applying boundary conditions to linear system.
Process 2: Applying boundary conditions to linear system.
Process 0: Applying boundary conditions to linear system.
Process 4: Applying boundary conditions to linear system.
Process 1: Applying boundary conditions to linear system.
Process 2: Applying boundary conditions to linear system.
Process 0: Applying boundary conditions to linear system.
Process 5: Applying boundary conditions to linear system.
Process 3: Applying boundary conditions to linear system.
Process 4: Applying boundary conditions to linear system.
Process 1: Applying boundary conditions to linear system.
Process 0: Applying boundary conditions to linear system.
Process 5: Applying boundary conditions to linear system.
Process 3: Applying boundary conditions to linear system.
Process 2: Applying boundary conditions to linear system.
Process 1: Solving linear system of size 13720 x 13720 (PETSc LU solver, mumps).
Process 4: Solving linear system of size 13720 x 13720 (PETSc LU solver, mumps).
Process 3: Solving linear system of size 13720 x 13720 (PETSc LU solver, mumps).
Process 2: Solving linear system of size 13720 x 13720 (PETSc LU solver, mumps).
Process 0: Solving linear system of size 13720 x 13720 (PETSc LU solver, mumps).
Process 5: Solving linear system of size 13720 x 13720 (PETSc LU solver, mumps).

and is stuck at this point forever.


-- 
You received this question notification because you are a member of
DOLFIN Team, which is an answer contact for DOLFIN.