TestGaussian.m

This is TestGaussian.m in view mode; [Download] [Up]
//
// BinTest
//
// This program tests the Random classes' abilities to generate
// random gaussian variables.
//


#import "Gaussian.h"
#import "StandardEngine.h"
#import "ElkinsEngine.h"
#import "R250Engine.h"
#import <math.h>
#import <stdio.h>
#import <stdlib.h>


#define DEFAULT_N	1000
#define DEFAULT_M	10
#define RANGE		3.5

#define MAX_BAR		40


//
// test_gaussian()
//

int test_gaussian(id myRand, int m, int n)
{
    int			max = 0;		// Size of largest bin.
    int			i;			// Generic loop variable.
    int			j;			// Generic loop variable.
    double		x;			// The random number.
    int			slot;			// The sorted position of a number.
    int			barsize;		// The length of a result bar.
    int			*binneg;		// The bin array (dynamically allocated).
    int			*binpos;		// The bin array (dynamically allocated).
    int			which;
    int			total = 0;
    
    //
    // Allocate and initialize slot array:
    //
    
    binneg = (int *)malloc(m * sizeof(int));
    binpos = (int *)malloc(m * sizeof(int));
    for(i = 0; i < m; i++) {
        binneg[i] = 0;
	binpos[i] = 0;
    }
    
    //
    // Make n random numbers:
    //
        
    for(i = 0; i < n; i++) {
        do {
	    x = [myRand gaussian];
//	    printf("Made gaussian %f\n", x);
	} while (x == 0.0);
	
	total++;
	which = (x > 0);
	slot = (int)floor(fabs(x) * m / RANGE);
	if((slot >= 0) && (slot < m) && (x != 0.0)) {		// If it is in range,
	    switch(which) {					//   count it.
	        case 0:
		    if((++binneg[slot]) > max) {
			max = binneg[slot];
		    }
		    break;
		case 1:
		    if((++binpos[slot]) > max) {
			max = binpos[slot];
		    }
		    break;
		default:
		    printf("Why is which = %d?\n", which);
		    break;
	    }
	}
    }
    
    //
    // Print the results:
    //
    
    for(i = m - 1; i >= 0; i--) {
        printf("%+1.6f to %+1.6f: %6d:", -((i + 1) * (RANGE / m)),
		-(i * (RANGE / m)), binneg[i]);
	
	barsize = (int)((double)binneg[i] / ((double)max / (double)MAX_BAR));
	for(j = 0; j < barsize; j++)
	    printf("*");
	printf("\n");
    }
    for(i = 0; i < m; i++) {
        printf("%+1.6f to %+1.6f: %6d:", i * (RANGE / m), (i + 1) * (RANGE / m), binpos[i]);
	
	barsize = (int)((double)binpos[i] / ((double)max / (double)MAX_BAR));
	for(j = 0; j < barsize; j++)
	    printf("*");
	printf("\n");
    }
    
    printf("Max   = %6d\n", max);
    printf("Total = %6d\n", total);
    
    //
    // Clean up:
    //
    
    free(binneg);
    free(binpos);
    
    //
    // Return to caller:
    //
    
    return 0;
}


//
// main()
//

int main(int argc, char *argv[])
{
    id			myRand;			// Random number generator.
    int			n;			// Quantity of numbers to generate.
    int			m;			// Number of bins in which to sort them.
    
    myRand = [[Gaussian alloc] initEngineClass:[ElkinsEngine class]];
    
    switch(argc) {
        case 2:
	    sscanf(argv[1], "%d", &n);
	    m = DEFAULT_M;
	    break;
	case 3:
	    sscanf(argv[1], "%d", &n);
	    sscanf(argv[2], "%d", &m);
	    break;
	default:
	    n = DEFAULT_N;
	    m = DEFAULT_M;
	    break;
    }
    
    printf("TestGaussian: Testing StandardEngine class:\n");
    myRand = [[Gaussian alloc] initEngineClass:[StandardEngine class]];
    test_gaussian(myRand, m, n);
    printf("\n\n");
    [myRand free];
    
    printf("TestGaussian: Testing ElkinsEngine class:\n");
    myRand = [[Gaussian alloc] initEngineClass:[ElkinsEngine class]];
    test_gaussian(myRand, m, n);
    printf("\n\n");
    [myRand free];
    
    printf("TestGaussian: Testing R250Engine class:\n");
    myRand = [[Gaussian alloc] initEngineClass:[R250Engine class]];
    test_gaussian(myRand, m, n);
    printf("\n\n");
    [myRand free];
    
    return 0;
}


//
// End of file.
//
These are the contents of the former NiCE NeXT User Group NeXTSTEP/OpenStep software archive, currently hosted by Netfuture.ch.