FICO
FICO Xpress Optimization Examples Repository
FICO Optimization Community FICO Xpress Optimization Home
Back to examples browserPrevious exampleNext example

Setops - Index sets

Description
A example showing the use of index sets and in particular, creating the union and intersection of index sets.


Source Files
By clicking on a file name, a preview is opened at the bottom of this page.
xbsetops.c[download]





xbsetops.c

/********************************************************
  BCL Example Problems
  ====================

  file xbsetops.c
  ```````````````
  Set operations.

  (c) 2008-2024 Fair Isaac Corporation
      author: S.Heipcke, Jan. 2000
********************************************************/

#include <stdio.h>
#include "xprb.h"

/* Two C functions are defined for union and intersection of (index) sets.
   Firstly, the new index set is created, then the elements are added,
   here by using function XPRBaddidxel. BCL also offers the possibility
   to attribute names to the new index sets.                             */

/* Define arrays with index names */
char *city_names[]={"rome", "bristol", "london", "paris", "liverpool"};
char *port_names[]={"plymouth", "bristol", "glasgow", "london", "calais", "liverpool"};

/***********************************************************************/

    /**** Create the union of two index sets ****/
void createunion(XPRBprob prob, XPRBidxset a, XPRBidxset b, XPRBidxset *c,char *name)
{
 int i;

 *c=XPRBnewidxset(prob,name,XPRBgetidxsetsize(a)+XPRBgetidxsetsize(b));
 for(i=0;i<XPRBgetidxsetsize(a);i++)
  XPRBaddidxel(*c,XPRBgetidxelname(a,i));
 for(i=0;i<XPRBgetidxsetsize(b);i++)
  XPRBaddidxel(*c,XPRBgetidxelname(b,i));
}

    /**** Create the intersection of two index sets ****/
void createinter(XPRBprob prob, XPRBidxset a, XPRBidxset b, XPRBidxset *c,char *name)
{
 int i;

 *c=XPRBnewidxset(prob,name,
   XPRBgetidxsetsize(a)<XPRBgetidxsetsize(b)?XPRBgetidxsetsize(a):XPRBgetidxsetsize(b));
 for(i=0;i<XPRBgetidxsetsize(a);i++)
  if(XPRBgetidxel(b,XPRBgetidxelname(a,i))>=0)
   XPRBaddidxel(*c,XPRBgetidxelname(a,i));
}

/***********************************************************************/

int main(int argc, char **argv)
{
 XPRBidxset ports, cities, both, places;
 int i;
 XPRBprob prob;

 prob=XPRBnewprob("Setops");     /* Initialize BCL */

        /* Create sets "cities" and "ports" and add the indices */
 cities=XPRBnewidxset(prob,"cities",5);
 for(i=0;i<5;i++) XPRBaddidxel(cities,city_names[i]);
 ports=XPRBnewidxset(prob,"ports",6);
 for(i=0;i<6;i++) XPRBaddidxel(ports,port_names[i]);

        /* Create the union of "cities" and "ports" and print it */
 createunion(prob,cities,ports,&places,"places");
 XPRBprintidxset(places);

        /* Create the intersection of "cities" and "ports", print it */
 createinter(prob,cities,ports,&both,"bothsets");
 XPRBprintidxset(both);

 return 0;
}

Back to examples browserPrevious exampleNext example