------- Additional Comments From bangerth at dealii dot org  2004-12-13 21:10 
-------
In case someone wondered about the mapping of the rank-3 array to the rank-2 
array in my previous testcase: here is something even simpler: 
========================== 
      subroutine choleskyzhp ()  
      integer    i,j,k  
      complex*16 a(500),x  
  
      do i = 1,5  
         do j = 1,5  
            do k = 1,j 
               x = a(i*5+j)*dconjg(a(j*5+j)) 
            enddo  
            a(i*5+k) = x  
         enddo  
      enddo  
  
      return  
      end  
  
  
      program test  
      call choleskyzhp 
      return  
      end  
============================ 
It fails in the same way as my previous testcase, i.e. in an movapd 
instruction. 
 
W. 

-- 


http://gcc.gnu.org/bugzilla/show_bug.cgi?id=14776

Reply via email to