Das deutsche Python-Forum

Code: Alles auswählen

r = 3000
treffer = 0
quadtreffer = r**2

for i in range(r):
    x2 = (i+0.5)**2
    for j in range(r):
        y2 = (j+0.5)**2
        #print("x2: ", x2, " y2: ", y2, '\n')
        if (x2 + y2 <= quadtreffer):
            treffer += 1

print("%.5f" %(4*treffer/quadtreffer))

Code: Alles auswählen

#include <iostream>
#include <math.h>
using namespace std;

int main(int argc, char **argv)
{
	int r = 3000, treffer = 0;
	double quadtreffer = pow(r,2);
	double x2 = 0.0, y2 = 0.0;

	for(int i=0; i<r; i++) {
		x2 = pow((i+0.5),2);
		for(int j=0; j<r; j++) {
			y2 = pow((j+0.5),2);
			//cout << "x2: " << x2 << " y2: " << y2 << endl;
			if((x2+y2) <= quadtreffer) {
				treffer++;
			}
		}
	}
	cout << (4*treffer/quadtreffer);
	return 0;
}

Code: Alles auswählen

def calculate_pi_a(radius=RADIUS):
    hit_zone = radius**2
    hit_count = 0
    for i in xrange(radius):
        x = (i + 0.5)**2
        for j in xrange(radius):
            y = (j + 0.5)**2
            if x + y <= hit_zone:
                hit_count += 1
    return 4.0 * hit_count / hit_zone

Code: Alles auswählen

#define RADIUS  3000
#define HIT_ZONE  ((double) RADIUS * RADIUS)

double calculate_pi_a(void)
{
    uint16_t i, j;
    double x, y;
    uint32_t hit_count = 0;

    for (i = 0; i < RADIUS; ++i) {
        x = i + 0.5;
        x *= x;
        for (j = 0; j < RADIUS; ++j) {
            y = j + 0.5;
            y *= y;
            if (x + y <= HIT_ZONE) {
                ++hit_count;
            }
        }
    }
    return 4 * hit_count / HIT_ZONE;
}

Code: Alles auswählen

def calculate_pi_b(radius=RADIUS):
    hit_zone = radius**2
    values = [(i + 0.5)**2 for i in xrange(radius)]
    hit_count = 0
    for x in values:
        for y in values:
            if x + y <= hit_zone:
                hit_count += 1
    return 4.0 * hit_count / hit_zone

Code: Alles auswählen

double calculate_pi_b(void)
{
    uint16_t i, j;
    double x, *values;
    uint32_t hit_count = 0;

    values = malloc(RADIUS * sizeof(*values));
    for (i = 0; i < RADIUS; ++i) {
        x = i + 0.5;
        values[i] = x * x;
    }

    for (i = 0; i < RADIUS; ++i) {
        for (j = 0; j < RADIUS; ++j) {
            if (values[i] + values[j] <= HIT_ZONE) {
                ++hit_count;
            }
        }
    }

    free(values);
    return 4 * hit_count / HIT_ZONE;
}

Code: Alles auswählen

def calculate_pi_numpy(radius=RADIUS):
    hit_zone = radius**2
    A = ((np.arange(radius) + 0.5)**2).reshape((1, -1)).repeat(radius, 0)
    A += A.copy().transpose()
    hit_count = (A <= hit_zone).sum()
    return 4.0 * hit_count / hit_zone

Code: Alles auswählen

double calculate_pi_c(void)
{
    uint16_t i, j;
    double x, *values;
    uint32_t hit_count = 0;

    values = malloc(RADIUS * sizeof(*values));
    for (i = 0; i < RADIUS; ++i) {
        x = i + 0.5;
        values[i] = x * x;
    }

    for (i = 0; i < RADIUS; ++i) {
        for (j = i + 1; j < RADIUS; ++j) {
            if (values[i] + values[j] <= HIT_ZONE) {
                ++hit_count;
            }
        }
    }
    hit_count *= 2;

    for (i = 0; i < RADIUS; ++i) {
        if (values[i] * 2 <= HIT_ZONE) {
            ++hit_count;
        }
    }

    free(values);
    return 4 * hit_count / HIT_ZONE;
}

Code: Alles auswählen

#define HIT_ZONE_FIXPOINT  ((RADIUS * RADIUS) << 2)

double calculate_pi_d(void)
{
    uint16_t i, j;
    uint32_t *values;   /* Fixed point values with two bits fractional part. */
    uint32_t hit_count = 0;

    values = malloc(RADIUS * sizeof(*values));
    /* 
     * `j` is a fixed point value with one bit fractional part here.
     * So it starts at 0.5 and is incremented by 1.0 each loop iteration.
     */
    for (i = 0, j = 1; i < RADIUS; ++i, j += 2) {
        values[i] = j * j;
    }

    for (i = 0; i < RADIUS; ++i) {
        for (j = i + 1; j < RADIUS; ++j) {
            if (values[i] + values[j] <= HIT_ZONE_FIXPOINT) {
                ++hit_count;
            }
        }
    }
    hit_count *= 2;

    for (i = 0; i < RADIUS; ++i) {
        if (values[i] * 2 <= HIT_ZONE_FIXPOINT) {
            ++hit_count;
        }
    }

    free(values);
    return 4.0 * hit_count / (HIT_ZONE_FIXPOINT >> 2);
}

Code: Alles auswählen

   10 REM@ £PROTOCOL
   20 REM@ £CONSTANT R,HZ
   30 REM@ £WORD I=FAST,J=FAST
   40 REM@ £FASTFOR
   50 TI$="000000":R=2999:HZ=(R+1)*(R+1):HC=0:DIM X(R)
   60 FOR I=0 TO R:X=I+0.5:X(I)=X*X:NEXT
   70 FOR I=0 TO R-1:FOR J=I+1 TO R:IF X(I)+X(J)<=HZ THEN HC=HC+1
   80 NEXT:NEXT:HC=HC*2
   90 FOR I=0 TO R:IF X(I)*2<=HZ THEN HC=HC+1
  100 NEXT
  110 PRINT 4*HC/HZ,TI$

Code: Alles auswählen

print("%.5f" %(4*float(treffer)/quadtreffer))