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
67
68
69
70
71
72
73
74
75
|
#include "rb_lapack.h"
extern integer ilaenv_(integer* ispec, char* name, char* opts, integer* n1, integer* n2, integer* n3, integer* n4);
static VALUE
rblapack_ilaenv(int argc, VALUE *argv, VALUE self){
VALUE rblapack_ispec;
integer ispec;
VALUE rblapack_name;
char *name;
VALUE rblapack_opts;
char *opts;
VALUE rblapack_n1;
integer n1;
VALUE rblapack_n2;
integer n2;
VALUE rblapack_n3;
integer n3;
VALUE rblapack_n4;
integer n4;
VALUE rblapack___out__;
integer __out__;
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 __out__ = NumRu::Lapack.ilaenv( ispec, name, opts, n1, n2, n3, n4, [:usage => usage, :help => help])\n\n\nFORTRAN MANUAL\n INTEGER FUNCTION ILAENV( ISPEC, NAME, OPTS, N1, N2, N3, N4 )\n\n* Purpose\n* =======\n*\n* ILAENV is called from the LAPACK routines to choose problem-dependent\n* parameters for the local environment. See ISPEC for a description of\n* the parameters.\n*\n* ILAENV returns an INTEGER\n* if ILAENV >= 0: ILAENV returns the value of the parameter specified by ISPEC\n* if ILAENV < 0: if ILAENV = -k, the k-th argument had an illegal value.\n*\n* This version provides a set of parameters which should give good,\n* but not optimal, performance on many of the currently available\n* computers. Users are encouraged to modify this subroutine to set\n* the tuning parameters for their particular machine using the option\n* and problem size information in the arguments.\n*\n* This routine will not function correctly if it is converted to all\n* lower case. Converting it to all upper case is allowed.\n*\n\n* Arguments\n* =========\n*\n* ISPEC (input) INTEGER\n* Specifies the parameter to be returned as the value of\n* ILAENV.\n* = 1: the optimal blocksize; if this value is 1, an unblocked\n* algorithm will give the best performance.\n* = 2: the minimum block size for which the block routine\n* should be used; if the usable block size is less than\n* this value, an unblocked routine should be used.\n* = 3: the crossover point (in a block routine, for N less\n* than this value, an unblocked routine should be used)\n* = 4: the number of shifts, used in the nonsymmetric\n* eigenvalue routines (DEPRECATED)\n* = 5: the minimum column dimension for blocking to be used;\n* rectangular blocks must have dimension at least k by m,\n* where k is given by ILAENV(2,...) and m by ILAENV(5,...)\n* = 6: the crossover point for the SVD (when reducing an m by n\n* matrix to bidiagonal form, if max(m,n)/min(m,n) exceeds\n* this value, a QR factorization is used first to reduce\n* the matrix to a triangular form.)\n* = 7: the number of processors\n* = 8: the crossover point for the multishift QR method\n* for nonsymmetric eigenvalue problems (DEPRECATED)\n* = 9: maximum size of the subproblems at the bottom of the\n* computation tree in the divide-and-conquer algorithm\n* (used by xGELSD and xGESDD)\n* =10: ieee NaN arithmetic can be trusted not to trap\n* =11: infinity arithmetic can be trusted not to trap\n* 12 <= ISPEC <= 16:\n* xHSEQR or one of its subroutines,\n* see IPARMQ for detailed explanation\n*\n* NAME (input) CHARACTER*(*)\n* The name of the calling subroutine, in either upper case or\n* lower case.\n*\n* OPTS (input) CHARACTER*(*)\n* The character options to the subroutine NAME, concatenated\n* into a single character string. For example, UPLO = 'U',\n* TRANS = 'T', and DIAG = 'N' for a triangular routine would\n* be specified as OPTS = 'UTN'.\n*\n* N1 (input) INTEGER\n* N2 (input) INTEGER\n* N3 (input) INTEGER\n* N4 (input) INTEGER\n* Problem dimensions for the subroutine NAME; these may not all\n* be required.\n*\n\n* Further Details\n* ===============\n*\n* The following conventions have been used when calling ILAENV from the\n* LAPACK routines:\n* 1) OPTS is a concatenation of all of the character options to\n* subroutine NAME, in the same order that they appear in the\n* argument list for NAME, even if they are not used in determining\n* the value of the parameter specified by ISPEC.\n* 2) The problem dimensions N1, N2, N3, N4 are specified in the order\n* that they appear in the argument list for NAME. N1 is used\n* first, N2 second, and so on, and unused problem dimensions are\n* passed a value of -1.\n* 3) The parameter value returned by ILAENV is checked for validity in\n* the calling subroutine. For example, ILAENV is used to retrieve\n* the optimal blocksize for STRTRI as follows:\n*\n* NB = ILAENV( 1, 'STRTRI', UPLO // DIAG, N, -1, -1, -1 )\n* IF( NB.LE.1 ) NB = MAX( 1, N )\n*\n* =====================================================================\n*\n* .. Local Scalars ..\n INTEGER I, IC, IZ, NB, NBMIN, NX\n LOGICAL CNAME, SNAME\n CHARACTER C1*1, C2*2, C4*2, C3*3, SUBNAM*6\n* ..\n* .. Intrinsic Functions ..\n INTRINSIC CHAR, ICHAR, INT, MIN, REAL\n* ..\n* .. External Functions ..\n INTEGER IEEECK, IPARMQ\n EXTERNAL IEEECK, IPARMQ\n* ..\n\n");
return Qnil;
}
if (rb_hash_aref(rblapack_options, sUsage) == Qtrue) {
printf("%s\n", "USAGE:\n __out__ = NumRu::Lapack.ilaenv( ispec, name, opts, n1, n2, n3, n4, [:usage => usage, :help => help])\n");
return Qnil;
}
} else
rblapack_options = Qnil;
if (argc != 7 && argc != 7)
rb_raise(rb_eArgError,"wrong number of arguments (%d for 7)", argc);
rblapack_ispec = argv[0];
rblapack_name = argv[1];
rblapack_opts = argv[2];
rblapack_n1 = argv[3];
rblapack_n2 = argv[4];
rblapack_n3 = argv[5];
rblapack_n4 = argv[6];
if (argc == 7) {
} else if (rblapack_options != Qnil) {
} else {
}
ispec = NUM2INT(rblapack_ispec);
opts = StringValueCStr(rblapack_opts);
n2 = NUM2INT(rblapack_n2);
n4 = NUM2INT(rblapack_n4);
name = StringValueCStr(rblapack_name);
n3 = NUM2INT(rblapack_n3);
n1 = NUM2INT(rblapack_n1);
__out__ = ilaenv_(&ispec, name, opts, &n1, &n2, &n3, &n4);
rblapack___out__ = INT2NUM(__out__);
return rblapack___out__;
}
void
init_lapack_ilaenv(VALUE mLapack, VALUE sH, VALUE sU, VALUE zero){
sHelp = sH;
sUsage = sU;
rblapack_ZERO = zero;
rb_define_module_function(mLapack, "ilaenv", rblapack_ilaenv, -1);
}
|
