petsc4Foam with compressible/lagrangian solver using a Pointwise mesh
Dear All, I am attempting to use petsc4Foam on the case described in this paper: https://www.researchgate.net/publication/346511359_Eulerian-Lagrangian_modelling_of_cough_droplets_irradiated_by_ultraviolet-C_light_in_relation_to_SARS-CoV-2_transmission. I have obtained the following error:
Initializing PETSc [0]PETSC ERROR: --------------------- Error Message ----------------------------
[0]PETSC ERROR: Argument out of range [0]PETSC ERROR: New nonzero at (0,4) caused a malloc Use MatSetOption(A, MAT_NEW_NONZERO_ALLOCATION_ERR, PETSC_FALSE) to turn off thi s check [0]PETSC ERROR: See https://urldefense.com/v3/https://www.mcs.anl.gov/petsc/documentation/faq.html;!!P1tgJ-3e!W6KyuR6BLWafSclZMXKFtNT43e2XJjhBf__3uWqb9dGMTaVFHSRP2FpSy4fPt35D$ for tro uble shooting. [0]PETSC ERROR: Petsc Release Version 3.13.3, Jul 01, 2020 [0]PETSC ERROR: Unknown Name on a named https://urldefense.com/v3/http://cresco6x077.portici.enea.it;!!P1tgJ-3e!W6KyuR6BLWafSclZMXKFtNT43e2XJjhBf__3uWqb9dGMTaVFHSRP2FpSyzVjLBlc$ by vdales S un Feb 7 11:04:02 2021 [0]PETSC ERROR: Configure options --CFLAGS="-O3 -xCORE-AVX512 -xHOST" --F90LAGS= "-O3 -xCORE-AVX512 -xHOST" --FFLAGS="-O3 -xHOST" --with-blaslapack-lib="-L/opt/i ntel/mkl/lib/intel64 -lmkl_intel_lp64 -lmkl_sequential -lmkl_core -lpthread -lm"  --with-debugging=no --with-dynamic-libraries --with-mpi-cc=/opt/intel/impi/2017 .3.196/intel64/bin --with-pic --with-shared-libraries --with-vendor-compilers=in tel --with-64-bit-indices=0 --with-precision=double --prefix=/gporq2/scratch_0/u sr/vdales/OpenFOAM/ThirdParty-v2006/platforms/linux64IccDPInt32/petsc-git PETSC_ ARCH=linux64IccDPInt32Opt
I have attempted to introduce
MatSetOption(Amat, MAT_NEW_NONZERO_ALLOCATION_ERR, PETSC_FALSE) ;
in computeMatAllocation. However, the solver crashes after 2/3 time-steps.
The speed-up in these steps is really significant. So, we are really intersted
on petsc4Foam.
Could you help me to fix this issue ?
Sincerely,
Valerio