24小时热门版块排行榜    

查看: 43617  |  回复: 864
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

小木虫于

金虫 (正式写手)


[交流] vasp 6.4.3 运行出错 LAPACK: Routine ZPOTRF failed! kpoint: 1 spin: 1

今天编译了vasp6.4.3, epyc 9654 cpu,ubuntu 24.04 server系统。 编译好后,一运行就出现这个错误:
Orbital orthonormalization failed in the inversion of matrix            |
|     LAPACK: Routine ZPOTRF failed! kpoint: 1 spin: 1  

不知道是否有人遇到这个问题,怎么解决的?不吝赐教,非常感谢

信息:结构没有任何问题,用vasp 5.4 运行一切正常,vasp6.4.3就不行
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小木虫于

金虫 (正式写手)


引用回帖:
17楼: Originally posted by rlafite at 2025-11-17 12:18:05
EPYC 9654 是 Zen4,大部分人遇到这个错误都是因为:

1. 链接到系统自带的 OpenBLAS/BLIS → 数值不稳定

Ubuntu 24.04 默认 LAPACK/OpenBLAS 非常不适合 VASP 6(尤其是大量 BLAS3 操作)。

解决:避免使用 ...

我是用spack安装的aocc和aocl,vasp的makefile.include 里面是直接连接到了spack安装的。makefile.include如下:
# Default precompiler options
CPP_OPTIONS = -DHOST=\"LinuxGNU\" \
              -DMPI -DMPI_BLOCK=8000 -Duse_collective \
              -DscaLAPACK \
              -DCACHE_SIZE=4000 \
              -Davoidalloc \
              -Dvasp6 \
              -Duse_bse_te \
              -Dtbdyn \
              -Dfock_dblbuf

CPP         = flang -E -ffree-form -C -w $*$(FUFFIX) >$*$(SUFFIX) $(CPP_OPTIONS)

FC          = mpif90
FCL         = mpif90

FREE        = -ffree-form -ffree-line-length-none
FFLAGS      = -w -fno-fortran-main -Mbackslash

OFLAG       = -O2
OFLAG_IN    = $(OFLAG)
DEBUG       = -O0

OBJECTS     = fftmpiw.o fftmpi_map.o fftw3d.o fft3dlib.o
OBJECTS_O1 += fftw3d.o fftmpi.o fftmpiw.o
OBJECTS_O2 += fft3dlib.o

# For what used to be vasp.5.lib
CPP_LIB     = $(CPP)
FC_LIB      = $(FC)
CC_LIB      = clang
CFLAGS_LIB  = -O
FFLAGS_LIB  = -O1
FREE_LIB    = $(FREE)

OBJECTS_LIB = linpack_double.o

# For the parser library
CXX_PARS    = clang++
LLIBS       = -lstdc++

##
## Customize as of this point! Of course you may change the preceding
## part of this file as well if you like, but it should rarely be
## necessary ...

# relevant target when cross-compiling for another architecture
VASP_TARGET_CPU ?= -march=native
FFLAGS     += $(VASP_TARGET_CPU)

# BLAS (mandatory)
AMDBLIS_ROOT ?= /public/spack/opt/spack/linux-zen4/amdblis-5.0-hi5do7p2j5yk3xxo5eaff76d57ozinlr
BLAS        = -L${AMDBLIS_ROOT}/lib -lblis

# LAPACK (mandatory)
AMDLIBFLAME_ROOT ?= /public/spack/opt/spack/linux-zen4/amdlibflame-5.0-72mkvewha6igrfhobf7uceemuy3xu4dp
LAPACK      = -L${AMDLIBFLAME_ROOT}/lib -lflame

# scaLAPACK (mandatory)
AMDSCALAPACK_ROOT ?= /public/spack/opt/spack/linux-zen4/amdscalapack-5.0-2ok2w3vt4dpnmouvy5ezfksheypcxg3b
SCALAPACK   = -L${AMDSCALAPACK_ROOT}/lib -lscalapack

LLIBS      += $(SCALAPACK) $(LAPACK) $(BLAS)

# FFTW (mandatory)
AMDFFTW_ROOT  ?= /public/spack/opt/spack/linux-zen4/amdfftw-5.0-fuzt5flfdtrmytkknf6dgpcomzpnbpiq
LLIBS      += -L$(AMDFFTW_ROOT)/lib -lfftw3
INCS       += -I$(AMDFFTW_ROOT)/include

# HDF5-support (optional but strongly recommended)
#CPP_OPTIONS+= -DVASP_HDF5
#HDF5_ROOT  ?= /path/to/your/hdf5/installation
#LLIBS      += -L$(HDF5_ROOT)/lib -lhdf5_fortran
#INCS       += -I$(HDF5_ROOT)/include

# For the VASP-2-Wannier90 interface (optional)
CPP_OPTIONS    += -DVASP2WANNIER90
WANNIER90_ROOT ?= /public/spack/opt/spack/linux-zen4/wannier90-3.1.0-qvmhejtugcrv7cz4xenk2y235lld2hrw
LLIBS          += -L$(WANNIER90_ROOT)/lib -lwannier
18楼2025-11-17 14:39:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 865 个回答

rlafite

至尊木虫 (正式写手)


★
小木虫于(金币+1): 谢谢参与
Install VASP 6.5.1
9楼2025-11-15 23:35:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小木虫于

金虫 (正式写手)


引用回帖:
9楼: Originally posted by rlafite at 2025-11-15 23:35:25
Install VASP 6.5.1

是vasp6.4.3有bug么?
14楼2025-11-16 22:27:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
简单回复
toppaper3楼
2025-11-14 22:05   回复  
小木虫于(金币+1): 谢谢参与
发自小木虫手机客户端
jiaoxg4楼
2025-11-15 09:25   回复  
小木虫于(金币+1): 谢谢参与
tfang6楼
2025-11-15 12:57   回复  
小木虫于(金币+1): 谢谢参与
发自小木虫手机客户端
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[公派出国] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +4 G2rhSvqvhNzl 2026-11-10 5/250 2026-11-13 11:37 by Equinoxhua
[考研] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +5 C6CRn66RYYXS 2026-11-10 6/300 2026-11-13 11:04 by Equinoxhua
[考博] 27博士申请(市政/环境) +3 ladf 2026-11-12 3/150 2026-11-13 08:19 by Equinoxhua
[公派出国] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 TgvkbZZl9tqB 2026-11-11 4/200 2026-11-12 12:51 by ECLTXnJnRSXf
[公派出国] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 TgvkbZZl9tqB 2026-11-11 4/200 2026-11-12 12:48 by ECLTXnJnRSXf
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 TgvkbZZl9tqB 2026-11-11 6/300 2026-11-12 12:42 by ECLTXnJnRSXf
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 Z7qfCsdIiD1T 2026-11-10 7/350 2026-11-12 11:31 by HajghKOMynmz
[考研] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +5 Z7qfCsdIiD1T 2026-11-10 8/400 2026-11-12 11:22 by HajghKOMynmz
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +5 czAiolzR0IJC 2026-11-10 7/350 2026-11-12 11:09 by HajghKOMynmz
[基金申请] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 czAiolzR0IJC 2026-11-10 8/400 2026-11-12 11:00 by HajghKOMynmz
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 G2rhSvqvhNzl 2026-11-10 6/300 2026-11-12 10:38 by HajghKOMynmz
[硕博家园] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 C6CRn66RYYXS 2026-11-10 6/300 2026-11-12 10:16 by HajghKOMynmz
[基金申请] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +6 zSbvVX9JFl9J 2026-11-09 9/450 2026-11-12 09:21 by HajghKOMynmz
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 C6CRn66RYYXS 2026-11-10 5/250 2026-11-12 06:03 by 9fOTqj074tIg
[教师之家] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +5 FChKX8meu1ds 2026-11-09 10/500 2026-11-12 05:19 by 9fOTqj074tIg
[教师之家] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +6 0TEMgxdT1eIw 2026-11-09 8/400 2026-11-12 05:00 by 9fOTqj074tIg
[硕博家园] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +6 0TEMgxdT1eIw 2026-11-09 9/450 2026-11-12 04:57 by 9fOTqj074tIg
[考博] 申请环境,微生物方向博士 +3 浅笑09 2026-11-08 3/150 2026-11-10 12:44 by 可淡不可忘
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +5 5sJFDKp2kwuI 2026-11-06 16/800 2026-11-09 11:47 by GwrQVVR6Oenx
[考博] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +5 5sJFDKp2kwuI 2026-11-06 15/750 2026-11-09 11:41 by GwrQVVR6Oenx
信息提示
请填处理意见