采用intel编译器并行编译执行杂化函数计算CP2K

2024-03-24 10:58

本文主要是介绍采用intel编译器并行编译执行杂化函数计算CP2K,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

在采用intel的ifort编译支持hybrid functional计算的CP2K时,尽管编译得到了cp2k.popt的可执行模块,但是在运行时常遇到如下的错误:
HFX_MEM_INFO would be printed out:
forrtl: severe (174): SIGSEGV, segmentation fault occurred
Image PC Routine Line Source
cp2k.popt 0000000001321BCE Unknown Unknown Unknown

经过多次的尝试,找到了一种解决办法,似乎是intel的ifort(至少11.1.072)还不支持ISO_C_BINDING。

下面将成功编译得到支持hybrid functional计算的CP2K的步骤给列出来,以作备份。
我用的编译器是intel的11.1.072/bin/intel64/ifort,并与openmpi/1.4.2绑定的mpif90,
用intel的编译器编译得到的fftw-3.2.2;
scalapck、lapack、blas库等用的是intel/11.1.072/mkl/lib/em64t带的。

1、下载libint-1.1.4.tar.gz ,它可从
http://libint.valeyev.net/
http://sourceforge.net/projects/libint/files/v1-releases/
下载到。

下载解压后运行该目录下的./configure。整个命令如下,
./configure --prefix='被安装到的目录' --with-libint-max-am=5 --with-libderiv-max-am1=4  F77=ifort 
这里“被安装的的目录”,根据自己的程序安装要求自行设置它的实际路径。
比如我实际的命令:
./configure --prefix='/home/zh/soft/libint-1.1.4' --with-libint-max-am=5 --with-libderiv-max-am1=4  F77=ifort 
然后就是
make
make install
这样就将编译好的libint 1.1.4版本安装到了所指定的目录。

2、下载CP2K。

http://cp2k.berlios.de/download.html
下载最新版本的CP2K。
我目前安装的版本是CP2K version 2.2.263 (Development Version)和CP2K version 2.1.397 (Development Version)。
下载解压后,进入cp2k/tools/hfx_tools/libint_tools目录下,参考该目录里面的“README_LIBINT”,运行如下的命令:
g++ -O3 -c libint_cpp_wrapper.cpp -I/'path_to_libint_include'
由于我将libint库安装到了/home/zh/soft/libint-1.1.4目录中,我实际运行的命令为:
g++ -O3 -c libint_cpp_wrapper.cpp -I/home/zh/soft/libint-1.1.4/include

然后进入CP2K的arch目录,即cp2k/arch,编辑里面的Linux-x86-64-intel.popt,它的里面的设置为:
# by default some intel compilers put temporaries on the stack
# this might lead to segmentation faults is the stack limit is set to low
# stack limits can be increased by sysadmins or e.g with ulimit -s 256000
# furthermore new ifort (10.0?) compilers support the option
# -heap-arrays 64
# add this to the compilation flags is the other options do not work
# The following settings worked for:
# - AMD64 Opteron
# - SUSE Linux Enterprise Server 10.0 (x86_64)
# - Intel(R) Fortran Compiler for Intel(R) EM64T-based applications, Version 10.0.025
# - AMD acml library version 3.6.0
# - MPICH2-1.0.5p4
# - SCALAPACK library          1.8.0
#
#
CC       = cc
CPP      =
FC       = mpif90
LD       = mpif90
AR       = ar -r
DFLAGS   = -D__INTEL -D__FFTSG -D__parallel -D__BLACS -D__SCALAPACK -D__FFTW3 -D__LIBINT -D__HAS_NO_ISO_C_BINDING
CPPFLAGS =
INTEL_INC = /usr/apps/isv/intel/11.1.072/include
LIBINT_INC = $(HOME)/soft/libint-1.1.4/include
FFTW_INC = /usr/apps/free/fftw/3.2.2_i/include
FCFLAGS  = $(DFLAGS) -I$(LIBINT_INC) -I$(INTEL_INC) -I$(FFTW_INC) -O2 -xW -heap-arrays 64 -funroll-loops -fpp -free
FCFLAGS2 = $(DFLAGS) -I$(LIBINT_INC) -I$(INTEL_INC) -O1 -xW  -heap-arrays 64 -fpp -free
LDFLAGS  = $(FCFLAGS) -I$(INTEL_INC)

FFTW_LIB = /usr/apps/free/fftw/3.2.2_i
INTEL_LIB = /usr/apps/isv/intel/11.1.072/mkl/lib/em64t
LIBINT = $(HOME)/soft/libint-1.1.4
LIBS     =-L$(INTEL_LIB)/ -lmkl_scalapack_lp64 -lmkl_blacs_openmpi_lp64 -lmkl_intel_lp64 -lmkl_sequential -lmkl_core \
           $(HOME)/soft/cp2k-2.1/tools/hfx_tools/libint_tools/libint_cpp_wrapper.o \
           $(LIBINT)/lib/libderiv.a \
           $(LIBINT)/lib/libint.a \
           -lstdc++ \
           $(FFTW_LIB)/lib/libfftw3.a

OBJECTS_ARCHITECTURE = machine_intel.o


graphcon.o: graphcon.F
        $(FC) -c $(FCFLAGS2) $<

qs_vxc_atom.o: qs_vxc_atom.F
        $(FC) -c $(FCFLAGS2) $<


编辑好了之后进入makefiles目录,然后运行如下命令:
make ARCH=Linux-x86-64-intel VERSION=popt
运行cp2k编译器来较慢,等10来分钟后,就可以得到可执行模块。

然后就是测试cp2k/tests/QS/regtest-hybrid目录里面的例子,比如CH3-hybrid-HSE06-lsd.inp,看能否正常运行。对于hybrid functional的计算,通常需要很大的内存。

这篇关于采用intel编译器并行编译执行杂化函数计算CP2K的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/841391

相关文章

使用C#代码计算数学表达式实例

《使用C#代码计算数学表达式实例》这段文字主要讲述了如何使用C#语言来计算数学表达式,该程序通过使用Dictionary保存变量,定义了运算符优先级,并实现了EvaluateExpression方法来... 目录C#代码计算数学表达式该方法很长,因此我将分段描述下面的代码片段显示了下一步以下代码显示该方法如

Oracle的to_date()函数详解

《Oracle的to_date()函数详解》Oracle的to_date()函数用于日期格式转换,需要注意Oracle中不区分大小写的MM和mm格式代码,应使用mi代替分钟,此外,Oracle还支持毫... 目录oracle的to_date()函数一.在使用Oracle的to_date函数来做日期转换二.日

如何用Java结合经纬度位置计算目标点的日出日落时间详解

《如何用Java结合经纬度位置计算目标点的日出日落时间详解》这篇文章主详细讲解了如何基于目标点的经纬度计算日出日落时间,提供了在线API和Java库两种计算方法,并通过实际案例展示了其应用,需要的朋友... 目录前言一、应用示例1、天安门升旗时间2、湖南省日出日落信息二、Java日出日落计算1、在线API2

如何使用 Bash 脚本中的time命令来统计命令执行时间(中英双语)

《如何使用Bash脚本中的time命令来统计命令执行时间(中英双语)》本文介绍了如何在Bash脚本中使用`time`命令来测量命令执行时间,包括`real`、`user`和`sys`三个时间指标,... 使用 Bash 脚本中的 time 命令来统计命令执行时间在日常的开发和运维过程中,性能监控和优化是不

C#如何优雅地取消进程的执行之Cancellation详解

《C#如何优雅地取消进程的执行之Cancellation详解》本文介绍了.NET框架中的取消协作模型,包括CancellationToken的使用、取消请求的发送和接收、以及如何处理取消事件... 目录概述与取消线程相关的类型代码举例操作取消vs对象取消监听并响应取消请求轮询监听通过回调注册进行监听使用Wa

PHP执行php.exe -v命令报错的解决方案

《PHP执行php.exe-v命令报错的解决方案》:本文主要介绍PHP执行php.exe-v命令报错的解决方案,文中通过图文讲解的非常详细,对大家的学习或工作有一定的帮助,需要的朋友可以参考下... 目录执行phpandroid.exe -v命令报错解决方案执行php.exe -v命令报错-PHP War

Oracle数据库执行计划的查看与分析技巧

《Oracle数据库执行计划的查看与分析技巧》在Oracle数据库中,执行计划能够帮助我们深入了解SQL语句在数据库内部的执行细节,进而优化查询性能、提升系统效率,执行计划是Oracle数据库优化器为... 目录一、什么是执行计划二、查看执行计划的方法(一)使用 EXPLAIN PLAN 命令(二)通过 S

C++11的函数包装器std::function使用示例

《C++11的函数包装器std::function使用示例》C++11引入的std::function是最常用的函数包装器,它可以存储任何可调用对象并提供统一的调用接口,以下是关于函数包装器的详细讲解... 目录一、std::function 的基本用法1. 基本语法二、如何使用 std::function

hdu1171(母函数或多重背包)

题意:把物品分成两份,使得价值最接近 可以用背包,或者是母函数来解,母函数(1 + x^v+x^2v+.....+x^num*v)(1 + x^v+x^2v+.....+x^num*v)(1 + x^v+x^2v+.....+x^num*v) 其中指数为价值,每一项的数目为(该物品数+1)个 代码如下: #include<iostream>#include<algorithm>

poj 1113 凸包+简单几何计算

题意: 给N个平面上的点,现在要在离点外L米处建城墙,使得城墙把所有点都包含进去且城墙的长度最短。 解析: 韬哥出的某次训练赛上A出的第一道计算几何,算是大水题吧。 用convexhull算法把凸包求出来,然后加加减减就A了。 计算见下图: 好久没玩画图了啊好开心。 代码: #include <iostream>#include <cstdio>#inclu