Actual source code: fnorm.F

  1: !
  2: !  "$Id: fnorm.F,v 1.8 2000/01/11 21:04:00 bsmith Exp $";
  3: !
  4: !    Fortran kernel for the Norm() vector routine
  5: !
 6:  #include include/finclude/petscdef.h
  7: !
  8:       subroutine FortranNormSqr(x,n,sum1)
  9:       implicit none
 10:       Scalar           x(*)
 11:       double precision sum1
 12:       integer n

 14:       integer i

 16:       do 10,i=1,n
 17:         sum1 = sum1 + PetscRealPart(x(i)*PetscConj(x(i)))
 18:  10   continue

 20:       return
 21:       end