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
