1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66
|
#include "rb_lapack.h"
extern VOID sladiv_(real* a, real* b, real* c, real* d, real* p, real* q);
static VALUE
rblapack_sladiv(int argc, VALUE *argv, VALUE self){
VALUE rblapack_a;
real a;
VALUE rblapack_b;
real b;
VALUE rblapack_c;
real c;
VALUE rblapack_d;
real d;
VALUE rblapack_p;
real p;
VALUE rblapack_q;
real q;
VALUE rblapack_options;
if (argc > 0 && TYPE(argv[argc-1]) == T_HASH) {
argc--;
rblapack_options = argv[argc];
if (rb_hash_aref(rblapack_options, sHelp) == Qtrue) {
printf("%s\n", "USAGE:\n p, q = NumRu::Lapack.sladiv( a, b, c, d, [:usage => usage, :help => help])\n\n\nFORTRAN MANUAL\n SUBROUTINE SLADIV( A, B, C, D, P, Q )\n\n* Purpose\n* =======\n*\n* SLADIV performs complex division in real arithmetic\n*\n* a + i*b\n* p + i*q = ---------\n* c + i*d\n*\n* The algorithm is due to Robert L. Smith and can be found\n* in D. Knuth, The art of Computer Programming, Vol.2, p.195\n*\n\n* Arguments\n* =========\n*\n* A (input) REAL\n* B (input) REAL\n* C (input) REAL\n* D (input) REAL\n* The scalars a, b, c, and d in the above expression.\n*\n* P (output) REAL\n* Q (output) REAL\n* The scalars p and q in the above expression.\n*\n\n* =====================================================================\n*\n* .. Local Scalars ..\n REAL E, F\n* ..\n* .. Intrinsic Functions ..\n INTRINSIC ABS\n* ..\n\n");
return Qnil;
}
if (rb_hash_aref(rblapack_options, sUsage) == Qtrue) {
printf("%s\n", "USAGE:\n p, q = NumRu::Lapack.sladiv( a, b, c, d, [:usage => usage, :help => help])\n");
return Qnil;
}
} else
rblapack_options = Qnil;
if (argc != 4 && argc != 4)
rb_raise(rb_eArgError,"wrong number of arguments (%d for 4)", argc);
rblapack_a = argv[0];
rblapack_b = argv[1];
rblapack_c = argv[2];
rblapack_d = argv[3];
if (argc == 4) {
} else if (rblapack_options != Qnil) {
} else {
}
a = (real)NUM2DBL(rblapack_a);
c = (real)NUM2DBL(rblapack_c);
b = (real)NUM2DBL(rblapack_b);
d = (real)NUM2DBL(rblapack_d);
sladiv_(&a, &b, &c, &d, &p, &q);
rblapack_p = rb_float_new((double)p);
rblapack_q = rb_float_new((double)q);
return rb_ary_new3(2, rblapack_p, rblapack_q);
}
void
init_lapack_sladiv(VALUE mLapack, VALUE sH, VALUE sU, VALUE zero){
sHelp = sH;
sUsage = sU;
rblapack_ZERO = zero;
rb_define_module_function(mLapack, "sladiv", rblapack_sladiv, -1);
}
|