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: > > 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 *** > > > > 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 > > > > -- > > 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/ > > -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, see https://groups.google.com/d/forum/dealii?hl=en --- 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+unsubscr...@googlegroups.com. To view this discussion on the web visit https://groups.google.com/d/msgid/dealii/16c13325-352c-4033-85f1-0e2e1011e98an%40googlegroups.com.