for (unsigned int d=0; d<dim; ++d)
{
Point<dim> shift;
-//TODO: unify the places where the finite differencing step length is used
+//TODO:[GK] unify the places where the finite differencing step length is used
shift (d) = 1.e-6;
// generate points and FEValues
const Tensor<1,dim>& left
= fe_internal.differences[d1+dim]->shape_grad(shape, q);
-//TODO: unify the places where the finite differencing step length is used
+//TODO:[GK] unify the places where the finite differencing step length is used
// compute the second
// derivative from a
// symmetric difference
// Step 2
sigma = q*t;
-//TODO: Find a really good breakdown criterion
-// The absolute one detects breakdown instead of convergence
+//TODO:[GK] Find a really good breakdown criterion. The absolute one detects breakdown instead of convergence
if (fabs(sigma/rho) < additional_data.breakdown)
return breakdown;
// Step 3
inverses[row]->gauss_jordan();
};
-//TODO: implement Tvmult
+//TODO:[WB,GK] implement Tvmult
template<typename number>
template<typename number2>
};
};
-//TODO: implement Tvmult
+//TODO:[WB,GK] implement Tvmult
template <typename number>
template <typename number2>
In the beginning the Universe was created. This has made a lot of
people very angry and has been widely regarded as a bad move.
Douglas Adams