Re: [petsc-users] Large MATMPIAIJ - 32bit integer overflow in nz value

2023-05-17 Thread Barry Smith
Yeah, this is silly. The check is just a "sanity-check" on the data in the file. We store redundant information in the matrix header in the file, header[3] is the total number of nonzeros in the matrix. When nz is too large, the correct value cannot fit in the header. Changing the file for

Re: [petsc-users] Large MATMPIAIJ - 32bit integer overflow in nz value

2023-05-17 Thread Matthew Knepley
On Wed, May 17, 2023 at 9:02 AM Fleischli Benno HSLU T&A < benno.fleisc...@hslu.ch> wrote: > Dear PETSc developers > > I am creating a very large parallel sparse matrix (MATMPIAIJ) with PETSc. > I write this matrix to disk. > The number of non-zeros exceeds the maximum number a 32-bit integer can

[petsc-users] Large MATMPIAIJ - 32bit integer overflow in nz value

2023-05-17 Thread
Dear PETSc developers I am creating a very large parallel sparse matrix (MATMPIAIJ) with PETSc. I write this matrix to disk. The number of non-zeros exceeds the maximum number a 32-bit integer can hold. When I read the matrix from disk i get an error because there was an overflow in the nz numbe