| 62 |
|
#include "config.h" |
| 63 |
|
#include <cmath> |
| 64 |
|
#include <cstdlib> |
| 65 |
+ |
#include <iostream> |
| 66 |
|
#include "math/Factorials.hpp" |
| 67 |
|
|
| 68 |
|
#ifndef NONBONDED_SLATERINTEGRALS_HPP |
| 95 |
|
{ |
| 96 |
|
RealType Term = 1.; |
| 97 |
|
RosenA_ = Term; |
| 98 |
< |
for (unsigned nu=1; nu<=n; nu++) |
| 98 |
> |
for (int nu=1; nu<=n; nu++) |
| 99 |
|
{ |
| 100 |
|
Term *= a / nu; |
| 101 |
|
RosenA_ += Term; |
| 174 |
|
* @note N. Rosen, Phys. Rev., 38 (1931), 255 |
| 175 |
|
*/ |
| 176 |
|
inline RealType RosenD(int m, int n, int p) |
| 176 |
– |
// [?ERROR?] can't return int |
| 177 |
|
{ |
| 178 |
– |
printf("RosenD\n"); |
| 179 |
– |
exit(0); |
| 178 |
|
if (m+n+p > maxFact) |
| 179 |
|
{ |
| 180 |
|
printf("Error, arguments exceed maximum factorial computed %d > %d\n", m+n+p, maxFact); |
| 234 |
|
else |
| 235 |
|
{ |
| 236 |
|
// Not implemented |
| 237 |
< |
printf("ERROR, sSTOCoulInt cannot compute from arguments\n"); |
| 238 |
< |
printf("a = %lf b = %lf m = %d n = %d R = %lf\n",a, b, m, n, R); |
| 239 |
< |
exit(0); |
| 237 |
> |
cerr << "ERROR, sSTOCoulInt cannot compute from arguments\n"; |
| 238 |
> |
cerr << "a = " << a << "\tb = " << b << "\tm = " << m <<"\tn = " << n << "\t R = " << R << "\n"; |
| 239 |
> |
//printf("ERROR, sSTOCoulInt cannot compute from arguments\n"); |
| 240 |
> |
//printf("a = %lf b = %lf m = %d n = %d R = %lf\n",a, b, m, n, R); |
| 241 |
> |
//exit(0); |
| 242 |
|
} |
| 243 |
|
} |
| 244 |
|
else |
