Skip to content
Snippets Groups Projects
Commit 6d308ef2 authored by lschulze's avatar lschulze
Browse files

Allow for bicubic and quadratic interpolation here

parent 6d9e0ef2
No related branches found
No related tags found
No related merge requests found
......@@ -137,16 +137,10 @@ void MonomialBasis<dim>::generateBasis(std::vector<unsigned int> m, unsigned int
grid_key_dx_iterator_sub_bc<dim> it(grid, start, stop, bc);
// Finally compute alpha_min
<<<<<<< Updated upstream
//unsigned char alphaMin = static_cast<unsigned char>(!(mSum % 2)); // if mSum is even, alpha_min must be 1
=======
unsigned char alphaMin = static_cast<unsigned char>(!(mSum % 2)); // if mSum is even, alpha_min must be 1
//std::cout<<"AlphaMin: "<<alphaMin<<std::endl;
alphaMin = 0; // we want to always have 1 in the basis
>>>>>>> Stashed changes
unsigned char alphaMin = 0; // we want to always have 1 in the basis
//std::cout<<"AlphaMin: "<<alphaMin<<std::endl;
while (it.isNext())
{
Point<dim, long int> p = it.get().get_k();
......@@ -179,12 +173,8 @@ void MonomialBasis<dim>::generateInterpolBasis(unsigned int orderLimit)
grid_key_dx_iterator_sub_bc<dim> it(grid, start, stop, bc);
// Finally compute alpha_min
<<<<<<< Updated upstream
//unsigned char alphaMin = static_cast<unsigned char>(!(mSum % 2)); // if mSum is even, alpha_min must be 1
//std::cout<<"AlphaMin: "<<alphaMin<<std::endl;
unsigned char alphaMin = 0; // we want to always have 1 in the basis
=======
>>>>>>> Stashed changes
//std::cout<<"AlphaMin: "<<alphaMin<<std::endl;
while (it.isNext())
{
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment