Mosel

Transcription

1

2 Mosel Mosel public Mosel Mosel Mosel Mosel Mosel i

3 ii A 231 A A A A B Mosel 248 B B

4 1 1.1 Xpress-Mosel Xpress-MP Mosel Xpress-Modeler Xpress-EMOSL Xpress-Mosel ODBC Mosel Xpress-MP Mosel Mosel Mosel Xpress-MP Xpress-IVE Mosel Mosel Mosel Mosel 1

5 2 1 Mosel Xpress-MP C Java Visual Basic Mosel Xpress-MP Essentials Xpress-MP 1.2 Mosel Mosel Mosel 1 Mosel Mosel Mosel Mosel Mosel 1.3 Mosel / Mosel / BInary Model(BIM) 2 / / BIM Mosel Mosel : VIrtual MAchine (VIMA) / ( ) ( )

6 : VIMA : mosel 2 : Mosel mmxprs Mosel Xpress-Optimizer maximize Mosel BIM Mosel / BIM Mosel Mosel Mosel Mosel 1.4 Mosel 2 3 Mosel ( ) Java ( Prolog Lisp) Mosel Pascal ( sum )

7 4 1 Mosel mpmodel Mosel Mosel 3 Mosel Mosel 4 Mosel Mosel Mosel 1.6 fixed width italic C 5 [...] {...} italic return CTRL CTRL UNIX PC ASCII 16 7C 8 124

8 2 Mosel Mosel Mosel Mosel 2.1! 5

9 6 2 MOSEL (! (! 2!)!) word Word and, array, as, boolean, break, case, declarations, div, do, mpvar, dynamic, elif, else, end, false, forall, forward, from, function, if, in, include, initialisations, initializations, integer, inter, is_binary, is_continuous, is_free, is_integer, is_partint, is_semcont, is_semint, is_sos1, is_sos2, linctr, max, min, mod, model, next, not, of, options, or, parameters, procedure, public, prod, range, real, repeat, set, string, sum, then, to, true, union, until, uses, while Mosel AND and And.. A+B! 1 -C+D! 2 A+B-! 3... C+D! 3

10 A+B ; -C+D!2 explterm(2.3 ) A+B;! 1 -C+D;! 2 A+B! C+D;! 3 word word todo todo [something] something {something...} something 2.2 Mosel model model name [Directives] [Parameters] [body] end-model model end-model model name Mosel directives parameters / body

11 8 2 MOSEL BIM 2.3 uses options 2 uses uses libname1 [,libname2...][;] options options optname1 [,optname2...] explterm noimplicit ( ) uses "mmsystem" options noimplicit,explterm 2.4 (XPRMrunmod(RUN) ) parameters ident1 = Expression1 [ident2 = Expression2...] end-parameters identi Expressioni XPRMrummod(RUN)

12 parameters size=12 F="myfile" end-parameters 2.5 Mosel include include filename filename.mos include filename a.mos model "Example for file inclusion" written ( From the main file ) include "b" end-model b.mos written ( From an included file ) b.mos a.mos model "Example for file inclusion" written ( From the main file ) written ( From an included file ) end-model uses options parameters 2.6 / 1 declarations end-declarations declarations Declare stat

13 10 2 MOSEL [Declare stat...] end-declarationss ident1 [, ident2..] : typename typename identi integer,real,string,boolean MP mpvar,linctr MP declarations i,j:integer str:string x,y,z:mpvar end-declarations integer real 1.7e e string boolean true false 0 false MP 2 mpvar linctr

14 ident1 [,ident2...]:set of typename typename identi ident1 [, ident2...]:range set of integer identi declarations sl: set of string rl: range end-declarations, ident1 [, ident2...] : [dynamic] array (list of sets) of typename list of sets / typename identi rs:set of real rs set of real declarations e: set of string t1: array (e, rs:set of real, range, integer) of real

15 12 2 MOSEL t2: array ({"i1","i2"}, 1..3) of integer end-declarations finalize dunamic Mosel mpvar mpvar create identifier = Expression identifier Expression declarations STR= my const string I1=12 R=1..10 S={2.3,5.6,7.01} end-declarations 2.7 array ident (Exp1 [, Exp2...])

16 array ident Expi i Exp1[,Exp2...] function ident function ident( Exp1 [, Exp2...] ) function ident expi i if Bool expr Exp1 Exp2 if(bool expr, Exp1, Exp2) Exp1 Exp2 Exp1 Exp2 Mosel integer -> real integer -> real,linctr real ->linctr mpvar ->linctr integer -> real,string,linctr real ->string,linctr mpvar ->linctr boolean -> string integer, real, string, boolean string boolean string string true MP!A=3.5,B=2 integer(a+b)!=5 string(a-b)!= 1.5 real(integer(a+b))!=5.5 ( ) ()

17 14 2 MOSEL Aggregate ident (Iterator1 [, Iterator2 ]) Expression Aggregate indent Expression ( ) Expression iterator Set expr ident1 [,ident2...]in Set expr[ Bool expr] ident = Expression [ Bool expr] 2 identi Set expr 3 identi 2 3 iterator Bool expr iterator (+ / + ) mod div mpvar sum integer real mpvar prod minimun maximum x*5.5+(2+z)^4+cos(12.4) sum(i in 1..10) (min(j in s) t(i)*(a(j) mod 2))

18 string string C a, b, f, n, r, t, v n C ddd1 ddd2 ddd3 C ddd1 ddd2 ddd3 C ddd1 ddd2 ddd3 C ddd1ddd2ddd3 string a1b2c3d5 e6 a1b2c3d5e6 a1b2c3d5 3d5 a1b2c {[ Exp1 [, Exp2 ]]} [{] Integer exp1.. Integer exp2 [}] 2 union inter {1,2,3}+{4,5,6}-{5..8}*{6,10}! = {1,2,3,4,5} { a, b, c }*{ b, c, d }! = { b, c } union(i in 1..4 i<>2) {i*3}! = {3,9,12} string ISO

19 16 2 MOSEL expr in Set expr Set Expr not in and ( ) or ( ) not and or 3<=x and y>=45 or t<>r and not r in {1..10} and(i in 1..10) 3<=x(i) mpvar Linear expr Linear expr1 Ctr cmp Linear expr2 Linear Expr SOS type mpvar ref mpvar type1 mpvar ref mpvar type2 Arith expr 2 3 is sos1 is sos2 2 1 is continuous is integer is binary is free 2 is partint is semicont is semint 3*y+sum(i in 1..10) x(i)*i >= z-t t is_integer sum(i in 1..10) x(i) is_sos1 y is_partint 5

20 *x>=5*y-10 3*x-5*y+10 ctl 3*x>=5*y-10 z-x+ctl z-2*x-5*y+10 x is type threshold x-threshold x is integer x-max INT [Exp1 [, Exp2...]] declarations T:array(1..2,1..3) of integer end-declarations T:=[1,2,3,4,5,6] writeln(t)! [1,2,3,4,5,6] writeln(getsize([2,3,4]))! Mosel 4

21 18 2 MOSEL ident ref := Expression ident ref + = Expression ident ref = Expression ident ref Expression ident ref ident ref c:= 3*x+y >= 5 c+= y! c 3*x+2*y-5 >= 0 c:= 3*x+y >= 5 c:= c + y! c 3*x+2*y-5 (c ) T:=[2,4,6,8] T(12):=[7,8,9,19]!A,S,SE A:=1!A S:={1,2,3}!S SE:={}! SE Mosel A(1):=1! A (1,..1) A(t):=2.5! t 1 10 f(t)>0! A noimplicit

22 <=x; x<=20 x is_integer procedure ident procedure ident (Exp1 [, Exp2 ]) procedure ident expi i?? writeln("hello!")! hello! initializations from Filename ident1 [ as Label1] [ident2 [ as Label2]...] end-initializations Filename identi string Labeli Label: value Label value integer,real,string,boolean []

23 20 2 MOSEL false true initialization to ident1 [as Label1] [ident2 [as Label2]...] end-initialization Filename 2 1 a.dat! Example of the use of initialization blocks t:[ (1 un) 10 (2 deux) 20 (3 trois) 30 ] t2:[ 10 (4) ] nb used : 0 model "Example initblk" declarations nb_used:integer s: set of string t: array(1..3,s) of real t2: array(1..5) of integer end-declarations initializations from a.dat t! t=[(1, un,10),(2, deux,20),(3, trois,30)] t2! t2=[10,0,0,30,40] 1 fname fname-

24 nb_used as "nb used"! nb_used=0 end-initializations nb_used+=1 t(2,"quatre"):=1000 initializations to t nb_used as "nb used" s end-initializations end-model a.dat! Example of the use of initialization blocks t:[(1 un ) 10 (2 deux ) 20 (2 quatre ) 1000 (3 trois ) 30] t2:[ 10 (4) ] nb used : 1 s : [ un deux trois quatre ] If if if Bool exp 1 then Statement list 1 [ elif Bool exp 2 then Statement list 2 ] [else Statement list E ] end if Bool exp 1 true Statement list 1 end-if elif if false else Statement list E end-if if c=1 then writeln( c=1 ) elif c=2 then writeln( c=2 ) else writeln( c<>1 and c<>2 )

25 22 2 MOSEL end-if Case Case case Expression 0 of Expression 1 : Statement 1 or Expression 1 : do Statement list 1 end-do [ Expression 2 : Statement 2 or Expression 2 : do Statement list 2 end-do ] [else Statement list E ] end case Expression 0 Expression i Statement i Statement list i end case else Statement list E end-if Expression i 2 Expression 0 case c of 1 : writeln( c=1 ) 2..5 : writeln( c in 2..5 ) 6,8,10: writeln( c in {6,8,10} ) else writeln( c in {7,9} or c >10 or c <1 ) end-case forall forall forall(iterator list) Statement forall(iterator list) do Statement list end-do Statement Statement list ( )

26 forall (i in 1..10,j in i<>j) do write( (,i,,,j, ) ) if isodd(i*j) then s+={i*j} end-if end-do While While while(bool expr) Statement while(bool expr) do Statement list end-do Statement Statement list, Bool expr true false while i:=1 while(i<=10) do write(,i) if isodd(i) then s+={i} end-if i+=1 end-do Repeat repeat repeat Statement1 [Statement2...] until Bool expr repeat until Bool expr true while i:=1 repeat write(,i) if isodd(i) then s+={i} end-if i+=1 until i>10

27 24 2 MOSEL break next break next break next break [n] next [n] n n-1! repeat!1: L1 forall (i in S) do!2: L2 while (C3) do!3: L3 break 3!4:L3 11 next!5:l3 (3 ) next 2!6:L3 i (2 ) end-do!7:l3 next 2!8:L2 L1 (11 ) break!9:l2 10 end-do!10:l2 until C1!11:L1 2.9 Mosel 2 ( writeln("hi!") ) ( round(12.3) ) ( ) ( ) procedure name proc [(list of args)] Proc body end-procedure

28 name proc list of args 2.6 ( A-12 t1:array(1..4) of real ) procedure myproc writeln("in myproc") end-procedure procedure withparams(a:array(r:range) of real, i,j:integer) writeln("i received: i=",i," j=",j) forall(n in r) writeln("a(",n,")=",a(n)) end-procedure declarations mytab:array(1..10) of real end-declarations myproc! myproc withparams(mytab,23,67)! withparams function name func [(list of args)]: Basic type Func body end-function Mosel integer,real,boolean,string returned function multiply_by_3(i:integer):integer returned:=i*3 end-function writeln("3*12=",multiply_by_3(12))!

29 26 2 MOSEL procedure alter(s:set of integer,i:integer) i+=1 s+={i} end-procedure gs:={1} gi:=5 alter(gs,gi) writeln(gs," ",gi)! displays: {1,6} 5 ( end-procedure end-function ) (mpvar) ( linctr) declarations! i,j:integer end-declarations procedure myproc declarations i:string! end-declarations i:="a string"! i j:=4 writeln("inside of myproc, i=",i," j=",j) end-procedure i:=45! i j:=10

30 myproc writeln("outside of myproc, i=",i," j=",j)! displays:! myproc, i=a string j=4! myproc, i=45 j=4 Mosel! 1 function random(limit:integer):integer returned:=round(.5+random*limit)! random end-function 1 2 ( ) 2 1 forwarf procedure Proc name[(list of prams)] forwarf procedure Func name[(list of prams)] :Basic type Proc name Func name Mosel forward function f2(x:integer):integer function f1(x:integer):integer returned:=x+if(x>0,f2(x-1),0)! f1 needs to know f2 end-function

31 28 2 MOSEL function f2(x:integer):integer returned:=x+if(x>0,f1(x-1),0)! f2 needs to know f1 end-function 2.10 public Mosel ( DISPLAY ) BIM -s public ( ) parameters public T="default"! T end-parameters declarations public a,b,c:integer! a,b,c d:real! d end-declarations forward public procedure myproc(i:integer)! myproc 2.11 / 2 (write, writeln,fflush) (read, readln, fskipline) Mosel fopen fclose ( ) fselect getfid model "test IO" def_out:=getfid(f_output)! ID

32 fopen("mylog.txt",f_output)! my_out:=getfid(f_output)! ID repeat fselect(def_out)! write("text? ")! text:= readln(text)! string fselect(my_out)! mylog.txt writeln(text)! string until text= fclose(f_output)! ( mylog.txt ) writeln("finished!")! end-model 2.12 Mosel ( ) Mosel ( ) mmive, mmodbc, mmsystem, mmxprs Mosel uses module name 6 7 mmsystem mmsystem OS mmsystem getsysstat mmsystem

33 30 2 MOSEL fdelete fmove getcwd getenv getfstat getsysstat gettime makedir removedir system mmxprs OS mmxprs Mosel Xpress-Optimizer Xpress-Optimizer 4 Mosel ( XPRS_LPSTATUS) getparam setparam XPRS_VERBOSE: / XPRS_LOADNAMES MPS / XPRS_PROBLEM : Mosel setparam("xprs_presolve",0); setparam("xprs_cutstrategy",0); setparam("xprs_extrarows",5000); setcallback Mosel Mosel ( ) / ) mmxprs : mmxprs

34 addcut addcuts clearmipdir clearmodcut delbasis delcuts dropcuts getcnlist getcplist getlb getprobstat getub loadbasis loadcuts loadprob maximize minimize savebasis setcallback setlb setmipdir setmodcut setub storecut storecuts MIP mmodbc Mosel ODBC ODBC ODBC / mydata pricelist

35 32 2 MOSEL articlenum colour price 1001 blue red black blue 3.99 Model model ODBC Example uses "mmodbc" declarations prices: array(range) of real end-declarations setparam("sqlverbose",true) SQLconnect("DSN=mydata") writeln("connection number: ",getparam("sqlconnection")) SQLexecute("select articlenum,price from pricelist",prices) SQLdisconnect end-model SQLVERBOSE ODBC SQLexecute price articlenum pricelist mmodbc : mmdobc SQLconnect SQLdisconnect SQLexecute SQLreadinteger SQLreadreal SQLreadstring ODBC ODBC SQL integer real string mmive mmive Xpress-MPIntegrated Visual Environment IVE 2 IVEaddtograph IVEinitgraph

36

37 3 Mosel 3.1 Mosel Mosel (.mos ) (.bim ) Mosel -h -v -s -c commands Mosel Commands ( OS ) Commands Mosel mosel -c "CLOAD -sg mymodel;run Mosel ( [] ) INFO [symbol]: ( ) SYSTEM [command]: OS QUIT: Mosel COMPILE[-sgep] filename[comment]: filename (BIM).mos.bim -e 34

38 3.1. MOSEL 35 -s BIM -g comment BIM LIST -p LOAD filename: filename BIM.bim model CLOAD[-sge] filename[comment]: filename COMPILE filename LOAD filename LIST: CLOAD LOAD name: model number: size: system: / user comment: (COMPILE,CLOAD ). DELETE,RUN.RESET SELECT number name: DELETE number name: (BIM ) ( ) RUN [parameters]: / param name = value dsoname

39 36 3 MOSEL ctrpar name dsoname.ctrpar name=value RESET: EXPORTPROB[-[ms] [filename[objective] :] filename LP MPS (-m ) -p -m MPS -s objective DISPLAY symbol: ( 0) ( 0) SYMBOLS [-cspo]: -c -s / -p / -o LSLIBS: (DSO) ( EXAMINE [-cspt] libname: libname / / -c -s / -t -p / FLUSHLIBS: CL CLOAD C COMPILE 1 / CLOAD mymodel RUN Xpress-MP

40 3.2. MOSEL Mosel Mosel C Mosel (BIM) BIM Mosel Mosel Mosel (Mosel ) 2 Mosel Xpress-MP Mosel xprm rt BIM xprm rt.h XRPMMmodel: XPRMdsolib: XPRMmpvar: XPRMlinctr: XPRMset: XPRMarray: XPRMproc: / XPRMinteger: (C int ) XPRMreal: (C double ) XPRMboolean: (C int: 0 = false, 1 = true) XPRMstring: (C const char * )

41 38 3 MOSEL Mosel (xprm mc) Mosel BIM xprm mc XPRMcompmod xprm rt xprm mc.h 3.3 Mosel Mosel XPRMinit Mosel Mosel XPRMfree (XPRMloadmod) (XPRMrunmod) (XPRMunloadmod) (XPRMgetmodinfo) Mosel (MPRMmodel ) XPRMloadmod BIM( ) XPRMmodel 2 (pp) Mosel BIM XPRMfindident XPRMgetnextparam XPRMgetnextproc ). s public (XPRMgetsettype) (XPRMgetsetsize) (XPRMgetelsetval, XPRMgetelsetndx, XPRMfirstsetndx,XPRMlastsetndx) Mosel (XPRMgetarrsize) (XPRMgetarrdim) (XPRMfirstarrentry, XPRMnextarrentry,XPRMlastarrentry) (XPRMgetarrval) 2 BIM Mosel

42 3.3. MOSEL 39 XPRMgetarrtype Mosel XPRM TYP INT, XPRM TYP REAL, XPRM TYP STRING, XPRM TYP BOOL, XPRM TYP MPVAR, XPRM TYP LINCTR, XPRM TYP SET, XPRM TYP ARR 2 XPRM GRP DYN: out of range error XPRM GRP GEN: ( ) XPRM GRP DYN XPRM GRP DYN XPRM GRP GEN ( ) Mosel XPRM GRP DYN Mosel 2 XPRMfirstarrtruentry,XPRMnextarrtruentry XPRMgetelsetndx main problem Mosel main problem XPRMgetobjval (XPRMgetvsol) (XPRMgetcsol) (XPRMgetrcost) (XPRMgetdual) (XPRMgetslack) XPRMexportprob mmxprs loadprob mmxprs maximize Mosel DSOs 2.12 Mosel 1 Mosel (XPRMgetnextdso) (XXRPXgetdsoinfo) (XPRMgetnextdsoparam)

43 40 3 MOSEL XPRMpreloaddso XPRMflushdso Mosel (xprm mc) XPRMcompmod Mosel xprm mc, xprm rt. XPRMcompmod XPRMinit Mosel 3.4 XPRMinit XPRMfree XPRMgetversion Mosel Mosel Mosel XPRMloadmod XPRMrunmod XPRMisrunmod XPRMstoprunmod XPRMunloadmod XPRMgetmodinfo XPRMgetnextmod XPRMfindmod Mosel Mosel XPRMfindident XPRMgetnextident XPRMgetnextparam XPRMgetnextproc XPRMgetprocinfo / / XPRMgetsetsize XPRMgetsettype XPRMfirstsetndx XPRMlastsetndx XPRMgetelsetval XPRMgetelsetndx

44 XPRMgetarrdim XPRMgetarrtype XPRMgetarrsize XPRMgetarrsets XPRMfirstarrentry XPRMlastarrentry XPRMnextarrentry XPRMnextarrentry trs XPRMfirstarrtruenty XPRMnextarrtruentry XPRMchkarrind XPRMcmpindices XPRMgetarrval 2 XPRMgetprobstat XPRMexportprob XPRMgetobjval XPRMgetobjval XPRMgetcsol XPRMgetrcost XPRMgetdual XPRMgetslack XPRMgetact XPRMgetvarnum XPRMgetctrnum XPRMautounloaddso XPRMfinddso XPRMflushdso XPRMgetdsoparam XPRMgetnextdso XPRMgetnextdsoconst XPRMgetnextdsoparam XPRMgetnextdsoproc XPRMgetdsoinfo XPRMpreloaddso DSO / DSO DSO / DSO DSO DSO DSO / DSO DSO

45 42 3 MOSEL XPRMcompmod Mosel BIM

46 4 Mosel 4.1 Mosel 2 Mosel Mosel Mosel Mosel MAX_INT integer MAX_REAL real M_E e M_PI π 4.2 / Mosel 43

47 44 4 MOSEL uses module name 1 2

48 abs integer real function abs(i: integer): integer function abs(r: real): real i r declarations number: real end-declarations writeln("enter a non-negative number.") readln(number) if(number <> abs(number)) then writeln("the number is negative.") end-if exp,ln,log,sqrt

49 46 4 MOSEL addcut procedure addcut(cuttype: integer, type: integer,linexp: linctr) cuttype type linexp CT_GEQ CT_LEQ CT_EQ mmxprs Boolean VERBOSE / addcuts, dropcuts, delcuts, loadcuts, storecut, storecuts, XPRSaddcuts ( )

50 addcuts procedure addcuts(cuttype: array(range) of integer, type: array(range) of integer, linexp: array(range) of linctr) cuttype type linexp CT_GEQ CT_LEQ CT_EQ mmxprs Boolean VERBOSE / 3 addcut, dropcuts, delcuts, loadcuts, storecut, storecuts, XPRSaddcuts( )

51 48 4 MOSEL arctan function arctan(r: real): real r 2 a b 2 declarations a,b: real end-declaration read(a,b) write("the triangle has angles ",arctan(a/b)*180/m_pi) writeln(", ",arctan(b/a)*180/m_pi," and 90 degrees") cos,sin

52 bittest function bittest(i: integer, mask: integer): integer i mask i 4 j 5 k 8 i := bittest(12,5) j := bittest(13,5) k := bittest(13,10)

53 50 4 MOSEL ceil function ceil(r: real): integer r d -1 d := ceil(-1.9) floor,round

54 clearmipdir MIP procedure clearmipdir mmxprs Boolean VERBOSE / setmipdir, XPRSgetdirs( ),XPRSloaddirs( )

55 52 4 MOSEL clearmodcut procedure clearmodcut mmxprs Boolean VERBOSE / declarations Items: set of integer ctr: array(items) of linctr end-declarations... forall(i in Items isodd(i)) setmodcut(ctr(i)) maximize(profit) writeln("profit: ",getobjval) clearmodcut setmodcut

56 cos function cos(r: real): real r 2 declarations b,c,a: real end-declarations read(b,c,a) write("the other side has length ") writeln(sqrt(b*b + c*c -2*b*c*cos(A*M_PI/180))) arctan,sin

57 54 4 MOSEL create procedure create(x: mpvar) x x(1) + x(3) + x(5) + x(7) declarations x: dynamic array(1..8) of mpvar end-declarations forall(i in 1..8 isodd(i)) create(x(i)) c := sum(i in 1..8) x(i) finalize(2 Mosel )

58 delbasis savebasis procedure delbasis(num: integer) procedure delbasis(name: string) num name mmxprs Boolean VERBOSE / 2 2 savebasis(2)... delbasis(2) loadbasis, savebasis, XPRSgetbasis ( ),XPRSloadbasis( )

59 56 4 MOSEL delcuts procedure delcuts(keepbasis: boolean, cuttype: integer, interpret: integer, delta: real, cuts: set of integer) procedure delcuts(keepbasis: boolean, cuttype: integer, interpret: integer, delta: real) keepbasis cuttype interpret false true cuttype 3 cuttype delta cuts mmxprs Boolean delta -MAX REAL cuttype VERBOSE / addcut, addcuts, dropcuts, loadcuts, storecut, storecuts, XPRSdelcuts ( )

60 dropcuts procedure dropcuts(cuttype: integer, interpret: integer, cuts: set of integer) procedure dropcuts(cuttype: integer, interpret: integer) cuttype interpret cuttype 3 cuttype cuts mmxprs Boolean cuttype VERBOSE / addcut, addcuts, delcuts, loadcuts, storecut, storecuts, XPRSdelcpcuts( )

61 58 4 MOSEL exlsts function exists(x): boolean x ( t(1)) declarations x: dynamic array(1..8) of mpvar end-declarations forall(i in 1..8 isodd(i)=false) create(x(i)) c:= sum(i in 1..8) x(i) forall(i in 1..8) if(exists(x(i)) = true) then writeln("x(",i,") exists") end-if create

62 exit procedure exit(code: integer) code 1 readln(number) if(number > 0) then writeln(log(number)) else exit(1) end-if exit 0 QUIT

63 60 4 MOSEL exp function exp(r: real): real r 2 exp declarations limit = 100 xp = 0 end-declarations read(number) forall(i in 0..limit) do xp := xp + pospow(number,limit-i)/factorial(limit-i) end-do writeln("calculated: ",xp,"; exp(",number,"): ", exp(number)) function factorial(a: integer): real if(a < 0) then returned := -1 else if(a = 1 or a = 0) then returned := 1 else returned := a*factorial(a-1) end-if end-if end-function function pospow(x:real, e:integer): real if(e<=0) then

64 returned := 1 else returned := x*pospow(x,e-1) end-if end-function abs,ln.log,sqrt

65 62 4 MOSEL exportprob procedure exportprob(options: integer, filename: string, obj: linctr) options EP MIN EP MAX EP MPS EP STRIP filename obj LP LP MPS LP prob1.lp Profit: declarations MinCost,Profit: linctr end-declarations exportprob(0,"",mincost) exportprob(ep_max+ep_strip,"prob1", Profit) LP.lp MPS.mps XPRMexportprob (EXPORTPROB), XPRBexportprob

66 fclose procedure fclose(stream: integer) stream F INPUT F OUTPUT Boolean IOCTRL IO / declarations profit: linctr variables: set of mpvar end-declarations... fopen("solution.dat",f_output) writeln("profit: ",getobjval) getvars(profit,variables) forall(x in variables) writeln(getsol(x)) fclose(f_output) fopen fopen, fselect, getfid, iseof, SQLconnect, SQLdisconnect.

67 64 4 MOSEL delete procedure fdelete(filename: string) filename mmsystem fdelete("logfile.log") fopen("logfile.log",f_output)... fclose(f_output) fmove, getfstat

68 fflush OS procedure fflush XPRMflushdso (FLUSHLIBS)

69 66 4 MOSEL finalize procedure finalize(s: set) s 3 declarations Set1: set of string x: array(set1) of mpvar end-declarations Set1 := {"first", "second", "fifth"} finalize(set1) declarations y: array(set1) of mpvar end-declarations create

70 floor function floor(r: real): integer r i 5 j -7 k 12 i := floor(5.6) j := floor(-6.7) k := floor(12.3) ceil,round

71 68 4 MOSEL fmove procedure fmove(namesrc: string, namedest: string) namesrc namedest / mmsystem forall(i in ) do file := "logfile" + i status := getfstat(file) if(bittest(status,sys_typ) = SYS_REG) then next else break end-if end-do fmove("logfile.log",file) fopen("logfile.log",f_output)... namesrc namedest 2 fdelete, getfstat

72 fopen procedure fopen(f: string, mode: integer) f mode F INPUT F OUTPUT F APPEND Boolean IOCTRL IO / fopen("data.dat",f_append) if(getparam(iostatus) <> 0) then exit(1) end-if writeln(getobjval) fclose(f_output) 1. write writeln read readln fskipline 2. IO IOCTRL false false IO LOSTATUS 0 getparam fclose, fselect, getfid, iseof, SQLconnect, SQLdisconnect

73 70 4 MOSEL fselect procedure fselect(stream: integer) stream ID mylog.txt ID def_out := getfid(f_output) fopen("mylog.txt",f_output)... my_out:=getfid(f_output) fselect(def_out) getfid fclose, fopen, getfid, iseof

74 fskipline procedure fskipline(filter: string) filter Boolean IOCTRL IO / fskipline("!#") fskipline("\n") n read,readln

75 72 4 MOSEL getact 0 function getact(c: linctr): real c function check(ctr: linctr, name: string): boolean if((gettype(ctr)=ct_geq or gettype(ctr)=ct_leq) and (getact(ctr)+getslack(ctr)+getcoeff(ctr)<>0)) then writeln("error in constraint ",name," ") returned := false end-if returned := true end-function maximize(profit) if(check(first,"first") and check(second,"second")) then writeln("profit: ",getobjval) end-if getdual, getslack, getsol, XPRMgetact

76 getcnlist procedure getcnlist(cuttype: integer, interpret: integer, cuts: set of integer) cuttype interpret cuttype 3 cuttype cuts mmxprs Boolean VERBOSE / cuts getcplist, XPRSgetcutlist( )

77 74 4 MOSEL getcoeff =-RHS function getcoeff(c: linctr): real function getcoeff(c: linctr, x: mpvar): real c x 3 getcoeff r 1.0 s 12.0 declarations x, y, z: mpvar end-declarations c := 4*x + y - z <= 12 r := getcoeff(c,z) s := getcoeff(c) getvars, setcoeff, XPRSgetrows ( )

78 getcplist procedure getcplist(cuttype: integer, interpret: integer, delta: real, cuts: set of integer, viol: array(range) of real) cuttype interpret cuttype 3 cuttype delta cuts viol mmxprs Boolean delta VERBOSE / cuts getcnlist, XPRSgetcpcutlist( ),XPRSgetcpcuts( )

79 76 4 MOSEL getcwd function getcwd: string mmsystem procedure printdir writeln("executing from ",getcwd) end-procedure Mosel getenv, getfstat, makedir, removedir

80 getdual 0 function getdual(c: linctr): real c first := 3*a + 2*b <= 400 second := a + 3*b <= 200 profit := a + 2*b maximize(profit) writeln("dual value of constraint first is ", getdual(first)); getobjval, getrcost, getslack, getsol, XPRMgetdual, XPRSgetsol ( )

81 78 4 MOSEL getenv OS function getenv(name: string): string name mmsystem PATH str := getenv("path") getcwd, getsysstat, system

82 getfid function getfid(stream: integer): integer stream F INPUT F OUTPUT def out file file def_out := getfid(f_output) fopen("file",f_output) writeln("objective function value: ",getobjval) my_out := getfid(f_output) fselect(def_out) fselect fclose, fopen, fselect, iseof

83 80 4 MOSEL getfirst 0 function getfirst(r: range): integer r r declarations r = 2..8 end-declarations... writeln("first element of r: ", getfirst(r),"\nlast element of r: ",getlast(r)) getlast

84 getfstat function getfstat(filename: string): integer filename mmsystem ftest fstat := getfstat("ftest") if bittest(fstat,sys_typ)=sys_dir then writeln("ftest is a directory") end-if if bittest(fstat,sys_write)=sys_write then writeln("ftest is writeable") end-if SYS TYP( ) SYS DIR SYS REG SYS OTH SYS MOD( ) SYS READ SYS WRITE SYS EXEC bittest

85 82 4 MOSEL getlast 0 function getlast(r: range): integer r r declarations r = 2..8 end-declarations... writeln("first element of r: ", getfirst(r),"\nlast element of r: ",getlast(r)) getfirst

86 getlb function getlb(x: mpvar): real x mmxprs Boolean VERBOSE / x declarations Items = x: array(items) of mpvar end-declarations... forall(i in Items) writeln(i,": ",getlb(x(i)),", ",getub(x(i))) setlb Mosel loadprob getub, setlb, setub, XPRSchgbounds ( ),XPRSgetlb ( )

87 84 4 MOSEL getobjval 0 function getobjval: real profit maximize(profit) writeln("maximum profit obtainable is ",getobjval) LP getdual, getprobstat, getrcost, getslack, getsol, LPOBJVAL ( ),MIPOBJVAL( )

88 getparam function getparam(name: string): integer string real boolean name SQL success := getparam("sqlsuccess") Mosel mmxprs.xprs VERBOSE setparam( 7 8 )

89 86 4 MOSEL getprobstat function getprobstat: integer mmxprs Boolean VERBOSE / getobjval maximize(profit) if(getprobstat = XPRS_OPT) then writeln("maximum profit is ",getobjval) else writeln("problem could not be solved to optimality") end-if getprobstat XPRS OPT XPRS UNF XPRS INF XPRS UNB getparam PRESOLVESTATE, LPSTATUS, MIPSTATUS ( ) XPRMgetprobstat, LPSTATUS ( ),MIPSTATUS ( )

90 getrcost 0 function getrcost(v: mpvar): real v declarations profit: linctr variables: set of mpvar end-declarations... if(getprobstat = XPRS_OPT) then getvars(profit,variables) forall(x in variables) writeln(getsol(x),"; ",getrcost(x)) end-if getdual, getobjval, getprobstat, getslack, getsol, XPRMgetrcost, XPRMgetsol( )

91 88 4 MOSEL getsize function getsize(a: array): integer function getsize(s: set): integer a s 8 2 getsize 2 declarations a: dynamic array(1..8) of real end-declarations a(1) := 4 a(5) := 7.2 l := getsize(a) XPRMgetarrsize, XPRMgetsetsize

92 getslack 0 function getslack(c: linctr): real c declarations ctr: set of linctr end-declarations... if(getprobstat = XPRS_OPT) then forall(c in ctr) if(getslack(c) <> 0) then writeln(getslack(c)) end-if end-if getdual, getobjval, getprobstat, getrcost, getsol, XPRMgetslack, XPRSgetsol ( )

93 90 4 MOSEL getsol 0 function getsol(v: mpvar): real function getsol(c: linctr): real v c declarations profit: linctr variables: set of mpvar end-declarations... if(getprobstat = XPRS_OPT) then getvars(profit,variables) forall(x in variables) writeln(getsol(x)) end-if getdual, getrcost, getobjval, getprobstat, XPRMgetcsol, XPRMgetvsol, XPRSgetsol ( )

94 getsysstat function getsysstat: integer mmsystem randomfile fdelete("randomfile") if(getsysstat <> 0) then writeln("randomfile could not be deleted.") end-if 0 getfstat, IOSTATUS, NBREAD, SQLSUCCESS

95 92 4 MOSEL gettime function gettime: real mmsystem starttime := gettime... write("time: ",gettime-starttime)

96 gettype function gettype(c: linctr): integer c declarations c: set of linctr end-declarations... if(getprobstat = XPRS_OPT) then forall(c in ctr) if(gettype(c) = CT_GEQ or gettype(c) = CT_LET) then writeln(getslack(c)) end-if end-if gettype CT EQ CT GEQ CT LEQ CT UNB CT SOS1 CT SOS2 = 1 2 CT CONT CT INT CT BIN CT PINT CT SEC CT SINT getdual, getslack, settype

97 94 4 MOSEL getub function getub(x: mpvar): real x mmxprs Boolean VERBOSE / x declarations Items = x: array(items) of mpvar end-declarations... forall(i in Items) writeln(i,": ",getlb(x(i)),", ",getub(x(i))) setub Mosel loadprob getlb, setlb, setub, XPRSchgbounds ( ),XPRSgetub ( )

98 getvars procedure getvars(c: linctr, s: set of mpvar) c s vset declarations c: linctr vset: set of mpvar end-declarations getvars(c,vset) forall(x in vset) writeln(getsol(x)) XPRMgetvarnum, COLS ( ), XPRSgetcols ( )

99 96 4 MOSEL initglobal procedure initglobal mmxprs XPRS NIG maxmize minimize maximize,minimize

100 iseof function iseof: boolean declarations d: integer end-declarations... fopen("datafile.dat",f_input) while(iseof = false) do readln(d) writeln(d) end-do fclose(f_input) fclose, fopen, fselect, getfid

101 98 4 MOSEL ishidden function ishidden(c: linctr): boolean c declarations ctr: set of linctr end-declarations... maximize(profit) writeln(getobjval) forall(c in ctr) if(ishidden(c)) then sethidden(c,false) end-if maximize(profit) writeln(getobjval) sethidden Mosel sethidden

102 isodd function isodd(i: integer): boolean i x(1) + x(3) + x(5) + x(7) declarations x: dynamic array(1..8) of mpvar end-declarations forall(i in 1..8 isodd(i)) create(x(i)) c := sum(i in 1..8) x(i)

103 100 4 MOSEL IVEaddtograph IVE procedure IVEaddtograph(xval: real, yval: real) xval yval x y mmive IVEinitgraph("Random graph","xcoord","value") forall (j in ) do IVEaddtograph(j,random) end-do 1. IVEinitgraph x IVEinitgraph

104 IVEinitgraph IVE procedure IVEinitgraph(name: string, x: string, y: string) name x y x y mmive Xpress-MP exp( 0.01x) cos(0.1x) model decay uses "mmive"; IVEinitgraph("Decaying oscillator","x","y") forall(i in ) IVEaddtograph(i,exp(-0.01*i)*cos(i*0.1)) end-model IVEaddtograph

105 102 4 MOSEL ln function ln(r: real): real r 1 i 0 j 1 k i := ln(1) j := ln(m_e) k := ln(0) 2 log exp, log, sqrt

106 loadbasis savebasis (delbasis ) procedure loadbasis(num: integer) procedure loadbasis(name: string) num name mmxprs Boolean VERBOSE / 2 declarations MinCost: linctr end-declarations savebasis(2)... loadprob(mincost) loadbasis(2) loadbasis maximize minimize loadprob delbasis, savebasis, XPRSloadbasis ( )

107 104 4 MOSEL loadcuts procedure loadcuts(cuttype: integer, interpret: integer, cuts: set of integer) procedure loadcuts(cuttype: integer, interpret: integer) cuttype interpret cuttype 3 cuttype cuts mmxprs Boolean cuttype VERBOSE / addcut, addcuts, delcuts, dropcuts, storecut, storecuts, XPRSloadcuts ( )

108 loadprob procedure loadprob(obj: linctr) procedure loadprob(force: boolean, obj: linctr) procedure loadprob(obj: linctr, extravar: set of mpvar) procedure loadprob(force: boolean, obj: linctr, extravar: set of mpvar) obj force extravar mmxprs Boolean LOADNAMES MPS / VERBOSE / 2 2 declarations MinCost: linctr end-declarations... savebasis(2)... loadprob(mincost) loasbasis(2) minimize maximize loadprob force extravar ( ) maximize, minimize, XPRSloadlp ( )

109 106 4 MOSEL log 10 function log(r: real): real r function logn(base,number: real): real if(number > 0 and base > 0) then returned := log(number)/log(base) else exit(1) end-if end-function exp,ln.sqrt

110 mkdir procedure makedir(dirname: string) dirname mmsystem makedir("logfiles") if(getsysstat = 0) then fopen("logfiles/log1.txt",f_output)... fclose(f_output) end-if getsysstat, removedir, system

111 108 4 MOSEL makesos (SOS) procedure makesos1(cs: linctr, s: set of mpvar, c: linctr) procedure makesos1(s: set of mpvar, c: linctr) procedure makesos2(cs: linctr, s: set of mpvar, c: linctr) procedure makesos2(s: set of mpvar, c: linctr) cs s c SOS s c cs

112 maximize procedure maximize(alg: integer, obj: linctr) procedure maximize(obj: linctr) alg XPRS BAR XPRS DUAL XPRS LIN XPRS TOP XPRS PRI XPRS GLB XPRS NIG obj LP LP iniglobal mmxprs Boolean LOADNAMES MPS / VERBOSE / Profit declarations Profit: linctr end-declarations maximize(xprs_dual+xprs_top,profit) 1. ( ) 2. XPRS LIN 3. XPRS TOP LP is LP initglobal, loadprob, minimize, XPRSmaxim (MAXIM) ( )

113 110 4 MOSEL maxlist function maxlist(i1: integer, i2: integer [, i3: integer ]): integer function maxlist(r1: real, r2: real [, r3: real ]): real i1,i2,.. r1,r2,.. maxlist r 7.0 r := maxlist(-1, 4.5, 2, 7, -0.3) minlist

114 minimize procedure minimize(alg: integer, obj: linctr) procedure minimize(obj: linctr) alg XPRS BAR XPRS DUAL XPRS LIN XPRS TOP XPRS PRI XPRS GLB XPRS NIG obj LP LP iniglobal mmxprs Boolean LOADNAMES MPS / VERBOSE / MinCost declarations MinCost: linctr end-declarations minimize(xprs_bar+xprs_lin, MinCost) 1. ( ) 2. XPRS LIN 3. XPRS TOP LP is LP

115 112 4 MOSEL initglobal, loadprob, maximize, XPRSminim (MINIM)( )

116 minlist function minlist(i1: integer, i2: integer [,i3: integer ]): integer function minlist(r1: real, r2: real [,r3: real ]): real i1,i2,.. r1,r2,.. maxlist r -1.0 r := minlist(-1, 4.5, 2, 7, -0.3) maxlist

117 114 4 MOSEL random [0, 1) function random: real i 1 10 setrandseed(round(gettime)) i := integer(round((10*random)+0.5)) setrandseed

118 read,readlin procedure read(e1: expr [,e2: expr ]) procedure readln procedure readln(e1: expr [,e2: expr ]) e1,e2,... Boolean IOCTRL IO / word toto(12 and 45)=word declarations i, j: integer s: string ts: array(range,range) of string end-declarations read(i,j,s) readln("toto(",i,"and",j,")=",ts(i,j)) 1. ei ei NBREAD 2. read Mosel 3. readlin readln readln NBREAD, SQLreadinteger, SQLreadreal, SQLreadstring, write, writeln

119 116 4 MOSEL removedir procedure removedir(dirname: string) dirname mmsystem dirlist declarations dirlist: set of string end-declarations... forall(dir in dirlist) if(bittest(getfstat(dir),sys_typ) = SYS_DIR) then removedir(dir) end-if getfstat, getsysstat, makedir

120 round function round(r: real): integer r c -2 d 2 c := round(-1.9) d := round(1.5) ceil,floor

121 118 4 MOSEL savebasis procedure savebasis(num: integer) procedure savebasis(name: string) num name mmxprs Boolean VERBOSE / 2 declarations MinCost: linctr end-declarations savebasis(2)... loadprob(mincost) loadbasis(2) delbasis loadbasis delbasis, loadbasis, XPRSwritebasis ( )

122 setcallback procedure setcallback(cbtype: integer, cb: string) cbtype cb. XPRS CB IL XPRS CB GL XPRS CB BL XPRS CB USN XPRS CB UPN XPRS CB UON XPRS CB UIN XPRS CB UIS XPRS CB UCN XPRS CB CMI XPRS CB CMS XPRS CB CM ( ) XPRS CB TCM / function cb:boolean.. XPRS CB IL, XPRS CB GL, XPRS CB BL, XPRS CB UON, XPRS CB CM, XPRS CB TCM function cb(node:integer):integer..xprs CB USN, XPRS CB UPN procedure cb..xprs CB UIN, XPRS CB UIS, XPRS CB CMI, XPRS CB CMS procedure cb(node:integer)..xprs CB UCN mmxprs Boolean VERBOSE / procedure printsol declarations objval: real end-declarations

123 120 4 MOSEL objval := getparam("xprs_lpobjval") writeln("solution value: ",objval) end-procedure setcallback(xprs_cb_uis,"printsol") Xpress-Optimizer Mosel / getparam 5.3

124 setcoeff procedure setcoeff(c: linctr, x: mpvar, r: real) procedure setcoeff(c: linctr, r: real) c x r c declarations x,y,z: mpvar end-declarations c := 4*x + y - z <= 12 setcoeff(c,y,2) setcoeff(c,8.1) 4*x + 2*y - z getcoeff, XPRSchgcoef ( ), XPRSchgmcoef ( )

125 122 4 MOSEL sethldden procedure sethidden(c: linctr, b: boolean) c b true false declarations x,y,z: mpvar end-declarations c := 4*x + y - z <= 12 sethidden(c,true) Mosel ishiddden ishidden

126 setlb procedure setlb(x: mpvar, r: real) x r mmxprs Boolean VERBOSE / b 50 profit loadprob(profit) setlb(b,50) maximize(profit) Mosel getlb, getub, loadprob, setub, XPRSchgbounds ( )

127 124 4 MOSEL setmlpdlr procedure setmipdir(x: mpvar, t: integer, r: real) procedure setmipdir(x: mpvar, t: integer) procedure setmipdir(c: linctr, t: integer, r: real) procedure setmipdir(c: linctr, t: integer) x c r t XPRS PR r XPRS UP XPRS DN XPRS PU r XPRS PO r mmxprs Boolean VERBOSE / clearmipdir, XPRSgetdirs ( ),XPRSloaddirs ( )

128 setmodcut procedure setmodcut(c: linctr) c mmxprs Boolean VERBOSE / declarations Items: set of integer ctr: array(items) of linctr end-declarations... forall(i in Items isodd(i)) setmodcut(ctr(i)) maximize(profit) writeln("profit: ",getobjval) clearmodcut clearmodcut, XPRSloadmodelcuts ( )

129 126 4 MOSEL setparam procedure setparam(name: string, val: integer string real boolean) name val 6 6 setparam("xprs_verbose",true) Mosel ( mmxprs.xprs LOADNAMES) getparam, 7

130 setrandseed procedure setrandseed(s: integer) s i 1 10 setrandseed(round(gettime)) i := integer(round((10*random)+0.5)) random random

131 128 4 MOSEL settype procedure settype(c: linctr, type: integer) c type ctr declarations ctr: linctr end-declarations... settype(ctr,ct_eq) maximize(profit) CT EQ = CT GEQ CT LEQ CT UNB CT SOS1 CT SOS2 1 2 CT CONT CT INT CT BIN CT PINT CT SEC CT SINT gettype, XPRSchgrowtype ( )

132 setub procedure setub(x: mpvar, r: real) x r mmxprs Boolean VERBOSE / loadprob(profit) setub(a,100) maximize(profit) Mosel loadprob getlb, getub, loadprob, setlb, XPRSchgbounds ( )

133 130 4 MOSEL sin ( ) function sin(r:real): real r IVE model sine uses "mmive"; IVEinitgraph("Oscillator","x","sin x") forall(i in ) IVEaddtograph(i,sin(i*0.1)) end-model srctan,cos

134 SQLconnect procedure SQLconnect(s: string) s mmodbs Boolean SQLVERBOSE ODBC / MySQL test yve DaSH SQLconnect("DSN=mysql;DB=test;UID=yves;PWD=DaSH") SQLCONNECTION SQLdisconnect, SQLexecute, SQLreadinteger, SQLreadreal, SQLreadstring

135 132 4 MOSEL SQLdisconnect procedure SQLdisconnect mmodbs Boolean SQLVERBOSE ODBC / test yve SQLconnect("DSN=mysql;DB=test;UID=yves;PWD=DaSH") i := SQLreadinteger("select articlenum from pricelist where colour=blue") SQLdisconnect SQLCONNECTION SQLconnect, SQLexecute, SQLreadinteger, SQLreadreal, SQLreadstring

136 SQLexecute SQL procedure SQLexecute(s: string) procedure SQLexecute(s: string, a: array) procedure SQLexecute(s: string, m: set) s a m SQL integer,real,string,boolean integer,real,string,boolean mmodbs Boolean SQLNDXCOL / SQLVERBOSE ODBC / Integer SQLCOLSIZE SQLCONNECTION ODBC 4 SQLexecute colour priicelist colour prices priicelist colour prices 2 ndx prices newtab declarations prices: array( ) of real colours: array( ) of string allcolours: set of string end-declarations... SQLexecute("select colour from pricelist",allcolours) SQLexecute("select articlenum,colour,price from pricelist",[colours,prices]) SQLexecute("create table newtab (ndx integer, price double)") SQLexecute("insert into table newtab (ndx, price) value (?,?)", prices)

137 134 4 MOSEL SQLconnect, SQLdisconnect, SQLreadinteger, SQLreadreal, SQLreadstring, SQLSUCCESS

138 SQLreadinteger function SQLreadinteger(s: string): integer s SQL mmodbs Boolean SQLNDXCOL / SQLVERBOSE ODBC / Integer SQLCOLSIZE SQLCONNECTION ODBC colour blue pricelist i := SQLreadinteger("select articlenum from pricelist where colour=blue") SQL read, readln, SQLexecute, SQLreadreal, SQLreadstring

139 136 4 MOSEL SQLreadreal function SQLreadreal(s: string): real s SQL mmodbs Boolean SQLNDXCOL / SQLVERBOSE ODBC / Integer SQLCOLSIZE SQLCONNECTION ODBC newtab 2 price r := SQLreadreal("select price from newtab where ndx=2") SQL read, readln, SQLexecute, SQLreadinteger, SQLreadstring

140 SQLreadstring function SQLreadstring(s: string): string s SQL mmodbs Boolean SQLNDXCOL / SQLVERBOSE ODBC / Integer SQLCOLSIZE SQLCONNECTION ODBC pricelist 1004 colour s := SQLreadstring("select colour from pricelist where articlenum=1004") SQL read, readln, SQLexecute, SQLreadinteger, SQLreadreal

141 138 4 MOSEL sqrt function sqrt(r: real): real r function isprime(number: integer): boolean if(number < 2) then returned := false end-if i := 1 repeat i+=1 until((number mod i = 0) or i > sqrt(number)) if(i > sqrt(number)) then returned := true else returned := false end-if end-function exp,ln.log

142 storecut function storecut(nodupl: integer, cuttype: integer, type: integer, linexp: linctr): integer nodupl cuttype type linexp CT_GEQ CT_LEQ CT_EQ mmxprs Boolean VERBOSE / loadcuts addcut, addcuts, dropcuts, delcuts, loadcuts, storecuts, XPRSstorecuts ( )

143 140 4 MOSEL storecuts procedure storecuts(nodupl: integer, cuttype: array(range) of integer, type: array(range) of integer, linexp: array(range) of linctr, ndx_a: array(range) of integer) procedure storecuts(nodupl: integer, cuttype: array(range) of integer, type: array(range) of integer, linexp: array(range) of linctr, ndx_s: set of integer) nodupl cuttype type linexp ndx a ndx b CT_GEQ CT_LEQ CT_EQ mmxprs Boolean VERBOSE / loadcuts ndx a ndx b 4 addcut, addcuts, dropcuts, delcuts, loadcuts, storecut, XPRSstorecuts ( )

144 strfmt function strfmt(str: string, len: integer): string function strfmt(i: integer, len: integer): string function strfmt(r: real, len: integer): string function strfmt(r: real, len: integer, dec: integer): string str i r len dec < 0 > text1 text2text3 writeln("text1", strfmt("text2",8), "text3") r := writeln(strfmt(r,0)," ", strfmt(r,4,2), strfmt(r,8,0)) 1. (write writeln ) 2. substr, write, writeln

145 142 4 MOSEL substr function substr(str: string, i1: integer, i2: integer): string str i1 i2 ample te write(substr("example text",3,10))

146 system procedure system(command: string) command mmsystem Newdir system("mkdir Newdir") getsysstat

147 144 4 MOSEL write,writeln procedure write(e1: expr [,e2: expr ]) procedure writeln procedure writeln(e1: expr [,e2: expr ]) e1,e2,.. Boolean IOCTRL IO / Set 1 A real:7.12, a Boolean:true Set1 := {"first", "second", "fifth"} write(set1) writeln b := true writeln("a real:", strfmt(7.1234, 4, 2), ", a Boolean:",b) writeln strfmt linctr,mpvar fselect, read, readln, strfmt

148 5 Xpress-Mosel ( ) 1 Mosel (xprm rt) (xprm mc) ( XPRMcompmod ) xprm rt.h xprm mc.h 5.1 Mosel 145

149 146 5 CLOAD COMPILE DELETE DISPLAY EXAMINE EXPRTPROB FLUSHLIBS INFO LIST LOAD LSLIBS QUIT RESRT RUN SELECT SYMBOLS SYSTEM BIM BIM DSO DSO Mosel OS 5.2 Mosel 1 2

150 XPRMautounloaddso void XPRMautounloaddso(int yesno); yesno 0 XPRMmodel mod; int nreturn;... XPRMautounloaddso(0);... for(mod = XPRMgetnextmod(NULL); mod!= NULL; mod = XPRMgetnextmod(mod)) XPRMrunmod(mod,&nReturn,NULL); XPRMflushdso(); XPRMflushdso XPRMflushdso (FLUSHLIBS)

151 148 5 XRPMchkarrind int XPRMchkarrind(XPRMarray array, int indices[]); array indices array n n- (4, 4) 2 A XPRMmodel model XPRMalltypes atarr; XPRMarray A; int indices[2];... XPRMfindident(model,"A",&atarr); A = atarr.array; indices[0] = 4; indices[1] = 4; if(xprmchkarrind(a,indices)!= 1) printf("index (%d,%d) lies within the range bounds of A\n",indices[0],indices[1]); 0 1 XPRMcmpindices

152 CLOAD (.bim ) Mosel CLOAD [-options] srcfile [userc] options g s BIM m.mod M.mod p e srcfile userc e.mos Mosel myprob.mos %> mosel -c CL myprob; RUN CLOAD COMPILE LOAD XPRMcompmod (COMPILE), XPRMloadmod (LOAD), XPRMrunmod (RUN)

153 150 5 XPRMcmpindices 2 int XPRMcmpindices(int transpose, int nbdim, int ind1[],int ind2[]) transpose nbdim ind1 ind2 0 ind1 ind2 nbdim nbdim A XPRMalltypes atarr; XPRMarray A; int *indices, *trueindices;... findident(model,"a",&atarr); A = atarr.array; indices = malloc(xprmgetarrdim(a)*sizeof(int)); trueindices = malloc(xprmgetarrdim(a)*sizeof(int)); XPRMfirstarrentry(A,indices); XPRMfirstarrtruentry(A,trueindices); if(xprmcmpindices(0, XPRMgetarrdim(A), trueindices, indices)!= 0) printf("dynamic array A has not had its first entry defined\n"); 2-1 ind1 ind ind2 ind1 XPRMchkarrind

154 XPRMcompmod COMPILE (.bim ) XPRMloadmod (LOAD) int XPRMcompmod(const char *options, const char *srcfile, const char *dstfile, const char *userc); COMPILE [-options] srcfile [userc] options g s BIM m.mod M.mod p e srcfile dstfile userc e.mos NULL srcfile.bim 1( ) Mosel myprob.mos XPRMmodel model; int nreturn; XPRMinit(); XPRMcompmod("","myprob.mos",NULL,"My Problem"); model = XPRMloadmod("myprob.bim",NULL); XPRMrunmod(model,&nReturn,NULL); XPRMfree(); 2( ) myprob.mos myprob.bim

155 152 5 COMPILE myprob LOAD myprob RUN QUIT 1. XPRMcompmod Mosel xprm mc.h ( ) 3 3. srcfile NULL 4. -e Mosel.mos.bim ( ) srcfile dstfile CLOAD, XPRMloadmod (LOAD), XPRMrunmod (RUN)

156 DELETE DELETE [number name] number name model Mosel simple.mos alterd.mos alterd simple CLOAD simple RUN CLOAD altered RUN DELETE 1. BIM LIST CLOAD, LIST, XPRMloadmod (LOAD), SELECT, XPRMunloadmod

157 154 5 DISPLAY DISPLAY symbol symbol Mosel myprob.mos (profit) CLOAD myprob RUN DISPLAY profit ( 0) ( 0) getobjval, getsol, SYMBOLS, XPRMgetcsol, XPRMgetobjval, XPRMgetvsol

158 EXAMINE / / EXAMINE [-options] libname options c s / p / t libname mmodbc / / EXAMINE mmodbc XPRMgetdsoinfo, XPRMgetdsoparam, XPRMgetnextdsoconst,XPRMgetnextdsoparam, XPRMgetnextdsoproc

159 156 5 XPRMexportprob EXPORTPROB MPS LP int XPRMexportprob(XPRMmodel model, const char *options,const char *fname, XPRMlinctr obj); EXPORTPROB [-options] [filename [obj]] model options m MPS LP ( ) p (LP ) s filename obj NULL NULL 1( ) LP XPRMmodel model; int nreturn;... XPRMrunmod(model,&nReturn,NULL); XPRMexportprob(model,"p",NULL,NULL) 2( ) MPS myprob.mat RUN EXPORTPROB -m myprob.mat NULL MPS.mat LP.lp ( sp ) Mosel exportprob, XPRBexportprob BCL

160 XPRMfinddso DSO XPRMdsolib XPRMfinddso(const char *libname); libname mmxprs XPRMdsolib lib; int libv; const char *name;... if((lib = XPRMfinddso("mmxprs"))!= NULL) { XPRMgetdsoinfo(lib,&name,NULL,&libv,NULL); printf("got module: %s version %d.%d.%d\n", name, libv/ , (libv/1000)-1000*(libv/ ), libv *(libv/1000)); } NULL XPRMgetnextdso

161 158 5 XPRMfindident int XPRMfindident(XPRMmodel model, const char *text,xprmalltypes *value); model text value a b XPRMalltypes atvar; XPRMmpvar a, b; in nreturn;... XPRMrunmod(simple,&nReturn,NULL); XPRMfindident(simple,"a",&atvar); a = atvar.mpvar; XPRMfindident(simple,"b",&atvar); b = atvar.mpvar; printf("a = %g, b = %g\n",xprmgetsol(a),xprmgetsol(b)); 1. XPRMfindident 0 XPRM TYP(t) XPRM STR(t) XPRM STR CONST XPRM STR REF XPRM STR SET XPRM STR ARR XPRM STR PROC

162 XPRM TYP NOT XPRM TYP INT XPRM TYP REAL XPRM TYP STRING XPRM TYP BOOL XPRM TYP MPVAR XPRM TYP LINCTR XPRM TYP SET XPRM TYP ARR XPRMgettypename 2. XPRMalltypes XPRMfindident value.integer value.real value.string value.boolean value.mpvar value.linctr value.set value.array value.proc XPRMgetnextident ( ) ( )

163 160 5 XPRMfindmod XPRMmodel XPRMfindmod(const char *name, int number); name number NULL 1 XPRMmodel myprob; XPRMinit(); XPRMloadprob("myprob.bim","My Problem");... myprob = XPRMfindmod("My Problem",-1); 1. XPRMfindmod NULL 2. (.bim ) (number = 1) (number =NULL) name XPRMloadmod

164 XPRMfirstarrentry int XPRMfirstarrentry(XPRMarray array, int indices[]); array indices n- n array A XPRMmodel model; XPRMalltypes atarr; XPRMarray A; int i,*indices;... XPRMfindident(model,"A",&atarr); A = atarr.array; indices = malloc(xprmgetarrdim(a)*sizeof(int)); XPRMfirstarrentry(A,indices); for(i=0;i<xprmgetarrdim(a);i++) printf("index of first entry of A in dimension %d is %d\n",i+1,indices[i]); XPRMfirstarrtruentry, XPRMfirstsetndx, XPRMlastarrentry, XPRMnextarrentry

165 162 5 XPRMfirstarrtruentry int XPRMfirstarrtruentry(XPRMarray array, int indices[]); array indices n- n array A XPRMmodel model; XPRMalltypes atarr; XPRMarray A; int i,*indices;... XPRMfindident(model,"A",&atarr); A = atarr.array; indices = malloc(xprmgetarrdim(a)*sizeof(int)); XPRMfirstarrtruentry(A,indices); for(i=0;i<xprmgetarrdim(a);i++) printf("index of first true entry of A in dimension %d is %d\n",i+1,indices[i]); XPRMfirstarrentry XPRMfirstarrentry, XPRMlastarrentry, XPRMnextarrtruentry

166 XPRMfirstsetndx int XPRMfirstsetndx(XPRMset set); set Items XPRMmodel model; XPRMalltypes atset, atelt; XPRMset items;... XPRMfindident(model,"Items",&atset); items = atset.set; if(xprm_typ(xprmgetsettype(items)) == XPRM_TYP_INT) for(i=firstsetndx(items);i<=xprmlastsetndx(items);i++) { XPRMgetelsetval(items,i,&atelt); printf("items[%d] = %d\n",j,atelt.integer); } 1. string 1 XPRMgetelsetndx, XPRMgetelsetval, XPRMgetsetsize,XPRMlastsetndx

167 164 5 XPRMflushdso FLUSHLIBS void XPRMflushdso(void); FLUSHLIBS 1( ) XPRMdsolib lib; const char *name; int libv;... XPRMflushdso(); for(lib=null;(lib = XPRMgetnextdso(lib))!= NULL;) { XPRMgetdsoinfo(lib,&name,NULL,&libv,NULL); printf("model uses %s, version %d.%d.%d\n", name, libv/ , (libv/1000)-1000*(libv/ ), libv-1000*(libv/1000)); } 2( ) RUN FLUSHLIBS LSLIBS fflush, XPRMautounloaddso