| 1 |
< |
/* |
| 1 |
> |
/* |
| 2 |
|
* Copyright (c) 2005 The University of Notre Dame. All Rights Reserved. |
| 3 |
|
* |
| 4 |
|
* The University of Notre Dame grants you ("Licensee") a |
| 41 |
|
|
| 42 |
|
#include <stdio.h> |
| 43 |
|
#include <cmath> |
| 44 |
< |
|
| 44 |
> |
#include <limits> |
| 45 |
|
#include "math/RealSphericalHarmonic.hpp" |
| 46 |
|
|
| 47 |
|
using namespace oopse; |
| 85 |
|
// check parameters |
| 86 |
|
if (m < 0 || m > l || fabs(x) > 1.0) { |
| 87 |
|
printf("LegendreP got a bad argument: l = %d\tm = %d\tx = %lf\n", l, m, x); |
| 88 |
< |
return NAN; |
| 88 |
> |
// return NAN; |
| 89 |
> |
return std::numeric_limits <double>:: quiet_NaN(); |
| 90 |
|
} |
| 91 |
|
|
| 92 |
|
double pmm = 1.0; |