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

linking with MKL for MDD

parent b2d0a21a
No related branches found
No related tags found
No related merge requests found
......@@ -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
......
......@@ -3,7 +3,6 @@
#include <assert.h>
#include <math.h>
#include <string.h>
#include<mkl_cblas.h>
typedef struct { /* complex number */
float r,i;
......
......@@ -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: ");
......
......@@ -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:
......
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