5 Apr
2016
5 Apr
'16
3 p.m.
Hi, I think I find a bug in function Nektar::MultiRegions::AssemblyMapCG::v_Assemble(const Array<OneD, const NekDouble> &loc, Array<OneD, NekDouble> &global) const, in file AssemblyMapCG.cpp. The fourth line of this function writes, "local = Array<OneD, NekDouble>(local.num_elements(),local.data());", which aims to make a copy of "loc" if "loc" and "global" have the same address. But this line makes a copy of "local", which hasn't been initialized yet. So this line should be substituted by "local = Array<OneD, NekDouble>(loc.num_elements(),loc.data());" Best regards, Ankang