void fxSerializeJSON(txMachine* the, txJSONSerializer* theSerializer)
{
txSlot* aSlot;
txInteger aFlag;
aSlot = fxGetInstance(the, mxThis);
theSerializer->offset = 0;
theSerializer->size = 1024;
theSerializer->buffer = c_malloc(1024);
if (!theSerializer->buffer)
mxUnknownError("out of memory");
if (mxArgc > 1) {
aSlot = mxArgv(1);
if (mxIsReference(aSlot)) {
aSlot = fxGetInstance(the, aSlot);
if (mxIsFunction(aSlot))
theSerializer->replacer = mxArgv(1);
else if (mxIsArray(aSlot))
mxSyntaxError("not yet implememented");
}
}
if (mxArgc > 2) {
aSlot = mxArgv(2);
if (mxIsReference(aSlot)) {
aSlot = fxGetInstance(the, aSlot);
if (mxIsNumber(aSlot) || mxIsString(aSlot))
aSlot = aSlot->next;
}
if ((aSlot->kind == XS_INTEGER_KIND) || (aSlot->kind == XS_NUMBER_KIND)) {
txInteger aCount = fxToInteger(the, aSlot), anIndex;
if (aCount < 0)
aCount = 0;
else if (aCount > 10)
aCount = 10;
for (anIndex = 0; anIndex < aCount; anIndex++)
theSerializer->indent[anIndex] = ' ';
}
else if (mxIsStringPrimitive(aSlot))
c_strncpy((char *)theSerializer->indent, aSlot->value.string, 10);
}
theSerializer->stack = the->stack;
mxPush(mxObjectPrototype);
fxNewObjectInstance(the);
aFlag = 0;
if (mxArgc > 0)
mxPushSlot(mxArgv(0));
else
mxPushUndefined();
mxPush(mxEmptyString);
fxSerializeJSONProperty(the, theSerializer, &aFlag);
the->stack++;
}
static Matrix *
SetNum (
Matrix * mat
)
{
Matrix * result = NULL;
int status;
if (mxIsString(mat)) {
mexSetTrapFlag(1);
status = mexCallMATLAB(1, & result, 1, & mat, "abs");
if (status == 1) {
result = NULL;
}
mexSetTrapFlag(0);
}
return (result);
}
void fx_String(txMachine* the)
{
if ((mxThis->kind == XS_REFERENCE_KIND) && ((mxThis->value.reference->flag & XS_SHARED_FLAG) == 0)) {
txSlot* anInstance = mxThis->value.reference;
if (mxIsString(anInstance)) {
if (mxArgc > 0) {
txString aString = fxToString(the, mxArgv(0));
anInstance->next->value.string = aString;
anInstance->next->next->value.integer = fxUnicodeLength(aString);
}
*mxResult = *mxThis;
}
return;
}
if (mxArgc > 0) {
*mxResult = *mxArgv(0);
fxToString(the, mxResult);
}
else {
mxResult->value.string = mxEmptyString.value.string;
mxResult->kind = mxEmptyString.kind;
}
}
void mexFunction(int nlhs, /* Num return vals on lhs */
Matrix *plhs[], /* Matrices on lhs */
int nrhs, /* Num args on rhs */
Matrix *prhs[] /* Matrices on rhs */
)
{
double *image,*filt, *temp, *result;
int x_fdim, y_fdim, x_idim, y_idim;
int x_rdim, y_rdim;
int x_start = 1;
int x_step = 1;
int y_start = 1;
int y_step = 1;
int x_stop, y_stop;
Matrix *arg;
double *mxMat;
char edges[15] = "reflect1";
if (nrhs<2) mexErrMsgTxt("requres at least 2 args.");
/* ARG 1: IMAGE */
arg = prhs[0];
if notDblMtx(arg) mexErrMsgTxt("IMAGE arg must be a non-sparse double float matrix.");
image = mxGetPr(arg);
x_idim = (int) mxGetM(arg); /* X is inner index! */
y_idim = (int) mxGetN(arg);
/* ARG 2: FILTER */
arg = prhs[1];
if notDblMtx(arg) mexErrMsgTxt("FILTER arg must be non-sparse double float matrix.");
filt = mxGetPr(arg);
x_fdim = (int) mxGetM(arg);
y_fdim = (int) mxGetN(arg);
if ((x_fdim > x_idim) || (y_fdim > y_idim))
{
mexPrintf("Filter: [%d %d], Image: [%d %d]\n",x_fdim,y_fdim,x_idim,y_idim);
mexErrMsgTxt("FILTER dimensions larger than IMAGE dimensions.");
}
/* ARG 3 (optional): EDGES */
if (nrhs>2)
{
if (!mxIsString(prhs[2]))
mexErrMsgTxt("EDGES arg must be a string.");
mxGetString(prhs[2],edges,15);
}
/* ARG 4 (optional): STEP */
if (nrhs>3)
{
arg = prhs[3];
if notDblMtx(arg) mexErrMsgTxt("STEP arg must be a double float matrix.");
if (mxGetM(arg) * mxGetN(arg) != 2)
mexErrMsgTxt("STEP arg must contain two elements.");
mxMat = mxGetPr(arg);
x_step = (int) mxMat[0];
y_step = (int) mxMat[1];
if ((x_step<1) || (y_step<1))
mexErrMsgTxt("STEP values must be greater than zero.");
}
void fx_JSON_stringify(txMachine* the)
{
txSlot* aSlot;
volatile txJSONSerializer aSerializer;
txFlag aFlag;
mxTry(the) {
aSlot = fxGetInstance(the, mxThis);
if (aSlot->flag & XS_SANDBOX_FLAG)
the->frame->flag |= XS_SANDBOX_FLAG;
else
the->frame->flag |= the->frame->next->flag & XS_SANDBOX_FLAG;
c_memset((txJSONSerializer*)&aSerializer, 0, sizeof(aSerializer));
aSerializer.offset = 0;
aSerializer.size = 1024;
aSerializer.buffer = c_malloc(1024);
if (!aSerializer.buffer)
fxThrowError(the, XS_RANGE_ERROR);
if (mxArgc > 1) {
aSlot = mxArgv(1);
if (mxIsReference(aSlot)) {
aSlot = fxGetInstance(the, aSlot);
if (mxIsFunction(aSlot))
aSerializer.replacer = mxArgv(1);
else if (mxIsArray(aSlot))
mxDebug0(the, XS_SYNTAX_ERROR, "not yet implememented");
}
}
if (mxArgc > 2) {
aSlot = mxArgv(2);
if (mxIsReference(aSlot)) {
aSlot = fxGetInstance(the, aSlot);
if (mxIsNumber(aSlot) || mxIsString(aSlot))
aSlot = aSlot->next;
}
if ((aSlot->kind == XS_INTEGER_KIND) || (aSlot->kind == XS_NUMBER_KIND)) {
txInteger aCount = fxToInteger(the, aSlot), anIndex;
if (aCount < 0)
aCount = 0;
else if (aCount > 10)
aCount = 10;
for (anIndex = 0; anIndex < aCount; anIndex++)
aSerializer.indent[anIndex] = ' ';
}
else if (aSlot->kind == XS_STRING_KIND)
c_strncpy((char *)aSerializer.indent, aSlot->value.string, 10);
}
aSerializer.stack = the->stack;
*(--the->stack) = mxObjectPrototype;
fxNewObjectInstance(the);
aFlag = 0;
if (mxArgc > 0)
*(--the->stack) = *mxArgv(0);
else
mxZeroSlot(--the->stack);
*(--the->stack) = mxEmptyString;
fxSerializeJSONProperty(the, (txJSONSerializer*)&aSerializer, &aFlag);
the->stack++;
if (aSerializer.offset) {
mxResult->value.string = (txString)fxNewChunk(the, aSerializer.offset + 1);
c_memcpy(mxResult->value.string, aSerializer.buffer, aSerializer.offset);
mxResult->value.string[aSerializer.offset] = 0;
mxResult->kind = XS_STRING_KIND;
}
c_free(aSerializer.buffer);
}
mxCatch(the) {
if (aSerializer.buffer)
c_free(aSerializer.buffer);
fxJump(the);
}
}
void
mexFunction(int nlhs, Matrix **plhs, int nrhs, Matrix **prhs)
{
FILE *nl;
char *buf1, buf[512], *what;
static fint n, nc, nz;
fint nerror;
real *J1, *W, *c, *f, *g, *v, *t, *x;
static real *J;
cgrad *cg, **cgp;
static size_t Jsize;
Jmp_buf err_jmp0;
ASL_pfgh *asl = (ASL_pfgh*)cur_ASL;
static fint nhnz;
static real *Hsp;
real *H, *He;
int *Ir, *Jc;
fint *hcs, *hr, i;
if (nrhs == 1 && mxIsString(prhs[0])) {
if (nlhs != 6)
usage();
if (mxGetString(prhs[0], buf1 = buf, sizeof(buf)))
mexErrMsgTxt("Expected 'stub' as argument\n");
at_end();
mexAtExit(at_end);
asl = (ASL_pfgh*)ASL_alloc(ASL_read_pfgh);
return_nofile = 1;
if (!(nl = jac0dim(buf1,strlen(buf)))) {
sprintf(msgbuf, "Can't open %.*s\n",
sizeof(msgbuf)-20, buf);
mexErrMsgTxt(msgbuf);
}
if (n_obj <= 0)
printf("Warning: objectve == 0\n");
n = n_var;
nc = n_con;
nz = nzc;
J = (real *)M1alloc(nz*sizeof(real));
X0 = mxGetPr(plhs[0] = mxCreateFull(n, 1, REAL));
LUv = mxGetPr(plhs[1] = mxCreateFull(n, 1, REAL));
Uvx = mxGetPr(plhs[2] = mxCreateFull(n, 1, REAL));
pi0 = mxGetPr(plhs[3] = mxCreateFull(nc, 1, REAL));
LUrhs = mxGetPr(plhs[4] = mxCreateFull(nc, 1, REAL));
Urhsx = mxGetPr(plhs[5] = mxCreateFull(nc, 1, REAL));
pfgh_read(nl, ASL_findgroups);
Jsize = nc*n*sizeof(real);
/* Arrange to compute the whole sparese Hessian */
/* of the Lagrangian function (both triangles). */
nhnz = sphsetup(0, 0, nc > 0, 0);
Hsp = (real *)M1alloc(nhnz*sizeof(real));
return;
}
if (!filename)
mexErrMsgTxt("spamfunc(\"stub\") has not been called\n");
nerror = -1;
err_jmp1 = &err_jmp0;
if (nlhs == 2) {
if (nrhs != 2)
usage();
x = sizechk(prhs[0],"x",n);
t = sizechk(prhs[1],"0 or 1", 1);
if (t[0] == 0.) {
f = mxGetPr(plhs[0] = mxCreateFull(1, 1, REAL));
c = mxGetPr(plhs[1] = mxCreateFull(nc, 1, REAL));
if (setjmp(err_jmp0.jb)) {
sprintf(msgbuf, "Trouble evaluating %s\n",
what);
mexErrMsgTxt(msgbuf);
}
what = "f";
*f = objval(0, x, &nerror);
what = "c";
conval(x, c, &nerror);
return;
}
g = mxGetPr(plhs[0] = mxCreateFull(n, 1, REAL));
J1 = mxGetPr(plhs[1] = mxCreateSparse(nc, n, nz, REAL));
what = "g";
objgrd(0, x, g, &nerror);
if (nc) {
what = "J";
jacval(x, J1, &nerror);
Ir = mxGetIr(plhs[1]);
memcpy(mxGetJc(plhs[1]), A_colstarts, (n+1)*sizeof(int));
cgp = Cgrad;
for(i = 0; i < nc; i++)
for(cg = *cgp++; cg; cg = cg->next)
Ir[cg->goff] = i;
}
return;
}
if (nlhs == 0 && nrhs == 3) {
/* eval2('solution message', x, v): x = primal, v = dual */
if (!mxIsString(prhs[0]))
usage();
x = sizechk(prhs[1],"x",n);
v = sizechk(prhs[2],"v",nc);
if (mxGetString(prhs[0], buf, sizeof(buf)))
mexErrMsgTxt(
"Expected 'solution message' as first argument\n");
write_sol(buf, x, v, 0);
return;
}
if (nlhs != 1 || nrhs != 1)
usage();
v = sizechk(prhs[0],"v",nc);
W = mxGetPr(plhs[0] = mxCreateSparse(n, n, nhnz, REAL));
what = "W";
sphes(H = Hsp, 0, 0, v);
/* Expand the Hessian lower triangle into the full Hessian... */
Ir = mxGetIr(plhs[0]);
Jc = mxGetJc(plhs[0]);
hcs = sputinfo->hcolstarts;
hr = sputinfo->hrownos;
for(i = 0; i <= n; i++)
Jc[i] = hcs[i];
He = H + hcs[n];
while(H < He) {
*W++ = *H++;
*Ir++ = *hr++;
}
}
