formula = form;
}
-//TODO: Discussion on an efficient implementation of Point additions.
+//TODO:[GK] Discussion on an efficient implementation of Point additions.
template <int dim>
double
using namespace std;
#endif
-//TODO:Optimize construction of vectors!
+//TODO:[GK] Optimize construction of vectors!
template <int dim>
void
}
-//TODO: Implement derivatives
+//TODO:[GK] Implement derivatives
Tensor<1,2>
LSingularityFunction::gradient (const Point<2> &p,
}
-//TODO: Implement derivatives
+//TODO:[GK] Implement derivatives
Tensor<1,2>
SlitSingularityFunction::gradient (const Point<2> &/*p*/,
const unsigned int n = GeometryInfo<dim>::faces_per_cell
* (sub ? GeometryInfo<dim>::subfaces_per_face : 1);
-//TODO: Can we keep it that simple for non-symmetric formulae?
+//TODO:[GK] Can we keep it that simple for non-symmetric formulae?
// Otherwise we'll have to include the weights in the project* functions.
// This may pose a problem anyway, since we do not want to switch endpoints
// in Gauss-Radau formulae.
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