Actual source code: fnorm.F

  1: !
  2: !  "$Id: fnorm.F,v 1.9 2001/08/07 03:05:24 balay 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:       PetscScalar      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