constrain_element_matrix changes the col_dofs vector size

[rebpoli]

I call

dof_map.constrain_element_matrix(matrix, row_dofs, col_dofs, false)

The col_dofs is modified and the size increases by 1.

This is crashing the code, because I use the col_dofs more than once.

Is that supposed to be?

Renato

[jwpeterson]

It is indeed possible for the size of col_dofs to be changed by a call to

dof_map.constrain_element_matrix(matrix, row_dofs, col_dofs, false)

This is because internally that function calls DofMap::build_constraint_matrix(), passing in col_dofs as an argument,
and within that function we have the following:

  // If we added any DOFS then we need to do this recursively.
  // It is possible that we just added a DOF that is also
  // constrained!
  //
  // Also, we need to handle the special case of an element having DOFs
  // constrained in terms of other, local DOFs
  
  ...
  
      // Add new dependency dofs to the end of the current dof set
      elem_dofs.insert(elem_dofs.end(),
                       dof_set.begin(), dof_set.end());
  

[rebpoli]

Got it. Thank you. Em seg., 29 de jan. de 2024 às 06:39, John W. Peterson < ***@***.***> escreveu:

…

It is indeed possible for the size of col_dofs to be changed by a call to dof_map.constrain_element_matrix(matrix, row_dofs, col_dofs, false) This is because internally that function calls DofMap::build_constraint_matrix(), passing in col_dofs as an argument, and within that function we have the following: // If we added any DOFS then we need to do this recursively. // It is possible that we just added a DOF that is also // constrained! // // Also, we need to handle the special case of an element having DOFs // constrained in terms of other, local DOFs ... // Add new dependency dofs to the end of the current dof set elem_dofs.insert(elem_dofs.end(), dof_set.begin(), dof_set.end()); — Reply to this email directly, view it on GitHub <#3773 (comment)>, or unsubscribe <https://github.com/notifications/unsubscribe-auth/AMOVERYASB4XJC4PAEMHVVTYQ6KB5AVCNFSM6AAAAABCNYPIJGVHI2DSMVQWIX3LMV43OSLTON2WKQ3PNVWWK3TUHMYTSMJUGYYTMMJWGI> . You are receiving this because you authored the thread.Message ID: ***@***.***>

[jwpeterson]

Closing because I don't think this is really an Issue so much as a clarification.