param *p;
{
const double hbar = GSL_CONST_MKSA_PLANCKS_CONSTANT_HBAR;
- /*
- DEBUG_PRINT2("% e", GSL_DBL_EPSILON);
- DEBUG_PRINT2("% e", DBL_MIN);
- DEBUG_PRINT2("% e", (2.0 * p->m / (hbar * hbar)));
- DEBUG_PRINT2("% e", (p->m * p->omega * p->omega));
- DEBUG_PRINT2("% e", (hbar * p->omega));
- */
DEBUG_PRINT2("% e", p->m);
DEBUG_PRINT2("% d", p->l);
DEBUG_PRINT2("% e", p->E);
dma_array_print(p->u);
dma_array_print(p->dudr);
*/
+ /*
+ DEBUG_PRINT2("% e", GSL_DBL_EPSILON);
+ DEBUG_PRINT2("% e", DBL_MIN);
+ DEBUG_PRINT2("% e", (2.0 * p->m / (hbar * hbar)));
+ DEBUG_PRINT2("% e", (p->m * p->omega * p->omega));
+ DEBUG_PRINT2("% e", (hbar * p->omega));
+ */
return;
}