static void prune_mat(gmx::ArrayRef<t_matrix> mat, gmx::ArrayRef<t_matrix> mat2, int skip)
{
- GMX_RELEASE_ASSERT(mat.size() == mat2.size(),
+ GMX_RELEASE_ASSERT(mat.size() == mat2.size() || mat2.empty(),
"Matrix pruning requires matrices of the same size");
for (gmx::index i = 0; i != gmx::ssize(mat); ++i)
{
int skip,
int mapoffset)
{
- GMX_RELEASE_ASSERT(mat.size() == mat2.size(),
+ GMX_RELEASE_ASSERT(mat.size() == mat2.size() || mat2.empty(),
"Combined matrix write requires matrices of the same size");
if (!mat2.empty())
{