Skip to content
Snippets Groups Projects
Commit e2ca12c1 authored by Jan Thorbecke's avatar Jan Thorbecke
Browse files

linking with MKL for MDD

parent aced96cd
No related branches found
No related tags found
No related merge requests found
...@@ -5,9 +5,14 @@ include ../Make_include ...@@ -5,9 +5,14 @@ include ../Make_include
######################################################################## ########################################################################
# define general include and system library # define general include and system library
ALLINC = -I. ALLINC = -I.
LIBS += -mkl -L$L -lgenfft $(LIBSM) #BLAS libs with Intel compiler
CFLAGS += -I$(MKLROOT)/include #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 #LIBS += -lblas -llapack -L$L -lgenfft $(LIBSM) -lc -lm
all: mdd all: mdd
......
...@@ -3,7 +3,6 @@ ...@@ -3,7 +3,6 @@
#include <assert.h> #include <assert.h>
#include <math.h> #include <math.h>
#include <string.h> #include <string.h>
#include<mkl_cblas.h>
typedef struct { /* complex number */ typedef struct { /* complex number */
float r,i; float r,i;
......
...@@ -5,9 +5,6 @@ ...@@ -5,9 +5,6 @@
#include "segy.h" #include "segy.h"
#include <assert.h> #include <assert.h>
extern FILE *fopen64 (__const char *__restrict __filename,
__const char *__restrict __modes);
typedef struct { /* complex number */ typedef struct { /* complex number */
float r,i; float r,i;
} complex; } complex;
...@@ -34,7 +31,7 @@ int readShotData(char *filename, float xmin, float dx, float *xrcv, float *xsrc, ...@@ -34,7 +31,7 @@ int readShotData(char *filename, float xmin, float dx, float *xrcv, float *xsrc,
/* Reading first header */ /* Reading first header */
if (filename == NULL) fp = stdin; if (filename == NULL) fp = stdin;
else fp = fopen64( filename, "r" ); else fp = fopen( filename, "r" );
if ( fp == NULL ) { if ( fp == NULL ) {
fprintf(stderr,"input file %s has an error\n", filename); fprintf(stderr,"input file %s has an error\n", filename);
perror("error in opening file: "); perror("error in opening file: ");
......
...@@ -88,6 +88,18 @@ To reproduce the Figures shown in the GEOPHYICS paper "Implementation of the Mar ...@@ -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. 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 MISC
---- ----
Other make commands which can be useful: Other make commands which can be useful:
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment