lit
    JINRLIB

DLSQQR -

E220

: S.Archer :

co popaoa acopa, popaa ao p popo cc m pa c n c oo a apao.
:

A * X = B , :

A(m*n) - aaa apa, cooca m pa c n ;
B(m*p) - p ;
X(n*p) - .

:

: - SUBROUTINE
: - DLSQQR
: - DDECOM, DPROC1, DPROD1, DPROD2, DSLV

:

CALL DLSQQR(A,X,B,M,N,IP,IM,IN,ERR), :

A(M,N) - (REAL*8) p acc, copa a o aa ap;
X(N,IP) - (REAL*8) p acc, copa a o ;
B(M,IP) - (REAL*8) p acc, copa a o IP pa pa ac;
M - (INTEGER) oco pa;
N - (INTEGER) oco c;
IP - (INTEGER) oco pa ac;
IM - (INTEGER) pa papoc acco A B opaop DIMENSION a popa, IM ≥ M;
IN - (INTEGER) pa papoc acc X opaop DIMENSION a popa, IN ≥ N;
ERR - oca pa. p opa popa oaa ERR=.TRUE. c p cc pa ao (pa ap n), o a o popa ERR=.FALSE.

:

  1. co, a. pao aopo a o.
    a apa. oca, "aocpo", 1976.

:

cc 6 pa c 2 c oo pao ac:
       . . .
       IMPLICIT REAL*8 (A-H,O-Z)
       DIMENSION A(6,2),B(6,1),X(2,1)
       LOGICAL ERR
       DATA A/7*1.0D0,2.0D0,3.0D0,4.0D0,5.0D0,6.0D0/,
      *     B/1.0D0,1.0D0,2.0D0,2.0D0,3.0D0,3.0D0/
       ERR=.TRUE.
       CALL DLSQQR(A,X,B,6,2,1,6,2,ERR)
       . . .
:
       X(1,1)=   .40D+00     X(2,1)=.4571428571D+00


home up e-mail