|
|
@ -1,8 +1,10 @@
|
|
|
|
from FFTMA import gen
|
|
|
|
from FFTMA import gen
|
|
|
|
import numpy as np
|
|
|
|
import numpy as np
|
|
|
|
import unittest
|
|
|
|
import unittest
|
|
|
|
|
|
|
|
import time
|
|
|
|
|
|
|
|
|
|
|
|
def generate(N):
|
|
|
|
def generate(N):
|
|
|
|
|
|
|
|
start_time = time.time()
|
|
|
|
nx, ny, nz = N,N,N
|
|
|
|
nx, ny, nz = N,N,N
|
|
|
|
dx, dy, dz = 1.0, 1.0, 1.0
|
|
|
|
dx, dy, dz = 1.0, 1.0, 1.0
|
|
|
|
seed= 1548762 #rdi(10000,99999)
|
|
|
|
seed= 1548762 #rdi(10000,99999)
|
|
|
@ -26,7 +28,10 @@ def generate(N):
|
|
|
|
variance=3.5682389
|
|
|
|
variance=3.5682389
|
|
|
|
typ=3
|
|
|
|
typ=3
|
|
|
|
|
|
|
|
|
|
|
|
return gen(nx, ny, nz, dx, dy, dz, seed, variograms, mean, variance, typ, 8)
|
|
|
|
k = gen(nx, ny, nz, dx, dy, dz, seed, variograms, mean, variance, typ, 8)
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
print(f"Generation with N = {N} time = {time.time() - start_time}s")
|
|
|
|
|
|
|
|
return k
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
def test(N):
|
|
|
|
def test(N):
|
|
|
@ -35,8 +40,6 @@ def test(N):
|
|
|
|
comparison = k == k_correct
|
|
|
|
comparison = k == k_correct
|
|
|
|
return comparison.all()
|
|
|
|
return comparison.all()
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
class TestGeneration(unittest.TestCase):
|
|
|
|
class TestGeneration(unittest.TestCase):
|
|
|
|
def test_8(self):
|
|
|
|
def test_8(self):
|
|
|
|
self.assertTrue(test(8))
|
|
|
|
self.assertTrue(test(8))
|
|
|
@ -46,3 +49,6 @@ class TestGeneration(unittest.TestCase):
|
|
|
|
|
|
|
|
|
|
|
|
def test_32(self):
|
|
|
|
def test_32(self):
|
|
|
|
self.assertTrue(test(32))
|
|
|
|
self.assertTrue(test(32))
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
def test_64(self):
|
|
|
|
|
|
|
|
self.assertTrue(test(64))
|