anchor.diff:              j = j+3
            end do
                     
  -         print *, 'Translational Velocity: ', n1, TwoNorm( LinVec_2 - LinVec_2_theta )
  +         print *, 'Translational Velocity: ', n1, TwoNorm( LinVec_2 - LinVec_2_approx )
                                 
         end do      
         
  
positive.code:                ! Mesh1_O%TranslationDisp = theta^S|_op + delta theta^S
              Mesh1_O%TranslationDisp(:,i) = Mesh1_O_op%TranslationDisp(:,i) + LinVec_1(j:j+2)
              
                 ! Mesh1_O%Orientation = theta^S|_op + delta theta^S
              call SmllRotTrans( 'orientation', LinVec_1_theta(j  )+LinVec_1_a(j  ) &
                                              , LinVec_1_theta(j+1)+LinVec_1_a(j+1) &
                                              , LinVec_1_theta(j+2)+LinVec_1_a(j+2) & 
                                              , Mesh1_O%Orientation(:,:,i), ErrStat=ErrStat, ErrMsg=ErrMsg)            
                                      
              j = j+3            
              !IF (ErrStat /= ErrID_None) CALL WrScr("*******"//TRIM(ErrMsg))
                          
  
negative.code:       do n1=1,15
           
           ! perturb x^S:
           call getRandomVector(LinVec_1)
           call getRandomVector(LinVec_1_a)
           call getRandomVector(LinVec_1_b)
           
           LinVec_2_approx = matmul( map_mod1_mod2%dm%tv, LinVec_1 ) + matmul( map_mod1_mod2%dm%mi, LinVec_1_a ) &
                          + matmul( map_mod1_mod2%dm%fx_p, LinVec_1_b )! approximate delta theta^D         
                        
           j=1
           do i=1,Mesh1_O%nnodes
              
                 ! Mesh1_O%Orientation = theta^S|_op + delta theta^S
              call SmllRotTrans( 'orientation', LinVec_1(j  ) &
                                              , LinVec_1(j+1) &
                                              , LinVec_1(j+2) & 
                                              , Orientation, ErrStat=ErrStat, ErrMsg=ErrMsg)            
              Mesh1_O%Orientation(:,:,i) = matmul( Mesh1_O_op%Orientation(:,:,i), Orientation )
  
                          
                 ! Mesh1_O%TranslationVel = theta^S|_op + delta theta^S
  