anchor.diff:        do n1=1,15
            
            ! perturb x^S:
  -         call getRandomVector(LinVec_1)
  -         
  -         LinVec_2_approx = matmul( map_mod1_mod2%dm%mi, LinVec_1 )  ! approximate delta theta^D         
  +         call random_number( LinVec_1 )
  +         LinVec_1 = (2*LinVec_1-1)/real(n1**3)  !note n1 from the loop this is called in
  +         
  +         LinVec_2_a = matmul( map_mod1_mod2%dm%mi, LinVec_1 )  ! approximate delta theta^D         
                         
            j=1
            do i=1,Mesh1_O%nnodes
  
positive.code:             j = j+3
           end do
                    
           print *, 'Translational Displacement: ', n1, TwoNorm( LinVec_2 - LinVec_2_approx )
                                
        end do      
        
  
negative.code:             j = j+3
           end do
                    
           print *, 'Rotational Velocity: ', n1, TwoNorm( LinVec_2 - LinVec_2_a )                              
        end do            
        
        
  