Hi All,
I figured it out. I was not calling localize from all the processor and
that could have been the problem. Nevermind the question.
Regards,
Prashant
On Sat, Jul 11, 2020 at 9:02 PM Prashant K. Jha <pjha....@gmail.com> wrote:
> Correction: For approach 1 and 2, I get error
>
> Assertion `this->verify(r.size())' failed.
>
> On Sat, Jul 11, 2020 at 8:55 PM Prashant K. Jha <pjha....@gmail.com>
> wrote:
>
> > Hi All,
> >
> > I have been struggling with a very simple problem. I create a system as
> > follows:
> >
> > EquationSystems es(mesh);
> > auto &taf_sys = es.add_system<TransientLinearImplicitSystem>("TAF");
> > taf_sys.add_variable("taf", FIRST);
> >
> > After I solve taf_sys, I want to localize the solution of taf_sys system
> > (ie "taf" field) on all processors in the taf_localize vector. I have
> tried
> > many approaches such as
> >
> > 1.
> > std::vector<double> taf_localize;
> > taf_sys.solution->localize(taf_localize);
> >
> > 2.
> > std::vector<double> taf_localize;
> > taf_sys.current_local_solution->localize(taf_localize);
> >
> > 3.
> > static int init = -1;
> > static std::unique_ptr<NumericVector<Number>> taf_localize =
> > NumericVector<Number>::build(comm());
> > if (init < 0) {
> > taf_localize->init (taf_sys.solution->size(), false, SERIAL);
> > init = 0;
> > }
> >
> > taf_sys.solution->localize(*taf_localize);
> >
> > They all fail with
> >
> > Assertion `(this->comm()).verify(std::to_string(642))' failed
> >
> > Can anyone help with this? Thank you.
> >
> > Regards,
> > Prashant
> >
> >
>
> _______________________________________________
> Libmesh-users mailing list
> Libmesh-users@lists.sourceforge.net
> https://lists.sourceforge.net/lists/listinfo/libmesh-users
>
_______________________________________________
Libmesh-users mailing list
Libmesh-users@lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/libmesh-users
