NAG FL Interface
f04lhf (real_​blkdiag_​fac_​solve)

1 Purpose

2 Specification

3 Description

4 References

5 Arguments

1: $\mathbf{trans}$ – Character(1) Input

On entry: specifies the equations to be solved.

$\mathbf{trans}=\text{'N'}$

Solve $AX=B$.

$\mathbf{trans}=\text{'T'}$

Solve $A^{\mathrm{T}}X=B$.

Constraint: $\mathbf{trans}=\text{'N'}$ or $\text{'T'}$.

2: $\mathbf{n}$ – Integer Input

On entry: $n$, the order of the matrix $A$.

Constraint: $\mathbf{n}>0$.

3: $\mathbf{nbloks}$ – Integer Input

On entry: the total number of blocks of the matrix $A$, as supplied to f04lhf.

Constraint: $0<\mathbf{nbloks}\le \mathbf{n}$.

4: $\mathbf{blkstr}(3,\mathbf{nbloks})$ – Integer array Input

On entry: information which describes the block structure of $A$, as supplied to f04lhf.

5: $\mathbf{a}\left(\mathbf{lena}\right)$ – Real (Kind=nag_wp) array Input

On entry: the elements in the factorization of $A$, as returned by f04lhf.

6: $\mathbf{lena}$ – Integer Input

On entry: the dimension of the array a as declared in the (sub)program from which f04lhf is called.

Constraint: $\mathbf{lena}\ge {\displaystyle \sum _{k=1}^{\mathbf{nbloks}}}\mathbf{blkstr}(1,k)\times \mathbf{blkstr}(2,k)$.

7: $\mathbf{pivot}\left(\mathbf{n}\right)$ – Integer array Input

On entry: details of the interchanges in the factorization, as returned by f04lhf.

8: $\mathbf{b}(\mathbf{ldb},\mathbf{ir})$ – Real (Kind=nag_wp) array Input/Output

On entry: the $n\times r$ right-hand side matrix $B$.

On exit: b is overwritten by the $n\times r$ solution matrix $X$.

9: $\mathbf{ldb}$ – Integer Input

On entry: the first dimension of the array b as declared in the (sub)program from which f04lhf is called.

Constraint: $\mathbf{ldb}\ge \mathbf{n}$.

10: $\mathbf{ir}$ – Integer Input

On entry: $r$, the number of right-hand sides.

Constraint: $\mathbf{ir}>0$.

11: $\mathbf{ifail}$ – Integer Input/Output

On entry: ifail must be set to $0$, $-1$ or $1$ to set behaviour on detection of an error; these values have no effect when no error is detected.

A value of $0$ causes the printing of an error message and program execution will be halted; otherwise program execution continues. A value of $-1$ means that an error message is printed while a value of $1$ means that it is not.

If halting is not appropriate, the value $-1$ or $1$ is recommended. If message printing is undesirable, then the value $1$ is recommended. Otherwise, the value $0$ is recommended. When the value $-\mathbf{1}$ or $\mathbf{1}$ is used it is essential to test the value of ifail on exit.

On exit: $\mathbf{ifail}=\mathbf{0}$ unless the routine detects an error or a warning has been flagged (see Section 6).

6 Error Indicators and Warnings

$\mathbf{ifail}=1$

On entry, $\mathbf{ir}=⟨\mathit{\text{value}}⟩$.
Constraint: $\mathbf{ir}\ge 1$.

On entry, $K=⟨\mathit{\text{value}}⟩$, $\mathbf{blkstr}(2,K)=⟨\mathit{\text{value}}⟩$, $\mathbf{blkstr}(3,K)=⟨\mathit{\text{value}}⟩$ and $\mathbf{blkstr}(1,K)=⟨\mathit{\text{value}}⟩$.
Constraint: $\mathbf{blkstr}(2,K)-\mathbf{blkstr}(3,K)\le \mathbf{blkstr}(1,K)$.

On entry, $K=⟨\mathit{\text{value}}⟩$, $\mathbf{blkstr}(2,K)=⟨\mathit{\text{value}}⟩$ and $\mathbf{blkstr}(1,K)=⟨\mathit{\text{value}}⟩$.
Constraint: $\mathbf{blkstr}(2,K)\ge \mathbf{blkstr}(1,K)$.

On entry, $K=⟨\mathit{\text{value}}⟩$, $\mathbf{blkstr}(3,K)=⟨\mathit{\text{value}}⟩$, $\mathbf{blkstr}(3,K-1)=⟨\mathit{\text{value}}⟩$ and $\mathbf{blkstr}(2,K)=⟨\mathit{\text{value}}⟩$.
Constraint: $\mathbf{blkstr}(3,K)+\mathbf{blkstr}(3,K-1)\ge \mathbf{blkstr}(2,K)$.

On entry, $K=⟨\mathit{\text{value}}⟩$ and $\mathbf{blkstr}(1,K)=⟨\mathit{\text{value}}⟩$.
Constraint: $\mathbf{blkstr}(1,K)\ge 1$.

On entry, $K=⟨\mathit{\text{value}}⟩$ and $\mathbf{blkstr}(2,K)=⟨\mathit{\text{value}}⟩$.
Constraint: $\mathbf{blkstr}(2,K)\ge 1$.

On entry, $K=⟨\mathit{\text{value}}⟩$ and $\mathbf{blkstr}(3,K)=⟨\mathit{\text{value}}⟩$.
Constraint: $\mathbf{blkstr}(3,K)\ge 0$.

On entry, $\mathbf{ldb}=⟨\mathit{\text{value}}⟩$ and $\mathbf{n}=⟨\mathit{\text{value}}⟩$.
Constraint: $\mathbf{ldb}\ge \mathbf{n}$.

On entry, lena is too small. $\mathbf{lena}=⟨\mathit{\text{value}}⟩$. Minimum possible dimension: $⟨\mathit{\text{value}}⟩$.

On entry, $\mathbf{n}=⟨\mathit{\text{value}}⟩$.
Constraint: $\mathbf{n}\ge 1$.

On entry, $\mathbf{n}=⟨\mathit{\text{value}}⟩$ and $\mathbf{nbloks}=⟨\mathit{\text{value}}⟩$.
Constraint: $\mathbf{n}\ge \mathbf{nbloks}$.

On entry, $\mathbf{nbloks}=⟨\mathit{\text{value}}⟩$.
Constraint: $\mathbf{nbloks}\ge 1$.

On entry, the following equality does not hold: $\mathbf{blkstr}(2,1)+\mathrm{sum}(\mathbf{blkstr}(2,k)-\mathbf{blkstr}(3,k-1):k=2,\mathbf{nbloks})=\mathbf{n}$.

On entry, the following equality does not hold: $\mathrm{sum}(\mathbf{blkstr}(1,k):k=1,\mathbf{nbloks})=\mathbf{n}$.

On entry, the following inequality was not satisfied for: $J=⟨\mathit{\text{value}}⟩$. $\mathrm{sum}(\mathbf{blkstr}(2,k)-\mathbf{blkstr}(3,k):k=1,J)\le$ $\mathrm{sum}(\mathbf{blkstr}(1,k):k=1,J)\le$ $\mathbf{blkstr}(2,1)+\mathrm{sum}(\mathbf{blkstr}(2,k)-\mathbf{blkstr}(3,k-1):k=2,J)$.

On entry, $\mathbf{trans}\ne \text{'N'}$ or $\text{'T'}$: $\mathbf{trans}=⟨\mathit{\text{value}}⟩$.

$\mathbf{ifail}=-99$

An unexpected error has been triggered by this routine. Please contact NAG.

See Section 7 in the Introduction to the NAG Library FL Interface for further information.

$\mathbf{ifail}=-399$

Your licence key may have expired or may not have been installed correctly.

See Section 8 in the Introduction to the NAG Library FL Interface for further information.

$\mathbf{ifail}=-999$

Dynamic memory allocation failed.

See Section 9 in the Introduction to the NAG Library FL Interface for further information.

7 Accuracy

8 Parallelism and Performance

9 Further Comments

10 Example

10.1 Program Text

10.2 Program Data

10.3 Program Results