Skip to content

Commit

Permalink
Cleaning code from checks (print)
Browse files Browse the repository at this point in the history
  • Loading branch information
RBergua authored Nov 23, 2023
1 parent 4c32cb0 commit 104c521
Showing 1 changed file with 0 additions and 11 deletions.
11 changes: 0 additions & 11 deletions modules/subdyn/src/FEM.f90
Original file line number Diff line number Diff line change
Expand Up @@ -1308,24 +1308,13 @@ SUBROUTINE ElemK_Spring(k11, k12, k13, k14, k15, k16, k22, k23, k24, k25, k26, k
K(12, 5) = K(11,6)
K(12, 11) = K(11,12)

! Temporary check. Looking at the spring element matrix (local coordinate system).
print*,'Spring element stiffness (local coordinate system)'
print*, K

! Temporary check. Looking at direction cosine matrix.
print*,'Direction cosine',DirCos

DC = 0.0_FEKi
DC( 1: 3, 1: 3) = DirCos
DC( 4: 6, 4: 6) = DirCos
DC( 7: 9, 7: 9) = DirCos
DC(10:12, 10:12) = DirCos

K = MATMUL( MATMUL(DC, K), TRANSPOSE(DC) ) ! TODO: change me if DirCos convention is transposed

! Temporary check. Looking at the spring element matrix (global coordinate system).
print*,'Spring element stiffness (global coordinate system)'
print*, K

END SUBROUTINE ElemK_Spring
!------------------------------------------------------------------------------------------------------
Expand Down

0 comments on commit 104c521

Please sign in to comment.