CPPLapack
 All Classes Files Functions Variables Friends Pages
Functions
_dsymatrix-dsymatrix.hpp File Reference

Go to the source code of this file.

Functions

_dsymatrix operator+ (const _dsymatrix &matA, const dsymatrix &matB)
 
_dsymatrix operator- (const _dsymatrix &matA, const dsymatrix &matB)
 
_dgematrix operator* (const _dsymatrix &matA, const dsymatrix &matB)
 

Function Documentation

_dsymatrix operator+ ( const _dsymatrix matA,
const dsymatrix matB 
)
inline

_dsymatrix+dsymatrix operator

Definition at line 3 of file _dsymatrix-dsymatrix.hpp.

References dsymatrix::darray, _dsymatrix::darray, i(), dsymatrix::n, and _dsymatrix::n.

4 {CPPL_VERBOSE_REPORT;
5 #ifdef CPPL_DEBUG
6  if(matA.n!=matB.n){
7  ERROR_REPORT;
8  std::cerr << "These two matrises can not make a summation." << std::endl
9  << "Your input was (" << matA.n << "x" << matA.n << ") + (" << matB.n << "x" << matB.n << ")." << std::endl;
10  exit(1);
11  }
12 #endif//CPPL_DEBUG
13 
14  for(CPPL_INT j=0; j<matA.n; j++){
15  for(CPPL_INT i=j; i<matA.n; i++){
16  matA.darray[j][i] += matB.darray[j][i];
17  }
18  }
19 
20  return matA;
21 }
_dgematrix i(const _dgbmatrix &mat)
double ** darray
array of pointers of column head addresses
Definition: _dsymatrix.hpp:13
CPPL_INT n
matrix column size
Definition: dsymatrix.hpp:10
CPPL_INT n
matrix column size
Definition: _dsymatrix.hpp:11
double ** darray
array of pointers of column head addresses
Definition: dsymatrix.hpp:12
_dsymatrix operator- ( const _dsymatrix matA,
const dsymatrix matB 
)
inline

_dsymatrix-dsymatrix operator

Definition at line 25 of file _dsymatrix-dsymatrix.hpp.

References dsymatrix::darray, _dsymatrix::darray, i(), dsymatrix::n, and _dsymatrix::n.

26 {CPPL_VERBOSE_REPORT;
27 #ifdef CPPL_DEBUG
28  if(matA.n!=matB.n){
29  ERROR_REPORT;
30  std::cerr << "These two matrises can not make a subtraction." << std::endl
31  << "Your input was (" << matA.n << "x" << matA.n << ") - (" << matB.n << "x" << matB.n << ")." << std::endl;
32  exit(1);
33  }
34 #endif//CPPL_DEBUG
35 
36  for(CPPL_INT j=0; j<matA.n; j++){
37  for(CPPL_INT i=j; i<matA.n; i++){
38  matA.darray[j][i] -= matB.darray[j][i];
39  }
40  }
41 
42  return matA;
43 }
_dgematrix i(const _dgbmatrix &mat)
double ** darray
array of pointers of column head addresses
Definition: _dsymatrix.hpp:13
CPPL_INT n
matrix column size
Definition: dsymatrix.hpp:10
CPPL_INT n
matrix column size
Definition: _dsymatrix.hpp:11
double ** darray
array of pointers of column head addresses
Definition: dsymatrix.hpp:12
_dgematrix operator* ( const _dsymatrix matA,
const dsymatrix matB 
)
inline

_dsymatrix*dsymatrix operator

Definition at line 47 of file _dsymatrix-dsymatrix.hpp.

References _(), dsymatrix::array, dgematrix::array, _dsymatrix::array, dsymatrix::complete(), _dsymatrix::destroy(), dsymatrix::m, dgematrix::m, dsymatrix::n, and _dsymatrix::n.

48 {CPPL_VERBOSE_REPORT;
49 #ifdef CPPL_DEBUG
50  if(matA.n!=matB.n){
51  ERROR_REPORT;
52  std::cerr << "These two matrises can not make a product." << std::endl
53  << "Your input was (" << matA.n << "x" << matA.n << ") * (" << matB.n << "x" << matB.n << ")." << std::endl;
54  exit(1);
55  }
56 #endif//CPPL_DEBUG
57 
58  matB.complete();
59 
60  dgematrix newmat(matA.n, matA.n);
61  char side ='l';
62  char uplo ='l';
63  double alpha =1.;
64  double beta =0.;
65 
66  dsymm_( &side, &uplo, &matA.n, &matB.n, &alpha, matA.array, &matA.n, matB.array, &matB.m, &beta, newmat.array, &newmat.m );
67 
68  matA.destroy();
69  return _(newmat);
70 }
void destroy() const
Real Double-precision General Dence Matrix Class.
Definition: dgematrix.hpp:3
double * array
1D array to store matrix data
Definition: _dsymatrix.hpp:12
CPPL_INT const & m
matrix row size
Definition: dsymatrix.hpp:9
CPPL_INT n
matrix column size
Definition: dsymatrix.hpp:10
void complete() const
double * array
1D array to store matrix data
Definition: dsymatrix.hpp:11
CPPL_INT n
matrix column size
Definition: _dsymatrix.hpp:11
_dcovector _(dcovector &vec)