1 /* roots/roots.h 2 * 3 * Copyright (C) 1996, 1997, 1998, 1999, 2000, 2007 Reid Priedhorsky, Brian Gough 4 * 5 * This program is free software; you can redistribute it and/or modify 6 * it under the terms of the GNU General Public License as published by 7 * the Free Software Foundation; either version 3 of the License, or (at 8 * your option) any later version. 9 * 10 * This program is distributed in the hope that it will be useful, but 11 * WITHOUT ANY WARRANTY; without even the implied warranty of 12 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU 13 * General Public License for more details. 14 * 15 * You should have received a copy of the GNU General Public License 16 * along with this program; if not, write to the Free Software 17 * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301, USA. 18 */ 19 20 /* roots.h -- declarations for internal root finding and RF support stuff. */ 21 22 #ifndef __ROOTS_H__ 23 #define __ROOTS_H__ 24 25 /* Call the pointed-to function with argument x, put its result in y, and 26 return an error if the function value is Inf/Nan. */ 27 28 #define SAFE_FUNC_CALL(f, x, yp) \ 29 do { \ 30 *yp = GSL_FN_EVAL(f,x); \ 31 if (!gsl_finite(*yp)) \ 32 GSL_ERROR("function value is not finite", GSL_EBADFUNC); \ 33 } while (0) 34 35 #endif /* __ROOTS_H__ */ 36 37 38