/********************************************************
  Xpress-BCL C# Example Problems
  ==============================

  file xbexpl2.cs
  ```````````````
  Transportation model from the Release 10,5 Supplement.

  (c) 2008 Fair Isaac Corporation
      authors: S.Heipcke, D.Brett.
********************************************************/

using System;
using System.Text;
using System.IO;
using BCL;


namespace Examples
{

    public class TestUGExpl2
    {

        const int MaxSuppliers = 100;         /* Max. number of suppliers */
        const int MaxCustomers = 1000;        /* Max. number of customers */
        const int MaxArcs = 10000;            /* Max. number of non-zero cost values */
        
        //Define XPRBDATAPATH to whichever folder you wish.
        const string XPRBDATAPATH = "../../data";

        const string DEMANDFILE = XPRBDATAPATH + "/trans/ex2dem1.dat";
        /* Demand data file (comma-separated format) */
        const string  AVAILFILE = XPRBDATAPATH + "/trans/ex2avail.dat";
        /* Supply data file (comma-separated format) */
        const string  COSTFILE = XPRBDATAPATH + "/trans/ex2cost.dat";  
        /* Cost data file (comma-separated format) */

        XPRBindexSet Suppliers;          /* Set of suppliers */
        XPRBindexSet Customers;          /* Set of customers */
        double[] AVAIL = new double[MaxSuppliers];      /* Availability of products */
        double[] DEMAND = new double[MaxCustomers];     /* Demand by customers */

        public struct CostStruct{
            public int suppl;
            public int custm;
            public double value;
        } ;                 /* Cost per supplier-customer pair */

        CostStruct[] COST = new CostStruct[MaxArcs];
        
        int NSuppl=0,NCustom=0, NArc=0;  /* Actual numbers of suppliers, customers, 
        and arcs */

        XPRBprob p = new XPRBprob("Trans");             /* Initialize a new problem in BCL */

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

        void modTrans()
        {
            XPRBexpr lobj;
            XPRBexpr[] av;
            XPRBexpr[] de;
            int s,c,a;
            XPRBvar[] x;

            /****VARIABLES****/
            x = new XPRBvar[NArc];
            if(x==null) System.Console.WriteLine("Allocating memory for variables failed.");
            for(a=0; a<NArc; x[a]=p.newVar("x"); lobj=new for(a=0; +=COST[a].value*x[a]; Define set objective function Create all constraints in a single loop Initialize the linear expressions av=new if(av==null) enough memory for AV de=new if(de==null) DE XPRBexprs: i=0; av[i]=new de[i]=new Add terms to one-by-one Terminate constraint definition for(s=0; <=AVAIL[s]); for(c=0;>= DEMAND[c]);

            /****SOLVING + OUTPUT****/
            p.exportProb(BCLconstant.XPRB_MPS,"trans");   /* Matrix generation & output to MPS file */

            p.setSense(BCLconstant.XPRB_MINIM);
            p.lpOptimize();
            System.Console.WriteLine("Objective: " + p.getObjVal());   /* Get objective value */

            for(a=0; a<NArc; Print out the solution values>0)
                {
                    System.Console.Write(Suppliers.getIndexName(COST[a].suppl) + " (" + AVAIL[COST[a].suppl] + ") -> ");
                    System.Console.Write(Customers.getIndexName(COST[a].custm) + " (" + DEMAND[COST[a].custm] + "): ");
                    System.Console.WriteLine(x[a].getSol());  
                }
        }

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

        /**** Read data from files ****/
        public void readData()
        {
            double value;
            FileStream file;
            StreamReader fileStreamIn;
            string name, name2;

            /* Create supplier and customer index sets */
            Suppliers=p.newIndexSet("suppl",MaxSuppliers);
            Customers=p.newIndexSet("custom",MaxCustomers);

            /* Read the demand data file */
            file = new FileStream(DEMANDFILE, FileMode.Open, FileAccess.Read);
            fileStreamIn = new StreamReader(file);
            object[] tempobj = new object[3];
            while (p.XPRBreadarrline(fileStreamIn, 200, "{T} , {g} ", out tempobj, 1) == 2)
            {
                name = (string)tempobj[0];
                value = (double)tempobj[1];
                DEMAND[(Customers + name)] = value;
            }
            fileStreamIn.Close();
            file.Close();
            NCustom = Customers.getSize();


            /* Read the supply data file */
            file = new FileStream(AVAILFILE, FileMode.Open, FileAccess.Read);
            fileStreamIn = new StreamReader(file);
            while (p.XPRBreadarrline(fileStreamIn, 200, "{T} , {g} ", out tempobj, 1) == 2)
            {
                name = (string)tempobj[0];
                value = (double)tempobj[1];
                AVAIL[Suppliers + name] = value;
            }
            fileStreamIn.Close();
            file.Close();
            NSuppl = Suppliers.getSize();

            /* Read the cost data file */
            NArc = 0;
            file = new FileStream(COSTFILE, FileMode.Open, FileAccess.Read);
            fileStreamIn = new StreamReader(file);
            while (p.XPRBreadarrline(fileStreamIn, 200, "{T} , {T} , {g} ", out tempobj, 1) == 3)
            {
                name = (string)tempobj[0];
                name2 = (string)tempobj[1];
                value = (double)tempobj[2];

                COST[NArc].suppl = Suppliers.getIndex(name);
                COST[NArc].custm = Customers.getIndex(name2);

                if (COST[NArc].custm < 0) 
                    System.Console.WriteLine("Cust(" + name2 + ")");
                if (COST[NArc].suppl < 0)
                    System.Console.WriteLine("Supp(" + name + ")");

                COST[NArc++].value = value;
            }
            fileStreamIn.Close();
            file.Close();

            System.Console.WriteLine("C: " + NCustom + "  S: " + NSuppl + "  A: " + NArc);
        }

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

        public static void Main()
        {
           XPRB.init();
           TestUGExpl2 TestInstance = new TestUGExpl2();

            TestInstance.readData();            /* Data input from file */
            TestInstance.modTrans();            /* Formulate and solve the problem */

            return;
        } 
    }
}