//array2d.C Oguz Yetkin 10/28/1998
//2-D dynamic array, use at your own pleasure/risk
//initially written for Physics499, recycled for Physics 691
//array2d implementation

#include<iostream.h>
#include<stdio.h>
//#include "ran2.h"
#include "array2d.h"

int debug=0;

double gauss(int gaussnum){
        double sum=0;
        for(int i=0;i<gaussnum;i++){
                sum+=(rand()/(double)RAND_MAX);
        }
        sum=sum-(gaussnum/2);
        return sum;
}





array2d::array2d(int h, int w){

        if(debug){
                cout<<"\nCONSTRUCTING ARRAY2d, memory being allocated!!!"<<endl;
        }

        width=w;
        height=h;
        error=0;
        //base=new double *[width];        //do these need to be reversed ???
        base=new double * [height];        
        for(int i=0;i<height;i++){
		base[i]=new double [width];
        }
}


array2d::~array2d(){


/* commented out 4/24/1998 --deletes are messing something up
        for(int i=0;i<height;i++){
                delete[]base[i];
        }       
  */

        delete[] base;
}


double  array2d::getvalue(int i, int j){

        if(i>=height||j>=width){
                error=1;
                cerr<<"\ngetvalue: illegal subscripts!"<<i<<" "<<j<<endl;
                cerr<<"\nthe array size is: "<<height<<" "<<width<<endl;
                return 1;
        }
        return base[i][j];

}



int array2d::putvalue(int i, int j, double  value){

        if(i>=height||j>=width){
                error=1;
                cerr<<"\nputvalue: illegal subscripts!"<<i<<" "<<j<<endl; 
                cerr<<"\nthe array is: "<<height<<" "<<width<<endl;
                return -1;
        }

        base[i][j]=value;
        return 1;               
}

void array2d::dump(){

        //print the array
        printf("\nthe array is %d by %d\n",height, width);      
        for(int i=0;i<height;i++){
                for(int j=0;j<width;j++){
                        //printf("%f ",base[i][j]);
                        printf("%.15f ",getvalue(i,j));
                }
                printf("%c",'\n');
        }
}


void array2d::dump(int beg, int end){   //restricts width


        for(int i=0;i<height;i++){
                for(int j=beg;j<end;j++){


                        printf("%f\t",getvalue(i,j));

                }
                printf("%c",'\n');

        }
}


void array2d::randomize(){

        for(int i=0;i<height;i++){


                for(int j=0;j<width;j++){





                        //base[i][j]=rand(0);





                                //putvalue(i,j,rand(0));


                                //putvalue(i,j,rand()); //for windows
                                  putvalue(i,j,gauss());        //changed 10/18/1998

               }
                //cout<<endl; 
        }       
}

//added 1/10/1999
void array2d::bzero(){

	for(int i=0;i<height;i++){ //the names 'height' and 'width' are probably meaningless

		for(int j=0;j<width;j++){

			putvalue(i,j,0);

		}

	}

}