免费国产精品自在自线-91精品国产色综合久久久浪潮-99热久久免费频精品-国产精品国模在线观看-久久亚洲国产精品成人?V秋霞-久久国产一级A片免费播放-亚洲国产欧洲综合97久久-久久国产白嫩美女呻吟高潮

ARTICLE DETAIL

資訊詳情

深耕商務(wù)建站與企業(yè)官網(wǎng)運(yùn)營(yíng)的一線實(shí)戰(zhàn)洞察。

Ubuntu 20.04 下使用 Intel oneAPI 編譯安裝 ALAMODE 的完整指南

Ubuntu 20.04 下使用 Intel oneAPI 編譯安裝 ALAMODE 的完整指南 ALAMODE 是基于第一性原理計(jì)算晶格動(dòng)力學(xué)和熱導(dǎo)率的開源工具。很多做聲子、熱導(dǎo)率、非諧效應(yīng)的材料計(jì)算研究者都會(huì)在 VASP 之外配套一個(gè)力常數(shù)和聲子計(jì)算程序ALAMODE 是這類工具里比較有針對(duì)性的一套。這篇文章要解決的是在 Ubuntu 20.04 上使用 Intel oneAPI 編譯器從零編譯安裝 ALAMODE并完成環(huán)境驗(yàn)證和基礎(chǔ)測(cè)試。先說結(jié)論ALAMODE 本身不依賴 GPU用 CPU 編譯和運(yùn)行。難點(diǎn)集中在依賴庫(kù)的拼裝上尤其是 FFTW、LAPACK 和 Eigen 三個(gè)底層庫(kù)。只要把 Intel 編譯器環(huán)境、MKL 數(shù)學(xué)庫(kù)和 Eigen 頭文件準(zhǔn)備到位剩下的configure make make install流程不算復(fù)雜。整個(gè)安裝過程適合在本地工作站或小型計(jì)算節(jié)點(diǎn)上完成通常不會(huì)出現(xiàn)“缺依賴導(dǎo)致無法繼續(xù)”的致命問題最多是在路徑配置上多花一點(diǎn)時(shí)間。本文會(huì)從 ALAMODE 的核心能力、適用場(chǎng)景、環(huán)境準(zhǔn)備、依賴庫(kù)編譯、源碼配置、功能測(cè)試、資源占用和常見問題幾個(gè)方面展開。每一節(jié)都會(huì)給出可復(fù)制的命令和配置示例方便你直接照著操作。如果你已經(jīng)裝過 Phonopy那么這篇文章的很多思路是相通的只是 ALAMODE 在非諧力常數(shù)和導(dǎo)熱率計(jì)算上更進(jìn)一步。1. ALAMODE 核心能力速覽能力項(xiàng)說明項(xiàng)目類型原子級(jí)材料計(jì)算開源軟件主要用于晶格動(dòng)力學(xué)與熱導(dǎo)率模擬開發(fā)語(yǔ)言Fortran 與 C/C 混合主要功能諧波/非諧力常數(shù)擬合、聲子色散、聲子態(tài)密度、聲子壽命、晶格熱導(dǎo)率計(jì)算硬件要求不需要 GPU純 CPU 計(jì)算內(nèi)存大小與超胞規(guī)模、原子數(shù)和振動(dòng)模式數(shù)相關(guān)支持操作系統(tǒng)Linux、macOS本文以 Ubuntu 20.04 為例編譯器支持Intel ifort/icc、GCC/gfortran本文采用 Intel oneAPI 編譯器關(guān)鍵依賴Intel MKL或 LAPACK/BLAS、FFTW、Eigen、MPI可選啟動(dòng)方式命令行工具安裝后生成alm和anphon兩個(gè)可執(zhí)行文件是否支持批量任務(wù)支持可通過 shell/Python 腳本批量處理多個(gè)輸入結(jié)構(gòu)也可以借助 MPI 并行適合場(chǎng)景基于 VASP 等第一性原理計(jì)算的聲子計(jì)算和晶格熱導(dǎo)率研究從表中可以看到ALAMODE 不屬于那種“下載即用”的軟件。它需要你提前準(zhǔn)備好一套科學(xué)計(jì)算編譯鏈。對(duì) Ubuntu 20.04 用戶來說最省心的方式是安裝 Intel oneAPI Base Toolkit 和 HPC Toolkit里面已經(jīng)包含了 ifort/icx 編譯器、Intel MPI 和 MKL。這樣 LAPACK 和 BLAS 可以直接從 MKL 獲取FFTW 也可以用 MKL 自帶的兼容接口只有 Eigen 需要單獨(dú)下載。2. 適用場(chǎng)景與使用邊界ALAMODE 主要適合做凝聚態(tài)物理和材料物理中與聲子相關(guān)的研究。典型使用場(chǎng)景包括從第一性原理計(jì)算得到的原子受力和能量出發(fā)擬合諧波力常數(shù)計(jì)算聲子色散和聲子態(tài)密度計(jì)算三階非諧力常數(shù)進(jìn)一步分析聲子-聲子散射、聲子壽命和模式群速度結(jié)合聲子玻爾茲曼輸運(yùn)方程計(jì)算材料的晶格熱導(dǎo)率模擬溫度變化對(duì)聲子頻率和色散關(guān)系的影響用于高溫相穩(wěn)定性分析。這些功能意味著 ALAMODE 通常需要和 VASP、QE 等第一性原理軟件配合使用。你需要先準(zhǔn)備好 DFPT 或有限位移法計(jì)算所需的超胞、位移構(gòu)型以及對(duì)應(yīng)的原子受力和總能數(shù)據(jù)。ALAMODE 本身不執(zhí)行電子結(jié)構(gòu)計(jì)算只負(fù)責(zé)從這些數(shù)據(jù)中提取力常數(shù)并做聲子相關(guān)的后處理。從邊界上看ALAMODE 不是萬(wàn)能的。它不擅長(zhǎng)處理強(qiáng)關(guān)聯(lián)電子體系或者包含明顯非諧局域模式的復(fù)雜體系。對(duì)于大超胞、高對(duì)稱性和更多位移構(gòu)型的任務(wù)計(jì)算量會(huì)快速增長(zhǎng)。此外ALAMODE 的編譯安裝對(duì) Fortran 編譯器版本和數(shù)學(xué)庫(kù)兼容性有一定要求如果你使用的編譯器是老版本 GCC 或者系統(tǒng)自帶的基礎(chǔ)庫(kù)在鏈接階段可能會(huì)遇到一些難以排查的錯(cuò)誤。使用 ALAMODE 處理實(shí)驗(yàn)或文獻(xiàn)數(shù)據(jù)時(shí)還要注意軟件許可和數(shù)據(jù)合規(guī)問題。ALAMODE 本身是開源許可證可以自由使用和修改但如果你用 VASP 生成輸入數(shù)據(jù)必須確保 VASP 的授權(quán)覆蓋你的研究場(chǎng)景。涉及未發(fā)表的結(jié)構(gòu)數(shù)據(jù)、合作方材料數(shù)據(jù)或商業(yè)項(xiàng)目時(shí)也要確認(rèn)數(shù)據(jù)使用邊界避免出現(xiàn)版權(quán)或保密問題。3. 環(huán)境準(zhǔn)備與前置條件在 Ubuntu 20.04 上編譯 ALAMODE建議先在干凈的系統(tǒng)環(huán)境下做一遍完整性檢查避免依賴沖突。下面是最小前置條件Ubuntu 20.04 操作系統(tǒng)內(nèi)核更新到最新補(bǔ)丁至少 8 GB 內(nèi)存磁盤剩余空間 5 GB 以上GCC、gfortran 和 make 已安裝Intel oneAPI Base Toolkit 和 HPC Toolkit 已安裝網(wǎng)絡(luò)可訪問 ALAMODE 源碼倉(cāng)庫(kù)和 Eigen 官網(wǎng)。如果還沒有安裝 Intel oneAPI可以通過以下命令檢查環(huán)境source /opt/intel/oneapi/setvars.sh which ifort which icc如果沒有找到ifort說明 HPC Toolkit 沒有安裝到位。需要注意的是Intel 已經(jīng)用ifx逐步替代ifort但 ALAMODE 的編譯腳本目前對(duì)ifort的支持更成熟。如果系統(tǒng)里只有ifx可以嘗試將FCifx傳入 configure但更穩(wěn)妥的路徑是安裝完整 HPC Toolkit它通常同時(shí)包含ifort和icx。系統(tǒng)基礎(chǔ)編譯工具也建議提前裝好sudo apt update sudo apt install -y build-essential gfortran wget git cmake接下來檢查是否已經(jīng)有mkl環(huán)境變量。通常安裝 Intel oneAPI 后setvars.sh會(huì)設(shè)置MKLROOTecho $MKLROOT如果輸出為空確認(rèn)setvars.sh是否正確 source。到這里基礎(chǔ)環(huán)境已經(jīng)準(zhǔn)備好下一步開始準(zhǔn)備 ALAMODE 的依賴庫(kù)。4. 依賴庫(kù)編譯準(zhǔn)備4.1 Intel MKL 的 LAPACK 支持ALAMODE 的alm和anphon核心功能需要大量線性代數(shù)運(yùn)算例如特征值分解、矩陣求逆、最小二乘擬合等。Intel MKL 提供了完整的 LAPACK 和 BLAS 實(shí)現(xiàn)性能上比 OpenBLAS 和系統(tǒng)自帶的 LAPACK 更好也能和 Intel 編譯器無縫配合。在 configure 階段可以使用--with-lapack$MKLROOT或者通過 LDFLAGS 顯式指定鏈接庫(kù)。MKL 的典型鏈接參數(shù)如下export LDFLAGS-L$MKLROOT/lib/intel64 -lmkl_intel_lp64 -lmkl_sequential -lmkl_core如果你的編譯器是ifort還需要確保lmkl_blacs_intelmpi_lp64等 MPI 相關(guān)庫(kù)存在。只做串行編譯時(shí)上面三件套已經(jīng)夠用。4.2 FFTW 的獲取方式ALAMODE 依賴 FFTW 進(jìn)行快速傅里葉變換這主要是為了在倒空間插值和聲子計(jì)算時(shí)提高效率。FFTW 有兩個(gè)常見來源獨(dú)立安裝的 FFTW或者使用 Intel MKL 自帶的 FFTW 兼容接口。獨(dú)立安裝 FFTW 的方式更通用尤其當(dāng)你可能還會(huì)把 ALAMODE 和 GCC 編譯器一起使用時(shí)。FFTW 3.x 的編譯很簡(jiǎn)單wget http://www.fftw.org/fftw-3.3.10.tar.gz tar -zxvf fftw-3.3.10.tar.gz cd fftw-3.3.10 ./configure --prefix$HOME/fftw --enable-shared --enable-single make -j$(nproc) make install這里--enable-single表示單精度 FFTW。ALAMODE 通常需要雙精度 FFTW如果你不確定可以去掉這個(gè)選項(xiàng)默認(rèn)就是雙精度。安裝完成后include目錄下會(huì)有fftw3.hlib目錄下會(huì)有l(wèi)ibfftw3.a或.so。如果不想額外編譯 FFTW也可以使用 MKL 的 FFTW 接口。在 configure 時(shí)可能需要配置--with-fftw$MKLROOT并確保頭文件路徑和庫(kù)文件路徑正確。不同版本的 MKL 對(duì)這個(gè)接口的封裝略有差異如果遇到找不到fftw3.h的問題建議直接安裝獨(dú)立 FFTW省時(shí)省力。4.3 Eigen 頭文件獲取Eigen 是一個(gè)純頭文件的 C 模板庫(kù)不需要編譯安裝過程只是解壓并放置到合適位置。ALAMODE 在擬合力常數(shù)時(shí)用到了 Eigen 的線性代數(shù)模板。從 Eigen 官網(wǎng)或 GitHub 下載最新穩(wěn)定版例如 3.4.0cd $HOME wget https://gitlab.com/libeigen/eigen/-/archive/3.4.0/eigen-3.4.0.tar.gz tar -zxvf eigen-3.4.0.tar.gz mv eigen-3.4.0 eigen這樣$HOME/eigen目錄下就能看到Eigen和unsupported兩個(gè)子目錄。configure 時(shí)可用--with-eigen$HOME/eigen指定路徑。Eigen 只需頭文件不參與編譯所以后面 ALAMODE 編譯時(shí)不會(huì)產(chǎn)生額外的.o文件。4.4 可選的 MPI 環(huán)境如果你的計(jì)算任務(wù)需要處理大超胞或多位移構(gòu)型建議安裝 Intel MPI。Intel oneAPI HPC Toolkit 通常會(huì)捆綁 Intel MPI安裝后運(yùn)行source /opt/intel/oneapi/setvars.sh which mpiifort如果函數(shù)正常返回說明 MPI 環(huán)境可用。ALAMODE 的 configure 腳本會(huì)自動(dòng)檢測(cè) MPI但如果沒有檢測(cè)到也不會(huì)影響串行編譯。你可以在 configure 參數(shù)中顯式禁用 MPI--without-mpi不過多數(shù)材料計(jì)算場(chǎng)景還是建議保留 MPI 支持因?yàn)闊釋?dǎo)率計(jì)算涉及大量聲子模式并行化能顯著縮短等待時(shí)間。5. ALAMODE 源碼下載與編譯配置5.1 下載源碼ALAMODE 的源碼托管在 GitHub 或官方網(wǎng)站上。建議下載最新 release 版本而不是直接 clone master 分支這樣能保證穩(wěn)定性和可復(fù)現(xiàn)性。以 tar.gz 包為例cd $HOME wget https://github.com/ttadano/alamode/releases/download/v1.4.0/alamode-v1.4.0.tar.gz tar -zxvf alamode-v1.4.0.tar.gz cd alamode-v1.4.0這里的版本號(hào)只是一個(gè)示例實(shí)際下載時(shí)請(qǐng)以官方倉(cāng)庫(kù)的 release 為準(zhǔn)。如果 GitHub 訪問速度慢也可以通過官方主頁(yè)下載。解壓后目錄下會(huì)有configure、Makefile.in、examples、src等文件。5.2 設(shè)置 Intel 編譯環(huán)境進(jìn)入源碼目錄后先 source Intel oneAPI 環(huán)境再顯式設(shè)置編譯器變量source /opt/intel/oneapi/setvars.sh export CCicc export CXXicpc export FCifort export F77ifort export FCFLAGS-O2 -xHost export CFLAGS-O2 -xHost如果你的 oneAPI 版本較新icx已經(jīng)替代iccifx也已經(jīng)出現(xiàn)。此時(shí)可以嘗試export CCicx export CXXicpx export FCifort但請(qǐng)注意ALAMODE 的某些 Fortran 代碼可能與ifx存在兼容性差異如果編譯報(bào)錯(cuò)回退到ifort是更穩(wěn)妥的方案。-xHost表示針對(duì)本機(jī) CPU 架構(gòu)優(yōu)化如果后續(xù)需要遷移到其他機(jī)器建議換成-marchcore-avx2等通用指令集。5.3 configure 生成 Makefile配置命令要根據(jù)上一節(jié)準(zhǔn)備的依賴路徑來寫。下面的示例假設(shè) FFTW 安裝到$HOME/fftwEigen 在$HOME/eigenMKL 使用系統(tǒng)默認(rèn)環(huán)境./configure --prefix$HOME/alamode \ --with-fftw$HOME/fftw \ --with-lapack$MKLROOT \ --with-eigen$HOME/eigen \ --with-mpiyes--with-mpiyes表示啟用 MPI。如果你沒有安裝 Intel MPI可以用--with-mpino關(guān)閉。configure 腳本會(huì)輸出檢查到的編譯器、庫(kù)路徑和功能狀態(tài)。關(guān)鍵看是否有checking for FFTW... yes、checking for LAPACK... yes、checking for Eigen... yes之類的結(jié)果。如果某個(gè)依賴顯示no需要回到上一步檢查路徑和頭文件。5.4 編譯與安裝configure 沒問題后直接執(zhí)行編譯make -j$(nproc) make installmake install會(huì)把可執(zhí)行文件安裝到--prefix指定的目錄下默認(rèn)產(chǎn)生$HOME/alamode/bin目錄其中包含alm和anphon兩個(gè)可執(zhí)行程序。安裝完成后將安裝目錄加入PATHexport PATH$HOME/alamode/bin:$PATH為了讓每次登錄自動(dòng)生效可以寫入~/.bashrc。這一步之后ALAMODE 就算安裝完成了。6. 功能測(cè)試與效果驗(yàn)證安裝完成后不要急著跑大體系先用程序自帶的幫助信息和示例數(shù)據(jù)做一次完整性驗(yàn)證。6.1 檢查可執(zhí)行文件which alm which anphon如果能正常返回路徑說明bin目錄已經(jīng)生效。再執(zhí)行alm --help anphon --help如果輸出幫助信息沒有出現(xiàn)error while loading shared libraries說明動(dòng)態(tài)庫(kù)鏈接正常。6.2 使用自帶示例驗(yàn)證流程ALAMODE 源碼包通常帶有examples目錄。比如常見的硅模型、金剛石模型等。你可以進(jìn)入某個(gè)示例目錄查看輸入文件格式然后嘗試運(yùn)行cd examples/si alm這里不一定會(huì)輸出一個(gè)完整的計(jì)算結(jié)果因?yàn)椴煌姹臼纠妮斎胛募煌?。更通用的?yàn)證方式是用mkalmt或anphon直接執(zhí)行一個(gè)簡(jiǎn)單任務(wù)。但為了不依賴特定示例文件你可以在自己已經(jīng)準(zhǔn)備好的第一性原理數(shù)據(jù)上測(cè)試。初次使用時(shí)建議選擇一個(gè)小超胞例如 2x2x2 的硅晶體這樣幾分鐘內(nèi)就能跑完便于快速確認(rèn)程序和依賴庫(kù)是否正常工作。6.3 判斷編譯成功的關(guān)鍵指標(biāo)編譯成功的判斷標(biāo)準(zhǔn)有幾個(gè)alm和anphon可以正常執(zhí)行沒有段錯(cuò)誤程序可以讀取你的力常數(shù)或位移輸入文件并給出力常數(shù)擬合報(bào)告聲子色散計(jì)算能輸出頻率數(shù)據(jù)數(shù)值在半導(dǎo)體的合理范圍比如硅的光學(xué)聲子頻率接近 16 THz日志文件末尾沒有出現(xiàn)NaN或Inf數(shù)值。如果輸出中有大量NaN通常不是編譯問題而是輸入數(shù)據(jù)質(zhì)量問題可能是位移構(gòu)型不完整或受力數(shù)據(jù)精度不夠。這種情況需要回到第一性原理計(jì)算階段檢查 VASP 設(shè)置。6.4 MPI 并行測(cè)試如果你的 ALAMODE 啟用了 MPI可以用 MPI 方式運(yùn)行示例mpirun -np 4 anphon觀察程序是否正常啟動(dòng)并利用多核計(jì)算。如果這步報(bào)錯(cuò)比如找不到libmpi.so說明 MPI 庫(kù)鏈接不完整需要檢查L(zhǎng)D_LIBRARY_PATH是否包含 Intel MPI 的 lib 目錄。7. 資源占用與性能觀察ALAMODE 是純 CPU 計(jì)算程序運(yùn)行時(shí)主要關(guān)注 CPU 使用率、內(nèi)存占用和并行效率。下面給出一些觀察方法和調(diào)優(yōu)思路。在終端用htop或top查看進(jìn)程信息。串行運(yùn)行時(shí)CPU 使用率應(yīng)該在一個(gè)核上接近 100%。并行運(yùn)行時(shí)多個(gè)核的 CPU 使用率都會(huì)升高。如果并行之后 CPU 使用率仍然只有 100%說明 MPI 沒有真正生效需要檢查命令行是否使用了mpirun以及程序啟動(dòng)時(shí)有沒有加載 MPI 庫(kù)。內(nèi)存方面ALAMODE 的內(nèi)存需求主要集中在力常數(shù)矩陣求解和聲子模式對(duì)角化。超胞原子數(shù)越多矩陣規(guī)模越大內(nèi)存占用越高。以常見 2x2x2 硅超胞為例內(nèi)存占用通常不超過 2 GB。如果你的超胞超過 100 個(gè)原子建議內(nèi)存至少 16 GB 以上。要臨時(shí)降低內(nèi)存峰值可以關(guān)閉并行環(huán)境變量中的超線程或者使用更緊湊的矩陣存儲(chǔ)模式。性能調(diào)優(yōu)可以參考以下幾點(diǎn)使用OMP_NUM_THREADS控制 OpenMP 線程數(shù)配合 MPI 做混合并行編譯時(shí)使用-O2或-O3優(yōu)化但不要盲目使用最高優(yōu)化級(jí)別某些優(yōu)化可能導(dǎo)致數(shù)值誤差使用 MKL 提供的多線程版本并設(shè)置MKL_NUM_THREADS參數(shù)在超算或服務(wù)器上盡量將計(jì)算任務(wù)綁定到物理核心避免 CPU 頻繁切換。8. 常見問題與排查方法在 Ubuntu 20.04 下編譯 ALAMODE最容易遇到的問題集中在編譯器環(huán)境、數(shù)學(xué)庫(kù)路徑和鏈接參數(shù)上。下表整理了常見現(xiàn)象、可能原因和解決方案。問題現(xiàn)象可能原因排查方式解決方案configure 報(bào)No Fortran compiler found未安裝 ifort或環(huán)境變量未設(shè)置執(zhí)行which ifort安裝 Intel oneAPI HPC Toolkit并source setvars.sh找不到fftw3.hFFTW 路徑配置錯(cuò)誤或未安裝 FFTW查看$HOME/fftw/include/fftw3.h是否存在重裝 FFTW或重新指定--with-fftw鏈接時(shí)報(bào)cannot find -lmkl_sequentialMKL 庫(kù)路徑未添加到 LDFLAGS檢查$MKLROOT/lib/intel64是否有對(duì)應(yīng)文件在 configure 前設(shè)置 LDFLAGSEigen 相關(guān)編譯報(bào)錯(cuò)Eigen 頭文件路徑錯(cuò)誤或版本過舊檢查Eigen/Core是否可訪問下載新版 Eigen并通過--with-eigen指定路徑alm啟動(dòng)時(shí)報(bào)libiomp5.so找不到Intel OpenMP 運(yùn)行庫(kù)未加載執(zhí)行l(wèi)dd alm/ldd anphon重新 source oneAPI 環(huán)境并設(shè)置LD_LIBRARY_PATH編譯過程中internal compiler error編譯器版本過舊或優(yōu)化級(jí)別過高查看編譯器版本重試降低優(yōu)化級(jí)別更新 Intel oneAPI 到最新版或改用-O1編譯MPI 并行程序報(bào)錯(cuò)MPI_Init失敗MPI 環(huán)境未正確加載或網(wǎng)絡(luò)接口不匹配運(yùn)行mpirun --version重新source setvars.sh或指定--mca參數(shù)計(jì)算結(jié)果全部為 NaN輸入數(shù)據(jù)或力常數(shù)擬合異常不一定是編譯問題檢查輸入文件格式和第一性原理輸出的受力精度增加 VASP 受力收斂標(biāo)準(zhǔn)重新生成輸入數(shù)據(jù)如果 configure 階段碰到LAPACK library not found可以手動(dòng)用 MKL 鏈接一個(gè)簡(jiǎn)單的 Fortran 測(cè)試程序來驗(yàn)證。比如寫一個(gè)只調(diào)用zheev的小程序用ifort配合 LDFLAGS 編譯能跑通就說明 MKL 環(huán)境正常。這樣能把編譯問題快速定位到 ALAMODE 本身還是系統(tǒng)數(shù)學(xué)庫(kù)上。9. 最佳實(shí)踐與使用建議編譯安裝完成后建議把一次完整的安裝流程記錄下來包括依賴版本、configure 參數(shù)、環(huán)境變量形成一份安裝筆記。后續(xù)換機(jī)器或者給同事復(fù)現(xiàn)時(shí)可以節(jié)省大量排查時(shí)間。在項(xiàng)目管理上推薦把源碼、依賴庫(kù)和計(jì)算結(jié)果分目錄存放。比如$HOME/ ├── alamode_src/ # ALAMODE 源碼 ├── fftw/ # 獨(dú)立編譯的 FFTW ├── eigen/ # Eigen 頭文件 ├── alamode_install/ # ALAMODE 安裝目錄 └── calc/ # 實(shí)際計(jì)算任務(wù)每個(gè)計(jì)算任務(wù)目錄下單獨(dú)放輸入文件、腳本和輸出結(jié)果避免多個(gè)項(xiàng)目混用同一個(gè)工作目錄。這樣即便要重新編譯 ALAMODE也不會(huì)影響正在跑的計(jì)算任務(wù)。使用 ALAMODE 時(shí)建議先小規(guī)模測(cè)試再上大規(guī)模任務(wù)。第一性原理計(jì)算和聲子計(jì)算都是重資源任務(wù)如果直接在幾百個(gè)原子的超胞上運(yùn)行后期可能會(huì)因?yàn)檩斎霐?shù)據(jù)問題浪費(fèi)大量機(jī)時(shí)。先跑一個(gè)小體系確認(rèn)程序、流程和結(jié)果合理性再擴(kuò)大范圍。對(duì)于需要長(zhǎng)期運(yùn)行的批量任務(wù)建議寫一個(gè)簡(jiǎn)單的 shell/python 包裝腳本自動(dòng)遍歷不同結(jié)構(gòu)、不同溫度或不同位移幅度并記錄每次運(yùn)行的日志。批量任務(wù)要設(shè)置超時(shí)和失敗重試機(jī)制避免某個(gè)計(jì)算任務(wù)卡住影響整個(gè)隊(duì)列。輸出文件夾也要包含清晰的命名比如結(jié)構(gòu)名、溫度、位移大小等。涉及第三方材料數(shù)據(jù)、未發(fā)表結(jié)構(gòu)或商業(yè)合作項(xiàng)目時(shí)務(wù)必確認(rèn)數(shù)據(jù)授權(quán)范圍。ALAMODE 的開源許可只覆蓋軟件本身不覆蓋你輸入的原子結(jié)構(gòu)和第一性原理計(jì)算結(jié)果。如果使用 VASP 等商業(yè)軟件生成數(shù)據(jù)還需要遵守對(duì)應(yīng)軟件的許可條款。10. 總結(jié)與下一步ALAMODE 在 Ubuntu 20.04 上的 Intel 版編譯安裝核心路徑是準(zhǔn)備好 Intel oneAPI 編譯環(huán)境編譯或配置好 FFTW、LAPACK/BLAS、Eigen 三個(gè)依賴庫(kù)然后運(yùn)行configure生成 Makefile最后make make install。整個(gè)過程不依賴 GPU主要考驗(yàn)的是編譯環(huán)境和路徑配置。最容易踩的坑有三個(gè)一是沒有 source Intel oneAPI 環(huán)境導(dǎo)致找不到 ifort二是 FFTW 路徑錯(cuò)誤導(dǎo)致 configure 檢測(cè)失敗三是 LAPACK 鏈接參數(shù)不完整導(dǎo)致編譯能通過但鏈接階段報(bào)錯(cuò)。如果你第一次安裝失敗優(yōu)先檢查這三處。安裝成功后先跑一個(gè)最小示例驗(yàn)證alm和anphon是否正常。如果小體系能輸出正確的聲子頻率說明整套流程已經(jīng)打通后面可以開始處理你自己的材料結(jié)構(gòu)。建議接下來的方向是先用簡(jiǎn)單體系做一次聲子色散計(jì)算逐步熟悉 ALAMODE 的輸入文件格式再嘗試計(jì)算三階非諧力常數(shù)和聲子壽命最后再做晶格熱導(dǎo)率分析。如果計(jì)算規(guī)模比較大可以嘗試 MPI 并行并對(duì)比不同并行方案對(duì)效率的影響。
返回列表
PREV
查看更多資訊
NEXT
返回資訊列表
五月丁香六月婷婷久久| 狠狠干综合网| 怡红院院久久| 色婷网站| 大香蕉伊人爱在线| 亚洲色婷婷| 五月婷综合性中心| 狠狠色狠狠干| 超碰操日| 成人AV中文字幕| 久久综合色情网站| 色婷婷成人做爰A片免费看网站 | 日本99在线| 99狠狠| 激情丁香久久| 日本啪啪天堂| 亚洲综合色成丁香五月色| 丁香五月婷在线观看| 婷婷五月天av小说| 综合色色色| 欧美色狠婷久| 亚洲国产另类av| 五月天激情久色| 久久丁香婷婷色情综合| se99高清无码| 丁香五月色激情| 久久久久这里只有精品| 亚洲午夜精品久久久久久人妖| 九九九激情网| 日本色色影片| 色情五月婷婷| 亚洲色热| 亚洲成人黄色网| 玖玖在线| 丁香九色不卡aaa| 女婷久久| 婷婷亚洲欧美丁香五月| 办公室少妇激情呻吟A片在线观看| 五月天色婷婷小说| 婷婷五月综合婷婷| 色情丁香五月婷婷精品| 9l视频自拍九色9l黑人| 97碰操| 六月婷婷激情| 日日夜夜婷婷| 丁香色婷婷| 五月天激情视频五月天| 婷婷五月激情基地| 99色色爰| 91九色熟女| 黄色录像网点| 国产精品涩涩涩视频网站| 狠狠狠狠狠干| 色综合婷婷| 婷婷激情五月天激情| 丁香久久AV| 五月婷婷综合网| 热99一二三| 991自拍视频| 色五月婷婷五月| www.激情| 久久视频婷婷| 中文成人在线| 色婷婷在线视频| 99热综合网| 伊人五月天婷婷| 一夜福利不卡| 99久久久| 9l视频自拍九色9l视频自拍九色9l社区| 深爱激情丁香| 久狠狠| 亚洲爱爱无码婷婷色五月| 亚洲欧洲色色| 色情免费视频播放| 天天日天天干天天插天天射| 天天舔天天摸天天透| 色综合久久88色综合天天人守婷| 九九视频精品这里只有| 五月婷婷色丁香| 青青草原中文字幕| 久久er99| 国产激情在线| 成人片久久网站| 五月婷婷婷婷婷婷艺术| 婷婷五月天激情在线| 91人人爽久久涩噜噜噜| 天天操比比| 婷婷色婷婷| 色五月婷婷综合| 天天综合网站| 九九综合色综合| 国产精品五月丁香| 婷婷伊人中文字幕| 超碰99热在线观看| 很很干五月天| 99热亚洲| 久久婷婷五月天丁香| 激情婷婷视频在线| 午夜色婷婷| 91色五月| 99热这里只有精品5| 亚洲精品视频在线| 超碰婷婷色| 激情五月天。| 99热国产国产| 99久久网站| 专区无日本视频高清8| 婷婷五月天网| 丁香啪啪| 久热播这里只有精品| 丁香五月中文字幕| 99在线免费视频| 五月丁香婷婷啪啪网| 成人丁香色| 色婷婷成人在线| 97视频久久| 欧美日本韩国亚洲| 色五月婷婷大| 99热这里有精品| 欧美97超碰| 男人天堂AV在线一区二区| 五月丁香色婷婷色| 六月婷婷五月天| 熟妇无码乱子成人精品| 成人五月天在线观看| 九九無妻| 国产操逼网站| 丁香五月综合在线视频| 五月做爱| 五月婷婷无码| 伊人9草在线观看| 五月天婷婷深深爱| 久草丁香婷婷1024| 天天操夜夜爽歪歪| 丁香五月婷婷亚洲另类| 欧美婷| 五月丁香六月婷婷色情| Caop在线| 色五月天堂| 狠狠综合色网| 9久热免费视频99| 成人精品一区日本无码网| 亚洲精品在线视频| 色五月色五天免费视频| 欧美成人AAA片一区国产精品| 99亚洲天堂| 青青草原福利在线| 婷婷干六月综合旧址| 五月婷婷综合网| www久久艹| 欧美槡BBBB槡BBB少妇| 丁香五月激情六月欧亚激情综合导航 | 久久九精品| 日日撸日日操| 国色天香成人网| 538午夜激情| 激情AV| 久久婷婷91| 五月天深爱激情网| 91oumei| 日日日日操| 亚洲午夜在线视频| 天天在线XXX| 色一情一乱一乱一区9| 亚洲va欧洲va国产va不卡| 999热这里只有美国精品| 天天射射夜| 久久超视频| 91成人看片| 五月色婷婷影院| 五月丁香久人妻中文| 五月丁香成人| 5月丁香综合网| 《亚洲操B久久免费在线观看,亚洲操B久久在线播放》在线播放 - 高清资源 - 97 | 97五月天| 999影院成人在线影院| 性爱五月婷| 五月丁香婷婷综合网| 五月天伊人久久久久| 色噜噜五月天| 亚洲中文AV网站| 亚洲激情免费视频观看| 五月丁香六月综合情在线观看| 狠狠草在线观看| 99人这里只有精品| 啪精品| 婷婷少妇激情| 色 五月俺去也| 79精品视频| 五月天激情图片网| 人人操AV| 991精品在线视频| 久操人妻| 狠狠爱丁香婷| 中文字幕在线日亚州9| 国产熟女一区二区三区五月婷| 玖玖婷婷五月天| 亚洲中文字幕网| 亚洲无AV在线中文字幕 | 丁香成人视频| 梁铮版蜘蛛女在线观看| 九九人人看| 婷婷天天日婷婷| 久草婷婷网| 狠狠插日日干撸| 狼人婷婷久久| 久久综合香蕉国产国产蜜臀AV| 久久丝袜婷婷| 色狠狠综合| 九九99男女视频在线观看| 五月天伊人久久久久| 超碰免费人人| 亚洲色图日韩网址| 亚洲欧美综合7777色亭亭| 玖玖在线| 99精品在线观看视频| 亚洲成人网站在线观看| 超碰99资源站| 成人AV在线网站| 日本美女上人| 婷婷五月丁香香蕉| 免费成人中文字幕| 人操人| 九九热精品| 啪啪五月天啪啪| 一二区成人电影| 亚洲综合婷婷| 五月综合在线| 成人国产欧美大片一区| 久久久中文| 婷婷色五月综合| 91精品电影18T| 日本九九视频| 婷婷五月天综合久久| 大香蕉福利导航| 天天爽天天爽| 日韩色色视频| 丁香激情网| 五月天婷婷av| 超碰女人天堂| 日本久久高清| 99爱视频免费| 婷婷在线播放| 99热这里都是精品| 婷婷五月天综合AV| 五月丁香婷草| 色色免费网站| 婷婷五月天基地| 丁香六月激情四射| 2025天天操| 天天草天天舔| 色婷婷先锋| 欧爱综合视频| 天天色综合图片| 欧美激情五月| 久久婷婷综合五月天| 午夜激情五月天| WWW.99热| 婷婷五月av| 丁香五月婷久久| 婷婷夜夜操| 美女妹子后射视频网站在线观看| 天天做好综合色| 国产色色色色| 99热亚洲| 日韩啪啪视频| 淫视馆av三区| 成人AV中文字幕| 六月婷婷香蕉| 色婷丁香五月| 这里只有精品9| 国产玖玖资源| 99精品久久久| 大陆极品少妇内射AAAAAA| 国产熟妇的荡欲午夜视频| 丝袜大香蕉| 极品少妇XXXX精品少妇偷拍| 99在线免费视频播放| 激情婷婷网| 99在线视频播放| 九色视频91疯狂| 久久这里都是精品免费| 激情五月婷婷色播网| 亚洲丁香五月在线观看| 97香蕉久久超级碰碰高清版| 五月丁香激情啪啪| 日本久碰| 天天综合久久| 色五月天婷婷| 99热| 91九色中文| 婷婷五月天综合网| 狠狠五月激情婷婷直播片| 婷婷基地爱| 亚洲网站999| 天天插综合| 成人午夜视频精品一区| 久99热在线观看| 五月天婷婷色在线视频免费观看 | 丁香六月婷婷| 丁香五月婷婷基地| 色五月婷婷自拍| 99婷婷五月天| 五月婷亚洲精品| 色色色色五月天| 《亚洲操B久久免费在线观看,亚洲操B久久在线播放》在线播放 - 高清资源 - 97 | 97干干干丁香| 激情九月婷婷| 日本成人噜噜噜| 99热国产| 色婷婷欧美| 九热久| 五月婷六月综合在线观看| 五月综合激情久久| 国产成人高清| 婷婷色色色| 五月香婷婷| 婷婷激情六月天视频| 色色色999| 久久伊人大香蕉| 国产综合A片| 五月丁香猫咪久久婷婷综合视频激情四射网入口 | 91九色在线视频| 婷婷5月天激情综合| 99色婷婷| 久久综合色五月| 六月丁香深深爱| 99色在线观看视频| 欧亚中文A V| 亚洲激情五月婷婷日日| 丁香六月啪啪| 国产精品视频免费看| 精品一区二区三区木瓜| 丁香五月天堂网| 五月天激情四射网站| 丁香五月婷婷综合网| 五月天激情小说网| 18av天堂| 高清无码网址| 婷婷五月播| 色碰碰| www网站在线观看| 五月激激网w'w'w| 日本久久综合| 中文字幕人妻在线| 亚洲欧美婷婷五月色综合| YW无码| 婷婷色5月天在线。| 大香蕉AV在线| site:pzdcoin.com| 99热66| 操操操www.com| 综合大香蕉| 91疯狂操操操操| 婷婷播5月| 夂夂夂夂夂夂夂夂夂夂夂夂夂夂夂夂夂夂夂亚洲亚洲亚洲亚洲亚洲亚洲亚洲亚洲色 | 久久婷婷五月综合啪| 日本nghangse中文字幕| 五月婷婷丁香综合| 欧美碰碰| 色色色五月婷| 色婷婷免费观看| 丁香蜜臀黄色婷婷五月天| 色吊丝永久访问网址| 五月天 无码| 在线综合91| 大陆肏屄视频| 日韩免费99| 亚洲第一综合| 五月激情综合网| 五月婷婷丁香综合| 欧美噜噜免费观看| 91丨九色丨国产| 天天天天色天天天天天干| 色婷婷亚洲综合av| 丁香五月婷婷激情网| 亚洲第一黄网| 第四色在线观看| 婷婷之六月丁香| 五月丁香激情综合网官网| 涩婷婷五月天| 婷婷性爱网| 五月丁香激情综合网| 色综合播放| AA片在线观看视频在线播放| 99精品综合视频| 激情综合网址| 七月婷婷色香综合网| 色五月丁香伊人五月| 欧美伊人9| 婷婷色丁香五月| site:jszngf.com| 欧美日韩欧美| 婷色五月天| 色丁香五月婷婷| 久草热8精品视频在线观看| 五月天色婷婷伊人网| 成人在线日韩欧美| 免费观看18视频网站| 伊人五月天婷婷| 日本黄色在线观看| 成人午夜天| 26uuu另类亚洲欧美日本一| 色综合视频| 99精品无码| 婷婷五月天激情影片| 天久久久久| 色婷婷亚洲婷婷| 五月开心深深爱激情综合| 狠狠色丁香久久| 亚洲乱码精品久久久久.. | 大色鬼综合| www,超碰| 91avse| 风流少妇A片一区二区蜜桃| 综合色五月| 丁香六月婷婷缴情欧美| 色婷婷狠狠久久YY| 91人人爱| 五月婷在线| 九九热这里只有精品6| 亭亭丁香aV| 亚洲色图啪啪| 亚洲狠狠终合停停终合| 五月天激情国产综合婷婷婷| 婷婷五月天丁香社区| 五月丁香亚洲五月| 五月天天丁香婷婷在线中| 亚洲无码激情| 青吴乐视频| 亚洲六月婷婷| 色婷婷久久| 天天爽天天摸天天爱| 管管補管管紱| 激情五月天啪啪视频| 久久人人九| 婷婷丁香色女人| 久色激情| 九九热在线精品视频| 色婷婷五月天天天做| 大香蕉五月婷婷| 成人 AV播放| 国产毛片精品一区二区色欲黄A片| 日本天堂爱爱| 六月丁香网| www.五月天。com| 俺也去婷婷五月天第五色| www.99日本| 91啦丨九色丨刺激中文| 久热99视频在线观看| 丁香五月性爱| 欧洲色色| 国产XXXX搡XXXXX搡麻豆| www.99热这里只有精品| 九九婷婷网五月天| 另类图片五月天| 色爱五月天| 久久久久亚洲AV无码网影音先锋| 丁香五月中文字幕色播| 五月天激情小说网| 久久综合中文字幕| 五月天激情图片| 99超级碰碰| 丁香婷婷激情五月天无毒不卡蜜桃| 激情五月婷婷色播网| 亚洲99手机免费看视频| 五月天婷婷激情春色小说| 天堂婷婷丁香六月网| 色另类五月天| 天天色天天色天天色天天色天天色| 国产丁香五月天婷婷| 丁香六月婷婷高清| 99热只有精| 色色色国产| 思思热AV| 婷婷大乡焦噜噜| 99热精品99| 99热新网址| 九九热在线99| 婷婷99狠狠| 色热久资源| 九九99视频精品| www.操逼comm| 丁香五月婷婷在线观看| 婷婷成人综合| 夜精品无码A片一区二区蜜桃| www.狠狠| 第四色色色色色丁香五月天| 日韩精品超碰在线观看| 五月天无码| 激情网第四色| 人人爱人人草| 亚洲另类电影| 五月婷久草| 激情五月色婷婷| 九九综合色| 婷丁香五月天| 五月天色综合| 91碰碰| 人人摸人人搞| 婷婷在线中文字幕| 色综合中文| 色综合色综合婷婷热| 成人草榴视频| 婷婷六月综合基地| 婷婷色无码| 另类小说激情五月天| 久久狼人天堂| 婷婷亚洲在线| 久久无码激情视频| 九九视频这里只有精品| 天天爱天天操| 超碰无码318604| 777久久综合视频| 婷婷开心激情五月激情网| 淫视馆AV在线| 五月婷护士| 丁香五月性爱| 色色亚洲五月天| 成人va在线播放| 亚洲色区17| 久热这里只有精品性色AV| 激情五月,婷婷五月,丁香五月| 香蕉综合在线| 99热在线网站| 日韩少妇内射免费播放| 9热视频在线观看| 甈你aaaaa| 97人人草| 色色色色色五月丁香| 婷婷五月天六月综合| 超碰资源在线| 去色色五月天| 天天狠狠色综合| 色噜噜狠狠色综合日日| 六月婷婷综合久久| 婷婷丁香五月色| 校花娇喘呻吟校长陈若雪视频| 成人 视频免费观看网站| 日本在线噜噜| 亚洲精品又粗又大又爽A片 | 色婷婷黄色网络| www.天天干| 久久九九激情五月天 | 色欲久久综合| 日韩成人网址| 色综合婷婷| 99热99思午夜精品| Www,五月天| 日韩黄在免| 综合激情啪啪| 六月丁香啪啪啪| 99久久久免费| 综合色五月天| 久色网址| 四色99久久| 美女黄频aⅴ视频| 丁香激情网| 亚洲欧洲一二| 亚洲va欧美va天堂v国产综合| 97人碰人操| 婷婷伊人激情婷婷| 国产精产国品一二三在观看| 米奇激情婷婷| 色色婷婷综合| 性生活久久人妻| 激情五月天婷婷播播久久综合91| 天天色天天噜| 色青青电影色五月| 79色色色色| 在线视频九色97| 五月丁香激情深爱婷婷| 黃色三级三级三级三级 qixing300.shrkbk.com www.jinbozs.com tianmiaosw.com | www。五月天激情| 久久婷婷色色| 欧美成人Va| 婷婷综合在线| 色综合综合色| 五月天丁香啪啪啪啪| 激情www| 婷婷色六月| 日本婷婷五月天| www久久艹| 日韩啪啪视频| 五月综合婷婷开心网| 日日躁夜夜躁狠狠久久AV| 婷婷六月天激情影院| 成人AV片播放| 噜噜噜久久| 综合另类视频| A片试看120分钟做受视频红杏 | 婷婷五月天成人在线视频| 九九视频这里只有精品| 欧美综合激情五月丁香| 丁香花大香蕉婷婷综合| 日本五月天激情| 婷婷综合色图| 狠狠CAO日日穞夜夜穞AV| 久久久久久综合五月婷婷| 激情五月狠狠喔| 日韩在线9| 99热精品在线播放观看| 9久热在线视频| 区美毛片子| 天天爽成人综合网站| 婷婷最新地址| 色九月婷婷综合| 第四色五月激情网| 色色丁香色五月| 99热久| 狠狠干五码| 在线免费视频caop| 99热最新国内| 91狠狠色| 日韩十国产极品久久| 色婷婷综合亚洲| 精品香蕉99久久久久网站| 狠狠狠狠狠狠狠狠| 欧美色色色色色色色色色色影视| 免费成人中文字幕| 色久综合| 天天综合五月| 激情综合网激情五月丁香| 婷久久| 日韩限制级大尺度黑料泄密大尺度视频一区二区在线观看 | 婷婷丁香综合在线| 天天综合网站| 色色五月天婷婷丁香| 开心五月激情站| 91视频综合网| 日韩欧美成人一区二区三区| 尔尔AV一区| 九九视频在线观看| 久色五月婷婷综合| 婷婷五月色播网| 婷婷中文字暮| 这里只有精品96| WWW,五月| 五月综合丁香婷婷| se99视频| 欧美日本一区二区三区| 狠狠色噜噜狠狠色噜噜噜999| 亚洲激情| 久久日婷婷| 五月色网| 在线中文av| 五月激情站| 五月综合影院| 综合激情综合啪啪| 亚州激情网站无码| 亚洲色精彩| 狠狠爱婷婷爱| 婷婷丁香色性爱| 99操九九网| 欧美成人AAA片一区国产精品| 色综合丁香| 丁香婷婷五月| 这里只有免费精品| 久99久热只有精品国产99| 就要爱综合| 九色视频91疯狂| 久久视频婷婷| 91九色国产熟女| www.激情com| 亚洲AV电影美洲AV电影| 五月色丁香婷婷综合| 日本波多野结衣视频| 偷拍丁香九月激情| 久久久久9久无码视频| 极品人妻VIDEOSSS人妻| 夫妇交换刺激做爰| 农村熟妇高潮精品A片| 婷婷色在线| 狠狠穞A片一區二區三區| www夜夜操| 色五月 五月婷婷| 婷婷综合五月| 国产成人精品123区免费视频| 久热99| 婷婷五月天伊人| 婷婷五月综合激情| 91久久精品国产91性色TV| 五月天久久色| 高清 码 免费看片短视频| 大香伊人婷婷| 偷拍91九色| 激情综合色网| 五月婷婷色男女| 国产五月婷| 日日夜夜狠狠干| 韩国激情五月天综合网| 爽极品色| 五月激情视频| 久久久久9久无码视频| 婷婷五月激情在线| 一起草无码| 五月丁香六月色婷婷综合五月天| 婷婷激情综合无月| 婷婷婷婷婷婷婷五月丁香| 99综合自拍| 五月天激情丁香| 五月网站| 婷婷九月在线| 色婷婷色五月另类综合| 99热只有这里才是精品| 91九色在线视频| 天天干天天爽| 嫩草视频。| 久久大香蕉视频| 五月丁香色婷婷婷基地| 丁香五月欧美激情| 色九九综合| 亚洲深喉AV| 丁香五月婷婷99| 九九热AV| 思思 热 99| 婷色五月| 白度黄视频| 日韩另类| 婷婷六月五月天综合| 婷婷在线综合| 丁香婷婷成年| 午夜大香蕉| www激情五月天| 婷婷色综合| 99视频啪啪| 久久婷婷青草五月天| 五月天综合图片| 婷婷丁香人妻天天爽| 99色亚洲| 日欧一片内射VA在线影院| 激情 婷婷| 视频这里只有精品16| 操操自拍| 亚洲成人丁香花| 精品九九在线观看视频| 五月天开心网| 婷婷在线日韩综合| 丁香五月婷婷啪| 天天摸天天透天天舔| 做爰丰满少妇1313| 色欲天天综合网| 亚洲精品午夜国产va久久成人| 丁香五月伊人| 91人妻人人操人人爽| 26uuu国产激情视频| 色丁香五月婷婷| 开心五月婷婷婷美女| A久网| 久热这里只有国产| 综合久久99| 日本一級黃色一級片| 538午夜激情| 大香蕉婷婷久久| 五月婷婷丁香综合,亚洲天堂| 日日操,日日爽| 99热这里只有精品22| 亚洲综合一区二区| 日本va欧美va欧美精品88| 激情婷婷丁香| 狠狠狠人妻| 99精品视频在线观看| 色噜噜狠狠色综合日日| 99九九99九九九视频精彩| 色爱综合五月| 婷婷五月天免费视频| 亚洲精品第一国产综合亚AV| 亚洲综合视频一下| 情欲禁地| www.伊人天堂偷偷婷婷| 乱精品一区字幕二区| 婷婷五月丁香五月丁香| 二色av| 99久久婷| 97色色综合| 激情五月综合网| 99精品在这里| 噜噜色五月| 色综合久| 97久久香草精品视频| 成人做爰A片免费看视频| 在线,国产,色,热视频| www.狠狠| WWW、日本色丁香、co m| 成人做爰A片免费看视频| 91精品熟女| 亚洲日韩一页精品发布| 天天操夜夜玩!| 五月丁香六月激情网站| 婷婷五月天堂一本在线| 丁香六月婷婷色XXXXX| 九九视频在线观看| 熟女网站久久| 国产成人亚洲综合A∨婷婷| 色天天综合天天综合频道。| 婷婷五月天网| 精品成人无码A片观看香草视频| 丰滿爆乳一区二区三区| 色色色色av色色色色| 大香AV| 狠狠色综合网| 色五月av伊人| 啪啪啪丁香五月| 五月天色区| 九九九热精品| 专区无日本视频高清8| 区区欧美你爱| 久久精品99久久久久久| 婷婷五月天影院| 五月天精品综合在线| 国产性爱色| 丁香五月色情| 亚洲字幕AV一区二区三区四区| 99热99热在线观看| 狠狠色综合图片| 婷婷久久免费| 99免费热在线精品| 五月婷婷激情啪啪| 色综合久久88色综合天天99| 在线视频色五月| 天堂久久久久天堂网| 五月天停婷基地| 5月激情天| 五月婷婷电影院| 久久总和99| 综合五月亭亭9| 天天综合91入口| 九九综合色| 日日干天天爽| 丁香五月天色| 色婷成人狠干| 中文字幕无码人妻少妇免费视频| 色五月婷婷色五月| 国产精品久久99| 久久综合香蕉国产国产蜜臀AV| 色五月天堂| 亚洲噜色| 91丨九色丨东北熟女| 久久91久久精品久久| 久婷婷| 久久色五月天综合网| 大香蕉av在线| 五月天婷婷涩涩| 大学生高潮无套内谢视频| 激情五月天com| 99re在线观看视频| 婷婷色啪| 五月天大香蕉视频| 亚洲va综合va国产va中文| 久久婷婷五月综合| 色五月五月天| 六月丁香婷婷五月天| 99热999| av一级棒av| 久久只有18视频| 五月婷婷六月丁香综合在线| 管管補管管紱| 婷婷综合在线| www.sebowuyue| 婷婷五月天综合久久日| 色五月婷婷天天干| 亚洲午夜视频| 婷婷久久伊人| 色色六月| 亚洲综合网在线| 亚洲色五月| 色999亚洲人成色| 色五月婷婷大香蕉| 999热在线视频| 99亚洲欧洲| 亚洲综合激| 成人色五月天婷婷| 成人欧美日韩| 91久久综合亚洲噜噜成人在线| 六月丁香AV| av免费在线看不卡无毒| 人妻丰满精品一区二区A片| 香蕉久久国产AV一区二区| 欧美日韩成人在线网| 99re这里只有精品视频了| 久久996re热这里只有精品无码| 日日噜人人人做人| 综合亚洲六月婷婷在线| 五月激情视频| av在线免费播放| 人妻久热| 欧美在线看| 欧美在线看| 97欧美在线| A级毛片高清免费不卡播放谢谢谢谢| 六月丁香啪| 日韩 mm 不卡| 午夜日日| 99九九玖玖| 久久网日本| 天天爽人人综合免费7799| 欧美一级色| 大香蕉婷婷久久| 97久久精品| 久久九九网| 99热这里只有精品在线观看| 天天综合精品| 日本欧美成人片AAAA| 久久精品A片777777| 激情综合丁香六| 97久久精品视频| 国精产品一区二区三区| 干一干xxxx| 九九九九国产| 五月色亭丁香| 六月婷婷狠狠做| 丁香五月婷婷欧美成人色图| 色99网站| 激情AV在线| 丁香五月大香蕉| 9色在线| 欧美日韩成卜| 成人短视频免费| 97性视频| 欧美群妇大交乱婬网| 五月丁香色五月| 亚洲亚洲永久无码777777| 五月天婷婷综合网| 五月婷婷婷| 黄色中文字目| 色噜久| 久久精品99| 天天噪夜夜爽| 六月婷婷色| 少妇综合网| 激情五月天婷婷播播久久综合91| 欧美月久久| 六月色婷婷色| 在线精品97| 99九九玖玖| 婷婷丁香成人在线视频| 五月丁香六月激情综合欧美| 日韩综合久久| 五月丁香婷婷成人网| 在线99精品| 国产乱码久久| 欧美成人精品老美女噜噜噜| 亚州操操| 欧亚成人A片一区二区| 99久久色| Se.婷婷五月天| 五月激情小说网| 在线只有精品| 日韩性爱无码| 九热在线这里有精品6| 日韩五月丁香| 五月天开心网| 久久综合色五月| 狠狠色丁香99| 激情久久肏屄视频| 丁香五月婷婷图片综合| 色激情综合| 五月激情久久综合网| 五月婷婷开心网| 操操熟女| 婷婷激情蜜桃玖玖丁香| 激情六月丁香| 色婷婷精| 人妻内射麻豆视频| 色在线99| WWW、99热| 欧美天天五月丁香免费观看| 日本爆乳片手机在线播放| 免费观看全黄做爰的视频| 五月婷婷人人人操| 久久性爱视频这里只有精品| 婷婷色在线| 亚洲成人影视在线观看| 欧美婷婷色| 无码地址| 色色五月天丁香| 色停停影院五月天| 森林影视大全,最好看的2019年视频 | 久久视频在线视频| 婷色五月天| 色色性爱视频| 五月婷婷色色爱| 丁香五月丐人妻| 久久久999精品| 五月天基地| 四月婷婷丁香| 免费观看全黄做爰的视频 | 丁香五月天激情综合| 丁香九月激情在线视频| 激情网第四色| 欧美交换配乱吟粗大25P| 日本久久99| 色播婷婷五月天| 五月天三级久久| 婷婷情爱五月天6| 国产黄色大片| 五月丁香激情欧洲啪啪| 51精品国自产在线| 深爱婷婷色| 男女啪啪做爰高潮无遮挡| 天堂AV三级| 激情综合五| 看婷婷五月天网| 四色99久久| 97操操| 久久久五月五丁香| 婷婷狠狠18禁久久| 99精品97| 中文字幕人妻在线| 国产在线黄色| www.yw尤物| 婷婷五月天成人在线视频| 久久丁香| 9er热在线精品视频| 天天爽天天干| 国产永久一二一起草| 婷婷伊人网| 综合亚洲五月天| 婷婷五月在线视频| 678五月丁香亚洲综合| 丰满少妇乱A片无码| 亚洲人妻一区二区 | www.99婷婷| 91热在线| 免费国产VA国产免费| 五月婷婷影院| 五月婷婷co.m| 九色综合五月天婷五月| 妻久久人久久| 碰97久久| 五月天婷婷无码| 色色色综合视频| 黄色热99| 天天做好综合色| 久婷久婷| 六月丁香啪啪啪| 日日操日日干| 久久aaaaa| 久婷婷| 婷婷成人综合免费视频| 色五月天婷婷婷婷婷婷婷婷婷婷婷婷婷婷婷婷婷婷婷婷婷婷婷婷婷婷婷婷 | 夜色.cnm| 91精品久久久久久久久| 伊人影音无码一区二区三区| 免费视频WWW在线观看网站| 婷婷五月中文字幕国产| se.久久视频在线观看| 色丁香五月婷婷| 337p午夜影院| 天堂网色婷婷| 99ri在线视频| 99热这里精| 综合久久久婷| www激情| 狠狠插狠狠| 免费播放片大片| 色久综合| 五月天六月色| 888精品福利地址| 偷拍99在线视频观看| 激情五月天色网站| 九玖欧洲亚洲| 精品久久9| 操操操97| 天堂久久婷婷| 日韩六十路91性交电影| 小骚穴电影| 少妇水多A片太爽了| 天天色粽合合合合合合合| 亚洲综合新99视频| 色婷亚洲五月丁香| 99色热| 79精品视频在线观看,| 丁香五月婷婷亚洲人| 无码人妻一区二区一牛影视| 五月天精品视频| 日日插日日干| 丁香丁香激情网| 日本久草福利| 婷婷五月天VI| 五月亭亭六月天| 直接看的AV| 色婷婷丁香五月观看| 任你搞网站| 天天色凹凸| 国产亚洲99久久精品| 五月婷婷色播| 色情五月天A片| 色丁香五月天射婷婷爱婷婷| 99热很操老逼| 强辱丰满人妻HD中文字幕| 婷婷丁香成人| 色中色综合| 深爱丁香激情| 激情五月天婷婷免费观看| 丁香五月天婷婷91| 六月丁花香啪啪激情欧美| 亭亭五月色男人| 99热精品在线观看| 丁香五月综合在线观看| 亚洲五月天天| 九九热超碰| 香蕉婷婷色五月| www.激情五月天.com| AA片在线观看视频在线播放| 色色五月婷婷| 婷婷午夜激情| www.99热精品| 毛片色五月| 丁香五月人妻熟女| 视频免费精品免费精品免费精品免费精品免费精品免费精品免费99 | 五月婷婷深深爱| 五月婷婷六月丁香激情| 天天天摸夜夜夜玩| 激情五月天之六月婷婷| 五月综合亚洲色| 成人午夜无码视频| 久久久99视频| 影音先锋毛片网站| 婷婷丁香五月天影院 | 丁香五月九九| 九九99视频精品| 九色99视频| 岛国在线观看91| 国产精品久久久久9999小说 | AV在线免费网站| 婷婷五月天小说| 狠狠操狠狠操| 五月丁香在线观看| 色五月激情| 色婷婷亚洲六月婷婷中文字幕| 日本激情综合| 这里只有精彩视| 岛国AV网站| 2025色婷婷| 来吧亚洲综合网| 色婷婷六月天| 五月停停色色丁香| 色色色婷婷五月| 午夜天堂一区人妻| 日韩成人综合网| 九月婷婷| www.色五月| 久久久天堂国产精品女人| 97超碰免费超级在线观看| 欧美婷婷丁香社区在线播放| 狠狠婷婷色| 九九99九九精品视频| 久久久噜噜噜www成人| 国产亚洲色婷婷久久99精品91 www.riverspirits.org www.hnnun.com www.changh | 怡红院视频| 婷婷另类小说| 色色五月天丁香婷婷| 日韩黄黄| www久| 国产做A爰片毛片A片美国| 婷婷97色| 激情性爱网站| 午夜爱插插| 开心六月丁香五月婷婷| 26uuu激情五月天| 丁香六月成人| 另类精品视频在线观看| 天天插天天插天天插天天插| 人妻体体内射精一区二区 | 五月丁香久久| 丁香九月激情| 风流少妇A片一区二区蜜桃 | 国产淫熟妇| 亚洲成人婷婷| 婷婷综合激情| 日韩一本操| 69色婷婷| WWW·天天操·视频?| 丁香激激情网| 九九大香蕉黄色影院| 天天综合网站| 午夜在线成人网站免费观看| 日韩在线观看亚洲| 七七九九色色| 九九久热| 99热最新| 久草 tingting| 婷婷丁香成人| 婷婷五月色综合| 色色丁香五月天| 国产AV影片| 91互操| 五月综合婷婷网| 欧美日韩成人在线网站| 狠狠狠狠狠干| 久久婷婷丁香五月宗合| www色婷婷久久综合久色| 玖玖热视频| 97色婷| 男人天堂 久久| 精品视频99看在线视频| 泰州成人视频| 国产特级毛片AAAAAAA高清| www.五月天色色色| 天天成人五月天| 99熟女视频| 日操夜操天天操不卡| 丁香五月激情综合| 日日爱激情| 色九月婷婷| 五月丁香网站| 色性五月天| 五月婷婷丁香六月| 日本3级片一区2区| 久久九九99亚洲国产久精综合| 五月色色色| 婷五月天| 五月丁香亚洲五月| www.狠狠| 2021日韩无码| 亚洲99在线| 夜夜操狠狠操| 午夜成人综合| 日韩色五月| 欧美五月停| 激情五月第四色| 99A级片| 五月久久| 思思热再线视频| 婷婷五月天综合网| 综合大香蕉| 久久婷婷五月天| 爱iii做iiii日日| 色综合久久伊伊婷婷五月| 成片免费播放| 丁香婷婷性久久| 激情五月视频在线婷婷| 97久久人人人干| 成人婷婷桔色| 久久综合五月| 97香蕉碰碰人妻国产欧美| 国产日韩欧美性爱| 久久国产高潮白浆免费观看99| 五月天久久www| 亚洲精品又粗又大又爽A片| 91久久久久久| 久久激情五月| 99在线观看视频| 久久伦乱| 五月丁香综合成人社区| 狠狠干在线| 天天肏天天肏| 亚洲综合热| 色综合播放| 欧美va欧美va差| 天天干天天av天天射| 丁香婷婷久久老熟女综合网| 国精产品一区一区三区免费视频| 99在线免费视| 噼里啪啦在线观看免费完整版视频| 婷婷在线播放av| 色色五月天婷婷| 婷婷五月激情网| 婷婷丁香五月亚洲| 快乐激情五月色婷婷| 五月天桃色深爱网| 華人性愛AV在線| 成人在线日韩| 五月天色导航| 69久久久|