Description

Check Convert's interval evaluator

Contributor: Steve Dirkse

Small Model of Type : GAMS

Category : GAMS Test library

Main file : convert10.gms

$title CONVERT test suite - check interval evaluator in Convert (CONVERT10,SEQ=668)

$onText
Check Convert's interval evaluator

Contributor: Steve Dirkse
$offText


variables
  x  / lo 4, up 4 /
  y  / lo 8, up 8 /
  z
  ;
equations
  f  'e1'
  g  'e2'
  h  'e3'
  ;

f .. sqr(x-2) + sqr(y-2) =e= 2*z;
g .. x + y =L= 1;
h .. exp(.25*x) + exp(.25*y) =G= -1;

model m / f, g, h /;

execute_loadpoint 'savepoint10';

$set JAC  jac10.gdx
$set JAC_ jac10_.gdx

$onEcho > convert.opt
dumpgdx %JAC%
gdxhessian 1
gdxquadratic 1
intervalEval yes
$offEcho
option nlp = convert;
m.optfile = 1;

execute 'rm -f %JAC%';
solve m using nlp min z;

execute '=test -e %JAC%';
abort$errorlevel 'we did not find Jacobian file %JAC% after Convert ran';
execute 'gdxdiff %JAC% %JAC_% relEps=2e-15 > %system.nullfile%';
abort$errorlevel 'Unexpected result diffing %JAC% with %JAC_%';