www.pudn.com > matlab7.x.rar > matsq_f.f

C=============================================================== 
C     matsq.f 
C 
C     Squares the input matrix 
C       
C     This is a MEX-file for MATLAB. 
C     Copyright (c) 1984-2000 The MathWorks, Inc. 
C     $Revision: 1.13 $ 
C=============================================================== 
 
 
 
C     The gateway routine 
      subroutine mexFunction(nlhs, plhs, nrhs, prhs) 
 
      integer mxGetM, mxGetN, mxGetPr 
      integer mxIsNumeric, mxCreateDoubleMatrix 
      integer plhs(*), prhs(*) 
      integer x_pr, y_pr 
      integer nlhs, nrhs 
      integer m, n, size 
      real*8  x(1000), y(1000) 
 
C     Check for proper number of arguments.  
      if(nrhs .ne. 1) then 
         call mexErrMsgTxt('One input required.') 
      elseif(nlhs .ne. 1) then 
         call mexErrMsgTxt('One output required.') 
      endif 
 
C     Get the size of the input array. 
      m = mxGetM(prhs(1)) 
      n = mxGetN(prhs(1)) 
      size = m*n 
 
C     Column * row should be smaller than 1000. 
      if(size.gt.1000) then 
         call mexErrMsgTxt('Row * column must be <= 1000.') 
      endif 
       
C     Check to ensure the array is numeric (not strings). 
      if(mxIsNumeric(prhs(1)) .eq. 0) then 
         call mexErrMsgTxt('Input must be a numeric array.') 
      endif 
 
C     Create matrix for the return argument. 
      plhs(1) = mxCreateDoubleMatrix(m, n, 0) 
      x_pr = mxGetPr(prhs(1)) 
      y_pr = mxGetPr(plhs(1)) 
 
C     Call the computational subroutine. 
      call matsq(%val(y_pr), %val(x_pr), m, n) 
 
      return 
      end 
 
 
C     Computational subroutine 
      subroutine matsq(y, x, m, n) 
      real*8 x(m,n), y(m,n) 
      integer m, n 
C      
      do 20 i=1,m 
         do 10 j=1,n 
            y(i,j)= x(i,j)**2 
 10      continue 
 20   continue 
      return 
      end