void mexStub4(int nlhs, mxArray* plhs[],
int nrhs, const mxArray* prhs[])
{
const char* mw_err_txt_ = 0;
int* in0_ =0; /* jt */
int* in1_ =0; /* l */
double* in2_ =0; /* val */
double* in3_ =0; /* tsh */
int in4_; /* N */
int in5_; /* j */
double in6_; /* gamma */
int out0_; /* keep */
if (mxGetM(prhs[0])*mxGetN(prhs[0]) != 0) {
in0_ = mxWrapGetArray_int(prhs[0], &mw_err_txt_);
if (mw_err_txt_)
goto mw_err_label;
} else
in0_ = NULL;
if (mxGetM(prhs[1])*mxGetN(prhs[1]) != 0) {
in1_ = mxWrapGetArray_int(prhs[1], &mw_err_txt_);
if (mw_err_txt_)
goto mw_err_label;
} else
in1_ = NULL;
if (mxGetM(prhs[2])*mxGetN(prhs[2]) != 0) {
in2_ = mxGetPr(prhs[2]);
} else
in2_ = NULL;
if (mxGetM(prhs[3])*mxGetN(prhs[3]) != 0) {
in3_ = mxGetPr(prhs[3]);
} else
in3_ = NULL;
in4_ = (int) mxWrapGetScalar(prhs[4], &mw_err_txt_);
if (mw_err_txt_)
goto mw_err_label;
in5_ = (int) mxWrapGetScalar(prhs[5], &mw_err_txt_);
if (mw_err_txt_)
goto mw_err_label;
in6_ = (double) mxWrapGetScalar(prhs[6], &mw_err_txt_);
if (mw_err_txt_)
goto mw_err_label;
if (mexprofrecord_)
mexprofrecord_[4]++;
out0_ = minisgammaloc(in0_, in1_, in2_, in3_, in4_, in5_, in6_);
plhs[0] = mxCreateDoubleMatrix(1, 1, mxREAL);
*mxGetPr(plhs[0]) = out0_;
mw_err_label:
if (in0_) mxFree(in0_);
if (in1_) mxFree(in1_);
if (mw_err_txt_)
mexErrMsgTxt(mw_err_txt_);
}
void mexStub1(int nlhs, mxArray* plhs[],
int nrhs, const mxArray* prhs[])
{
const char* mw_err_txt_ = 0;
int in0_; /* 1 */
int in1_; /* 2 */
int out0_; /* a */
in0_ = (int) mxWrapGetScalar(prhs[0], &mw_err_txt_);
if (mw_err_txt_)
goto mw_err_label;
in1_ = (int) mxWrapGetScalar(prhs[1], &mw_err_txt_);
if (mw_err_txt_)
goto mw_err_label;
if (mexprofrecord_)
mexprofrecord_[1]++;
MWF77_foo(&out0_, &in0_, &in1_);
plhs[0] = mxCreateDoubleMatrix(1, 1, mxREAL);
*mxGetPr(plhs[0]) = out0_;
mw_err_label:
if (mw_err_txt_)
mexErrMsgTxt(mw_err_txt_);
}
void mexStub6(int nlhs, mxArray* plhs[],
int nrhs, const mxArray* prhs[])
{
const char* mw_err_txt_ = 0;
double* in0_ =0; /* W */
int in1_; /* M */
int in2_; /* T */
int in3_; /* K */
int in4_; /* fac */
double* in5_ =0; /* Y */
int in6_; /* Nt */
double in7_; /* tpad */
double in8_; /* eta */
int in9_; /* skip */
double in10_; /* gamma */
double* in11_ =0; /* nlps */
int in12_; /* wantR */
int in13_; /* verb */
double* out1_=0; /* t */
int* out2_=0; /* l */
double* out3_=0; /* a */
int* out4_=0; /* Ns */
double* out5_=0; /* R */
mwSize dim14_; /* maxNs */
mwSize dim15_; /* maxNs */
mwSize dim16_; /* maxNs */
mwSize dim17_; /* 1 */
mwSize dim18_; /* M */
mwSize dim19_; /* Nt */
dim14_ = (mwSize) mxWrapGetScalar(prhs[14], &mw_err_txt_);
dim15_ = (mwSize) mxWrapGetScalar(prhs[15], &mw_err_txt_);
dim16_ = (mwSize) mxWrapGetScalar(prhs[16], &mw_err_txt_);
dim17_ = (mwSize) mxWrapGetScalar(prhs[17], &mw_err_txt_);
dim18_ = (mwSize) mxWrapGetScalar(prhs[18], &mw_err_txt_);
dim19_ = (mwSize) mxWrapGetScalar(prhs[19], &mw_err_txt_);
if (mxGetM(prhs[0])*mxGetN(prhs[0]) != 0) {
in0_ = mxGetPr(prhs[0]);
} else
in0_ = NULL;
in1_ = (int) mxWrapGetScalar(prhs[1], &mw_err_txt_);
if (mw_err_txt_)
goto mw_err_label;
in2_ = (int) mxWrapGetScalar(prhs[2], &mw_err_txt_);
if (mw_err_txt_)
goto mw_err_label;
in3_ = (int) mxWrapGetScalar(prhs[3], &mw_err_txt_);
if (mw_err_txt_)
goto mw_err_label;
in4_ = (int) mxWrapGetScalar(prhs[4], &mw_err_txt_);
if (mw_err_txt_)
goto mw_err_label;
if (mxGetM(prhs[5])*mxGetN(prhs[5]) != 0) {
in5_ = mxWrapGetArray_double(prhs[5], &mw_err_txt_);
if (mw_err_txt_)
goto mw_err_label;
} else
in5_ = NULL;
in6_ = (int) mxWrapGetScalar(prhs[6], &mw_err_txt_);
if (mw_err_txt_)
goto mw_err_label;
in7_ = (double) mxWrapGetScalar(prhs[7], &mw_err_txt_);
if (mw_err_txt_)
goto mw_err_label;
in8_ = (double) mxWrapGetScalar(prhs[8], &mw_err_txt_);
if (mw_err_txt_)
goto mw_err_label;
in9_ = (int) mxWrapGetScalar(prhs[9], &mw_err_txt_);
if (mw_err_txt_)
goto mw_err_label;
in10_ = (double) mxWrapGetScalar(prhs[10], &mw_err_txt_);
if (mw_err_txt_)
goto mw_err_label;
if (mxGetM(prhs[11])*mxGetN(prhs[11]) != 0) {
in11_ = mxGetPr(prhs[11]);
} else
in11_ = NULL;
in12_ = (int) mxWrapGetScalar(prhs[12], &mw_err_txt_);
if (mw_err_txt_)
goto mw_err_label;
in13_ = (int) mxWrapGetScalar(prhs[13], &mw_err_txt_);
if (mw_err_txt_)
goto mw_err_label;
out1_ = (double*) mxMalloc(dim14_*sizeof(double));
out2_ = (int*) mxMalloc(dim15_*sizeof(int));
out3_ = (double*) mxMalloc(dim16_*sizeof(double));
out4_ = (int*) mxMalloc(dim17_*sizeof(int));
out5_ = (double*) mxMalloc(dim18_*dim19_*sizeof(double));
if (mexprofrecord_)
mexprofrecord_[6]++;
multiglutton(in0_, in1_, in2_, in3_, in4_, in5_, in6_, in7_, in8_, in9_, in10_, in11_, out1_, out2_, out3_, out4_, out5_, in12_, in13_);
plhs[0] = mxCreateDoubleMatrix(mxGetM(prhs[5]), mxGetN(prhs[5]), mxREAL);
mxWrapCopy_double(plhs[0], in5_, mxGetM(prhs[5])*mxGetN(prhs[5]));
plhs[1] = mxCreateDoubleMatrix(dim14_, 1, mxREAL);
mxWrapCopy_double(plhs[1], out1_, dim14_);
plhs[2] = mxCreateDoubleMatrix(dim15_, 1, mxREAL);
mxWrapCopy_int(plhs[2], out2_, dim15_);
plhs[3] = mxCreateDoubleMatrix(dim16_, 1, mxREAL);
mxWrapCopy_double(plhs[3], out3_, dim16_);
plhs[4] = mxCreateDoubleMatrix(dim17_, 1, mxREAL);
mxWrapCopy_int(plhs[4], out4_, dim17_);
plhs[5] = mxCreateDoubleMatrix(dim18_, dim19_, mxREAL);
mxWrapCopy_double(plhs[5], out5_, dim18_*dim19_);
mw_err_label:
if (in5_) mxFree(in5_);
if (out1_) mxFree(out1_);
if (out2_) mxFree(out2_);
if (out3_) mxFree(out3_);
if (out4_) mxFree(out4_);
if (out5_) mxFree(out5_);
if (mw_err_txt_)
mexErrMsgTxt(mw_err_txt_);
}
void mexStub3(int nlhs, mxArray* plhs[],
int nrhs, const mxArray* prhs[])
{
const char* mw_err_txt_ = 0;
double* in0_ =0; /* S */
int in1_; /* K */
int in2_; /* Nsh */
double* in3_ =0; /* nlps */
int* in4_ =0; /* jt */
int* in5_ =0; /* l */
double* in6_ =0; /* s */
int* out3_=0; /* Nmin */
mwSize dim7_; /* 1 */
dim7_ = (mwSize) mxWrapGetScalar(prhs[7], &mw_err_txt_);
if (mxGetM(prhs[0])*mxGetN(prhs[0]) != 0) {
in0_ = mxGetPr(prhs[0]);
} else
in0_ = NULL;
in1_ = (int) mxWrapGetScalar(prhs[1], &mw_err_txt_);
if (mw_err_txt_)
goto mw_err_label;
in2_ = (int) mxWrapGetScalar(prhs[2], &mw_err_txt_);
if (mw_err_txt_)
goto mw_err_label;
if (mxGetM(prhs[3])*mxGetN(prhs[3]) != 0) {
in3_ = mxGetPr(prhs[3]);
} else
in3_ = NULL;
if (mxGetM(prhs[4])*mxGetN(prhs[4]) != 0) {
in4_ = mxWrapGetArray_int(prhs[4], &mw_err_txt_);
if (mw_err_txt_)
goto mw_err_label;
} else
in4_ = NULL;
if (mxGetM(prhs[5])*mxGetN(prhs[5]) != 0) {
in5_ = mxWrapGetArray_int(prhs[5], &mw_err_txt_);
if (mw_err_txt_)
goto mw_err_label;
} else
in5_ = NULL;
if (mxGetM(prhs[6])*mxGetN(prhs[6]) != 0) {
in6_ = mxWrapGetArray_double(prhs[6], &mw_err_txt_);
if (mw_err_txt_)
goto mw_err_label;
} else
in6_ = NULL;
out3_ = (int*) mxMalloc(dim7_*sizeof(int));
if (mexprofrecord_)
mexprofrecord_[3]++;
locvalidmins(in0_, in1_, in2_, in3_, in4_, in5_, in6_, out3_);
plhs[0] = mxCreateDoubleMatrix(mxGetM(prhs[4]), mxGetN(prhs[4]), mxREAL);
mxWrapCopy_int(plhs[0], in4_, mxGetM(prhs[4])*mxGetN(prhs[4]));
plhs[1] = mxCreateDoubleMatrix(mxGetM(prhs[5]), mxGetN(prhs[5]), mxREAL);
mxWrapCopy_int(plhs[1], in5_, mxGetM(prhs[5])*mxGetN(prhs[5]));
plhs[2] = mxCreateDoubleMatrix(mxGetM(prhs[6]), mxGetN(prhs[6]), mxREAL);
mxWrapCopy_double(plhs[2], in6_, mxGetM(prhs[6])*mxGetN(prhs[6]));
plhs[3] = mxCreateDoubleMatrix(dim7_, 1, mxREAL);
mxWrapCopy_int(plhs[3], out3_, dim7_);
mw_err_label:
if (in4_) mxFree(in4_);
if (in5_) mxFree(in5_);
if (in6_) mxFree(in6_);
if (out3_) mxFree(out3_);
if (mw_err_txt_)
mexErrMsgTxt(mw_err_txt_);
}
void mexStub2(int nlhs, mxArray* plhs[],
int nrhs, const mxArray* prhs[])
{
const char* mw_err_txt_ = 0;
double* in0_ =0; /* W */
int in1_; /* M */
int in2_; /* T */
int in3_; /* K */
int in4_; /* fac */
double* in5_ =0; /* Y */
int in6_; /* Nt */
double* in7_ =0; /* tsh */
int in8_; /* Nsh */
double in9_; /* eta */
int* in10_ =0; /* shflags */
double* out0_=0; /* S */
mwSize dim11_; /* K */
mwSize dim12_; /* Nsh */
dim11_ = (mwSize) mxWrapGetScalar(prhs[11], &mw_err_txt_);
dim12_ = (mwSize) mxWrapGetScalar(prhs[12], &mw_err_txt_);
if (mxGetM(prhs[0])*mxGetN(prhs[0]) != 0) {
in0_ = mxGetPr(prhs[0]);
} else
in0_ = NULL;
in1_ = (int) mxWrapGetScalar(prhs[1], &mw_err_txt_);
if (mw_err_txt_)
goto mw_err_label;
in2_ = (int) mxWrapGetScalar(prhs[2], &mw_err_txt_);
if (mw_err_txt_)
goto mw_err_label;
in3_ = (int) mxWrapGetScalar(prhs[3], &mw_err_txt_);
if (mw_err_txt_)
goto mw_err_label;
in4_ = (int) mxWrapGetScalar(prhs[4], &mw_err_txt_);
if (mw_err_txt_)
goto mw_err_label;
if (mxGetM(prhs[5])*mxGetN(prhs[5]) != 0) {
in5_ = mxGetPr(prhs[5]);
} else
in5_ = NULL;
in6_ = (int) mxWrapGetScalar(prhs[6], &mw_err_txt_);
if (mw_err_txt_)
goto mw_err_label;
if (mxGetM(prhs[7])*mxGetN(prhs[7]) != 0) {
in7_ = mxGetPr(prhs[7]);
} else
in7_ = NULL;
in8_ = (int) mxWrapGetScalar(prhs[8], &mw_err_txt_);
if (mw_err_txt_)
goto mw_err_label;
in9_ = (double) mxWrapGetScalar(prhs[9], &mw_err_txt_);
if (mw_err_txt_)
goto mw_err_label;
if (mxGetM(prhs[10])*mxGetN(prhs[10]) != 0) {
in10_ = mxWrapGetArray_int(prhs[10], &mw_err_txt_);
if (mw_err_txt_)
goto mw_err_label;
} else
in10_ = NULL;
out0_ = (double*) mxMalloc(dim11_*dim12_*sizeof(double));
if (mexprofrecord_)
mexprofrecord_[2]++;
fillscore(in0_, in1_, in2_, in3_, in4_, in5_, in6_, in7_, in8_, out0_, in9_, in10_);
plhs[0] = mxCreateDoubleMatrix(dim11_, dim12_, mxREAL);
mxWrapCopy_double(plhs[0], out0_, dim11_*dim12_);
mw_err_label:
if (out0_) mxFree(out0_);
if (in10_) mxFree(in10_);
if (mw_err_txt_)
mexErrMsgTxt(mw_err_txt_);
}
void mexStub3(int nlhs, mxArray* plhs[],
int nrhs, const mxArray* prhs[])
{
const char* mw_err_txt_ = 0;
double* in0_ =0; /* beta */
int* in1_ =0; /* nterms */
int* in2_ =0; /* m1 */
int* in3_ =0; /* m2 */
dcomplex* in4_ =0; /* mpole */
int* in5_ =0; /* lmp */
dcomplex* in6_ =0; /* mpout */
int* in7_ =0; /* lmpn */
mwSize dim8_; /* 1 */
mwSize dim9_; /* 1 */
mwSize dim10_; /* 1 */
mwSize dim11_; /* 1 */
mwSize dim12_; /* 1 */
mwSize dim13_; /* 1 */
dim8_ = (mwSize) mxWrapGetScalar(prhs[8], &mw_err_txt_);
dim9_ = (mwSize) mxWrapGetScalar(prhs[9], &mw_err_txt_);
dim10_ = (mwSize) mxWrapGetScalar(prhs[10], &mw_err_txt_);
dim11_ = (mwSize) mxWrapGetScalar(prhs[11], &mw_err_txt_);
dim12_ = (mwSize) mxWrapGetScalar(prhs[12], &mw_err_txt_);
dim13_ = (mwSize) mxWrapGetScalar(prhs[13], &mw_err_txt_);
if (mxGetM(prhs[0])*mxGetN(prhs[0]) != dim8_) {
mw_err_txt_ = "Bad argument size: beta"; goto mw_err_label;
}
if (mxGetM(prhs[1])*mxGetN(prhs[1]) != dim9_) {
mw_err_txt_ = "Bad argument size: nterms"; goto mw_err_label;
}
if (mxGetM(prhs[2])*mxGetN(prhs[2]) != dim10_) {
mw_err_txt_ = "Bad argument size: m1"; goto mw_err_label;
}
if (mxGetM(prhs[3])*mxGetN(prhs[3]) != dim11_) {
mw_err_txt_ = "Bad argument size: m2"; goto mw_err_label;
}
if (mxGetM(prhs[5])*mxGetN(prhs[5]) != dim12_) {
mw_err_txt_ = "Bad argument size: lmp"; goto mw_err_label;
}
if (mxGetM(prhs[7])*mxGetN(prhs[7]) != dim13_) {
mw_err_txt_ = "Bad argument size: lmpn"; goto mw_err_label;
}
if (mxGetM(prhs[0])*mxGetN(prhs[0]) != 0) {
in0_ = mxGetPr(prhs[0]);
} else
in0_ = NULL;
if (mxGetM(prhs[1])*mxGetN(prhs[1]) != 0) {
in1_ = mxWrapGetArray_int(prhs[1], &mw_err_txt_);
if (mw_err_txt_)
goto mw_err_label;
} else
in1_ = NULL;
if (mxGetM(prhs[2])*mxGetN(prhs[2]) != 0) {
in2_ = mxWrapGetArray_int(prhs[2], &mw_err_txt_);
if (mw_err_txt_)
goto mw_err_label;
} else
in2_ = NULL;
if (mxGetM(prhs[3])*mxGetN(prhs[3]) != 0) {
in3_ = mxWrapGetArray_int(prhs[3], &mw_err_txt_);
if (mw_err_txt_)
goto mw_err_label;
} else
in3_ = NULL;
if (mxGetM(prhs[4])*mxGetN(prhs[4]) != 0) {
in4_ = mxWrapGetArray_dcomplex(prhs[4], &mw_err_txt_);
if (mw_err_txt_)
goto mw_err_label;
} else
in4_ = NULL;
if (mxGetM(prhs[5])*mxGetN(prhs[5]) != 0) {
in5_ = mxWrapGetArray_int(prhs[5], &mw_err_txt_);
if (mw_err_txt_)
goto mw_err_label;
} else
in5_ = NULL;
if (mxGetM(prhs[6])*mxGetN(prhs[6]) != 0) {
in6_ = mxWrapGetArray_dcomplex(prhs[6], &mw_err_txt_);
if (mw_err_txt_)
goto mw_err_label;
} else
in6_ = NULL;
if (mxGetM(prhs[7])*mxGetN(prhs[7]) != 0) {
in7_ = mxWrapGetArray_int(prhs[7], &mw_err_txt_);
if (mw_err_txt_)
goto mw_err_label;
} else
in7_ = NULL;
if (mexprofrecord_)
mexprofrecord_[3]++;
MWF77_rotviarecur3f90(in0_, in1_, in2_, in3_, in4_, in5_, in6_, in7_);
plhs[0] = mxCreateDoubleMatrix(mxGetM(prhs[4]), mxGetN(prhs[4]), mxCOMPLEX);
mxWrapCopy_dcomplex(plhs[0], in4_, mxGetM(prhs[4])*mxGetN(prhs[4]));
plhs[1] = mxCreateDoubleMatrix(mxGetM(prhs[6]), mxGetN(prhs[6]), mxCOMPLEX);
mxWrapCopy_dcomplex(plhs[1], in6_, mxGetM(prhs[6])*mxGetN(prhs[6]));
mw_err_label:
if (in1_) mxFree(in1_);
if (in2_) mxFree(in2_);
if (in3_) mxFree(in3_);
if (in4_) mxFree(in4_);
if (in5_) mxFree(in5_);
if (in6_) mxFree(in6_);
if (in7_) mxFree(in7_);
if (mw_err_txt_)
mexErrMsgTxt(mw_err_txt_);
}
