| author | francesco <francesco.bbt@gmail.com> | 2011年04月24日 12:44:25 +0200 |
|---|---|---|
| committer | francesco <francesco.bbt@gmail.com> | 2011年04月24日 12:44:25 +0200 |
| commit | 90a9dc178acca62749df1393eacffd264350cae8 (patch) | |
| tree | ec2e1518ec60fbd14393a11b6a2d42404737664b /sf_implement.h | |
| parent | 2b1733183834a93f4cfcff03de15b372bc69b4ac (diff) | |
| download | gsl-shell-90a9dc178acca62749df1393eacffd264350cae8.tar.gz | |
| -rw-r--r-- | sf_implement.h | 33 |
diff --git a/sf_implement.h b/sf_implement.h index 65030f06..34e12824 100644 --- a/sf_implement.h +++ b/sf_implement.h @@ -313,22 +313,6 @@ int GSH_LUA_NAME(hypergU) (lua_State *L) return push_gsl_result (L, &res); } -int GSH_LUA_NAME(hyperg2F1conj) (lua_State *L) -{ - Complex a = luaL_checkcomplex(L, 1); - double c = luaL_checknumber(L, 2); - double x = luaL_checknumber (L, 3); - gsl_sf_result res; - int status; - - status = gsl_sf_hyperg_2F1_conj_e (creal(a), cimag(a), c, x, &res); - - if (status != GSL_SUCCESS) - return luaL_error (L, "hyperg2F1conj: %s", gsl_strerror (status)); - - return push_gsl_result (L, &res); -} - int GSH_LUA_NAME(laguerre) (lua_State *L) { int i = luaL_checkinteger (L, 1); @@ -517,20 +501,3 @@ int GSH_LUA_NAME(legendreQ) (lua_State *L) return push_gsl_result (L, &res); } - -int GSH_LUA_NAME(cdilog) (lua_State *L) -{ - Complex z = luaL_checkcomplex (L, 1); - gsl_sf_result rr, ri; - double zr = creal(z), zi = cimag(z); - double r = sqrt(zr*zr+zi*zi), th = atan2(cimag(z), creal(z)); - int status; - - status = gsl_sf_complex_dilog_e (r, th, &rr, &ri); - - if (status != GSL_SUCCESS) - return luaL_error (L, "cdilog: %s", gsl_strerror (status)); - - lua_pushcomplex (L, rr.val + I * ri.val); - return 1; -} |