|
SyFi
0.3
|
Go to the source code of this file.
Functions | |
| int | main () |
| int main | ( | ) |
Definition at line 8 of file robust.cpp.
References SyFi::compare_archives(), SyFi::Robust::compute_basis_functions(), SyFi::div(), SyFi::initSyFi(), SyFi::istr(), SyFi::StandardFE::N(), and SyFi::StandardFE::nbf().
{
initSyFi(2);
ReferenceTriangle triange;
Robust fe(triange);
fe.compute_basis_functions();
for (unsigned int i=0; i< fe.nbf(); i++) {
cout <<"fe.N("<<i<<")= "<<fe.N(i)<<endl;
cout <<"div(fe.N("<<i<<"))= "<<(div(fe.N(i))).expand()<<endl;
}
// regression test
archive ar;
for (unsigned int i=0; i< fe.nbf(); i++) {
ar.archive_ex(fe.N(i), istr("N",i).c_str());
ar.archive_ex((div(fe.N(i))).expand(), istr("divN",i).c_str());
}
ofstream vfile("robust.gar.v");
vfile << ar; vfile.close();
if(!compare_archives("robust.gar.v", "robust.gar.r")) {
cerr << "Failure!" << endl;
return -1;
}
return 0;
}