/******************************************************* Mosel Example Problems ====================== file runrtparqueued.c ````````````````````` Running several instances of a model from another Mosel model. - Queuing submodels for parallel execution in a distributed architecture (one or several models per node) - Before running this model, you need to set up the list NodeList with machine names/addresses of your local network. All nodes that are used need to have the same version of Xpress installed and suitably licensed, and the server "xprmsrv" must have been started on these machines. The maximum number of models per node in array MaxMod needs to be adapted to the number of executions licensed on the corresponding nodes. All files are local to the root node, no write access is required at remote nodes. (c) 2012 Fair Isaac Corporation author: S. Heipcke & Y. Colombani, Apr. 2012 *******************************************************/ #include <stdio.h>
#include <stdlib.h>
#include <string.h>
#include "xprd.h"

#define J 10         /* Number of jobs to run */
#define NUMPAR 2     /* Number of parallel model executions */
                     /* (preferrably <= no. of processors) */

#define MAXQSIZE 20  /* Maximum number of elements in a queue (>=J) */

                     /* Use the name or IP address of a machine in
                        your local network, or "" for current node */
const char *NodeList[]={"localhost","localhost"};
#define NodeListSize (sizeof(NodeList)/sizeof(NodeList[0]))
#define nbNodes (NodeListSize<J?NodeListSize:J) int const char typedef struct void XPRDcontext XPRDmosel XPRDmodel s_queue double lastId=0; xprd=XPRDinit(); Create an XPRD context Setting up remote Mosel instances for(n=0;n<nbNodes;n++) mosInst[n]=XPRDconnect(xprd,NodeList[n], if(mosInst[n]==NULL) to connect MaxMod[n]=NUMPAR; Adapt this setting number of processors and licences per node Compile the model file on first XPRDcompmod(mosInst[0],="" ,="rmt:rtparams.mos" Loading nct=0; for(m=0;(m<MaxMod[n]) Load bim modPar[nct]=XPRDloadmod(mosInst[n],>lastId) lastId=XPRDgetnumber(modPar[nct]);
   nct++;
  }

 jobid=malloc(sizeof(int)*(lastId+1));
 modid=malloc(sizeof(int)*(lastId+1));
 modNode=malloc(sizeof(char *)*(lastId+1));

 for(j=0;j<nct;j++) i=XPRDgetnumber(modPar[j]); modid[i]=j; Store the model ID modNode[i]=strdup(XPRDsysinfo(XPRDgetmosel(modPar[j]),XPRD_SYS_NODE, for(i=0;i<J;i++) Define list of jobs JobSize=queue_size(&JobList); number Start initial lot runs for(j=0;j<nct;j++) Run all remaining Wait for termination if(eventClass==XPRD_EVENT_END) We are only interested in="end" events Keep track job>0)     /* Start a new run if queue not empty */
    start_next_job(&JobList,sender);
  }
 }

 for(n=0;n<nbNodes;n++) Disconnect remote instances Terminate XPRD Cleaning up return Start next job in a queue void char int i=queue_get(jobList); Retrieve and remove first entry from on sprintf(params,="PARAM1=%d,PARAM2=%g,PARAM3='a string',PARAM4=%s" (i%2==0)?"true":"false"); jobid[XPRDgetnumber(model)]=i; Initialize Add an element to>size>=MAXQSIZE) return -1;
 else
 {
  q->buf[q->head++]=e;
  if(q->head>=MAXQSIZE) q->head=0;
  return ++q->size;
 }
}

/**********************************/
/* Remove an element from a queue */
/**********************************/
int queue_get(s_queue *q)
{
 int rts;

 if(q->size<1) return -1;
 else
 {
  rts=q->buf[q->tail++];
  if(q->tail>=MAXQSIZE) q->tail=0;
  q->size--;
  return rts;
 }
}

/***************************/
/* Get the size of a queue */
/***************************/
int queue_size(s_queue *q)
{ return q->size; }