本文整理汇总了C++中MatrixTemplate::isSquare方法的典型用法代码示例。如果您正苦于以下问题：C++ MatrixTemplate::isSquare方法的具体用法？C++ MatrixTemplate::isSquare怎么用？C++ MatrixTemplate::isSquare使用的例子？那么, 这里精选的方法代码示例或许可以为您提供帮助。您也可以进一步了解该方法所在类MatrixTemplate的用法示例。
在下文中一共展示了MatrixTemplate::isSquare方法的4个代码示例，这些例子默认根据受欢迎程度排序。您可以为喜欢或者感觉有用的代码点赞，您的评价将有助于系统推荐出更棒的C++代码示例。
示例1: LtBackSubstitute
bool LtBackSubstitute(const MatrixTemplate<T>& a, const VectorTemplate<T>& b, VectorTemplate<T>& x)
{
Assert(a.isSquare());
Assert(a.n == b.n);
if(x.empty()) x.resize(a.n);
else Assert(a.n == x.n);
int n=a.n;
T aii,sum;
for(int i=n-1; i>=0; i--) {
aii=a(i,i);
sum=b[i];
for(int j=i+1; j<n; j++)
sum-=a(j,i)*x[j];
if(aii == 0) {
if(!FuzzyZero(sum,(T)kBackSubZeroTolerance)) {
//LOG4CXX_ERROR(KrisLibrary::logger(),"LtBackSubstitute: dividing by zero: "<<sum<<"/"<<aii);
return false;
}
x[i]=0;
}
else
x[i]=sum/aii;
}
return true;
}示例2: LBackSubstitute
bool LBackSubstitute(const MatrixTemplate<T>& a, const VectorTemplate<T>& b, VectorTemplate<T>& x)
{
Assert(a.isSquare());
Assert(a.n == b.n);
if(x.empty()) x.resize(a.n);
else Assert(a.n == x.n);
int n=a.n;
T aii,sum;
for(int i=0; i<n; i++) {
aii=a(i,i);
sum=b[i];
for(int j=0; j<i; j++)
sum-=a(i,j)*x[j];
if(aii == 0) {
if(!FuzzyZero(sum,(T)kBackSubZeroTolerance)) {
//cerr<<"LBackSubstitute: dividing by zero: "<<sum<<"/"<<aii<<endl;
//cerr<<MatrixPrinter(a)<<endl;
return false;
}
x[i]=0;
}
else
x[i]=sum/aii;
}
return true;
}示例3: U1BackSubstitute
void U1BackSubstitute(const MatrixTemplate<T>& a, const VectorTemplate<T>& b, VectorTemplate<T>& x)
{
Assert(a.isSquare());
Assert(a.n == b.n);
if(x.empty()) x.resize(a.n);
else Assert(a.n == x.n);
int n=a.n;
T sum;
for(int i=n-1; i>=0; i--) {
sum=b(i);
for(int j=i+1; j<n; j++)
sum-=a(i,j)*x[j];
x[i]=sum;
}
}示例4: L1BackSubstitute
void L1BackSubstitute(const MatrixTemplate<T>& a, const VectorTemplate<T>& b, VectorTemplate<T>& x)
{
Assert(a.isSquare());
Assert(a.n == b.n);
if(x.empty()) x.resize(a.n);
else Assert(a.n == x.n);
int n=a.n;
T sum;
for(int i=0; i<n; i++) {
sum=b[i];
for(int j=0; j<i; j++)
sum-=a(i,j)*x[j];
x[i]=sum;
}
}