##******************************************************************************
##  Copyright(C) 2006-2013 Intel Corporation. All Rights Reserved.
##  
##  The source code, information  and  material ("Material") contained herein is
##  owned  by Intel Corporation or its suppliers or licensors, and title to such
##  Material remains  with Intel Corporation  or its suppliers or licensors. The
##  Material  contains proprietary information  of  Intel or  its  suppliers and
##  licensors. The  Material is protected by worldwide copyright laws and treaty
##  provisions. No  part  of  the  Material  may  be  used,  copied, reproduced,
##  modified, published, uploaded, posted, transmitted, distributed or disclosed
##  in any way  without Intel's  prior  express written  permission. No  license
##  under  any patent, copyright  or  other intellectual property rights  in the
##  Material  is  granted  to  or  conferred  upon  you,  either  expressly,  by
##  implication, inducement,  estoppel or  otherwise.  Any  license  under  such
##  intellectual  property  rights must  be express  and  approved  by  Intel in
##  writing.
##  
##  *Third Party trademarks are the property of their respective owners.
##  
##  Unless otherwise  agreed  by Intel  in writing, you may not remove  or alter
##  this  notice or  any other notice embedded  in Materials by Intel or Intel's
##  suppliers or licensors in any way.
##
##******************************************************************************
##  Content:
##      Intel(R) Math Kernel Library Poisson Solver examples creation
##      and run
##
##******************************************************************************

help:
	@echo "Usage: make {libia32|soia32|libintel64|sointel64} [function=name]"
	@echo "[compiler=compiler_name] [interface=interface_name] [threading=threading_name]"
	@echo "[parallel=parallel_name] [omp=omp_name][gnu_path=gnu_lib_path]"
	@echo
	@echo "name           - function name. Please see pdepoissonf.lst file"
	@echo
	@echo "compiler_name  - can be gnu, pgi, or intel. Default value is intel."
	@echo "                 Intel (R) Fortran Compiler as default."
	@echo "                 If compiler=gnu, then GNU gfortran compiler will be used."
	@echo "                 Note: GNU g77 compiler is not supported."
	@echo "                 If compiler=pgi, then PGI Fortran 95 compiler will be used."
	@echo
	@echo "interface_name - can be lp64 or ilp64 for intel64. Default value is lp64."
	@echo
	@echo "threading_name - can be parallel or sequential. Default value is parallel."
	@echo "parallel_name  - can be intel, pgi (only if compiler=pgi), or gnu (only if compiler=gnu)."
	@echo "                 Default value is intel."
	@echo "omp_name       - can be iomp5 if parallel=intel, or"
	@echo "                 can be iomp5 or gomp if parallel=gnu, or"
	@echo "                 can be pgmp if parallel=pgi."
	@echo "                 Default value is iomp5."
	@echo "gnu_lib_path   - If you are using gnu threading layer, specify path to GNU libraries,"
	@echo "                 such as libgomp and libgfortran, with gnu_path variable."
	@echo "                 Default value is /usr/lib."

##------------------------------------------------------------------------------
## examples of usage:
##
## make libia32 function=d_Helmholtz_2D_f  - build by Intel(R) Fortran Compiler (as default) and
##                                           run d_Helmholtz_2D_f example for 32-bit applications,
##                                           static linking
##
## make libintel64                         - build by Intel(R) Fortran Compiler (as default) and run all
##                                           examples of MKL PDE POISSON for Intel(R) 64 processor family
##                                           applications, static linking
##
##------------------------------------------------------------------------------

include pdepoissonf.lst

ifndef function
function = $(PDE_POISSON_F)
endif

ifeq (,$(filter gnu pgi,$(compiler)))
   override compiler=intel
   override parallel=intel
endif

ifndef interface
interface=lp64
endif

ifndef threading
threading=parallel
endif

ifndef $(gnu_path)
   gnu_path=/usr/lib
endif

RES = $(addsuffix .res ,$(function))
MOD = $(addsuffix .mod ,$(modules))

ifndef MKLROOT
MKLROOT = ../..
endif
MKL_PATH = "$(MKLROOT)/lib/$(_IA)"
CMPLR_PATH = "$(MKLROOT)/../compiler/lib/$(_IA)"

ifndef $(gnu_path)
   gnu_path=/usr/lib
endif

ifeq ($(compiler),gnu)
   FC=gfortran
   ifeq ($(FC),gfortran)
      IFACE_COMP_PART=gf
   else
      IFACE_COMP_PART=intel
   endif
else
   ifeq ($(compiler),pgi)
      FC=pgf95
   else
      override FC=ifort
   endif
   IFACE_COMP_PART=intel
endif

ifeq (,$(filter pgi,$(parallel)))
   override parallel=intel
   override omp=iomp5
else
   ifeq ($(parallel),gnu)
      ifneq ($(omp),gomp)
         override omp=iomp5
      endif
   else
      override omp=pgmp
   endif
endif

ifeq ($(parallel),pgi)
   IFACE_THREADING_PART=pgi
else
   IFACE_THREADING_PART=intel
endif


ifeq ($(interface),ilp64)
   IFACE_LIB=$(MKL_PATH)/libmkl_$(IFACE_COMP_PART)_ilp64.$(EXT)
   FOPTS = -i8
   ifeq ($(compiler),gnu)
      override FC=gfortran
      FOPTS=-fdefault-integer-8
   endif
else
   IFACE_LIB=$(MKL_PATH)/libmkl_$(IFACE_COMP_PART)_lp64.$(EXT)
   FOPTS =
endif

ifeq ($(_IA),ia32)
   ifeq ($(compiler),intel)
       SPEC_OPT=-mia32
#This option is required by Intel(R) 11.0 compiler to produce workable binaries for Pentium(R) III.
#If you don't need it, you can remove this option.
   endif
IFACE_LIB=$(MKL_PATH)/libmkl_$(IFACE_COMP_PART).$(EXT)
FOPTS =
endif

ifeq ($(threading),sequential)
   THREADING_LIB=$(MKL_PATH)/libmkl_sequential.$(EXT)
   OMP_LIB =
   GFORTRAN_LIB=
else
   THREADING_LIB=$(MKL_PATH)/libmkl_$(IFACE_THREADING_PART)_thread.$(EXT)
   ifeq ($(parallel),pgi)
      FOPTS += -mp
   endif
   ifeq ($(omp),gomp)
      OMP_LIB = -L$(gnu_path) -l$(omp)
   else
      ifeq ($(omp),pgmp)
        OMP_LIB =
      else
        OMP_LIB = -L$(CMPLR_PATH) -l$(omp)
      endif
   endif
endif

CORE_LIB=$(MKL_PATH)/libmkl_core.$(EXT)

MKL_LIBS=-Wl,--start-group $(IFACE_LIB) $(THREADING_LIB) $(CORE_LIB) -Wl,--end-group $(OMP_LIB)

ifeq ($(_IA),ia32)
RES_DIR=_results/$(compiler)_$(threading)_$(_IA)_$(RES_EXT)$Z
else
RES_DIR=_results/$(compiler)_$(interface)_$(threading)_$(_IA)_$(RES_EXT)$Z
endif

libia32 lib32: clean
	$(MAKE) $(MOD) $(RES) _IA=ia32 EXT=a RES_EXT=lib
soia32 so32: clean
	$(MAKE) $(MOD) $(RES) _IA=ia32 EXT=so RES_EXT=so
libintel64 libem64t: clean
	$(MAKE) $(MOD) $(RES) _IA=intel64 EXT=a RES_EXT=lib
sointel64 soem64t: clean
	$(MAKE) $(MOD) $(RES) _IA=intel64 EXT=so RES_EXT=so
clean:
	rm -f *.mod *.o

#-------------------------------------------------------------------------------

vpath %.f90 source
vpath %.f90 $(MKLROOT)/include

$(MOD): %.mod: %.f90
	mkdir -p ./$(RES_DIR)
	$(FC) $(SPEC_OPT) $(FOPTS) -w -c $<

$(RES): %.res: %.f90
	mkdir -p ./$(RES_DIR)
	$(FC) $(SPEC_OPT) $(FOPTS) -w $< -L$(MKL_PATH) $(MKL_LIBS) -lpthread -ldl -o $(RES_DIR)/$*.out
	export LD_LIBRARY_PATH=$(MKL_PATH):$(LD_LIBRARY_PATH):$(CMPLR_PATH); $(RES_DIR)/$*.out >$(RES_DIR)/$@
#-------------------------------------------------------------------------------
