====================== *** PetIBM - Start *** ====================== directory: . Parsing file ./cartesianMesh.yaml... done. Parsing file ./flowDescription.yaml... done. Parsing file ./simulationParameters.yaml... done. --------------------------------------- Cartesian grid --------------------------------------- number of cells: 20 x 20 --------------------------------------- --------------------------------------- Flow --------------------------------------- dimensions: 2 viscosity: 0.01 initial velocity field: 0. 0. boundary conditions (component, type, value): ->location: xMinus (left) 0 DIRICHLET 0. 1 DIRICHLET 0. ->location: xPlus (right) 0 DIRICHLET 0. 1 DIRICHLET 0. ->location: yMinus (bottom) 0 DIRICHLET 0. 1 DIRICHLET 0. ->location: yPlus (top) 0 DIRICHLET 1. 1 DIRICHLET 0. --------------------------------------- --------------------------------------- Time-stepping --------------------------------------- formulation: Navier-Stokes solver (Perot, 1993) convection: Euler-explicit diffusion: Euler-implicit time-increment: 0.0005 starting time-step: 0 number of time-steps: 1 saving-interval: 1 --------------------------------------- ---------------------------------------- KSP info: Velocity system ---------------------------------------- KSP Object:(velocity_) 2 MPI processes type: bcgs maximum iterations=10000 tolerances: relative=1e-08, absolute=0., divergence=10000. left preconditioning using nonzero initial guess using DEFAULT norm type for convergence test PC Object:(velocity_) 2 MPI processes type: jacobi PC has not been set up so information may be incomplete linear system matrix = precond matrix: Mat Object: 2 MPI processes type: mpiaij rows=760, cols=760 total: nonzeros=3644, allocated nonzeros=3800 total number of mallocs used during MatSetValues calls =0 not using I-node (on process 0) routines ---------------------------------------- KSP info: Poisson system ---------------------------------------- KSP Object:(poisson_) 2 MPI processes type: cg maximum iterations=20000 tolerances: relative=1e-08, absolute=0., divergence=10000. left preconditioning using nonzero initial guess using DEFAULT norm type for convergence test PC Object:(poisson_) 2 MPI processes type: gamg PC has not been set up so information may be incomplete MG: type is MULTIPLICATIVE, levels=0 cycles=unknown Cycles per PCApply=0 Using Galerkin computed coarse grid matrices GAMG specific options Threshold for dropping small values from graph 0. AGG specific options Symmetric graph false linear system matrix = precond matrix: Mat Object: 2 MPI processes type: mpiaij rows=400, cols=400 total: nonzeros=1920, allocated nonzeros=1920 total number of mallocs used during MatSetValues calls =0 not using I-node (on process 0) routines [time-step 0] Writing grid into file... done. Residual norms for velocity_ solve. 0 KSP preconditioned resid norm 1.148589971947e-03 true resid norm 2.324746103222e+00 ||r(i)||/||b|| 1.000000000000e+00 1 KSP preconditioned resid norm 8.875846864181e-09 true resid norm 1.793849713607e-05 ||r(i)||/||b|| 7.716325284387e-06 2 KSP preconditioned resid norm 1.418061760018e-13 true resid norm 2.860743141913e-10 ||r(i)||/||b|| 1.230561538720e-10 Linear velocity_ solve converged due to CONVERGED_RTOL iterations 2 KSP Object:(velocity_) 2 MPI processes type: bcgs maximum iterations=10000 tolerances: relative=1e-08, absolute=0., divergence=10000. left preconditioning using nonzero initial guess using PRECONDITIONED norm type for convergence test PC Object:(velocity_) 2 MPI processes type: jacobi linear system matrix = precond matrix: Mat Object: 2 MPI processes type: mpiaij rows=760, cols=760 total: nonzeros=3644, allocated nonzeros=3800 total number of mallocs used during MatSetValues calls =0 not using I-node (on process 0) routines ** On entry to DGESVD parameter number 6 had an illegal value