LatticeStuff.cpp

#include "Headers.h"
#include "Parameters.h" // >> contains all the req parameters

void create_Bmat_bimodal(std::vector <int>& Bmat)
{
	long i = 0, var = 0, count1 = 0, count2 = 0;
	std::uniform_int_distribution<int>  distr(0, 1000);
	std::random_device  rand_dev;
	std::mt19937        generator(rand_dev());
	std::cout << "B mtrix ->" << tab;
	for (i = 0; i < N; i++)
	{
		//Bmat[i] = distr(generator);// floor(rand() / 1000);
		var = distr(generator);
		if (var < 500)
		{
			Bmat[i] = w;
			count1++;
		}
		else
		{
			Bmat[i] = -w;
			count2++;
		}
		//std::cout << Bmat[i] << d;
	}	std::cout << std::endl << "counts : " << count1 << d << count2 << std::endl;


}

// http://www.cplusplus.com/reference/random/normal_distribution/
void create_Bmat_gaussian(std::vector <int>& Bmat)
{
	double var = 0;
	long i = 0,  count1 = 0, count2 = 0;
	std::normal_distribution<double>  distr(0, 1);
	std::random_device  rand_dev;
	std::mt19937        generator(rand_dev());
	std::cout << "B mtrix ->" << tab;
	for (i = 0; i < N; i++)
	{
		//Bmat[i] = distr(generator);// floor(rand() / 1000);
		var = (distr(generator));
		Bmat[i] = var;
		if (var < 0)
		{
			//Bmat[i] = w;
			count1++;
		}
		else
		{
			//Bmat[i] = -w;
			count2++;
		}
		//std::cout << Bmat[i] << d;
	}	std::cout << std::endl << "counts : " << count1 << d << count2 << std::endl;


}

void create_Exmat(std::vector < std::vector <int> >& Exmat, std::vector <int> & latt)
{
	int i = 0;
	for (i = 0; i < N; i++)
	{
		if (((i + 1) % VER != 0))
		{
			Exmat[i][i + 1] = 1 * latt[i] * latt[i + 1];
			//cout << Exmat[i][i + 1];
		}
		if (i + VER < N)
		{
			Exmat[i][i + VER] = 1 * latt[i] * latt[i + VER];
		}
		//else 
	}
}


void create_CapacityMat(std::vector < std::vector <int> >& CapacityMat, std::vector < std::vector <int> >& Exmat)
{
	for (int i = 1; i <= N; i++)
	{
		//cout << "\n";
		for (int j = 1; j <= N; j++)
		{
			if (i < j)
			{
				CapacityMat[i][j] = 4 * J * Exmat[i - 1][j - 1] * 10;

			}
			else
				CapacityMat[i][j] = 0;
			//cout << CapacityMat[i][j] << " ";
		}
	}
}


void create_Wmat(std::vector < int >& Wmat,std::vector < std::vector <int> >& CapacityMat, std::vector < int >& Bmat, int del )
{
	int i = 0, j = 0, cap = 0;
	for (i = 1; i <= N; i++)
	{
		for (j = 1; j <= N; j++)
		{
			cap += CapacityMat[i][j] - CapacityMat[j][i];
			//cout << cap<<"\t"; 
		}
		Wmat[i - 1] = -2 * Bmat[i - 1] * del - cap / 2;
		cap = 0;
		//cout <<Wmat[i-1]<< endl;
	}
}

// making reisdual graph by overwriting the flow mat,flow[i][j] = CapacityMat[i][j] - flow[i][j];
void create_Augumented_CapacityMat(std::vector < int >& Wmat, std::vector < std::vector <int> >& CapacityMat)
{
	int i = 0;
	for (i = 1; i <= N; i++)
	{
		if (Wmat[i - 1] > 0)
		{
			CapacityMat[0][i] = 0;
			CapacityMat[i][N + 1] = Wmat[i - 1];
		}

		else
		{
			CapacityMat[0][i] = -Wmat[i - 1];
			CapacityMat[i][N + 1] = 0;
		}

	}
}

void create_Residual_graph(std::vector < std::vector <int> >& CapacityMat, std::vector < std::vector <int> >& flow)
{
	int i = 0, j = 0;
	for (i = 0; i < V; i++)
	{
		for (int j = 0; j < V; j++)
		{
			flow[i][j] = CapacityMat[i][j] - flow[i][j];
			//cout<< CapacityMat[i][j]-flow[i][j];
		}
		//cout<<endl;
	}
}

void printMatrix(std::vector <std::vector <int> >  & M, int len) {
	int i, j;
	for (i = 0; i < len; i++) {
		for (j = 0; j < len; j++)
			std::cout << "\t" << M[i][j];
		std::cout << "\n";
	}
}

void printMatrix(std::vector<int>  & M, int len) {
	int i;
	for (i = 0; i < len; i++) {
		std::cout << "\t" << M[i];

	}std::cout << "\n";
}