rm files
Oli
3 years ago
parent
562b948508
commit
6b97a9dfeb
@ -1,34 +0,0 @@
|
||||
from FFTMA import gen
|
||||
from refine import refine as ref
|
||||
import numpy as np
|
||||
from random import randint as rdi
|
||||
|
||||
N = 128
|
||||
|
||||
for i in range(1):
|
||||
nx, ny, nz = N, N, N
|
||||
dx, dy, dz = 1.0, 1.0, 1.0
|
||||
seed = 1548762 # rdi(10000,99999)
|
||||
var = 1
|
||||
vario = 2
|
||||
alpha = 1
|
||||
lcx = 2
|
||||
lcy = 4
|
||||
lcz = 16
|
||||
ap1x = 1
|
||||
ap1y = 0
|
||||
ap1z = 0
|
||||
ap2x = 0
|
||||
ap2y = 1
|
||||
ap2z = 0
|
||||
|
||||
v1 = (var, vario, alpha, lcx, lcy, lcz, ap1x, ap1y, ap1z, ap2x, ap2y, ap2z)
|
||||
variograms = [v1]
|
||||
|
||||
mean = 15.3245987
|
||||
variance = 3.5682389
|
||||
typ = 1
|
||||
|
||||
k = gen(nx, ny, nz, dx, dy, dz, seed, variograms, mean, variance, typ)
|
||||
|
||||
np.save("out" + str(i) + ".npy", ref(k, 4, 4, 4))
|
||||
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
@ -1,73 +0,0 @@
|
||||
#include <stdlib.h>
|
||||
#include <string.h>
|
||||
#include "genlib.h"
|
||||
#include "geostat.h"
|
||||
|
||||
#ifndef _TOOLSFFTMA_H
|
||||
#define _TOOLSFFTMA_H
|
||||
|
||||
/* Create december, the 16th 2002 */
|
||||
/* Modified december, the 9th 2002 */
|
||||
|
||||
/* List of functions: */
|
||||
/* ------------------ */
|
||||
/* kgeneration, FFTMA2, prebuild_gwn, build_real, clean_real */
|
||||
|
||||
|
||||
/*FAST FOURIER TRANSFORM Pressure Simulation */
|
||||
/*Turns a Gaussian white noise vector into a */
|
||||
/*spatially correlated vector */
|
||||
/*input: */
|
||||
/*variogram: structure defining the variogram */
|
||||
/* model */
|
||||
/*grid: structure defining the grid */
|
||||
/*n: vector with the number of cells along the */
|
||||
/* X, Y and Z axes for the underlying grid */
|
||||
/* i = [0 1 2] */
|
||||
/* --> 0 0 0 : n will be computed and */
|
||||
/* updated as output */
|
||||
/* --> nx ny nz: these dimensions are used */
|
||||
/*realin: structure defining a realization - */
|
||||
/* must be a Gaussian white noise */
|
||||
/*gradient: macroscopic gradient pression vector */
|
||||
/*output: */
|
||||
/*realout: structure defining a realization - */
|
||||
/*realout2: structure defining a pressure field */
|
||||
/*realout3: structure defining a xvelocity field */
|
||||
/*realout4: structure defining a yvelocity field */
|
||||
/*realout5: structure defining a zvelocity field */
|
||||
|
||||
void FFTMA2(struct vario_mod variogram,struct grid_mod grid,int n[3],struct realization_mod *realin,struct realization_mod *realout);
|
||||
|
||||
/* prebuild_gwn */
|
||||
/* Produce a first construction in real space of the Gaussian white noise */
|
||||
/* grid: structure defining the grid */
|
||||
/*n: vector with the number of cells along the */
|
||||
/* X, Y and Z axes for the underlying grid */
|
||||
/* i = [0 1 2] */
|
||||
/* --> 0 0 0 : n will be computed and */
|
||||
/* updated as output */
|
||||
/* --> nx ny nz: these dimensions are used */
|
||||
/*realin: structure defining a realization - */
|
||||
/* must be a Gaussian white noise */
|
||||
/*realization: structure defining a realization*/
|
||||
|
||||
void prebuild_gwn(struct grid_mod grid,int n[3],struct realization_mod *realin,double *realization,int solver);
|
||||
|
||||
/* build_real */
|
||||
/* build a realization in the spectral domain */
|
||||
/*n: vector with the number of cells along the */
|
||||
/* X, Y and Z axes for the underlying grid */
|
||||
/* i = [0 1 2] */
|
||||
/* --> 0 0 0 : n will be computed and */
|
||||
/* updated as output */
|
||||
/* --> nx ny nz: these dimensions are used */
|
||||
/*covar: vector defining the covariance in spectral domain */
|
||||
/*realization: vector defining the real part */
|
||||
/*ireal: vector defining the i-part */
|
||||
|
||||
void build_real(int n[3],int NTOT,double *covar,double *realization,double *ireal);
|
||||
|
||||
void clean_real(struct realization_mod *realin,int n[3],struct grid_mod grid,double *vectorresult,struct realization_mod *realout);
|
||||
|
||||
#endif // define _TOOLSFFTMA_H
|
||||
@ -1,23 +0,0 @@
|
||||
#include <string.h>
|
||||
#include "geostat.h"
|
||||
#include "pressure.h"
|
||||
|
||||
|
||||
#ifndef _TOOLSFFTPSIM_H
|
||||
#define _TOOLSFFTPSIM_H
|
||||
|
||||
/* Create december, the 16th 2002 */
|
||||
|
||||
/* List of functions: */
|
||||
/* ------------------ */
|
||||
/* FFTPSim, FFTPressure, clean_real2 */
|
||||
|
||||
/*Functions */
|
||||
/*----------*/
|
||||
|
||||
void FFTPSim(struct vario_mod variogram,struct statistic_mod stat,struct grid_mod grid,int n[3],struct realization_mod *realin,struct pressure_mod gradient,struct realization_mod *realout,struct realization_mod *realout2,struct realization_mod *realout3,struct realization_mod *realout4,struct realization_mod *realout5);
|
||||
|
||||
void normAxes(double *vec, double *normed);
|
||||
void waveVectorCompute1D(int n,double *vec);
|
||||
|
||||
#endif // define _TOOLSFFTPSIM_H
|
||||
@ -1,24 +0,0 @@
|
||||
#include <stdlib.h>
|
||||
#include <string.h>
|
||||
#include "genlib.h"
|
||||
#include "geostat.h"
|
||||
#include "pressure.h"
|
||||
|
||||
#ifndef _TOOLSIO_H
|
||||
#define _TOOLSIO_H
|
||||
|
||||
/* Create december, the 16th 2002 */
|
||||
/* Modified december, the 9th 2002 */
|
||||
|
||||
/* List of functions: */
|
||||
/* ------------------ */
|
||||
/* testmemory */
|
||||
|
||||
|
||||
/*FUNCTIONS*/
|
||||
/*----------*/
|
||||
|
||||
/* Allocation test */
|
||||
void testmemory(double *realint);
|
||||
|
||||
#endif // define _TOOLSIO_H
|
||||
Loading…
Reference in New Issue