Hi,
I am trying to develop an exponential filter for my code and for that I need
to construct the Vandermonde matrix using legendre polynomials. I've taken a
very simple code to find the output of legendre polynomials at various
orders as a test case and I seem to be getting the wrong results. I am using
version 7.0.0
Code:
#include <iostream>
#include <base/polynomial.h>
using namespace dealii;
int main ()
{
degree = 1;
Polynomials::Legendre test(degree);
std::cout << test.value(1) << std::endl;
}
Output:
For degree = 0
1
For degree = 1
1.73205
For degree = 2
2.23607
Shouldn't the value be 1 for all the three cases or am I missing a very
simple point here?
--
Thanks,
gk
_______________________________________________
dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
