Dear Prof. Wolfgang, I appreciate your answer.
The exact same steps I used for elliptic problem with FE_Q and solution that is written in vector. It was working perfectly. New things here, we are using fe system and block vector, others should be the same. Best, Najwa On Mon, 2 Dec 2024 at 11:42 PM Wolfgang Bangerth <bange...@colostate.edu> wrote: > > Najwa: > I assume you have taken a look at step-83, the tutorial program that > demonstrates serialization? > https://dealii.org/developer/doxygen/deal.II/step_83.html > > Separately, the way to debug these sorts of problems is to create as > small a program as possible that still shows the issue. It doesn't have > to do anything useful, for example you can replace assemble_system() and > solve() by just setting the vector to ones. Can you come up with such a > minimal program? > > Best > W. > > > On 12/2/24 06:34, Najwa Alshehri wrote: > > *** Caution: EXTERNAL Sender *** > > > > Dear All, > > > > Thank you for your always support. > > > > I am trying to serialize a Stokes solution in 2d ( dim for velocity +1 > > for pressure) to use it as a reference solution. Please stop me here if > > there is a better way to do this. > > > > Unofortunatly, for some reason that I do not know, the saved_solution is > > all zeros or code breaks. > > > > To give more details I am solving using : > > , fe(FE_Q<dim>(degree + 1), dim, FE_Q<dim>(degree), 1) > > with BlockVector solution. This solution is not zeros ( I checked that). > > > > Then I open a file and save in order the following: > > <<triangulation > > <<fe > > <<dof_handler > > <<solution > > using > > boost::archive::text_iarchive > > This finishes with no issues. I even checked if the file is unreadable > > or empty, it is not. > > > > When I deserialize to check the saved solution is matching my computed > > solution, I use same fe > > FESystem<dim> saved_fe(FE_Q<dim>(degree + 1), dim, FE_Q<dim>(degree), 1); > > and deserialize in the following order > > >>saved_triangulation > > >>safe_fe > > reint and distribute_dofe for the saved_dof_handler using the saved > > triangulation and fe. > > >>saved_dof_handler > > >>saved solution > > > > > > After a deep investigation, it seems like " >>saved_dof_handler" is > > where the code beaks with error > > "---------------------------------------------------- > > Exception on processing: > > input stream error > > Aborting! > > ---------------------------------------------------- " > > which means that saved_dof_handler cannot be read form the saved file. > > However, I checked after distributing dofs and "saved_dof_handler" has > > the correct n_dofs(). > > > > In this code I was just checking, so dof_handler is not really needed > > here, but in the real code it is need to write a function out of the > > solution. For the sake of debugging, I removed " >>saved_dof_handler". > > Removing it gives all zeros solution vector. This is suspicious, is it > > connected to the saved_dof_handler?? > > > > I think I have checked many other things to find the source of the issue > > with no success. Conclusion, saved triangulation and saved fe seems to > > be correct. > > > > Any one knows what is happening here? I would appreciate any guidance > > one what could be a source for the issue. > > > > Best, > > Najwa > > > > -- > > 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=05%7C02%7CWolfgang.Bangerth% > 40colostate.edu > %7Cd23d2bd47fb74e7a13d008dd12d617b6%7Cafb58802ff7a4bb1ab21367ff2ecfc8b%7C0%7C0%7C638687432905944661%7CUnknown%7CTWFpbGZsb3d8eyJFbXB0eU1hcGkiOnRydWUsIlYiOiIwLjAuMDAwMCIsIlAiOiJXaW4zMiIsIkFOIjoiTWFpbCIsIldUIjoyfQ%3D%3D%7C0%7C%7C%7C&sdata=aL1SluRAwCThs3i1ORCDkJqgy6G5yYqZIzrus4R1hTk%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=05%7C02%7CWolfgang.Bangerth% > 40colostate.edu > %7Cd23d2bd47fb74e7a13d008dd12d617b6%7Cafb58802ff7a4bb1ab21367ff2ecfc8b%7C0%7C0%7C638687432905970544%7CUnknown%7CTWFpbGZsb3d8eyJFbXB0eU1hcGkiOnRydWUsIlYiOiIwLjAuMDAwMCIsIlAiOiJXaW4zMiIsIkFOIjoiTWFpbCIsIldUIjoyfQ%3D%3D%7C0%7C%7C%7C&sdata=2qVckhsHD78WJek30A49bAR3RHjH4Zt1JpXJchZe8es%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+unsubscr...@googlegroups.com > > <mailto:dealii+unsubscr...@googlegroups.com>. > > To view this discussion visit https://groups.google.com/d/msgid/dealii/ > > b2e37661-b458-4dbb-a849-52c6f1744ec2n%40googlegroups.com <https:// > > nam10.safelinks.protection.outlook.com/? > > url=https%3A%2F%2Fgroups.google.com%2Fd%2Fmsgid%2Fdealii%2Fb2e37661- > > b458-4dbb- > > a849-52c6f1744ec2n%2540googlegroups.com > %3Futm_medium%3Demail%26utm_source%3Dfooter&data=05%7C02%7CWolfgang.Bangerth% > 40colostate.edu > %7Cd23d2bd47fb74e7a13d008dd12d617b6%7Cafb58802ff7a4bb1ab21367ff2ecfc8b%7C0%7C0%7C638687432905991612%7CUnknown%7CTWFpbGZsb3d8eyJFbXB0eU1hcGkiOnRydWUsIlYiOiIwLjAuMDAwMCIsIlAiOiJXaW4zMiIsIkFOIjoiTWFpbCIsIldUIjoyfQ%3D%3D%7C0%7C%7C%7C&sdata=bgq9j1DMj%2BiUfEMR5WCYiQKw0IYC9OkrnVjtLUdS3GU%3D&reserved=0>. > > -- > 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 a topic in the > Google Groups "deal.II User Group" group. > To unsubscribe from this topic, visit > https://groups.google.com/d/topic/dealii/LYNR9WAf-lI/unsubscribe. > To unsubscribe from this group and all its topics, send an email to > dealii+unsubscr...@googlegroups.com. > To view this discussion visit > https://groups.google.com/d/msgid/dealii/bd4632fd-eebb-4fc5-bc8c-6ed96086276d%40colostate.edu > . > -- 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 visit https://groups.google.com/d/msgid/dealii/CABQM8dR%3DMXCPUHAEdkPxnRdXfztgZ2Z6j6zCYoq9N7RJAxha4g%40mail.gmail.com.
