PETSc error 63

[Мария Бронзова]

Deal all,

I have been having troubles with running my script with PETSc. I have adjusted my code for parallel computing based on the example from step-40. But I am getting the following error:

----------------------------------------------------
Exception on processing:
--------------------------------------------------------
An error occurred in line <1483> of file </home/mariia/dealii-9.3.1/include/deal.II/lac/petsc_matrix_base.h> in function
    void dealii::PETScWrappers::MatrixBase::add(dealii::PETScWrappers::MatrixBase::size_type, dealii::PETScWrappers::MatrixBase::size_type, const size_type*, const PetscScalar*, bool, bool)
The violated condition was:
    ierr == 0
Additional information:
    deal.II encountered an error while calling a PETSc function.
    The description of the error provided by PETSc is "Argument out of
    range".
    The numerical value of the original error code is 63.
--------------------------------------------------------

I was trying to understand, where an argument gets out of range, but I didn't have any success so far. 

I created a very simple version of my script without any complications in it for an easier debugging, I attach it in this mail.

If you have any hints, I would be extremely grateful!

Thank you!

Regards,

Mariia

[Wolfgang Bangerth]

Mariia:
the error suggests that you are trying to write into an element of a matrix
that is outside the size of the matrix. You need to find out why that is so.
Run the program in a debugger and then inspect why the arguments you give to
that function are out of range.

Best
W.


On 2/3/22 07:03, Мария Бронзова wrote:
> *** Caution: EXTERNAL Sender ***

> --
> The deal.II project is located at http://www.dealii.org/
> <https://nam10.safelinks.protection.outlook.com/?url=http%3A%2F%2Fwww.dealii.org%2F&data=04%7C01%7CWolfgang.Bangerth%40colostate.edu%7Cf25b5705e9144246e7ca08d9e71deb29%7Cafb58802ff7a4bb1ab21367ff2ecfc8b%7C0%7C0%7C637794939344043914%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000&sdata=PzTqlU7N2wlcnwWHSk9iHykPHlibfbQwrQ2DuWiqtoY%3D&reserved=0>
> For mailing list/forum options, see
> https://groups.google.com/d/forum/dealii?hl=en
> <https://nam10.safelinks.protection.outlook.com/?url=https%3A%2F%2Fgroups.google.com%2Fd%2Fforum%2Fdealii%3Fhl%3Den&data=04%7C01%7CWolfgang.Bangerth%40colostate.edu%7Cf25b5705e9144246e7ca08d9e71deb29%7Cafb58802ff7a4bb1ab21367ff2ecfc8b%7C0%7C0%7C637794939344043914%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000&sdata=rp7fMX9v%2B%2F9cINjVEBNqBpt4TTOmGQz%2FNN7nWT8U0hI%3D&reserved=0>
> ---
> You received this message because you are subscribed to the Google Groups
> "deal.II User Group" group.
> To unsubscribe from this group and stop receiving emails from it, send an
> email to dealii+un...@googlegroups.com
> <mailto:dealii+un...@googlegroups.com>.
> To view this discussion on the web visit
> https://groups.google.com/d/msgid/dealii/7c5e3dfd-12e8-47be-9a71-b78674e4bd6an%40googlegroups.com
> <https://nam10.safelinks.protection.outlook.com/?url=https%3A%2F%2Fgroups.google.com%2Fd%2Fmsgid%2Fdealii%2F7c5e3dfd-12e8-47be-9a71-b78674e4bd6an%2540googlegroups.com%3Futm_medium%3Demail%26utm_source%3Dfooter&data=04%7C01%7CWolfgang.Bangerth%40colostate.edu%7Cf25b5705e9144246e7ca08d9e71deb29%7Cafb58802ff7a4bb1ab21367ff2ecfc8b%7C0%7C0%7C637794939344043914%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000&sdata=yeXssjQAv1vFbozrqfGAqEYqhaVpoOphDWwP8IlWsi8%3D&reserved=0>.


--
------------------------------------------------------------------------
Wolfgang Bangerth email: bang...@colostate.edu
www: http://www.math.colostate.edu/~bangerth/

[Мария Бронзова]

Thank you a lot for the reply, I found my mistake, it was a dull one: I removed the sparsity pattern creation line by mistake! Now it is running but the output file is unreadable, somewhat like:

<?xml version="1.0" ?> <!-- # vtk DataFile Version 3.0 #This file was generated by the deal.II library on 12:29:29 at 2022/2/4 --> <VTKFile type="UnstructuredGrid" version="0.1" compressor="vtkZLibDataCompressor" byte_order="LittleEndian"> <UnstructuredGrid> <Piece NumberOfPoints="8" NumberOfCells="1" > <Points> <DataArray type="Float32" NumberOfComponents="3" format="binary"> AQAAAGAAAABgAAAALgAAAA==eNo1ioEJADAIw/L/T9MH3E+rrBWCUQI04s4cUXHvyj+Oe0F6/rRv0q8/x+QrmQ== </DataArray> </Points> <Cells> <DataArray type="Int32" Name="connectivity" format="binary"> AQAAACAAAAAgAAAAHQAAAA==eNpjYGBgYARiZiBmAmIWIGYFYnYgZgNiAAF4AB0= </DataArray> <DataArray type="Int32" Name="offsets" format="binary"> AQAAAAQAAAAEAAAADAAAAA==eNrjYGBgAAAAJAAJ </DataArray> <DataArray type="UInt8" Name="types" format="binary"> AQAAAAEAAAABAAAACQAAAA==eNrjAQAADQAN </DataArray> </Cells> <PointData Scalars="scalars"> <DataArray type="Float32" Name="u1_re" format="binary"> AQAAACAAAAAgAAAAHAAAAA==eNpjYEAF8sKCm9IThYwk+QU3daYKGQEAIVQEHw== </DataArray> <DataArray type="Float32" Name="u1_im" format="binary"> AQAAACAAAAAgAAAAHAAAAA==eNpjYEAFloKdhmuZezZaSXQamvzt3ggALskGaw== </DataArray> 

Would you possibly have a hint on why that could be happening? 

Thank you!

Kind regards,

Mariia

четверг, 3 февраля 2022 г. в 15:50:15 UTC+1, Wolfgang Bangerth:

[Wolfgang Bangerth]

On 2/4/22 04:42, Мария Бронзова wrote:
>
> Would you possibly have a hint on why that could be happening?

This is how VTU files are supposed to look like. It's the format Visit and
Paraview understand. What were you hoping to see?

Best
W.

[Мария Бронзова]

Oh, I see. I just got used to the output formats I used to get before I tried to introduce parallel computing. Wasn't expecting to see letters instead of numbers. Thought I might have missed some command in the code.

Thank you for the clarification!

Regards,

Mariia