*>
*> \param[in,out] AFP
*> \verbatim
-*> AFP is DOUBLE PRECISION array, dimension
-*> (N*(N+1)/2)
+*> AFP is DOUBLE PRECISION array, dimension (N*(N+1)/2)
*> If FACT = 'F', then AFP is an input argument and on entry
*> contains the triangular factor U or L from the Cholesky
*> factorization A = U**T*U or A = L*L**T, in the same storage
*>
*> \param[in,out] AP
*> \verbatim
-*> AP is DOUBLE PRECISION array, dimension
-*> (N*(N+1)/2)
+*> AP is DOUBLE PRECISION array, dimension (N*(N+1)/2)
*> On entry, the upper or lower triangle of the symmetric matrix
*> A, packed columnwise in a linear array. The j-th column of A
*> is stored in the array AP as follows:
*>
*> \param[in,out] AFP
*> \verbatim
-*> AFP is DOUBLE PRECISION array, dimension
-*> (N*(N+1)/2)
+*> AFP is DOUBLE PRECISION array, dimension (N*(N+1)/2)
*> If FACT = 'F', then AFP is an input argument and on entry
*> contains the block diagonal matrix D and the multipliers used
*> to obtain the factor U or L from the factorization
*>
*> \param[in,out] AFP
*> \verbatim
-*> AFP is REAL array, dimension
-*> (N*(N+1)/2)
+*> AFP is REAL array, dimension (N*(N+1)/2)
*> If FACT = 'F', then AFP is an input argument and on entry
*> contains the triangular factor U or L from the Cholesky
*> factorization A = U**T*U or A = L*L**T, in the same storage
*>
*> \param[in,out] AP
*> \verbatim
-*> AP is REAL array, dimension
-*> (N*(N+1)/2)
+*> AP is REAL array, dimension (N*(N+1)/2)
*> On entry, the upper or lower triangle of the symmetric matrix
*> A, packed columnwise in a linear array. The j-th column of A
*> is stored in the array AP as follows:
*>
*> \param[in,out] AFP
*> \verbatim
-*> AFP is REAL array, dimension
-*> (N*(N+1)/2)
+*> AFP is REAL array, dimension (N*(N+1)/2)
*> If FACT = 'F', then AFP is an input argument and on entry
*> contains the block diagonal matrix D and the multipliers used
*> to obtain the factor U or L from the factorization