/*
 * PSEUDO-RANDOM NUMBER GENERATOR
 *
 * Peter Signell, 3/18/2002
 */
#include"lib351.h"

void main(){

    // output initialization
	ofstream dataOutStr;
    outFileMaker("m354p1c.out", &dataOutStr);
    
    // program variables definition
 	int intN, intNumInBlock, intIndex, intBeginRow, intNumLeftOver;
	float fltRanNum, fltM, fltX, fltC;
	float fltNormdRanNum[1001];

	// cycle through the runs
  	bool goout = false;
  	while (!goout){

    	// data input
    	cout<<"Please, supply values :"<<endl;
    	cout<<"\nnumber of pseudorandom numbers desired, N ="<<endl;
    	cin>>intN;
    	cout<<"\nPower Residue modulus, M ="<<endl;
    	cin>>fltM;
    	cout<<"\nfirst number of the series (the 'seed'), X ="<<endl;
    	cin>>fltX;
    	cout<<"Power Residue constant, C ="<<endl;
    	cin>>fltC;

		// print the header
		float fltMaxNum = fltM - 1.0;
		dataOutStr<<"\f"<<endl;
		dataOutStr<<"Pseudorandom Number Generator (Power Residue Method)"<<endl<<endl;
        dataOutStr<<"number of numbers, N = "<<intN<<endl;
		dataOutStr<<"Power Residue modulus, M = "<<fltM<<endl;
        dataOutStr<<"seed, X ="<<fltX<<endl;
		dataOutStr<<"Power Residue constant, C ="<<fltC<<endl<<endl;

    	// get the pseudorandom numbers as an array
    	fltRanNum = fltX;
    	for (intIndex=1; intIndex<=intN; intIndex++) {
    		fltRanNum = fmodf(fltC*fltRanNum,fltM);
    		fltNormdRanNum[intIndex] = fltRanNum/fltMaxNum;
    	}
    	
    	// output
    	intNumInBlock = (intN/6)*6;
    	intNumLeftOver = intN-intNumInBlock;
    	for (intBeginRow=1; intBeginRow<=intNumInBlock; intBeginRow=intBeginRow+6) {
    		for (intIndex=intBeginRow; intIndex<=intBeginRow+5; intIndex++) {
    			dataOutStr<<setiosflags(ios::fixed)<<setw(8)<<setprecision(5)<<fltNormdRanNum[intIndex];
    		}
    		dataOutStr<<endl;
    	}
    	for (intIndex=intNumInBlock+1; intIndex<=intN; intIndex++) {
    		dataOutStr<<setiosflags(ios::fixed)<<setw(8)<<setprecision(5)<<fltNormdRanNum[intIndex];
    	}
    	dataOutStr<<endl;

		// another run?
		goout = query();

	}
    // close the output file
	fileClose("m354p1c.out", &dataOutStr);
};