On Fri, Apr 8, 2022 at 8:25 AM 강다영 <dydy1...@pusan.ac.kr> wrote:

> Dear all,
>
> I have a question about the overflow error I met during the SCM training in
> Libmesh.
>
> For the training, I used the mesh model with a degree of freedom of
> 318,021.
> And I got an error message as below:
>
> [0]PETSC ERROR: --------------------- Error Message
> --------------------------------------------------------------
> [0]PETSC ERROR: No support for this operation for this object type
> [0]PETSC ERROR: Product of two integer 317065 317065 overflow, you must
> ./configure PETSc with --with-64-bit-indices for the case you are running
> [0]PETSC ERROR: See https://www.mcs.anl.gov/petsc/documentation/faq.html
> for trouble shooting.
> [0]PETSC ERROR: Petsc Release Version 3.12.2, unknown
> [0]PETSC ERROR: example-opt on a arch-linux2-c-debug named X570-AORUS-ELITE
> by dc308-ubuntu Thu Apr  7 12:46:34 2022
> [0]PETSC ERROR: Configure options --with-cc=mpicc --with-cxx=mpicxx
> --with-fc=mpif90 --download-fblaslapack --download-parmetis
> --download-metis --download-scalapack --download-hdf5 --download-superlu
> --download-mumps
> [0]PETSC ERROR: #1 PetscIntMultError() line 2167 in
> /home/dc308-ubuntu/Packages/petsc/include/petscsys.h
> [0]PETSC ERROR: #2 BVCreate_Svec() line 461 in
> /home/dc308-ubuntu/Packages/slepc/src/sys/classes/bv/impls/svec/svec.c
> [0]PETSC ERROR: #3 BVSetSizesFromVec() line 186 in
> /home/dc308-ubuntu/Packages/slepc/src/sys/classes/bv/interface/bvbasic.c
> [0]PETSC ERROR: #4 EPSAllocateSolution() line 570 in
> /home/dc308-ubuntu/Packages/slepc/src/eps/interface/epssetup.c
> [0]PETSC ERROR: #5 EPSSetUp_LAPACK() line 36 in
> /home/dc308-ubuntu/Packages/slepc/src/eps/impls/lapack/lapack.c
> [0]PETSC ERROR: #6 EPSSetUp() line 173 in
> /home/dc308-ubuntu/Packages/slepc/src/eps/interface/epssetup.c
> [0]PETSC ERROR: #7 EPSSolve() line 136 in
> /home/dc308-ubuntu/Packages/slepc/src/eps/interface/epssolve.c
>
> In the 3rd line of the error message, I noticed that the number "317065"
> matches with (the number of degree of freedom of model - the number of
> constrained degree of freedom of the model).
> Therefore, I decreased the degree of freedom to 7,503 to deal with this
> error.
> With this model, it was possible to finish the SCM training without an
> error.
> So I guess a large degree of freedom caused this overflow error.
>
> However, I'd still like to fix this overflow error to use the original
> model with a degree of freedom of 318,021 (for more accurate solutions).
>
> Could you possibly give some information on how to finish SCM training for
> a model with a large degree of freedom without the overflow error?
> And I'd appreciate it if you could let me know the maximum degree of
> freedom that will not cause the overflow error.
>

Hi Dayoung,

I'm speculating here, but it looks like you are using a LAPACK-based (hence
dense) solver, so that will limit the size of problem that you can solve.
This solver apparently needs to build a single array of size NDOFS**2, so
the largest size problem you will be able to run is likely approximately
NODFS = sqrt(INT_MAX) = sqrt(2147483647) = 46340. I don't know exactly how
to fix this issue, but it would probably involve using a different solver
at the very least. Finally, our mailing list is not really followed by many
people, so you may also have more luck by starting a GitHub Discussion (
https://github.com/libMesh/libmesh/discussions) on this topic.

-- 
John

_______________________________________________
Libmesh-users mailing list
Libmesh-users@lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/libmesh-users

Reply via email to