diff --git a/MDD/Makefile b/MDD/Makefile index 7cb5f2d377810cec4cc7708a78c7459523454edb..bbc45bb37e86e1ecb9d6ce403d76b2e7b471c4ef 100644 --- a/MDD/Makefile +++ b/MDD/Makefile @@ -5,9 +5,14 @@ include ../Make_include ######################################################################## # define general include and system library ALLINC = -I. -LIBS += -mkl -L$L -lgenfft $(LIBSM) -CFLAGS += -I$(MKLROOT)/include - +#BLAS libs with Intel compiler +#LIBS += -mkl -L$L -lgenfft $(LIBSM) +#General BLAS library +LIBS += -L$L -lgenfft $(LIBSM) +#General BLAS library +LIBS += $(BLAS) + +#CFLAGS += -I$(MKLROOT)/include #LIBS += -lblas -llapack -L$L -lgenfft $(LIBSM) -lc -lm all: mdd diff --git a/MDD/deconvolve.c b/MDD/deconvolve.c index 2ef6d49937a86782f604b88ddcb5623eb792029c..6b88d68d438a081ef7328e2aaa6c5dea1408ab27 100644 --- a/MDD/deconvolve.c +++ b/MDD/deconvolve.c @@ -3,7 +3,6 @@ #include <assert.h> #include <math.h> #include <string.h> -#include<mkl_cblas.h> typedef struct { /* complex number */ float r,i; diff --git a/MDD/readShotData.c b/MDD/readShotData.c index ce774364c552931aba49bda7fb25ae8aa952bbfd..20f953e2ac140a6032ed3c71198e482ee06bb69f 100644 --- a/MDD/readShotData.c +++ b/MDD/readShotData.c @@ -5,9 +5,6 @@ #include "segy.h" #include <assert.h> -extern FILE *fopen64 (__const char *__restrict __filename, - __const char *__restrict __modes); - typedef struct { /* complex number */ float r,i; } complex; @@ -34,7 +31,7 @@ int readShotData(char *filename, float xmin, float dx, float *xrcv, float *xsrc, /* Reading first header */ if (filename == NULL) fp = stdin; - else fp = fopen64( filename, "r" ); + else fp = fopen( filename, "r" ); if ( fp == NULL ) { fprintf(stderr,"input file %s has an error\n", filename); perror("error in opening file: "); diff --git a/README b/README index 6b9f9f5bd6bd80c034c062022d7a334f3266a46e..16f6b372e69c416c32efe670171e40dc75afd046 100644 --- a/README +++ b/README @@ -88,6 +88,18 @@ To reproduce the Figures shown in the GEOPHYICS paper "Implementation of the Mar To reproduce the Figures shown in the Scientific Reports paper "Virtual acoustics in inhomogeneous media with single-sided access" the scripts in marchenko/demo/ScientificReports directory can be used. The README in this directory gives more instructions and guidelines. +MDD +--- +The MDD kernels depend on BLAS and LAPACK calls. Free downloads of these libraries can be found on + +https://www.netlib.org/blas/index.html +https://www.netlib.org/lapack/index.html + +If you are running on Intel processors you can download (for free) Intel's highly optimised MKL package: + +https://software.intel.com/en-us/mkl/choose-download + + MISC ---- Other make commands which can be useful: