/******************************************************** Xpress-BCL Java Example Problems ================================ file xbcoco2.java ````````````````` Coco Problem Phase 2. Use parameters, data tables and subscripted variables to separate the model structure from the data. Read data tables in from text data files. (c) 2008-2024 Fair Isaac Corporation author: S.Heipcke, Jan. 2000, rev. Mar. 2011 ********************************************************/ import java.io.*; import com.dashoptimization.*; public class xbcoco2 { static final int NP = 2; /* Number of products (p) */ static final int NF = 2; /* factories (f) */ static final int NR = 2; /* raw materials (r) */ static final String REVFILE = System.getProperty("XPRBDATA") + "/coco/rev.dat"; static final String CMAKEFILE = System.getProperty("XPRBDATA") + "/coco/cmake.dat"; static final String CBUYFILE = System.getProperty("XPRBDATA") + "/coco/cbuy.dat"; static final String REQFILE = System.getProperty("XPRBDATA") + "/coco/req.dat"; static final String MXSELLFILE= System.getProperty("XPRBDATA") + "/coco/maxsell.dat"; static final String MXMAKEFILE= System.getProperty("XPRBDATA") + "/coco/mxmake.dat"; /****TABLES****/ static double[] REV; /* Unit selling price of product p */ static double[][] CMAK; /* Unit cost to make product p at factory f */ static double[] CBUY; /* Unit cost to buy raw material r */ static double[][] REQ; /* Requirement by unit of prod. p for raw mat. r */ static double[] MXSELL; /* Max. amount of p that can be sold */ static double[] MXMAKE; /* Max. amount factory f can make over all prod.s */ static double[][] PROFIT; /* Profit contribution of product p at factory f */ /***********************************************************************/ static void modCoco2() { XPRBvar[][] make; XPRBexpr lobj, lc; int p,f; try (XPRBprob pb = new XPRBprob("Coco2")) { /* Initialize BCL and create a new problem */ /****VARIABLES****/ make = new XPRBvar[NP][NF]; for(p=0;p<NP;p++)         /* Amount of prod. p to make at factory f */
                for(f=0;f<NF;f++)
                    make[p][f] = pb.newVar("make_"+ (p+1) +"f"+ (f+1));

            /****OBJECTIVE****/
            lobj = new XPRBexpr();
            for(p=0;p<NP;p++)         /* Objective: maximize total profit */
                for(f=0;f<NF;f++)  lobj.add(make[p][f].mul(PROFIT[p][f]));
            pb.setObj(lobj);

            /****CONSTRAINTS****/

            for(p=0;p<NP;p++) {          /* Limit on the amount of product p to be sold */
                lc = new XPRBexpr();
                for(f=0;f<NF;f++)  lc .add(make[p][f]);
                pb.newCtr("MxSell", lc.lEql(MXSELL[p]) );
            }

            for(f=0;f<NF;f++) {          /* Capacity limit at factory f */
                lc = new XPRBexpr();
                for(p=0;p<NP;p++) lc.add(make[p][f]);
                pb.newCtr("MxMake", lc.lEql(MXMAKE[f]) );
            }

            /****SOLVING + OUTPUT****/
            pb.setSense(XPRB.MAXIM);  /* Choose the sense of the optimization */
            pb.lpOptimize("");        /* Solve the LP-problem */
            System.out.println("Objective: " + pb.getObjVal()); /* Get objective value */

            for(p=0;p<NP;p++)         /* Print the solution values */
                for(f=0;f<NF;f++)
                    System.out.print(make[p][f].getName()+ ":"+ make[p][f].getSol() +" ");
            System.out.println();
        }
    }

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

    /**** Initialize the stream tokenizer ****/
    static StreamTokenizer initST(FileReader file) {
        StreamTokenizer st=null;

        st= new StreamTokenizer(file);   /* Initialize the stream tokenizer */
        st.commentChar('!');             /* Use the character '!' for comments */
        st.eolIsSignificant(true);       /* Return end-of-line character */
        st.ordinaryChar(',');            /* Use ',' as separator */
        st.parseNumbers();               /* Read numbers as numbers (not strings)*/
        return st;
    }

    /**** Read single line data file ****/
    static void readOneLineFile(String filename, double[] data) throws IOException {
        FileReader datafile=null;
        StreamTokenizer st;
        int i=0;

        datafile = new FileReader(filename);
        st = initST(datafile);
        do {
            do {
                st.nextToken();
            } while(st.ttype==st.TT_EOL);    /* Skip empty lines */
            while(st.ttype == st.TT_NUMBER) {
                data[i++] = st.nval;
                if(st.nextToken() != ',') break;
                st.nextToken();
            }
        } while( st.ttype == st.TT_EOL );
        datafile.close();
    }

    /**** Read multi line data file ****/
    static void readMultiLineFile(String filename, double[][] data) throws IOException {
        FileReader datafile=null;
        StreamTokenizer st;
        int i=0, j=0;

        datafile = new FileReader(filename);
        st = initST(datafile);
        do {
            do {
                st.nextToken();
            } while(st.ttype==st.TT_EOL);    /* Skip empty lines */
            while(st.ttype == st.TT_NUMBER) {
                data[j][i++] = st.nval;
                if(st.nextToken() != ',') { j++; i=0; break;}
                st.nextToken();
            }
        } while( st.ttype == st.TT_EOL );
        datafile.close();
    }

    /**** Read data from files ****/
    static void readData() throws IOException {
        int p,f,r;

        REV = new double[NP];
        CMAK = new double[NP][NF];
        CBUY = new double[NR];
        REQ = new double[NP][NR];
        MXSELL = new double[NP];
        MXMAKE = new double[NF];
        PROFIT = new double[NP][NF];

        /* Read the revenu data file */
        readOneLineFile(REVFILE,REV);

        /* Read the production cost data file */
        readMultiLineFile(CMAKEFILE,CMAK);

        /* Read the raw material cost data file */
        readOneLineFile(CBUYFILE,CBUY);

        /* Read the resource requirement data file */
        readMultiLineFile(REQFILE,REQ);

        /* Read the max. sales quantities data file */
        readOneLineFile(MXSELLFILE,MXSELL);

        /* Read the production capacities data file */
        readOneLineFile(MXMAKEFILE,MXMAKE);

        /* Calculate the table PROFIT */
        for(p=0;p<NP;p++)
            for(f=0;f<NF;f++) {
                PROFIT[p][f] = REV[p] - CMAK[p][f];
                for(r=0;r<NR;r++) PROFIT[p][f] -= (REQ[p][r]*CBUY[r]);
            }
    }

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

    public static void main(String[] args) {
        try {
            readData();           /* Data input from file */
        }
        catch(IOException e) {
            System.err.println(e.getMessage());
            System.exit(1);
        }
        modCoco2();            /* Model and solve the problem */
    }
}