changeoperator.m

student teacher and proferssor llove this project,this a very important exmple

   nif{2}=[nif{2},'         out=1.0',r];
   nif{2}=[nif{2},'        else',r];
   nif{2}=[nif{2},'         out=a',r];
   nif{2}=[nif{2},'         do i=1,int(b-1)',r];
   nif{2}=[nif{2},'          out=matmul(out,a)',r];
   nif{2}=[nif{2},'         enddo',r];
   nif{2}=[nif{2},'        endif',r];
   nif{2}=[nif{2},'       end function mxpow_r2r',r];
   nif{2}=[nif{2},'       function mxpow_r2i(a,b) result(out)',r];
   nif{2}=[nif{2},'        real, dimension(:,:) :: a',r];
   nif{2}=[nif{2},'        integer :: b',r];
   nif{2}=[nif{2},'        real, dimension(size(a,dim=1),size(a,dim=2)) :: out',r];
   nif{2}=[nif{2},'        integer a_m,ia_n, b_m, b_n, i, j',r];
   nif{2}=[nif{2},'        a_m=size(a,dim=1);a_n=size(a,dim=2);',r];
   nif{2}=[nif{2},'        if (b==0) then',r];
   nif{2}=[nif{2},'         out=1.0',r];
   nif{2}=[nif{2},'        else',r];
   nif{2}=[nif{2},'         out=a',r];
   nif{2}=[nif{2},'         do i=1,int(b-1)',r];
   nif{2}=[nif{2},'          out=matmul(out,a)',r];
   nif{2}=[nif{2},'         enddo',r];
   nif{2}=[nif{2},'        endif',r];
   nif{2}=[nif{2},'       end function mxpow_r2i',r];
   nif{2}=[nif{2},'       function mxpow_c2r(a,b) result(out)',r];
   nif{2}=[nif{2},'        complex, dimension(:,:) :: a',r];
   nif{2}=[nif{2},'        real :: b',r];
   nif{2}=[nif{2},'        complex, dimension(size(a,dim=1),size(a,dim=2)) :: out',r];
   nif{2}=[nif{2},'        integer a_m, a_n, b_m, b_n, i, j',r];
   nif{2}=[nif{2},'        a_m=size(a,dim=1);a_n=size(a,dim=2);',r];
   nif{2}=[nif{2},'        if (b==0.0) then',r];
   nif{2}=[nif{2},'         out=(1.0,0.0)',r];
   nif{2}=[nif{2},'        else',r];
   nif{2}=[nif{2},'         out=a',r];
   nif{2}=[nif{2},'         do i=1,int(b-1)',r];
   nif{2}=[nif{2},'          out=matmul(out,a)',r];
   nif{2}=[nif{2},'         enddo',r];
   nif{2}=[nif{2},'        endif',r];
   nif{2}=[nif{2},'       end function mxpow_c2r',r];
   nif{2}=[nif{2},'       function mxpow_c2i(a,b) result(out)',r];
   nif{2}=[nif{2},'        complex, dimension(:,:) :: a',r];
   nif{2}=[nif{2},'        integer :: b',r];
   nif{2}=[nif{2},'        complex, dimension(size(a,dim=1),size(a,dim=2)) :: out',r];
   nif{2}=[nif{2},'        integer a_m,ia_n, b_m, b_n, i, j',r];
   nif{2}=[nif{2},'        a_m=size(a,dim=1);a_n=size(a,dim=2);',r];
   nif{2}=[nif{2},'        if (b==0) then',r];
   nif{2}=[nif{2},'         out=(1.0,0.0)',r];
   nif{2}=[nif{2},'        else',r];
   nif{2}=[nif{2},'         out=a',r];
   nif{2}=[nif{2},'         do i=1,int(b-1)',r];
   nif{2}=[nif{2},'          out=matmul(out,a)',r];
   nif{2}=[nif{2},'         enddo',r];
   nif{2}=[nif{2},'        endif',r];
   nif{2}=[nif{2},'       end function mxpow_c2i',r];
   nif{2}=[nif{2},'       function mxpow_r2r2(a,b) result(out)',r];
   nif{2}=[nif{2},'        real, dimension(:,:) :: a',r];
   nif{2}=[nif{2},'        real, dimension(:,:) :: b',r];
   nif{2}=[nif{2},'        real, dimension(size(a,dim=1),size(a,dim=2)) :: out',r];
   nif{2}=[nif{2},'        integer a_m, a_n, b_m, b_n, i, j',r];
   nif{2}=[nif{2},'        a_m=size(a,dim=1);a_n=size(a,dim=2);',r];
   nif{2}=[nif{2},'        if (b(1,1)==0.0) then',r];
   nif{2}=[nif{2},'         out=1.0',r];
   nif{2}=[nif{2},'        else',r];
   nif{2}=[nif{2},'         out=a',r];
   nif{2}=[nif{2},'         do i=1,int(b(1,1)-1)',r];
   nif{2}=[nif{2},'          out=matmul(out,a)',r];
   nif{2}=[nif{2},'         enddo',r];
   nif{2}=[nif{2},'        endif',r];
   nif{2}=[nif{2},'       end function mxpow_r2r2',r];
   nif{2}=[nif{2},'       function mxpow_r2i2(a,b) result(out)',r];
   nif{2}=[nif{2},'        real, dimension(:,:) :: a',r];
   nif{2}=[nif{2},'        integer, dimension(:,:) :: b',r];
   nif{2}=[nif{2},'        real, dimension(size(a,dim=1),size(a,dim=2)) :: out',r];
   nif{2}=[nif{2},'        integer a_m,ia_n, b_m, b_n, i, j',r];
   nif{2}=[nif{2},'        a_m=size(a,dim=1);a_n=size(a,dim=2);',r];
   nif{2}=[nif{2},'        if (b(1,1)==0) then',r];
   nif{2}=[nif{2},'         out=1.0',r];
   nif{2}=[nif{2},'        else',r];
   nif{2}=[nif{2},'         out=a',r];
   nif{2}=[nif{2},'         do i=1,int(b(1,1)-1)',r];
   nif{2}=[nif{2},'          out=matmul(out,a)',r];
   nif{2}=[nif{2},'         enddo',r];
   nif{2}=[nif{2},'        endif',r];
   nif{2}=[nif{2},'       end function mxpow_r2i2',r];
   nif{2}=[nif{2},'       function mxpow_c2r2(a,b) result(out)',r];
   nif{2}=[nif{2},'        complex, dimension(:,:) :: a',r];
   nif{2}=[nif{2},'        real, dimension(:,:) :: b',r];
   nif{2}=[nif{2},'        complex, dimension(size(a,dim=1),size(a,dim=2)) :: out',r];
   nif{2}=[nif{2},'        integer a_m, a_n, b_m, b_n, i, j',r];
   nif{2}=[nif{2},'        a_m=size(a,dim=1);a_n=size(a,dim=2);',r];
   nif{2}=[nif{2},'        if (b(1,1)==0.0) then',r];
   nif{2}=[nif{2},'         out=(1.0,0.0)',r];
   nif{2}=[nif{2},'        else',r];
   nif{2}=[nif{2},'         out=a',r];
   nif{2}=[nif{2},'         do i=1,int(b(1,1)-1)',r];
   nif{2}=[nif{2},'          out=matmul(out,a)',r];
   nif{2}=[nif{2},'         enddo',r];
   nif{2}=[nif{2},'        endif',r];
   nif{2}=[nif{2},'       end function mxpow_c2r2',r];
   nif{2}=[nif{2},'       function mxpow_c2i2(a,b) result(out)',r];
   nif{2}=[nif{2},'        complex, dimension(:,:) :: a',r];
   nif{2}=[nif{2},'        integer, dimension(:,:) :: b',r];
   nif{2}=[nif{2},'        complex, dimension(size(a,dim=1),size(a,dim=2)) :: out',r];
   nif{2}=[nif{2},'        integer a_m,ia_n, b_m, b_n, i, j',r];
   nif{2}=[nif{2},'        a_m=size(a,dim=1);a_n=size(a,dim=2);',r];
   nif{2}=[nif{2},'        if (b(1,1)==0) then',r];
   nif{2}=[nif{2},'         out=(1.0,0.0)',r];
   nif{2}=[nif{2},'        else',r];
   nif{2}=[nif{2},'         out=a',r];
   nif{2}=[nif{2},'         do i=1,int(b(1,1)-1)',r];
   nif{2}=[nif{2},'          out=matmul(out,a)',r];
   nif{2}=[nif{2},'         enddo',r];
   nif{2}=[nif{2},'        endif',r];
   nif{2}=[nif{2},'       end function mxpow_c2i2',r];
   nif{2}=[nif{2},'       function mxpow_i2r(a,b) result(out)',r];
   nif{2}=[nif{2},'        integer, dimension(:,:) :: a',r];
   nif{2}=[nif{2},'        real :: b',r];
   nif{2}=[nif{2},'        real, dimension(size(a,dim=1),size(a,dim=2)) :: out',r];
   nif{2}=[nif{2},'        integer a_m, a_n, b_m, b_n, i, j',r];
   nif{2}=[nif{2},'        a_m=size(a,dim=1);a_n=size(a,dim=2);',r];
   nif{2}=[nif{2},'        if (b==0.0) then',r];
   nif{2}=[nif{2},'         out=1.0',r];
   nif{2}=[nif{2},'        else',r];
   nif{2}=[nif{2},'         out=a',r];
   nif{2}=[nif{2},'         do i=1,int(b-1)',r];
   nif{2}=[nif{2},'          out=matmul(out,a)',r];
   nif{2}=[nif{2},'         enddo',r];
   nif{2}=[nif{2},'        endif',r];
   nif{2}=[nif{2},'       end function mxpow_i2r',r];
   nif{2}=[nif{2},'       function mxpow_i2i(a,b) result(out)',r];
   nif{2}=[nif{2},'        integer, dimension(:,:) :: a',r];
   nif{2}=[nif{2},'        integer :: b',r];
   nif{2}=[nif{2},'        integer, dimension(size(a,dim=1),size(a,dim=2)) :: out',r];
   nif{2}=[nif{2},'        integer a_m,ia_n, b_m, b_n, i, j',r];
   nif{2}=[nif{2},'        a_m=size(a,dim=1);a_n=size(a,dim=2);',r];
   nif{2}=[nif{2},'        if (b==0) then',r];
   nif{2}=[nif{2},'         out=1',r];
   nif{2}=[nif{2},'        else',r];
   nif{2}=[nif{2},'         out=a',r];
   nif{2}=[nif{2},'         do i=1,int(b-1)',r];
   nif{2}=[nif{2},'          out=matmul(out,a)',r];
   nif{2}=[nif{2},'         enddo',r];
   nif{2}=[nif{2},'        endif',r];
   nif{2}=[nif{2},'       end function mxpow_i2i',r];
   nif{2}=[nif{2},'       function mxpow_i2r2(a,b) result(out)',r];
   nif{2}=[nif{2},'        integer, dimension(:,:) :: a',r];
   nif{2}=[nif{2},'        real, dimension(:,:) :: b',r];
   nif{2}=[nif{2},'        real, dimension(size(a,dim=1),size(a,dim=2)) :: out',r];
   nif{2}=[nif{2},'        integer a_m, a_n, b_m, b_n, i, j',r];
   nif{2}=[nif{2},'        a_m=size(a,dim=1);a_n=size(a,dim=2);',r];
   nif{2}=[nif{2},'        if (b(1,1)==0.0) then',r];
   nif{2}=[nif{2},'         out=1.0',r];
   nif{2}=[nif{2},'        else',r];
   nif{2}=[nif{2},'         out=a',r];
   nif{2}=[nif{2},'         do i=1,int(b(1,1)-1)',r];
   nif{2}=[nif{2},'          out=matmul(out,a)',r];
   nif{2}=[nif{2},'         enddo',r];
   nif{2}=[nif{2},'        endif',r];
   nif{2}=[nif{2},'       end function mxpow_i2r2',r];
   nif{2}=[nif{2},'       function mxpow_i2i2(a,b) result(out)',r];
   nif{2}=[nif{2},'        integer, dimension(:,:) :: a',r];
   nif{2}=[nif{2},'        integer, dimension(:,:) :: b',r];
   nif{2}=[nif{2},'        integer, dimension(size(a,dim=1),size(a,dim=2)) :: out',r];
   nif{2}=[nif{2},'        integer a_m,ia_n, b_m, b_n, i, j',r];
   nif{2}=[nif{2},'        a_m=size(a,dim=1);a_n=size(a,dim=2);',r];
   nif{2}=[nif{2},'        if (b(1,1)==0) then',r];
   nif{2}=[nif{2},'         out=1',r];
   nif{2}=[nif{2},'        else',r];
   nif{2}=[nif{2},'         out=a',r];
   nif{2}=[nif{2},'         do i=1,int(b(1,1)-1)',r];
   nif{2}=[nif{2},'          out=matmul(out,a)',r];
   nif{2}=[nif{2},'         enddo',r];
   nif{2}=[nif{2},'        endif',r];
   nif{2}=[nif{2},'       end function mxpow_i2i2',r];
  case '.^'
   %Never needed!
  case '*'
   %Already taken care of by matmul
  case '/'
   %Never needed!
  case '.*'
   %Already correct
  case './'
   %Already correct
  case '+'
   %Never needed!
  case '-'
   %Never needed!
  case '\'
   %Should we do this?
 end
end