#include <gpc.h>

Inheritance diagram for GPC_LAPACK_pinv:

Collaboration diagram for GPC_LAPACK_pinv:

Public Member Functions
	GPC_LAPACK_pinv (integer n, integer m)

	~GPC_LAPACK_pinv (void)

integer	Inv (integer ndima, integer nrowa, integer ncola, doublereal *a)

Public Member Functions inherited from GPCInv
	GPCInv (void)

virtual	~GPCInv (void)

Protected Attributes
integer	m

integer	n

integer	iMin

integer	iMax

integer	iWork

doublereal *	pdS

doublereal *	pdU

doublereal *	pdVt

doublereal *	pdWork

Protected Attributes inherited from GPCInv
doublereal *	pdBase

Detailed Description

Definition at line 75 of file gpc.h.

Constructor & Destructor Documentation

GPC_LAPACK_pinv::GPC_LAPACK_pinv	(	integer	n,
		integer	m
	)

Definition at line 600 of file gpc.cc.

References ASSERT, iMax, iMin, iWork, n, GPCInv::pdBase, pdS, pdU, pdVt, pdWork, and SAFENEWARR.

 : GPCInv(), 
 m(m), n(n), 
 iMin(0), 
 iMax(0), 
 iWork(0),
 pdS(NULL), 
 pdU(NULL), 
 pdVt(NULL), 
 pdWork(NULL)
 {
         ASSERT(m > 0);
         ASSERT(n > 0);   
    
         iMin = std::min(m, n);  /* usati per vari parametri */
         iMax = std::max(m, n);       /* serve solo qui, non serve tenerlo */
    
         iWork = std::max(3*iMin+iMax, 5*iMin-4); /* aumentare per maggiore eff. */
    
         /* aree di lavoro per SVD:
          * s:    (min(m,n))
          * u:    (ldu*min(m,n))
          *       ldu >= m
          * vt:   (ldvt*n)
          *       ldvt >= min(m,n)
          * work: lwork
          *       lwork >= max(3*min(m,n)+max(m,n),5*min(m,n)-4)
          */
      
         integer i = iMin        /* s */
                 +m*iMin         /* u */
                 +iMin*n         /* vt */
 #ifndef OPTIMIZE_WORK
                 +iWork;         /* work */
 #else /* OPTIMIZE_WORK */
                 ;
 
         /* per ottimizzare le dimensioni di work: viene allocato a parte */
         SAFENEWARR(pdWork, doublereal, iWork);
 #endif /* OPTIMIZE_WORK */
      
         SAFENEWARR(pdBase, doublereal, i);
    
         pdS = pdBase;
         pdU = pdS+iMin;
         pdVt = pdU+m*iMin;
    
 #ifndef OPTIMIZE_WORK
         pdWork = pdVt+iMin*n;
 #endif /* !OPTIMIZE_WORK */
 }

GPC_LAPACK_pinv::~GPC_LAPACK_pinv ( void )

Definition at line 652 of file gpc.cc.

References NO_OP, pdWork, and SAFEDELETEARR.

 {
 #ifndef OPTIMIZE_WORK
         NO_OP;                  /* pdBase viene distrutto da GPCInv */
 #else /* OPTIMIZE_WORK */
         SAFEDELETEARR(pdWork);
 #endif /* OPTIMIZE_WORK */
 }

Member Function Documentation

integer GPC_LAPACK_pinv::Inv	(	integer	ndima,
		integer	nrowa,
		integer	ncola,
		doublereal *	a
	)

virtual

Implements GPCInv.

Definition at line 662 of file gpc.cc.

References ASSERT, gpc_pinv(), iMin, iWork, m, n, pdS, pdU, pdVt, pdWork, SAFEDELETEARR, and SAFENEWARR.

 {
         ASSERT(nrowa == m);
         ASSERT(ncola == n);
    
         integer info = gpc_pinv(ndima, nrowa, ncola, a,
                                 pdS, m, pdU, iMin, pdVt, iWork, pdWork);
    
         if (info != 0) {
                 return info;
         }
 
 #ifdef OPTIMIZE_WORK
         if (((integer *)pdWork)[0] > iWork) {
                 iWork = ((integer *)pdWork)[0];
                 SAFEDELETEARR(pdWork);
                 SAFENEWARR(pdWork, doublereal, iWork);
         }
 #endif /* OPTIMIZE_WORK */
    
         return integer(0);
 }