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; |