eNorm Function

private function eNorm(n, x)

Computes the Euclidean norm of a vector.

The Euclidean norm is computed by accumulating the sum of squares in three different sums. The sums of squares for the small and large components are scaled so that no overflows occur. Non-destructive underflows are permitted. Underflows and overflows do not occur in the computation of the unscaled sum of squares for the intermediate components.

The definitions of small, intermediate and large components depend on two constants, RDWARF and RGIANT. The main restrictions on these constants are that RDWARF^2 not underflow and RGIANT^2 not overflow.

Arguments

Type IntentOptional Attributes Name
integer :: n

Length of the vector
real(kind=dp) :: x(n)

Vector whose norm is desired

Return Value real(kind=dp)

Euclidean norm of the vector


Called by

proc~~enorm~~CalledByGraph proc~enorm eNorm proc~dogleg dogleg proc~dogleg->proc~enorm proc~hybrd hybrd proc~hybrd->proc~enorm proc~hybrd->proc~dogleg proc~qrfac qrFac proc~hybrd->proc~qrfac proc~qrfac->proc~enorm proc~fsolve fsolve proc~fsolve->proc~hybrd proc~solve_junction solve_junction proc~solve_junction->proc~fsolve proc~junction_resolution_from_and_to_ports Junction_resolution_from_and_to_ports proc~junction_resolution_from_and_to_ports->proc~solve_junction proc~main_loop main_loop proc~main_loop->proc~junction_resolution_from_and_to_ports program~reims_p reims_p program~reims_p->proc~main_loop