?? spdfd.h
字號(hào):
// LAPACK++ (V. 1.1)// (C) 1992-1996 All Rights Reserved.#ifndef _LA_SPD_FACT_DOUBLE_H_#define _LA_SPD_FACT_DOUBLE_H_#include LA_SPD_MAT_DOUBLE_H#include "lapack.h"class LaSpdFactDouble{ int size_; int gdim_; LaSpdMatDouble S_; public: LaSpdFactDouble(); LaSpdFactDouble(int,int); LaSpdFactDouble(const LaSpdFactDouble&); ~LaSpdFactDouble(); LaSpdFactDouble& ref(LaSpdFactDouble&); LaSpdFactDouble& ref(LaSpdMatDouble&); LaSpdFactDouble& copy(const LaSpdFactDouble&); LaSpdFactDouble& copy(const LaSpdMatDouble&); LaSpdMatDouble& S() { return S_; } int size() { return size_; } int gdim() { return gdim_; }};inline LaSpdFactDouble::LaSpdFactDouble(): S_(), size_(0), gdim_(0){}inline LaSpdFactDouble::LaSpdFactDouble(int m,int n):S_(m,n), size_(n), gdim_(m){}inline LaSpdFactDouble::LaSpdFactDouble(const LaSpdFactDouble &X){ size_ = X.size_; gdim_ = X.gdim_; S_.copy(X.S_);}inline LaSpdFactDouble::~LaSpdFactDouble(){}inline LaSpdFactDouble& LaSpdFactDouble::ref(LaSpdFactDouble &X){ size_ = X.size_; gdim_ = X.gdim_; S_.ref(X.S_); return *this;}inline LaSpdFactDouble& LaSpdFactDouble::ref(LaSpdMatDouble &X){ size_ = X.size(1); gdim_ = X.gdim(0); S_.ref(X); return *this;} inline LaSpdFactDouble& LaSpdFactDouble::copy(const LaSpdFactDouble &X){ size_ = X.size_; gdim_ = X.gdim_; S_.copy(X.S_); return *this;}inline LaSpdFactDouble& LaSpdFactDouble::copy(const LaSpdMatDouble &X){ size_ = X.size(1); gdim_ = X.gdim(0); S_.copy(X); return *this;}inline void LaSpdMatFactorize(LaSpdMatDouble &A, LaSpdFactDouble &AF){ char uplo = 'L'; integer N = A.size(1), lda = A.gdim(0), info = 0; AF.copy(A); F77NAME(dpotrf)(&uplo, &N, &(AF.S()(0,0)), &lda, &info);}inline void LaLinearSolve(LaSpdFactDouble &AF, LaGenMatDouble &X, LaGenMatDouble &B){ char uplo = 'L'; integer N = AF.size(), nrhs = X.size(1), lda = AF.gdim(), ldb = B.size(0), info = 0; X.inject(B); F77NAME(dpotrs)(&uplo, &N, &nrhs, &(AF.S()(0,0)), &lda, &X(0,0), &ldb, &info);}#endif // _LA_SPD_FACT_DOUBLE_H_
?? 快捷鍵說(shuō)明
復(fù)制代碼
Ctrl + C
搜索代碼
Ctrl + F
全屏模式
F11
切換主題
Ctrl + Shift + D
顯示快捷鍵
?
增大字號(hào)
Ctrl + =
減小字號(hào)
Ctrl + -