Commit dd093a0e authored by Giovanni La Mura's avatar Giovanni La Mura
Browse files

Fix lurking bug in AT matrix copy constructor and MPI constructor

parent 699591ef

src/libnptm/Commons.cpp

+2 −2
Original line number Diff line number Diff line
@@ -1646,7 +1646,7 @@ ParticleDescriptorInclusion::ParticleDescriptorInclusion(const ParticleDescripto 
  vec_at = new dcomplex[_nlemt * _ndm];

  for (int vi = 0; vi < _nlemt * _ndm; vi++) vec_at[vi] = rhs.vec_at[vi];

  at = new dcomplex*[_nlemt];

  for (int ai = 0; ai < _nlemt; ai++) at[ai] = vec_at + (ai * _nlemt);

  for (int ai = 0; ai < _nlemt; ai++) at[ai] = vec_at + (ai * _ndm);

}



#ifdef MPI_VERSION
@@ -1727,7 +1727,7 @@ ParticleDescriptorInclusion::ParticleDescriptorInclusion(const mixMPI *mpidata) 
  vec_at = new dcomplex[_nlemt * _ndm];

  MPI_Bcast(vec_at, _nlemt * _ndm, MPI_C_DOUBLE_COMPLEX, 0, MPI_COMM_WORLD);

  at = new dcomplex*[_nlemt];

  for (int ai = 0; ai < _nlemt; ai++) at[ai] = vec_at + (ai * _nlemt);

  for (int ai = 0; ai < _nlemt; ai++) at[ai] = vec_at + (ai * _ndm);

}

#endif // MPI_VERSION

// >>> End of ParticleDescriptorInclusion class implementation. <<< //