jess
/
fast_sparse_interpolation
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

changed main.cpp

This commit is contained in:
David Holzmüller 2019-04-07 22:26:09 +02:00
parent 71e618290f
commit 245764aef0
1 changed files with 5 additions and 5 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

10
test/main.cpp
View File

@ -112,8 +112,8 @@ class Timer {
// using namespace fsi; // using namespace fsi;
int main() { int main() {
constexpr size_t d = 2; constexpr size_t d = 8;
size_t bound = 2; size_t bound = 24;
fsi::TemplateBoundedSumIterator<d> it(bound); fsi::TemplateBoundedSumIterator<d> it(bound);
// fsi::BoundedSumIterator it(d, bound); // fsi::BoundedSumIterator it(d, bound);
std::vector<MonomialFunctions> phi(d); std::vector<MonomialFunctions> phi(d);
@ -122,7 +122,7 @@ int main() {
auto rhs = evaluateFunction(it, f, x); auto rhs = evaluateFunction(it, f, x);
auto op = createInterpolationOperator(it, phi, x); auto op = createInterpolationOperator(it, phi, x);
std::cout << rhs << "\n"; // std::cout << rhs << "\n";
Timer timer; Timer timer;
op.prepareSolve(); op.prepareSolve();
@ -133,12 +133,12 @@ int main() {
auto c = op.solve(rhs); auto c = op.solve(rhs);
// auto c = fsi::interpolate(f, it, phi, x); // auto c = fsi::interpolate(f, it, phi, x);
std::cout << "Time for solve(): " << timer.elapsed() << " s\n"; std::cout << "Time for solve(): " << timer.elapsed() << " s\n";
std::cout << c << "\n"; // std::cout << c << "\n";
timer.reset(); timer.reset();
auto b = op.apply(c); auto b = op.apply(c);
std::cout << "Time for apply(): " << timer.elapsed() << " s\n"; std::cout << "Time for apply(): " << timer.elapsed() << " s\n";
std::cout << b << "\n"; // std::cout << b << "\n";
std::cout << "Number of points: " << it.numValues() << "\n"; std::cout << "Number of points: " << it.numValues() << "\n";