Next: ../STGIRONS/oo_1lapack++

Up: Presentation

Previous: ../STGIRONS/oo_sml

A Simple Example

F77 LAPACK / BLAS

        DCOPY(N, AXMB, 1, B, 1)			
        DGEMV('No Transpose', N, N, 1.0, A, N, X, 1, -1.0, AXMB, 1) 
        D0 100 I=1,N 
100     R(I) = DASUM(N, A(I), N) 
        INDEX = IDAMAX(N, R, 1) 
        RES = DASUM(N,AXMB,1) / ( N* R(INDEX) * DASUM(N,X) 
      &         * DLMACH('e'))

LAPACK++

   res =           Norm_Inf(A*x-b) / 
          (N * Norm_Inf(A) * Norm_Inf(x) * Mach_deps() );