anchor.diff:              !!
               !! Note: \f$ \underline{\underline{R}}\underline{\underline{R}}_0 \f$ is used to go from the material basis into the inertial basis
               !!       and the transpose for the other direction.
  +!FIXME: note that p%uu0 is not the same as p%QuadPt (which is used in mesh orientation)  
            rot0_T0(:) = p%uu0(4:6,idx_qp,nelem)                ! initial rotation parameter set at gauss point at T=0
         
            CALL BD_CrvCompose(cc,rot_c_xi,rot0_T0,FLAG_R1R2)   ! cc = rot_c_xi composed with rot0_T0
  
positive.code:    m%elf(:,:,:)      = 0.0_BDKi
     m%elm(:,:,:,:)    = 0.0_BDKi
  
     DO idx_qp=1,p%nqp
  
  
  
negative.code:    ENDDO
  
  
  write(51,*) 'p%uu0'
     DO nelem = 1,p%elem_total
         DO idx_qp = 1,p%nqp
              !> ### Calculate the the initial displacement fields in an element
  