init commit

src/Model.h

@@ -0,0 +1,26 @@
+#pragma once
+
+#include <memory>
+
+struct Facility
+{
+    int m_Capacity;
+    int m_Cost;
+    float m_X;
+    float m_Y;
+};
+
+struct Customer
+{
+    int m_Demand;
+    float m_X;
+    float m_Y;
+};
+
+
+struct Cost
+{
+    Customer m_Customer;
+    Facility m_Facility;
+    int m_Cost;
+};

src/Parser.cpp

@@ -0,0 +1,65 @@
+#include "Parser.h"
+
+#include <fstream>
+#include <iostream>
+#include <sstream>
+
+FileData ParseFile(const std::string& a_FileName) {
+	std::ifstream in(a_FileName);
+	if (!in) {
+		std::cerr << "Impossible d'ouvrir " << a_FileName << "!\n";
+		return {};
+	}
+
+	FileData fileData;
+
+	std::string line;
+	while (getline(in, line))
+	{
+		if (line == "FACILITIES [capacity cost x y]") {
+			std::cout << "Reading facilities ...\n";
+			break;
+		}
+	}
+
+	// facilities
+	while (getline(in, line))
+	{
+		if (line == "CUSTOMERS [demand x y]") {
+			std::cout << "Reading customers ...\n";
+			break;
+		}
+		Facility facility;
+		std::stringstream stream(line);
+		stream >> facility.m_Capacity >> facility.m_Cost >> facility.m_X >> facility.m_Y;
+		fileData.m_Facilities.push_back(facility);
+	}
+
+	// customers
+	while (getline(in, line))
+	{
+		if (line == "COSTS CUSTOMERS/FACILITIES [id_cus id_fac cost]") {
+			std::cout << "Reading costs ...\n";
+			break;
+		}
+
+		Customer customer;
+		std::stringstream stream(line);
+		stream >> customer.m_Demand >> customer.m_X >> customer.m_Y;
+		fileData.m_Customers.push_back(customer);
+	}
+
+	// costs
+	while (getline(in, line))
+	{
+		Cost cost;
+		int facilityID, customerID;
+		std::stringstream stream(line);
+		stream >> customerID >> facilityID >> cost.m_Cost;
+		cost.m_Customer = fileData.m_Customers[customerID - 1];
+		cost.m_Facility = fileData.m_Facilities[facilityID - 1];
+		fileData.m_Costs.push_back(cost);
+	}
+
+	return fileData;
+}

src/Parser.h

@@ -0,0 +1,13 @@
+#pragma once
+
+#include "Model.h"
+#include <vector>
+
+struct FileData
+{
+	std::vector<Facility> m_Facilities;
+	std::vector<Customer> m_Customers;
+	std::vector<Cost> m_Costs;
+};
+
+FileData ParseFile(const std::string& a_FileName);

src/main.cpp

@@ -0,0 +1,56 @@
+#include "Parser.h"
+#include <fmt/format.h>
+#include <iostream>
+#include <kiwi/kiwi.h>
+
+using namespace kiwi;
+
+int main(int argc, char** argv) {
+	auto fileData = ParseFile("RhoneCities_100.flpb");
+
+	std::vector<Variable> y;
+	std::vector<std::vector<Variable>> x(fileData.m_Facilities.size());
+
+	std::cout << "Adding variables ...\n";
+	Solver solver;
+
+	for (std::size_t j = 0; j < fileData.m_Facilities.size(); j++) {
+		y.emplace_back(fmt::format("y_{}", j));
+		for (std::size_t i = 0; i < fileData.m_Customers.size(); i++) {
+			x[i].emplace_back(fmt::format("x_{}{}", i, j));
+		}
+	}
+
+	std::cout << "Adding constraints 1...\n";
+	for (std::size_t j = 0; j < fileData.m_Facilities.size(); j++) {
+		Expression sum;
+
+		for (std::size_t i = 0; i < fileData.m_Customers.size(); i++) {
+			sum = sum + fileData.m_Customers[i].m_Demand * x[i][j];
+		}
+
+		solver.addConstraint({sum <= fileData.m_Facilities[j].m_Capacity * y[j]});
+	}
+
+	std::cout << "Adding constraints 2...\n";
+	for (std::size_t i = 0; i < fileData.m_Customers.size(); i++) {
+		Expression sum;
+
+		for (std::size_t j = 0; j < fileData.m_Facilities.size(); j++) {
+			sum = sum + x[i][j];
+		}
+
+		solver.addConstraint({sum == 1});
+	}
+
+	std::cout << "Updating variables...\n";
+	solver.updateVariables();
+
+	for (std::size_t j = 0; j < fileData.m_Facilities.size(); j++)
+	{
+		std::cout << y[j].name() << "=" << y[j].value() << "\n";
+	}
+	
+
+	return 0;
+}