亭亭五月天在线观看,亭亭五月天在线观看,国产最新av一区二区,国产 高清 中文字幕,99re热久久亚洲综合精品成人,熟妇 一区二区三区,一级做a爰片性色毛片武则天,美女的骚穴视频播放,国产美女午夜免费视频

24小時(shí)熱門版塊排行榜    

查看: 3659  |  回復(fù): 11

mengfc

金蟲 (正式寫手)

[求助] 編譯vasp出錯(cuò)

用ifort+mkl+mpich編譯vasp
提示error #10236: File not found:  '/opt/libs/libgoto/libgoto_p4_512-r0.6.so'

makefile如下,請(qǐng)大家?guī)兔纯?br /> .SUFFIXES: .inc .f .f90 .F
#-----------------------------------------------------------------------
# Makefile for Intel Fortran compiler for P4 systems
#
# The makefile was tested only under Linux on Intel platforms
# (Suse 5.3- Suse 9.0)
# the followin compiler versions have been tested
# 5.0, 6.0, 7.0 and 7.1 (some 8.0 versions seem to fail compiling the code)
# presently we recommend version 7.1 or 7.0, since these
# releases have been used to compile the present code versions
#
# it might be required to change some of library pathes, since
# LINUX installation vary a lot
# Hence check ***ALL**** options in this makefile very carefully
#-----------------------------------------------------------------------
#
# BLAS must be installed on the machine
# there are several options:
# 1) very slow but works:
#   retrieve the lapackage from ftp.netlib.org
#   and compile the blas routines (BLAS/SRC directory)
#   please use g77 or f77 for the compilation. When I tried to
#   use pgf77 or pgf90 for BLAS, VASP hang up when calling
#   ZHEEV  (however this was with lapack 1.1 now I use lapack 2.0)
# 2) most desirable: get an optimized BLAS
#
# the two most reliable packages around are presently:
# 3a) Intels own optimised BLAS (PIII, P4, Itanium)
#     http://developer.intel.com/software/products/mkl/
#   this is really excellent when you use Intel CPU's
#
# 3b) or obtain the atlas based BLAS routines
#     http://math-atlas.sourceforge.net/
#   you certainly need atlas on the Athlon, since the  mkl
#   routines are not optimal on the Athlon.
#   If you want to use atlas based BLAS, check the lines around LIB=
#
# 3c) mindblowing fast SSE2 (4 GFlops on P4, 2.53 GHz)
#   Kazushige Goto's BLAS
#   http://www.cs.utexas.edu/users/kgoto/signup_first.html
#
#-----------------------------------------------------------------------

# all CPP processed fortran files have the extension .f90
SUFFIX=.f90

#-----------------------------------------------------------------------
# fortran compiler and linker
#-----------------------------------------------------------------------
FC=ifort
# fortran linker
FCL=$(FC)


#-----------------------------------------------------------------------
# whereis CPP ?? (I need CPP, can't use gcc with proper options)
# that's the location of gcc for SUSE 5.3
#
#  CPP_   =  /usr/lib/gcc-lib/i486-linux/2.7.2/cpp -P -C
#
# that's probably the right line for some Red Hat distribution:
#
#  CPP_   =  /usr/lib/gcc-lib/i386-redhat-linux/2.7.2.3/cpp -P -C
#
#  SUSE X.X, maybe some Red Hat distributions:

CPP_ =  ./preprocess <$*.F | /usr/bin/cpp -P -C -traditional >$*$(SUFFIX)

#-----------------------------------------------------------------------
# possible options for CPP:
# NGXhalf             charge density   reduced in X direction
# wNGXhalf            gamma point only reduced in X direction
# avoidalloc          avoid ALLOCATE if possible
# IFC                 work around some IFC bugs
# CACHE_SIZE          1000 for PII,PIII, 5000 for Athlon, 8000-12000 P4
# RPROMU_DGEMV        use DGEMV instead of DGEMM in RPRO (depends on used BLAS)
# RACCMU_DGEMV        use DGEMV instead of DGEMM in RACC (depends on used BLAS)
#-----------------------------------------------------------------------

CPP     = $(CPP_)  -DHOST=\"LinuxIFC\" \
          -Dkind8 -DNGXhalf -DCACHE_SIZE=12000 -DPGF90 -Davoidalloc \
#          -DRPROMU_DGEMV  -DRACCMU_DGEMV

#-----------------------------------------------------------------------
# general fortran flags  (there must a trailing blank on this line)
#-----------------------------------------------------------------------

FFLAGS =  -FR -lowercase -assume byterecl

#-----------------------------------------------------------------------
# optimization
# we have tested whether higher optimisation improves performance
# -axK  SSE1 optimization,  but also generate code executable on all mach.
#       xK improves performance somewhat on XP, and a is required in order
#       to run the code on older Athlons as well
# -xW   SSE2 optimization
# -axW  SSE2 optimization,  but also generate code executable on all mach.
# -tpp6 P3 optimization
# -tpp7 P4 optimization
#-----------------------------------------------------------------------

OFLAG=-O3 -xW

OFLAG_HIGH = $(OFLAG)
OBJ_HIGH =

OBJ_NOOPT =
DEBUG  = -FR -O0
INLINE = $(OFLAG)


#-----------------------------------------------------------------------
# the following lines specify the position of BLAS  and LAPACK
# on P4, VASP works fastest with the libgoto library
# so that's what I recommend
#-----------------------------------------------------------------------

# Atlas based libraries
#ATLASHOME= $(HOME)/archives/BLAS_OPT/ATLAS/lib/Linux_P4SSE2/
#BLAS=   -L$(ATLASHOME)  -lf77blas -latlas

# use specific libraries (default library path might point to other libraries)
#BLAS= $(ATLASHOME)/libf77blas.a $(ATLASHOME)/libatlas.a

# use the mkl Intel libraries for p4 (www.intel.com)
# mkl.5.1
# set -DRPROMU_DGEMV  -DRACCMU_DGEMV in the CPP lines
#BLAS=-L/opt/intel/mkl/lib/32 -lmkl_p4  -lpthread

# mkl.5.2 requires also to -lguide library
# set -DRPROMU_DGEMV  -DRACCMU_DGEMV in the CPP lines
#BLAS=-L/opt/intel/mkl/lib/32 -lmkl_p4 -lguide -lpthread

# even faster Kazushige Goto's BLAS
# http://www.cs.utexas.edu/users/kgoto/signup_first.html
BLAS=  /opt/libs/libgoto/libgoto_p4_512-r0.6.so

# LAPACK, simplest use vasp.4.lib/lapack_double
LAPACK= ../vasp.4.lib/lapack_double.o

# use atlas optimized part of lapack
#LAPACK= ../vasp.4.lib/lapack_atlas.o  -llapack -lcblas

# use the mkl Intel lapack
#LAPACK= -lmkl_lapack

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

LIB  = -L../vasp.4.lib -ldmy \
     ../vasp.4.lib/linpack_double.o $(LAPACK) \
     $(BLAS)

# options for linking (for compiler version 6.X, 7.1) nothing is required
LINK    =
# compiler version 7.0 generates some vector statments which are located
# in the svml library, add the LIBPATH and the library (just in case)
#LINK    =  -L/opt/intel/compiler70/ia32/lib/ -lsvml

#-----------------------------------------------------------------------
# fft libraries:
# VASP.4.6 can use fftw.3.0.X (http://www.fftw.org)
# since this version is faster on P4 machines, we recommend to use it
#-----------------------------------------------------------------------

FFT3D   = fft3dfurth.o fft3dlib.o
#FFT3D   = fftw3d.o fft3dlib.o   /opt/libs/fftw-3.0.1/lib/libfftw3.a


#=======================================================================
# MPI section, uncomment the following lines
#
# one comment for users of mpich or lam:
# You must *not* compile mpi with g77/f77, because f77/g77            
# appends *two* underscores to symbols that contain already an        
# underscore (i.e. MPI_SEND becomes mpi_send__).  The pgf90/ifc
# compilers however append only one underscore.
# Precompiled mpi version will also not work !!!
#
# We found that mpich.1.2.1 and lam-6.5.X to lam-7.0.4 are stable
# mpich.1.2.1 was configured with
#  ./configure -prefix=/usr/local/mpich_nodvdbg -fc="pgf77 -Mx,119,0x200000"  \
# -f90="pgf90 " \
# --without-romio --without-mpe -opt=-O \
#
# lam was configured with the line
#  ./configure  -prefix /opt/libs/lam-7.0.4 --with-cflags=-O -with-fc=ifc \
# --with-f77flags=-O --without-romio
#
# please note that you might be able to use a lam or mpich version
# compiled with f77/g77, but then you need to add the following
# options: -Msecond_underscore (compilation) and -g77libs (linking)
#
# !!! Please do not send me any queries on how to install MPI, I will
# certainly not answer them !!!!
#=======================================================================
#-----------------------------------------------------------------------
# fortran linker for mpi: if you use LAM and compiled it with the options
# suggested above,  you can use the following line
#-----------------------------------------------------------------------

FC=mpif90
FCL=$(FC)

#-----------------------------------------------------------------------
# additional options for CPP in parallel version (see also above):
# NGZhalf               charge density   reduced in Z direction
# wNGZhalf              gamma point only reduced in Z direction
# scaLAPACK             use scaLAPACK (usually slower on 100 Mbit Net)
#-----------------------------------------------------------------------

CPP    = $(CPP_) -DMPI  -DHOST=\"LinuxIFC\" -DIFC \
     -Dkind8 -DNGZhalf -DCACHE_SIZE=4000 -DPGF90 -Davoidalloc \
     -DMPI_BLOCK=500  \
##    -DRPROMU_DGEMV  -DRACCMU_DGEMV

#-----------------------------------------------------------------------
# location of SCALAPACK
# if you do not use SCALAPACK simply uncomment the line SCA
#-----------------------------------------------------------------------

BLACS=$(HOME)/archives/SCALAPACK/BLACS/
SCA_=$(HOME)/archives/SCALAPACK/SCALAPACK

SCA= $(SCA_)/libscalapack.a  \
$(BLACS)/LIB/blacsF77init_MPI-LINUX-0.a $(BLACS)/LIB/blacs_MPI-LINUX-0.a $(BLACS)/LIB/blacsF77init_MPI-LINUX-0.a

SCA=

#-----------------------------------------------------------------------
# libraries for mpi
#-----------------------------------------------------------------------

LIB     = -L../vasp.4.lib -ldmy  \
      ../vasp.4.lib/linpack_double.o $(LAPACK) \
      $(SCA) $(BLAS)

# FFT: fftmpi.o with fft3dlib of Juergen Furthmueller
FFT3D   = fftmpi.o fftmpi_map.o fft3dlib.o

# fftw.3.0.1 is slighly faster and should be used if available
#FFT3D   = fftmpiw.o fftmpi_map.o fft3dlib.o   /opt/libs/fftw-3.0.1/lib/libfftw3.a

#-----------------------------------------------------------------------
# general rules and compile lines
#-----------------------------------------------------------------------
BASIC=   symmetry.o symlib.o   lattlib.o  random.o   

SOURCE=  base.o     mpi.o      smart_allocate.o      xml.o  \
         constant.o jacobi.o   main_mpi.o  scala.o   \
         asa.o      lattice.o  poscar.o   ini.o      setex.o     radial.o  \
         pseudo.o   mgrid.o    mkpoints.o wave.o      wave_mpi.o  $(BASIC) \
         nonl.o     nonlr.o    dfast.o    choleski2.o    \
         mix.o      charge.o   xcgrad.o   xcspin.o    potex1.o   potex2.o  \
         metagga.o  constrmag.o pot.o      cl_shift.o force.o    dos.o      elf.o      \
         tet.o      hamil.o    steep.o    \
         chain.o    dyna.o     relativistic.o LDApU.o sphpro.o  paw.o   us.o \
         ebs.o      wavpre.o   wavpre_noio.o broyden.o \
         dynbr.o    rmm-diis.o reader.o   writer.o   tutor.o xml_writer.o \
         brent.o    stufak.o   fileio.o   opergrid.o stepver.o  \
         dipol.o    xclib.o    chgloc.o   subrot.o   optreal.o   davidson.o \
         edtest.o   electron.o shm.o      pardens.o  paircorrection.o \
         optics.o   constr_cell_relax.o   stm.o    finite_diff.o \
         elpol.o    setlocalpp.o

INC=

vasp: $(SOURCE) $(FFT3D) $(INC) main.o
        rm -f vasp
        $(FCL) -o vasp $(LINK) main.o  $(SOURCE)   $(FFT3D) $(LIB)
makeparam: $(SOURCE) $(FFT3D) makeparam.o main.F $(INC)
        $(FCL) -o makeparam  $(LINK) makeparam.o $(SOURCE) $(FFT3D) $(LIB)
zgemmtest: zgemmtest.o base.o random.o $(INC)
        $(FCL) -o zgemmtest $(LINK) zgemmtest.o random.o base.o $(LIB)
dgemmtest: dgemmtest.o base.o random.o $(INC)
        $(FCL) -o dgemmtest $(LINK) dgemmtest.o random.o base.o $(LIB)
ffttest: base.o smart_allocate.o mpi.o mgrid.o random.o ffttest.o $(FFT3D) $(INC)
        $(FCL) -o ffttest $(LINK) ffttest.o mpi.o mgrid.o random.o smart_allocate.o base.o $(FFT3D) $(LIB)
kpoints: $(SOURCE) $(FFT3D) makekpoints.o main.F $(INC)
        $(FCL) -o kpoints $(LINK) makekpoints.o $(SOURCE) $(FFT3D) $(LIB)

clean:       
        -rm -f *.g *.f *.o *.L *.mod ; touch *.F

main.o: main$(SUFFIX)
        $(FC) $(FFLAGS)$(DEBUG)  $(INCS) -c main$(SUFFIX)
xcgrad.o: xcgrad$(SUFFIX)
        $(FC) $(FFLAGS) $(INLINE)  $(INCS) -c xcgrad$(SUFFIX)
xcspin.o: xcspin$(SUFFIX)
        $(FC) $(FFLAGS) $(INLINE)  $(INCS) -c xcspin$(SUFFIX)

makeparam.o: makeparam$(SUFFIX)
        $(FC) $(FFLAGS)$(DEBUG)  $(INCS) -c makeparam$(SUFFIX)

makeparam$(SUFFIX): makeparam.F main.F
#
# MIND: I do not have a full dependency list for the include
# and MODULES: here are only the minimal basic dependencies
# if one strucuture is changed then touch_dep must be called
# with the corresponding name of the structure
#
base.o: base.inc base.F
mgrid.o: mgrid.inc mgrid.F
constant.o: constant.inc constant.F
lattice.o: lattice.inc lattice.F
setex.o: setexm.inc setex.F
pseudo.o: pseudo.inc pseudo.F
poscar.o: poscar.inc poscar.F
mkpoints.o: mkpoints.inc mkpoints.F
wave.o: wave.inc wave.F
nonl.o: nonl.inc nonl.F
nonlr.o: nonlr.inc nonlr.F

$(OBJ_HIGH):
        $(CPP)
        $(FC) $(FFLAGS) $(OFLAG_HIGH) $(INCS) -c $*$(SUFFIX)
$(OBJ_NOOPT):
        $(CPP)
        $(FC) $(FFLAGS) $(INCS) -c $*$(SUFFIX)

fft3dlib_f77.o: fft3dlib_f77.F
        $(CPP)
        $(F77) $(FFLAGS_F77) -c $*$(SUFFIX)

.F.o:
        $(CPP)
        $(FC) $(FFLAGS) $(OFLAG) $(INCS) -c $*$(SUFFIX)
.F$(SUFFIX):
        $(CPP)
$(SUFFIX).o:
        $(FC) $(FFLAGS) $(OFLAG) $(INCS) -c $*$(SUFFIX)

# special rules
#-----------------------------------------------------------------------
# these special rules are cummulative (that is once failed
#   in one compiler version, stays in the list forever)
# -tpp5|6|7 P, PII-PIII, PIV
# -xW use SIMD (does not pay of on PII, since fft3d uses double prec)
# all other options do no affect the code performance since -O1 is used
#-----------------------------------------------------------------------

fft3dlib.o : fft3dlib.F
        $(CPP)
        $(FC) -FR -lowercase -O1 -tpp7 -xW -prefetch- -unroll0  -vec_report3 -c $*$(SUFFIX)
fft3dfurth.o : fft3dfurth.F
        $(CPP)
        $(FC) -FR -lowercase -O1 -c $*$(SUFFIX)

radial.o : radial.F
        $(CPP)
        $(FC) -FR -lowercase -O1 -c $*$(SUFFIX)

symlib.o : symlib.F
        $(CPP)
        $(FC) -FR -lowercase -O1 -c $*$(SUFFIX)

symmetry.o : symmetry.F
        $(CPP)
        $(FC) -FR -lowercase -O1 -c $*$(SUFFIX)

dynbr.o : dynbr.F
        $(CPP)
        $(FC) -FR -lowercase -O1 -c $*$(SUFFIX)

broyden.o : broyden.F
        $(CPP)
        $(FC) -FR -lowercase -O2 -c $*$(SUFFIX)

us.o : us.F
        $(CPP)
        $(FC) -FR -lowercase -O1 -c $*$(SUFFIX)

wave.o : wave.F
        $(CPP)
        $(FC) -FR -lowercase -O0 -c $*$(SUFFIX)

LDApU.o : LDApU.F
        $(CPP)
        $(FC) -FR -lowercase -O2 -c $*$(SUFFIX)
回復(fù)此樓

» 猜你喜歡

» 本主題相關(guān)價(jià)值貼推薦,對(duì)您同樣有幫助:

已閱   回復(fù)此樓   關(guān)注TA 給TA發(fā)消息 送TA紅花 TA的回帖

saikun

木蟲 (正式寫手)

【答案】應(yīng)助回帖


youzhizhe(金幣+1): 謝謝交流。 2011-05-02 09:14:12
mengfc(金幣+10): 2011-05-02 09:33:49
mengfc(金幣+5): 謝謝回復(fù),請(qǐng)看十樓的makefile有沒有問題? 2011-05-02 16:02:46
把goto-blas那句改成用mkl或先編譯goto-blas就可以了。

但我覺得mkl比較快(xeon上) 本來intel上有一篇指南的,但好象沒了。只能在google cache里找到,你看看。

=======================
Using Intel® MKL in VASP
October 5, 2010 12:50 AM PDT

Step 1 - Overview

This guide is intended to help current VASP* (Vienna Ab-Initio Package Simulation*) users get better benchmark performance by utilizing Intel® Math Kernel Library (Intel® MKL).

VASP is a package for performing ab-initio quantum-mechanical molecular dynamics (MD) using pseudo potentials and a plane wave basis set. The approach implemented in VAMP/VASP is based on a finite-temperature local-density approximation (with the free energy as variational quantity) and an exact evaluation of the instantaneous electronic ground state at each MD-step using efficient matrix diagonalization schemes and an efficient Pulay mixing. These techniques avoid all problems occurring in the original Car-Parrinello method which is based on the simultaneous integration of electronic and ionic equations of motion. The interaction between ions and electrons is described using ultrasoft Vanderbilt pseudopotentials (US-PP) or the projector augmented wave method (PAW). Both techniques allow a considerable reduction of the necessary number of plane-waves per atom for transition metals and first row elements. Forces and stress can be easily calculated with VAMP/VASP and used to relax atoms into their instantaneous groundstate. [Ref: VASP*]

Version Information

This application note was created to help users who benchmark clusters using VASP to also incorporate the latest version of Intel® MKL on Linux* platforms on Intel® Xeon® *processor-based systems. Specifically, we'll address Intel® MKL version 10.2.

More information on VASP can be found from http://cms.mpi.univie.ac.at/vasp/*
Step 2 - Downloading VASP Source Code

VASP is not public-domain or shareware, and will be distributed only after a license contract has been signed. Enquiries must be sent to Juerrgen Hafner. The enquiry should contain a short description of the short term research aims (less than half a page).

Prerequisites:

    Intel® MKL contains highly optimized FFT and also the wrappers for FFTW, which can be obtained using the following options:
        Download a FREE evaluation version of the Intel® MKL product.
        Download the FREE non-commercial* version of the Intel® MKL product.
        Purchase a copy of Intel® MKL.
    All of these can be obtained at: Intel® Math Kernel Library Product Page.
    Obtain Intel® C++ and Fortran compilers.
    Obtain Intel® MPI.
    Obtain Open Source MPI (MPICH2)

Step 3 - Configuration

    Use the following commands to extract the VASP files:

    $tar -xvzf vasp.tgz

This will create a directory called vasp. Make sure that MPI, C++ and Fortran compilers are installed and in the PATH. Also set LD_LIBRARY_PATH to your compiler (C++ and Fortran), MPI, and Intel® MKL libraries. Note: This application note is written specifically for use with the Intel® compilers and MPI.


Step 4 - Building VASP

Build FFTW wrapper library for Intel® MKL

    Change directory to /interfaces/fftw3xf.
    Assuming default installation--

    $cd /opt/intel/mkl/10.2.0.013/interfaces/fftw3xf

Build the fftw 3.x Fortran wrapper library for Intel® 64 (previously known as em64t) platform using the Intel compilers.

    $make libem64t compiler=intel
     
This will create the fftw wrapper library "libfftw3xf_intel.a" in lib/em64t directory of your Intel® MKL installation.

Build libdmy.a

    Change directory to vasp.x.x.lib.
   
Run the following command from vasp.x.x.lib, using the Makefile for Linux* using the Intel compiler.

    #make -f makefile.linux_ifc_P4

After a successful compilation, libdmy.a will be built in the same directory.

Build VASP

    Change directory to vasp.x.x.
    Edit the makefile.linux_ifc_P4 to link with Intel® MKL libraries.
    Go to the section where BLAS, LAPACK and FFT libraries are specified and change as below, assuming the default installation /opt/intel/mkl/10.2.0.013 for MKL 10.2.


    BLAS=-L/opt/intel/mkl/10.2.0.013/lib/em64t -lmkl_intel_lp64 -lmkl_sequential -lmkl_core -lpthread  
     
    LAPACK=-L/opt/intel/mkl/10.2.0.013/lib/em64t -lmkl_intel_lp64  -lmkl_sequential -lmkl_core -lpthread  
      
    FFT3D= fftmpiw.o fftmpi_map.o fft3dlib.o /opt/intel/mkl/10.2.0.013/lib/em64t/libfftw3xf_intel.a
     
      
    Modify Fortran compiler flags line to include fftw3.f from Intel® MKL include directory as below:


    FFLAGS= -I/opt/intel/mkl/10.2.0.013/include/fftw -FR -lower_case
     

This will eliminate compilation errors pertaining to "cannot open include file `fftw3.f' when Intel® MKL FFT is used rather than VASP's own FFTs.

    To use Intel® MKL ScaLAPACK, modify additional options for CPP for parallel version adding -DscaLAPACK to the CPP line and specify Intel® MKL ScaLAPACK as below:


    SCA=/opt/intel/mkl/10.2.0.013/lib/em64t/libmkl_scalapack_lp64.a /opt/intel/mkl/10.2.0.13/lib/em64t/libmkl_blacs_intelmpi_lp64.a

    To use Intel® MPI specify Fortran compiler as FC=mpiifort.
    Run the following command to build vasp:

    $make -f makefile.linux_ifc_P4
        
This will create the VASP executable in the current directory.

Step 5 - Running VASP

    Run vasp by executing the following command with your workloads:

    $mpirun -n 4 ./vasp
2樓2011-05-01 23:34:04
已閱   回復(fù)此樓   關(guān)注TA 給TA發(fā)消息 送TA紅花 TA的回帖

hakuna

木蟲 (知名作家)


youzhizhe(金幣+1): 謝謝交流。 2011-05-02 09:14:38
BLAS=  /opt/libs/libgoto/libgoto_p4_512-r0.6.so
是這個(gè)數(shù)學(xué)庫的問題,你確認(rèn)安裝有這個(gè)數(shù)學(xué)庫嗎?
建議你用intel的數(shù)學(xué)庫,如果非要用libgoto,還是自己編一個(gè)新版本的吧,上面的這個(gè)太老了,應(yīng)該是五六年前的東西了,自己編譯一個(gè)新的也不麻煩.....
3樓2011-05-01 23:36:23
已閱   回復(fù)此樓   關(guān)注TA 給TA發(fā)消息 送TA紅花 TA的回帖

mengfc

金蟲 (正式寫手)

引用回帖:
Originally posted by saikun at 2011-05-01 23:34:04:
把goto-blas那句改成用mkl或先編譯goto-blas就可以了。

但我覺得mkl比較快(xeon上) 本來intel上有一篇指南的,但好象沒了。只能在google cache里找到,你看看。

=======================
Using Intel& ...

你好,謝謝你的回帖。我按著這上面又改了一下maekfile,還是不行。把錯(cuò)誤信息貼在后面,請(qǐng)指點(diǎn)一下。

問題一:在vasp.4.lib文件夾中make時(shí),提示
                gcc -E -P -C preclib.F >preclib.f
ifort -O0 -FI -FR  -c preclib.f
cc -O -c timing_.c
cc -O -c derrf_.c
cc -O -c dclock_.c
gcc -E -P -C diolib.F >diolib.f
ifort -O0 -FI -FR  -c diolib.f
gcc -E -P -C dlexlib.F >dlexlib.f
ifort -O0 -FI -FR  -c dlexlib.f
gcc -E -P -C drdatab.F >drdatab.f
ifort -O0 -FI -FR  -c drdatab.f
ifort -O0 -FI  -c lapack_double.f
ifort -O0 -FI  -c linpack_double.f
ifort -O0 -FI  -c lapack_atlas.f
rm libdmy.a
rm: cannot remove `libdmy.a': No such file or directory
make: [libdmy.a] Error 1 (ignored)
ar vq libdmy.a preclib.o timing_.o derrf_.o dclock_.o  diolib.o dlexlib.o drdatab.o
ar: creating libdmy.a
a - preclib.o
a - timing_.o
a - derrf_.o
a - dclock_.o
a - diolib.o
a - dlexlib.o
a - drdatab.o

這是否正確編譯了庫文件?

問題二:在vasp.4.6 文件夾中make時(shí),出錯(cuò):
elpol.f90(693): warning #7866: The statement following this DEC loop optimization directive must be an iterative do-stmt, a vector assignment, an OMP pdo-directive, or an OMP parallel-do-directive.
!DIR$ IVDEP
------^
ifort: warning #10147: no action performed for specified file(s)
./preprocess setlocalpp.f90 -DMPI  -DHOST=\"LinuxIFC\" -DIFC -Dkind8 -DNGZhalf -DCACHE_SIZE=4000 -DPGF90 -Davoidalloc -DMPI_BLOCK=500
mpif90 /opt/intel/composerxe-2011.3.174/mkl/include/fftw -FR -lowercase   -O3 -xW   -c setlocalpp.f90
ifort: command line remark #10010: option '-lowercase' is deprecated and will be removed in a future release. See '-help deprecated'
ifort: command line remark #10279: option '-xW' is deprecated and will be removed in a future release. See '-help deprecated'
ifort: warning #10147: no action performed for specified file(s)
make: *** No rule to make target `FFT3D', needed by `vasp'.  Stop.

makefile 如下,不知應(yīng)如何修改?


.SUFFIXES: .inc .f .f90 .F
#-----------------------------------------------------------------------
# Makefile for Intel Fortran compiler for P4 systems
#
# The makefile was tested only under Linux on Intel platforms
# (Suse 5.3- Suse 9.0)
# the followin compiler versions have been tested
# 5.0, 6.0, 7.0 and 7.1 (some 8.0 versions seem to fail compiling the code)
# presently we recommend version 7.1 or 7.0, since these
# releases have been used to compile the present code versions
#
# it might be required to change some of library pathes, since
# LINUX installation vary a lot
# Hence check ***ALL**** options in this makefile very carefully
#-----------------------------------------------------------------------
#
# BLAS must be installed on the machine
# there are several options:
# 1) very slow but works:
#   retrieve the lapackage from ftp.netlib.org
#   and compile the blas routines (BLAS/SRC directory)
#   please use g77 or f77 for the compilation. When I tried to
#   use pgf77 or pgf90 for BLAS, VASP hang up when calling
#   ZHEEV  (however this was with lapack 1.1 now I use lapack 2.0)
# 2) most desirable: get an optimized BLAS
#
# the two most reliable packages around are presently:
# 3a) Intels own optimised BLAS (PIII, P4, Itanium)
#     http://developer.intel.com/software/products/mkl/
#   this is really excellent when you use Intel CPU's
#
# 3b) or obtain the atlas based BLAS routines
#     http://math-atlas.sourceforge.net/
#   you certainly need atlas on the Athlon, since the  mkl
#   routines are not optimal on the Athlon.
#   If you want to use atlas based BLAS, check the lines around LIB=
#
# 3c) mindblowing fast SSE2 (4 GFlops on P4, 2.53 GHz)
#   Kazushige Goto's BLAS
#   http://www.cs.utexas.edu/users/kgoto/signup_first.html
#
#-----------------------------------------------------------------------

# all CPP processed fortran files have the extension .f90
SUFFIX=.f90

#-----------------------------------------------------------------------
# fortran compiler and linker
#-----------------------------------------------------------------------
FC=ifort
# fortran linker
FCL=$(FC)


#-----------------------------------------------------------------------
# whereis CPP ?? (I need CPP, can't use gcc with proper options)
# that's the location of gcc for SUSE 5.3
#
#  CPP_   =  /usr/lib/gcc-lib/i486-linux/2.7.2/cpp -P -C
#
# that's probably the right line for some Red Hat distribution:
#
#  CPP_   =  /usr/lib/gcc-lib/i386-redhat-linux/2.7.2.3/cpp -P -C
#
#  SUSE X.X, maybe some Red Hat distributions:

CPP_ =  ./preprocess <$*.F | /usr/bin/cpp -P -C -traditional >$*$(SUFFIX)

#-----------------------------------------------------------------------
# possible options for CPP:
# NGXhalf             charge density   reduced in X direction
# wNGXhalf            gamma point only reduced in X direction
# avoidalloc          avoid ALLOCATE if possible
# IFC                 work around some IFC bugs
# CACHE_SIZE          1000 for PII,PIII, 5000 for Athlon, 8000-12000 P4
# RPROMU_DGEMV        use DGEMV instead of DGEMM in RPRO (depends on used BLAS)
# RACCMU_DGEMV        use DGEMV instead of DGEMM in RACC (depends on used BLAS)
#-----------------------------------------------------------------------

CPP     = $(CPP_)  -DHOST=\"LinuxIFC\" \
          -Dkind8 -DNGXhalf -DCACHE_SIZE=12000 -DPGF90 -Davoidalloc \
#          -DRPROMU_DGEMV  -DRACCMU_DGEMV

#-----------------------------------------------------------------------
# general fortran flags  (there must a trailing blank on this line)
#-----------------------------------------------------------------------

FFLAGS =  /opt/intel/composerxe-2011.3.174/mkl/include/fftw -FR -lowercase  

#-----------------------------------------------------------------------
# optimization
# we have tested whether higher optimisation improves performance
# -axK  SSE1 optimization,  but also generate code executable on all mach.
#       xK improves performance somewhat on XP, and a is required in order
#       to run the code on older Athlons as well
# -xW   SSE2 optimization
# -axW  SSE2 optimization,  but also generate code executable on all mach.
# -tpp6 P3 optimization
# -tpp7 P4 optimization
#-----------------------------------------------------------------------

OFLAG=-O3 -xW

OFLAG_HIGH = $(OFLAG)
OBJ_HIGH =

OBJ_NOOPT =
DEBUG  = -FR -O0
INLINE = $(OFLAG)


#-----------------------------------------------------------------------
# the following lines specify the position of BLAS  and LAPACK
# on P4, VASP works fastest with the libgoto library
# so that's what I recommend
#-----------------------------------------------------------------------

# Atlas based libraries
#ATLASHOME= $(HOME)/archives/BLAS_OPT/ATLAS/lib/Linux_P4SSE2/
#BLAS=   -L$(ATLASHOME)  -lf77blas -latlas

# use specific libraries (default library path might point to other libraries)
#BLAS= $(ATLASHOME)/libf77blas.a $(ATLASHOME)/libatlas.a

# use the mkl Intel libraries for p4 (www.intel.com)
# mkl.5.1
# set -DRPROMU_DGEMV  -DRACCMU_DGEMV in the CPP lines
#BLAS=-L/opt/intel/mkl/lib/32 -lmkl_p4  -lpthread

# mkl.5.2 requires also to -lguide library
# set -DRPROMU_DGEMV  -DRACCMU_DGEMV in the CPP lines
BLAS=-L/opt/intel/composerxe-2011.3.174/mkl/lib/ia32 -lmkl_intel-lmkl_sequential -lmkl_core -lpthread

# even faster Kazushige Goto's BLAS
# http://www.cs.utexas.edu/users/kgoto/signup_first.html
#BLAS=  /opt/libs/libgoto/libgoto_p4_512-r0.6.so

# LAPACK, simplest use vasp.4.lib/lapack_double
#LAPACK= ../vasp.4.lib/lapack_double.o

# use atlas optimized part of lapack
#LAPACK= ../vasp.4.lib/lapack_atlas.o  -llapack -lcblas

# use the mkl Intel lapack
LAPACK=-L/opt/intel/composerxe-2011.3.174/mkl/lib/ia32 -lmkl_intel-lmkl_sequential -lmkl_core -lpthread

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

LIB  = -L../vasp.4.lib -ldmy \
     ../vasp.4.lib/linpack_double.o $(LAPACK) \
     $(BLAS)

# options for linking (for compiler version 6.X, 7.1) nothing is required
LINK    =
# compiler version 7.0 generates some vector statments which are located
# in the svml library, add the LIBPATH and the library (just in case)
#LINK    =  -L/opt/intel/compiler70/ia32/lib/ -lsvml

#-----------------------------------------------------------------------
# fft libraries:
# VASP.4.6 can use fftw.3.0.X (http://www.fftw.org)
# since this version is faster on P4 machines, we recommend to use it
#-----------------------------------------------------------------------

#FFT3D   = fft3dfurth.o fft3dlib.o
FFT3D   = fftmpiw.o fftmpi_map.o fft3dlib.o   /opt/intel/composerxe-2011.3.174/mkl/lib/ia32/libfftw3xf_intel.a


#=======================================================================
# MPI section, uncomment the following lines
#
# one comment for users of mpich or lam:
# You must *not* compile mpi with g77/f77, because f77/g77            
# appends *two* underscores to symbols that contain already an        
# underscore (i.e. MPI_SEND becomes mpi_send__).  The pgf90/ifc
# compilers however append only one underscore.
# Precompiled mpi version will also not work !!!
#
# We found that mpich.1.2.1 and lam-6.5.X to lam-7.0.4 are stable
# mpich.1.2.1 was configured with
#  ./configure -prefix=/usr/local/mpich_nodvdbg -fc="pgf77 -Mx,119,0x200000"  \
# -f90="pgf90 " \
# --without-romio --without-mpe -opt=-O \
#
# lam was configured with the line
#  ./configure  -prefix /opt/libs/lam-7.0.4 --with-cflags=-O -with-fc=ifc \
# --with-f77flags=-O --without-romio
#
# please note that you might be able to use a lam or mpich version
# compiled with f77/g77, but then you need to add the following
# options: -Msecond_underscore (compilation) and -g77libs (linking)
#
# !!! Please do not send me any queries on how to install MPI, I will
# certainly not answer them !!!!
#=======================================================================
#-----------------------------------------------------------------------
# fortran linker for mpi: if you use LAM and compiled it with the options
# suggested above,  you can use the following line
#-----------------------------------------------------------------------

FC=mpif90
FCL=$(FC)

#-----------------------------------------------------------------------
# additional options for CPP in parallel version (see also above):
# NGZhalf               charge density   reduced in Z direction
# wNGZhalf              gamma point only reduced in Z direction
# scaLAPACK             use scaLAPACK (usually slower on 100 Mbit Net)
#-----------------------------------------------------------------------

CPP    = $(CPP_) -DMPI  -DHOST=\"LinuxIFC\" -DIFC \
     -Dkind8 -DNGZhalf -DCACHE_SIZE=4000 -DPGF90 -Davoidalloc \
     -DMPI_BLOCK=500  \
##    -DRPROMU_DGEMV  -DRACCMU_DGEMV

#-----------------------------------------------------------------------
# location of SCALAPACK
# if you do not use SCALAPACK simply uncomment the line SCA
#-----------------------------------------------------------------------

BLACS=$(HOME)/archives/SCALAPACK/BLACS/
SCA_=$(HOME)/archives/SCALAPACK/SCALAPACK

SCA= $(SCA_)/libscalapack.a  \
$(BLACS)/LIB/blacsF77init_MPI-LINUX-0.a $(BLACS)/LIB/blacs_MPI-LINUX-0.a $(BLACS)/LIB/blacsF77init_MPI-LINUX-0.a

SCA=/opt/intel/composerxe-2011.3.174/mkl/lib/ia32/libmkl_scalapack_core.a /opt/intel/composerxe-2011.3.174/mkl/lib/ia32/libmkl_blacs_intelmpi.a

#-----------------------------------------------------------------------
# libraries for mpi
#-----------------------------------------------------------------------

#LIB     = -L../vasp.4.lib -ldmy  \
#      ../vasp.4.lib/linpack_double.o $(LAPACK) \
#      $(SCA) $(BLAS)

# FFT: fftmpi.o with fft3dlib of Juergen Furthmueller
FFT3D   = FFT3D   = fftmpiw.o fftmpi_map.o fft3dlib.o   /opt/intel/composerxe-2011.3.174/mkl/lib/ia32/libfftw3xf_intel.a

# fftw.3.0.1 is slighly faster and should be used if available
#FFT3D   = fftmpiw.o fftmpi_map.o fft3dlib.o   /opt/libs/fftw-3.0.1/lib/libfftw3.a

#-----------------------------------------------------------------------
# general rules and compile lines
#-----------------------------------------------------------------------
BASIC=   symmetry.o symlib.o   lattlib.o  random.o   

SOURCE=  base.o     mpi.o      smart_allocate.o      xml.o  \
         constant.o jacobi.o   main_mpi.o  scala.o   \
         asa.o      lattice.o  poscar.o   ini.o      setex.o     radial.o  \
         pseudo.o   mgrid.o    mkpoints.o wave.o      wave_mpi.o  $(BASIC) \
         nonl.o     nonlr.o    dfast.o    choleski2.o    \
         mix.o      charge.o   xcgrad.o   xcspin.o    potex1.o   potex2.o  \
         metagga.o  constrmag.o pot.o      cl_shift.o force.o    dos.o      elf.o      \
         tet.o      hamil.o    steep.o    \
         chain.o    dyna.o     relativistic.o LDApU.o sphpro.o  paw.o   us.o \
         ebs.o      wavpre.o   wavpre_noio.o broyden.o \
         dynbr.o    rmm-diis.o reader.o   writer.o   tutor.o xml_writer.o \
         brent.o    stufak.o   fileio.o   opergrid.o stepver.o  \
         dipol.o    xclib.o    chgloc.o   subrot.o   optreal.o   davidson.o \
         edtest.o   electron.o shm.o      pardens.o  paircorrection.o \
         optics.o   constr_cell_relax.o   stm.o    finite_diff.o \
         elpol.o    setlocalpp.o

INC=

vasp: $(SOURCE) $(FFT3D) $(INC) main.o
        rm -f vasp
        $(FCL) -o vasp $(LINK) main.o  $(SOURCE)   $(FFT3D) $(LIB)
makeparam: $(SOURCE) $(FFT3D) makeparam.o main.F $(INC)
        $(FCL) -o makeparam  $(LINK) makeparam.o $(SOURCE) $(FFT3D) $(LIB)
zgemmtest: zgemmtest.o base.o random.o $(INC)
        $(FCL) -o zgemmtest $(LINK) zgemmtest.o random.o base.o $(LIB)
dgemmtest: dgemmtest.o base.o random.o $(INC)
        $(FCL) -o dgemmtest $(LINK) dgemmtest.o random.o base.o $(LIB)
ffttest: base.o smart_allocate.o mpi.o mgrid.o random.o ffttest.o $(FFT3D) $(INC)
        $(FCL) -o ffttest $(LINK) ffttest.o mpi.o mgrid.o random.o smart_allocate.o base.o $(FFT3D) $(LIB)
kpoints: $(SOURCE) $(FFT3D) makekpoints.o main.F $(INC)
        $(FCL) -o kpoints $(LINK) makekpoints.o $(SOURCE) $(FFT3D) $(LIB)

clean:       
        -rm -f *.g *.f *.o *.L *.mod ; touch *.F

main.o: main$(SUFFIX)
        $(FC) $(FFLAGS)$(DEBUG)  $(INCS) -c main$(SUFFIX)
xcgrad.o: xcgrad$(SUFFIX)
        $(FC) $(FFLAGS) $(INLINE)  $(INCS) -c xcgrad$(SUFFIX)
xcspin.o: xcspin$(SUFFIX)
        $(FC) $(FFLAGS) $(INLINE)  $(INCS) -c xcspin$(SUFFIX)

makeparam.o: makeparam$(SUFFIX)
        $(FC) $(FFLAGS)$(DEBUG)  $(INCS) -c makeparam$(SUFFIX)

makeparam$(SUFFIX): makeparam.F main.F
#
# MIND: I do not have a full dependency list for the include
# and MODULES: here are only the minimal basic dependencies
# if one strucuture is changed then touch_dep must be called
# with the corresponding name of the structure
#
base.o: base.inc base.F
mgrid.o: mgrid.inc mgrid.F
constant.o: constant.inc constant.F
lattice.o: lattice.inc lattice.F
setex.o: setexm.inc setex.F
pseudo.o: pseudo.inc pseudo.F
poscar.o: poscar.inc poscar.F
mkpoints.o: mkpoints.inc mkpoints.F
wave.o: wave.inc wave.F
nonl.o: nonl.inc nonl.F
nonlr.o: nonlr.inc nonlr.F

$(OBJ_HIGH):
        $(CPP)
        $(FC) $(FFLAGS) $(OFLAG_HIGH) $(INCS) -c $*$(SUFFIX)
$(OBJ_NOOPT):
        $(CPP)
        $(FC) $(FFLAGS) $(INCS) -c $*$(SUFFIX)

fft3dlib_f77.o: fft3dlib_f77.F
        $(CPP)
        $(F77) $(FFLAGS_F77) -c $*$(SUFFIX)

.F.o:
        $(CPP)
        $(FC) $(FFLAGS) $(OFLAG) $(INCS) -c $*$(SUFFIX)
.F$(SUFFIX):
        $(CPP)
$(SUFFIX).o:
        $(FC) $(FFLAGS) $(OFLAG) $(INCS) -c $*$(SUFFIX)

# special rules
#-----------------------------------------------------------------------
# these special rules are cummulative (that is once failed
#   in one compiler version, stays in the list forever)
# -tpp5|6|7 P, PII-PIII, PIV
# -xW use SIMD (does not pay of on PII, since fft3d uses double prec)
# all other options do no affect the code performance since -O1 is used
#-----------------------------------------------------------------------

fft3dlib.o : fft3dlib.F
        $(CPP)
        $(FC) -FR -lowercase -O1 -tpp7 -xW -prefetch- -unroll0 -vec_report3 -c $*$(SUFFIX)
fft3dfurth.o : fft3dfurth.F
        $(CPP)
        $(FC) -FR -lowercase -O1 -c $*$(SUFFIX)

radial.o : radial.F
        $(CPP)
        $(FC) -FR -lowercase -O1 -c $*$(SUFFIX)

symlib.o : symlib.F
        $(CPP)
        $(FC) -FR -lowercase -O1 -c $*$(SUFFIX)

symmetry.o : symmetry.F
        $(CPP)
        $(FC) -FR -lowercase -O1 -c $*$(SUFFIX)

dynbr.o : dynbr.F
        $(CPP)
        $(FC) -FR -lowercase -O1 -c $*$(SUFFIX)

broyden.o : broyden.F
        $(CPP)
        $(FC) -FR -lowercase -O2 -c $*$(SUFFIX)

us.o : us.F
        $(CPP)
        $(FC) -FR -lowercase -O1 -c $*$(SUFFIX)

wave.o : wave.F
        $(CPP)
        $(FC) -FR -lowercase -O0 -c $*$(SUFFIX)

LDApU.o : LDApU.F
        $(CPP)
        $(FC) -FR -lowercase -O2 -c $*$(SUFFIX)
4樓2011-05-02 09:31:56
已閱   回復(fù)此樓   關(guān)注TA 給TA發(fā)消息 送TA紅花 TA的回帖

enola

捐助貴賓 (正式寫手)

【答案】應(yīng)助回帖

mengfc(金幣+5): 謝謝參與 2011-05-03 08:16:20
FFT3D   = FFT3D   = fftmpiw.o fftmpi_map.o fft3dlib.o   /opt/intel/composerxe-2011.3.174/mkl/lib/ia32/libfftw3xf_intel.a
是不是這里有錯(cuò)啊。。。寫了兩個(gè)FFT3D啊
timeflies..
5樓2011-05-02 11:30:37
已閱   回復(fù)此樓   關(guān)注TA 給TA發(fā)消息 送TA紅花 TA的回帖

enola

捐助貴賓 (正式寫手)

SCA=
#BLACS=
#SCA_=
不用加入的。。。
timeflies..
6樓2011-05-02 11:31:55
已閱   回復(fù)此樓   關(guān)注TA 給TA發(fā)消息 送TA紅花 TA的回帖

賀儀

鐵桿木蟲 (著名寫手)

【答案】應(yīng)助回帖

mengfc(金幣+5): 2011-05-02 15:19:56
# even faster Kazushige Goto's BLAS
# http://www.cs.utexas.edu/users/kgoto/signup_first.html
#BLAS=  /opt/libs/libgoto/libgoto_p4_512-r0.6.so

# LAPACK, simplest use vasp.4.lib/lapack_double
#LAPACK= ../vasp.4.lib/lapack_double.o

# use atlas optimized part of lapack
#LAPACK= ../vasp.4.lib/lapack_atlas.o  -llapack -lcblas

# use the mkl Intel lapack
LAPACK=-L/opt/intel/composerxe-2011.3.174/mkl/lib/ia32 -lmkl_intel-lmkl_sequential -lmkl_core -lpthread

還是用自帶的吧:
# even faster Kazushige Goto's BLAS
# http://www.cs.utexas.edu/users/kgoto/signup_first.html
#BLAS=  /opt/libs/libgoto/libgoto_p4_512-r0.6.so

# LAPACK, simplest use vasp.4.lib/lapack_double
LAPACK= ../vasp.4.lib/lapack_double.o

# use atlas optimized part of lapack
#LAPACK= ../vasp.4.lib/lapack_atlas.o  -llapack -lcblas

# use the mkl Intel lapack
#LAPACK=-L/opt/intel/composerxe-2011.3.174/mkl/lib/ia32 -lmkl_intel-lmkl_sequential #-lmkl_core -lpthread
嘟啦啦嘟啦啦嘟
7樓2011-05-02 12:03:42
已閱   回復(fù)此樓   關(guān)注TA 給TA發(fā)消息 送TA紅花 TA的回帖

小木蟲于

金蟲 (正式寫手)

【答案】應(yīng)助回帖

mengfc(金幣+5): 謝謝回復(fù),請(qǐng)看十樓的makefile有沒有問題? 2011-05-02 16:01:52
intel fortran compiler 12 編譯能通過么?我試了不行啊,不知哪里錯(cuò)了。
8樓2011-05-02 13:32:46
已閱   回復(fù)此樓   關(guān)注TA 給TA發(fā)消息 送TA紅花 TA的回帖

mengfc

金蟲 (正式寫手)

引用回帖:
Originally posted by 賀儀 at 2011-05-02 12:03:42:
# even faster Kazushige Goto's BLAS
# http://www.cs.utexas.edu/users/kgoto/signup_first.html
#BLAS=  /opt/libs/libgoto/libgoto_p4_512-r0.6.so

# LAPACK, simplest use vasp.4.lib/lapac ...

謝謝你,越裝越?jīng)]信心了。。。。
/opt/intel/composerxe-2011.3.174/mkl/lib/ia32/libmkl_intel.so: undefined reference to `mkl_lapack_dgtcon'
/opt/intel/composerxe-2011.3.174/mkl/lib/ia32/libmkl_intel.so: undefined reference to `mkl_lapack_zcposv'
/opt/intel/composerxe-2011.3.174/mkl/lib/ia32/libmkl_intel.so: undefined reference to `mkl_lapack_ssyevd'
/opt/intel/composerxe-2011.3.174/mkl/lib/ia32/libmkl_intel.so: undefined reference to `mkl_spblas_mkl_cspblas_ccoosymv'
/opt/intel/composerxe-2011.3.174/mkl/lib/ia32/libmkl_intel.so: undefined reference to `mkl_pds_dss_factor_complex_d'
make: *** [vasp] Error 1

這個(gè)錯(cuò)誤是什么原因呢?
9樓2011-05-02 15:20:32
已閱   回復(fù)此樓   關(guān)注TA 給TA發(fā)消息 送TA紅花 TA的回帖

mengfc

金蟲 (正式寫手)

引用回帖:
Originally posted by enola at 2011-05-02 11:31:55:
SCA=
#BLACS=
#SCA_=
不用加入的。。。

[mengfc@localhost vasp.4.6]$ ./vasp
Error reading item 'IMAGES' from file INCAR.

這樣算編譯成功了嗎?

另附我的makefile,請(qǐng)幫忙看看可以嗎?
.SUFFIXES: .inc .f .f90 .F
#-----------------------------------------------------------------------
# Makefile for Intel Fortran compiler for P4 systems
#
# The makefile was tested only under Linux on Intel platforms
# (Suse 5.3- Suse 9.0)
# the followin compiler versions have been tested
# 5.0, 6.0, 7.0 and 7.1 (some 8.0 versions seem to fail compiling the code)
# presently we recommend version 7.1 or 7.0, since these
# releases have been used to compile the present code versions
#
# it might be required to change some of library pathes, since
# LINUX installation vary a lot
# Hence check ***ALL**** options in this makefile very carefully
#-----------------------------------------------------------------------
#
# BLAS must be installed on the machine
# there are several options:
# 1) very slow but works:
#   retrieve the lapackage from ftp.netlib.org
#   and compile the blas routines (BLAS/SRC directory)
#   please use g77 or f77 for the compilation. When I tried to
#   use pgf77 or pgf90 for BLAS, VASP hang up when calling
#   ZHEEV  (however this was with lapack 1.1 now I use lapack 2.0)
# 2) most desirable: get an optimized BLAS
#
# the two most reliable packages around are presently:
# 3a) Intels own optimised BLAS (PIII, P4, Itanium)
#     http://developer.intel.com/software/products/mkl/
#   this is really excellent when you use Intel CPU's
#
# 3b) or obtain the atlas based BLAS routines
#     http://math-atlas.sourceforge.net/
#   you certainly need atlas on the Athlon, since the  mkl
#   routines are not optimal on the Athlon.
#   If you want to use atlas based BLAS, check the lines around LIB=
#
# 3c) mindblowing fast SSE2 (4 GFlops on P4, 2.53 GHz)
#   Kazushige Goto's BLAS
#   http://www.cs.utexas.edu/users/kgoto/signup_first.html
#
#-----------------------------------------------------------------------

# all CPP processed fortran files have the extension .f90
SUFFIX=.f90

#-----------------------------------------------------------------------
# fortran compiler and linker
#-----------------------------------------------------------------------
FC=ifort
# fortran linker
FCL=$(FC)


#-----------------------------------------------------------------------
# whereis CPP ?? (I need CPP, can't use gcc with proper options)
# that's the location of gcc for SUSE 5.3
#
#  CPP_   =  /usr/lib/gcc-lib/i486-linux/2.7.2/cpp -P -C
#
# that's probably the right line for some Red Hat distribution:
#
#  CPP_   =  /usr/lib/gcc-lib/i386-redhat-linux/2.7.2.3/cpp -P -C
#
#  SUSE X.X, maybe some Red Hat distributions:

CPP_ =  ./preprocess <$*.F | /usr/bin/cpp -P -C -traditional >$*$(SUFFIX)

#-----------------------------------------------------------------------
# possible options for CPP:
# NGXhalf             charge density   reduced in X direction
# wNGXhalf            gamma point only reduced in X direction
# avoidalloc          avoid ALLOCATE if possible
# IFC                 work around some IFC bugs
# CACHE_SIZE          1000 for PII,PIII, 5000 for Athlon, 8000-12000 P4
# RPROMU_DGEMV        use DGEMV instead of DGEMM in RPRO (depends on used BLAS)
# RACCMU_DGEMV        use DGEMV instead of DGEMM in RACC (depends on used BLAS)
#-----------------------------------------------------------------------

CPP     = $(CPP_)  -DHOST=\"LinuxIFC\" \
          -Dkind8 -DNGXhalf -DCACHE_SIZE=12000 -DPGF90 -Davoidalloc \
          -DRPROMU_DGEMV  -DRACCMU_DGEMV

#-----------------------------------------------------------------------
# general fortran flags  (there must a trailing blank on this line)
#-----------------------------------------------------------------------

FFLAGS =  -FR -lowercase  -heap-arrays  64 -assume byterecl

#-----------------------------------------------------------------------
# optimization
# we have tested whether higher optimisation improves performance
# -axK  SSE1 optimization,  but also generate code executable on all mach.
#       xK improves performance somewhat on XP, and a is required in order
#       to run the code on older Athlons as well
# -xW   SSE2 optimization
# -axW  SSE2 optimization,  but also generate code executable on all mach.
# -tpp6 P3 optimization
# -tpp7 P4 optimization
#-----------------------------------------------------------------------

OFLAG=-O3 -xW

OFLAG_HIGH = $(OFLAG)
OBJ_HIGH =

OBJ_NOOPT =
DEBUG  = -FR -O0
INLINE = $(OFLAG)


#-----------------------------------------------------------------------
# the following lines specify the position of BLAS  and LAPACK
# on P4, VASP works fastest with the libgoto library
# so that's what I recommend
#-----------------------------------------------------------------------

# Atlas based libraries
#ATLASHOME= $(HOME)/archives/BLAS_OPT/ATLAS/lib/Linux_P4SSE2/
#BLAS=   -L$(ATLASHOME)  -lf77blas -latlas

# use specific libraries (default library path might point to other libraries)
#BLAS= $(ATLASHOME)/libf77blas.a $(ATLASHOME)/libatlas.a

# use the mkl Intel libraries for p4 (www.intel.com)
# mkl.5.1
# set -DRPROMU_DGEMV  -DRACCMU_DGEMV in the CPP lines
#BLAS=-L/opt/intel/mkl/lib/32 -lmkl_p4  -lpthread

# mkl.5.2 requires also to -lguide library
# set -DRPROMU_DGEMV  -DRACCMU_DGEMV in the CPP lines
BLAS=-L/opt/intel/composerxe-2011.3.174/mkl/lib/ia32 -lmkl_intel -lmkl_sequential -lmkl_core -lpthread

# even faster Kazushige Goto's BLAS
# http://www.cs.utexas.edu/users/kgoto/signup_first.html
#BLAS=  /opt/libs/libgoto/libgoto_p4_512-r0.6.so

# LAPACK, simplest use vasp.4.lib/lapack_double
LAPACK= ../vasp.4.lib/lapack_double.o

# use atlas optimized part of lapack
#LAPACK= ../vasp.4.lib/lapack_atlas.o  -llapack -lcblas

# use the mkl Intel lapack
#LAPACK= -lmkl_lapack

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

LIB  = -L../vasp.4.lib -ldmy \
     ../vasp.4.lib/linpack_double.o $(LAPACK) \
     $(BLAS)

# options for linking (for compiler version 6.X, 7.1) nothing is required
LINK    =
# compiler version 7.0 generates some vector statments which are located
# in the svml library, add the LIBPATH and the library (just in case)
#LINK    =  -L/opt/intel/compiler70/ia32/lib/ -lsvml

#-----------------------------------------------------------------------
# fft libraries:
# VASP.4.6 can use fftw.3.0.X (http://www.fftw.org)
# since this version is faster on P4 machines, we recommend to use it
#-----------------------------------------------------------------------

#FFT3D   = fft3dfurth.o fft3dlib.o
#FFT3D   = fftw3d.o fft3dlib.o   /opt/libs/fftw-3.0.1/lib/libfftw3.a


#=======================================================================
# MPI section, uncomment the following lines
#
# one comment for users of mpich or lam:
# You must *not* compile mpi with g77/f77, because f77/g77            
# appends *two* underscores to symbols that contain already an        
# underscore (i.e. MPI_SEND becomes mpi_send__).  The pgf90/ifc
# compilers however append only one underscore.
# Precompiled mpi version will also not work !!!
#
# We found that mpich.1.2.1 and lam-6.5.X to lam-7.0.4 are stable
# mpich.1.2.1 was configured with
#  ./configure -prefix=/usr/local/mpich_nodvdbg -fc="pgf77 -Mx,119,0x200000"  \
# -f90="pgf90 " \
# --without-romio --without-mpe -opt=-O \
#
# lam was configured with the line
#  ./configure  -prefix /opt/libs/lam-7.0.4 --with-cflags=-O -with-fc=ifc \
# --with-f77flags=-O --without-romio
#
# please note that you might be able to use a lam or mpich version
# compiled with f77/g77, but then you need to add the following
# options: -Msecond_underscore (compilation) and -g77libs (linking)
#
# !!! Please do not send me any queries on how to install MPI, I will
# certainly not answer them !!!!
#=======================================================================
#-----------------------------------------------------------------------
# fortran linker for mpi: if you use LAM and compiled it with the options
# suggested above,  you can use the following line
#-----------------------------------------------------------------------

FC=mpif90
FCL=$(FC)

#-----------------------------------------------------------------------
# additional options for CPP in parallel version (see also above):
# NGZhalf               charge density   reduced in Z direction
# wNGZhalf              gamma point only reduced in Z direction
# scaLAPACK             use scaLAPACK (usually slower on 100 Mbit Net)
#-----------------------------------------------------------------------

CPP    = $(CPP_) -DMPI  -DHOST=\"LinuxIFC\" -DIFC \
     -Dkind8 -DNGZhalf -DCACHE_SIZE=4000 -DPGF90 -Davoidalloc \
     -DMPI_BLOCK=500  \
    -DRPROMU_DGEMV  -DRACCMU_DGEMV

#-----------------------------------------------------------------------
# location of SCALAPACK
# if you do not use SCALAPACK simply uncomment the line SCA
#-----------------------------------------------------------------------

BLACS=$(HOME)/archives/SCALAPACK/BLACS/
SCA_=$(HOME)/archives/SCALAPACK/SCALAPACK

SCA= $(SCA_)/libscalapack.a  \
$(BLACS)/LIB/blacsF77init_MPI-LINUX-0.a $(BLACS)/LIB/blacs_MPI-LINUX-0.a $(BLACS)/LIB/blacsF77init_MPI-LINUX-0.a

SCA=/opt/intel/composerxe-2011.3.174/mkl/lib/ia32/libmkl_scalapack_core.a  /opt/intel/composerxe-2011.3.174/mkl/lib/ia32/libmkl_blacs_intelmpi.a

#-----------------------------------------------------------------------
# libraries for mpi
#-----------------------------------------------------------------------

#LIB     = -L../vasp.4.lib -ldmy  \
#      ../vasp.4.lib/linpack_double.o $(LAPACK) \
#      $(SCA) $(BLAS)

# FFT: fftmpi.o with fft3dlib of Juergen Furthmueller
FFT3D   = fftmpi.o fftmpi_map.o fft3dlib.o

# fftw.3.0.1 is slighly faster and should be used if available
#FFT3D   = fftmpiw.o fftmpi_map.o fft3dlib.o  /opt/intel/composerxe-2011.3.174/mkl/interfaces/fftw3xf/libfftw3xf_intel.a  


#-----------------------------------------------------------------------
# general rules and compile lines
#-----------------------------------------------------------------------
BASIC=   symmetry.o symlib.o   lattlib.o  random.o   

SOURCE=  base.o     mpi.o      smart_allocate.o      xml.o  \
         constant.o jacobi.o   main_mpi.o  scala.o   \
         asa.o      lattice.o  poscar.o   ini.o      setex.o     radial.o  \
         pseudo.o   mgrid.o    mkpoints.o wave.o      wave_mpi.o  $(BASIC) \
         nonl.o     nonlr.o    dfast.o    choleski2.o    \
         mix.o      charge.o   xcgrad.o   xcspin.o    potex1.o   potex2.o  \
         metagga.o  constrmag.o pot.o      cl_shift.o force.o    dos.o      elf.o      \
         tet.o      hamil.o    steep.o    \
         chain.o    dyna.o     relativistic.o LDApU.o sphpro.o  paw.o   us.o \
         ebs.o      wavpre.o   wavpre_noio.o broyden.o \
         dynbr.o    rmm-diis.o reader.o   writer.o   tutor.o xml_writer.o \
         brent.o    stufak.o   fileio.o   opergrid.o stepver.o  \
         dipol.o    xclib.o    chgloc.o   subrot.o   optreal.o   davidson.o \
         edtest.o   electron.o shm.o      pardens.o  paircorrection.o \
         optics.o   constr_cell_relax.o   stm.o    finite_diff.o \
         elpol.o    setlocalpp.o

INC=

vasp: $(SOURCE) $(FFT3D) $(INC) main.o
        rm -f vasp
        $(FCL) -o vasp $(LINK) main.o  $(SOURCE)   $(FFT3D) $(LIB)
makeparam: $(SOURCE) $(FFT3D) makeparam.o main.F $(INC)
        $(FCL) -o makeparam  $(LINK) makeparam.o $(SOURCE) $(FFT3D) $(LIB)
zgemmtest: zgemmtest.o base.o random.o $(INC)
        $(FCL) -o zgemmtest $(LINK) zgemmtest.o random.o base.o $(LIB)
dgemmtest: dgemmtest.o base.o random.o $(INC)
        $(FCL) -o dgemmtest $(LINK) dgemmtest.o random.o base.o $(LIB)
ffttest: base.o smart_allocate.o mpi.o mgrid.o random.o ffttest.o $(FFT3D) $(INC)
        $(FCL) -o ffttest $(LINK) ffttest.o mpi.o mgrid.o random.o smart_allocate.o base.o $(FFT3D) $(LIB)
kpoints: $(SOURCE) $(FFT3D) makekpoints.o main.F $(INC)
        $(FCL) -o kpoints $(LINK) makekpoints.o $(SOURCE) $(FFT3D) $(LIB)

clean:       
        -rm -f *.g *.f *.o *.L *.mod ; touch *.F

main.o: main$(SUFFIX)
        $(FC) $(FFLAGS)$(DEBUG)  $(INCS) -c main$(SUFFIX)
xcgrad.o: xcgrad$(SUFFIX)
        $(FC) $(FFLAGS) $(INLINE)  $(INCS) -c xcgrad$(SUFFIX)
xcspin.o: xcspin$(SUFFIX)
        $(FC) $(FFLAGS) $(INLINE)  $(INCS) -c xcspin$(SUFFIX)

makeparam.o: makeparam$(SUFFIX)
        $(FC) $(FFLAGS)$(DEBUG)  $(INCS) -c makeparam$(SUFFIX)

makeparam$(SUFFIX): makeparam.F main.F
#
# MIND: I do not have a full dependency list for the include
# and MODULES: here are only the minimal basic dependencies
# if one strucuture is changed then touch_dep must be called
# with the corresponding name of the structure
#
base.o: base.inc base.F
mgrid.o: mgrid.inc mgrid.F
constant.o: constant.inc constant.F
lattice.o: lattice.inc lattice.F
setex.o: setexm.inc setex.F
pseudo.o: pseudo.inc pseudo.F
poscar.o: poscar.inc poscar.F
mkpoints.o: mkpoints.inc mkpoints.F
wave.o: wave.inc wave.F
nonl.o: nonl.inc nonl.F
nonlr.o: nonlr.inc nonlr.F

$(OBJ_HIGH):
        $(CPP)
        $(FC) $(FFLAGS) $(OFLAG_HIGH) $(INCS) -c $*$(SUFFIX)
$(OBJ_NOOPT):
        $(CPP)
        $(FC) $(FFLAGS) $(INCS) -c $*$(SUFFIX)

fft3dlib_f77.o: fft3dlib_f77.F
        $(CPP)
        $(F77) $(FFLAGS_F77) -c $*$(SUFFIX)

.F.o:
        $(CPP)
        $(FC) $(FFLAGS) $(OFLAG) $(INCS) -c $*$(SUFFIX)
.F$(SUFFIX):
        $(CPP)
$(SUFFIX).o:
        $(FC) $(FFLAGS) $(OFLAG) $(INCS) -c $*$(SUFFIX)

# special rules
#-----------------------------------------------------------------------
# these special rules are cummulative (that is once failed
#   in one compiler version, stays in the list forever)
# -tpp5|6|7 P, PII-PIII, PIV
# -xW use SIMD (does not pay of on PII, since fft3d uses double prec)
# all other options do no affect the code performance since -O1 is used
#-----------------------------------------------------------------------

fft3dlib.o : fft3dlib.F
        $(CPP)
        $(FC) -FR -lowercase -O1 -xW -unroll0 -vec_report3 -c $*$(SUFFIX)
fft3dfurth.o : fft3dfurth.F
        $(CPP)
        $(FC) -FR -lowercase -O1 -c $*$(SUFFIX)

radial.o : radial.F
        $(CPP)
        $(FC) -FR -lowercase -O1 -c $*$(SUFFIX)

symlib.o : symlib.F
        $(CPP)
        $(FC) -FR -lowercase -O1 -c $*$(SUFFIX)

symmetry.o : symmetry.F
        $(CPP)
        $(FC) -FR -lowercase -O1 -c $*$(SUFFIX)

dynbr.o : dynbr.F
        $(CPP)
        $(FC) -FR -lowercase -O1 -c $*$(SUFFIX)

broyden.o : broyden.F
        $(CPP)
        $(FC) -FR -lowercase -O2 -c $*$(SUFFIX)

us.o : us.F
        $(CPP)
        $(FC) -FR -lowercase -O1 -c $*$(SUFFIX)

wave.o : wave.F
        $(CPP)
        $(FC) -FR -lowercase -O0 -c $*$(SUFFIX)

LDApU.o : LDApU.F
        $(CPP)
        $(FC) -FR -lowercase -O2 -c $*$(SUFFIX)
10樓2011-05-02 16:00:49
已閱   回復(fù)此樓   關(guān)注TA 給TA發(fā)消息 送TA紅花 TA的回帖
相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主 mengfc 的主題更新
最具人氣熱帖推薦 [查看全部] 作者 回/看 最后發(fā)表
[考研] 271材料工程求調(diào)劑 +5 .6lL 2026-03-18 5/250 2026-03-19 03:07 by 無懈可擊111
[考研] 304求調(diào)劑 +6 司空. 2026-03-18 6/300 2026-03-18 23:03 by 星空星月
[考研] 材料專業(yè)求調(diào)劑 +5 hanamiko 2026-03-18 5/250 2026-03-18 20:19 by 星空星月
[考研] 26調(diào)劑/材料/英一數(shù)二/總分289/已過A區(qū)線 +7 步川酷紫123 2026-03-13 7/350 2026-03-18 17:12 by 盡舜堯1
[考研] 286求調(diào)劑 +6 lemonzzn 2026-03-16 9/450 2026-03-18 15:31 by vgtyfty
[考研] 314求調(diào)劑 +8 無懈可擊的巨人 2026-03-12 8/400 2026-03-18 14:50 by haxia
[考研] 一志愿西南交大,求調(diào)劑 +4 材化逐夢(mèng)人 2026-03-18 4/200 2026-03-18 14:22 by 007_lilei
[考研] 收復(fù)試調(diào)劑生 +4 雨后秋荷 2026-03-18 4/200 2026-03-18 14:16 by elevennnne
[考研] 0703化學(xué)調(diào)劑 +4 pupcoco 2026-03-17 7/350 2026-03-18 12:14 by djl2006
[考研] 生物學(xué)071000 329分求調(diào)劑 +3 我愛生物生物愛?/a> 2026-03-17 3/150 2026-03-18 10:12 by macy2011
[考研] 278求調(diào)劑 +5 煙火先于春 2026-03-17 5/250 2026-03-18 08:43 by 星空星月
[考研] 334求調(diào)劑 +3 志存高遠(yuǎn)意在機(jī)?/a> 2026-03-16 3/150 2026-03-18 08:34 by lm4875102
[考研] 268求調(diào)劑 +7 好運(yùn)連綿不絕 2026-03-12 8/400 2026-03-17 20:28 by xilongliang
[考研] 294求調(diào)劑 +3 Zys010410@ 2026-03-13 4/200 2026-03-15 10:59 by zhq0425
[考研] 288求調(diào)劑 +4 奇點(diǎn)0314 2026-03-14 4/200 2026-03-14 23:04 by JourneyLucky
[考研] 材料與化工 323 英一+數(shù)二+物化,一志愿:哈工大 本人本科雙一流 +4 自由的_飛翔 2026-03-13 5/250 2026-03-14 19:39 by hmn_wj
[考研] 材料工程調(diào)劑 +9 咪咪空空 2026-03-12 9/450 2026-03-13 22:05 by 星空星月
[碩博家園] 085600 260分求調(diào)劑 +3 天空還下雨么 2026-03-13 5/250 2026-03-13 18:46 by 天空還下雨么
[考研] 材料專碩350 求調(diào)劑 +4 王金科 2026-03-12 4/200 2026-03-13 16:02 by ruiyingmiao
[考研] 土木第一志愿276求調(diào)劑,科研和技能十分豐富,求新興方向的導(dǎo)師收留 +3 土木小天才 2026-03-12 3/150 2026-03-13 15:01 by JourneyLucky
信息提示
請(qǐng)?zhí)钐幚硪庖?/div>
www,日韩av,com| 亚洲男人的天堂最新网址| www,日韩av,com| 国产成人在线观看视频播放| 在线视频国产精品欧美| 抽插小穴啊啊啊视频| 日本少妇人妻中文在线| xxnxx国产美女| 国产igao激情在线视频入口 | 中文字幕精品人妻久久久久| 亚洲一区二区在线激情| 久久午夜免费鲁丝片| 精品不卡一区二区三区| 伊人综合在线视频免费观看| 欧美黄色性视频网站| 中文字幕丰满子伦无码专区 | 人妻系列中文字幕大乳丰满人妻 | 久久视频 在线播放| 亚洲春色av中文字幕| 午夜精品久久久久久久久久蜜桃| 青青青青青爽视频在线| 女人的天堂 av在线| 九热精品视频在线观看| 91九色pony蝌蚪| 松本菜奈实最新av在线| 亚洲成人偷拍自拍在线| 欧美激情视频第一页| 4438x亚洲最大的成人| 老色鬼精品视频在线观看播放| 九九热视频1这里只有精品| 操死你美女在线视频| 天天插天天干天天狠| 福利小视频免费在线| 免费在线观看视频啪啪| 女人的天堂 av在线| 精品不卡一区二区三区| 十八禁黄色免费污污污亚洲| 亚洲天堂色综合久久| 999精品视频免费在线观看| av无限看熟女人妻另类av| 嗯~嗯~啊啊啊~高潮了软件| 久久人妻诱惑我视频| 亚洲成年人精品国产| 亚洲无人区乱码中文字幕一区| 在线视频国产精品欧美| 亚洲欧美精品日韩偷拍| 女人高潮潮呻吟喷水网站| 国产福利一区二区三区在线观看| 国产人妻熟女ⅹxx丝袜| 欧美啪啪一区二区三区| 两个奶被揉得又硬又翘怎么回事 | 中文字幕麻绳捆绑的人妻| 核xp工厂精品久久亚洲| 一区二区三区内射美女| 性高潮视频在线观看日韩| 夜夜操夜夜爱夜夜摸| 国产视频成人一区二区| 日本东京热视频欧美视频| 欧美操大黑鸡巴视频在线观看| 日本韩国欧美在线视频| 黑人爆操女人免费视频| 婷婷综合缴情亚洲五月伊人| 香港日本台湾经典三级| 天天日天天干天天日天天干天天| 中文字幕熟女人妻丝袜丝在线| 乌克兰美女操逼高清内射视频| 日韩成人精品久久久免费看| 天天干夜夜撸天天操| 中文字幕丰满子伦无码专区| 人妻系列级片在线观看视频| 涩涩黄片在线免费观看| 男人的天堂在线2025| av毛片在线观看网址| 91色老久久精品偷偷蜜臀| 99在线视频精品观看高| 大尺度久久久久久久| avtt中文字幕手机版| 国产清纯一区二区在线观看 | 国产午夜在线播放视频| 亚洲欧美小说中文字幕| 91精品国产欧美在线| 一区二区三区观看在线| 亚洲欧美日韩中文视频| 午夜92福利1000| 在线观看网站伊人网| 老熟女 露脸 嗷嗷叫| 四季av人妻一区二区三区| 国产激情在线观看一区二区三区| 999国产精品视频免费看| 正在播放麻豆精品一区二区| 天天操天天舔天天做| 大尺度av毛片在线网址| 国产黑色丝袜 在线日韩欧美| 制服丝袜中文字幕熟女人妻| 久久视频 在线播放| av中文字幕国产精品| 婷婷一区二区三区五月丁| 放荡人妻极品少妇全集| 五月婷婷伊人久久中文字幕| 亚洲欧美国产人成在线| 久久久久久a女人处女| 精品免费一区二区三区四区视频| 夫亡人妻被强干中文字幕| 亚洲熟女乱一区二区精品成人| 最新日韩av电影在线播放| 精品不卡一区二区三区| 老司国产精品视频免费观看| 大香蕉在线欧美在线视频 | 老司机在线视频福利观看| 另类欧美激情校园春色| 国内自拍第一区二区三区| 中文字幕熟女乱一区二区| 天天干天天色综合久久| 快色视频在线观看免费| tushy一区二区三区视频| 欧美日韩成人高清中文网| 人妻视频网站快射视频网站| 青青操天堂在线观看视频| 超碰在线观看97资源| 欧美色区国产日韩亚洲区| 在线观看中文字幕少妇av| 人妻熟女 亚洲 一页二页| 久操资源在线免费播放| 日本成人福利电影网| 最新福利二区三区视频| 狂操鸡巴小骚逼视频免费观看| 国产最新av在线免费观看| 色网站在线观看免费| 日本香港韩国三级黄色| 中文字幕中文字幕在线中…一区| 99久久人人爽亚洲精品美女| 天天操天天搞天天操| 亚洲熟妇丰满多毛xxxx网站| 亚洲成人自拍图片网站| 高清国产美女a一级毛片| 国产精品剧情在线亚洲| 蜜桃臀av在线一区二区| 亚洲蜜桃久久久久久| 伊人精品成人综合网| 一区二区九日韩美女| av男人站在线观看| 成人av在线视频免费| 99女福利女女视频在线播放| 亚洲成人三级黄色片| 好看的日本中文字幕在线观看二区| 2020国产成人精品视频| 4438全国成人免费视频| 偷拍熟女大胆免费视频| 视频免费在线观看网站| 中文乱码字幕人妻熟女人妻| 日本福利网站一区二区| 青青青在线视频观看97| 小妹妹爱大棒棒免费观看视频 | 免费观看在线中文字幕视频| 欧美不卡一二三区精品| 公侵犯人妻中文字幕巨| 欧美一级日韩一级亚洲一级va| 国产一级一国产一级毛片| 国产精品免费看一区二区三区| 国产 亚洲 欧美 自拍| 国产肥胖熟女又色又爽免费视频| 人妻熟女 亚洲 一页二页| 成人大片男人的天堂| 亚洲一区视频中文字幕在线播放| 熟妇人妻丰满久久久久久久| 免费啪啪啪网站在线观看| 伊人情人成综合视频| 一区二区三区四区 在线播放| 欧美精品999不卡| 国产精品午夜无码AV体验区| 亚洲日本欧美韩国另类综合| 白白色在线免费视频发布视频 | 麻豆白洁少妇在线播放| 亚洲熟女一区二区六区| 最新日韩av电影在线播放| 亚洲情色777中文字幕| 欧美日韩高清片在线观看| 韩日一级人添人人澡人人妻精品| 亚洲高清免费在线观看视频| 亚洲美女a级黄色在线播放| 亚洲一区二区三区四区入口| 日本福利片在线播放| 开心激情五月天作爱片| 69视频在线精品国自产拍| 精产国品一二三产品区别91| 大香蕉尹人在线最新| 大香蕉尹人在线最新| 国语对白性爱三级片免费看| 日本欧美国产在线一区| 九色91操最新在线观看网址| 亚洲欧美日韩电影一区| 一区二区三区av免费天天看| 国产精品福利久久久久| 狠狠干狠狠操免费视频| 在线有码人妻自拍视频| 日韩av电影中文在线免费观看| av天堂hezyo| 日韩无码国产一区二区| 熟妇人妻av无码中文字幕| 欧美区一区二区三视频| 最新福利二区三区视频| 国语精品视频自产自拍| 成人午夜高清福利视频| 午夜福利片无码10000| 5566熟女人妻人妻| 最新免费在线观看污视频| 亚洲av综合av一去二区三区| 国产福利三级在线观看| 日韩成人精品久久久免费看| 制服丝袜 中文字幕 日韩| 国产自拍偷拍在线精品| 一级做性色a爱片久久片| 国产亚洲精品啪啪视频| 精品国产污污污免费入口| 放荡人妻极品少妇全集| a级黄片免费观看| 中国特黄色性生活片| 日本高清在线观看不卡视频| 男人电影天堂在线观看| 十八禁黄色免费污污污亚洲| 亚洲综合熟女乱中文| 爱搞视频在线观看视频91| 亚洲一区二区偷拍女厕所| 夫妻黄色一级性生活片| 欧美精品熟妇免费在线| 欧美性感美女热舞视频| 婷婷综合缴情亚洲五月伊人| 天天在线播放日韩av| 亚洲精品综合欧美精品综合| 可在线免费观看av| 亚洲成人自拍av在线| 2020年亚洲男人天堂网| 亚洲一区二区中文字幕久久 | 亚洲一区视频中文字幕在线播放| 91性高湖久久久久久久久久| 夜色17s精品人妻熟女av| 国产探花自拍亚洲av| 日日躁夜夜躁狠狠操| 大尺度久久久久久久| 凹凸视频一区二区在线观看| 麻豆国产91制片厂| 人妻系列级片在线观看视频| 成年男女免费视频网站无毒| 亚洲中文字幕最新地址| 少妇精品视频一区二区免费看| 老司机伊人99久久精品| 91污污在线观看视频| 一区二区三区不卡免费视频网站 | 婷婷六月天在线视频| 999精品视频免费在线观看| 免费在线观看视频啪啪| 最近中文字幕免费视频一| 亚洲午夜高清在线观看| 国产,亚洲,欧美综合| 欧美一区日韩二区三区四区| 97香蕉久久国产超碰| 天天操天天日天天插天天舔| 国产精品网站亚洲发布| 亚洲精品1卡2卡3卡| 91香蕉国产亚洲一二三区| 老熟女xxxⅹhd老熟女性| 人妻系列在线免费视频| 亚洲自拍偷拍一区二区中文字幕 | 亚洲成人激情在线综合| 精品国产人伦一区二区三区| av成人三级高清日韩| 一级做性色a爱片久久片| 精品欧美乱码久久久| 欧美精品激情在线不卡| 成人免费视频现网站99在线观看 | 伊人精品久久一区二区| 天天早上头和脸出汗是怎么办| 河北全程露脸对白自拍| 美国伦理片午夜理论片| 360偷拍蜜桃臀69式| 精品国产av虐杀两警花| 亚洲天堂av最新在线| av一区二区三区蜜桃| 视频在线 一区二区| 成人免费电影二区三区| 男女啪啪啪网站在线观看免费| 在线观看黄页网站视频网站| 人妻女侠被擒受辱记| 大成色亚洲一二三区| 大香焦一道本一区二区三区| 91精品视频在线观看视频| 青娱乐这里只有精品| 欧美区日本区国产区| 豆豆专区操逼性视频在线| 先锋人妻啪啪中文字幕| 亚洲欧美成人午夜一区二区| 国产高清在线观看av| 国产肥胖熟女又色又爽免费视频| a级黄片免费观看| 伦理在线观看未删减中文字幕 | 午夜精品久久秘?18免费观看| 天天爽天天操天天插| 丰满少妇_区二区三区| 中文字幕在线观看av观看| 亚洲天堂色综合久久| 啊不行啊操逼好爽大鸡吧视频| 中文字幕欧美一区二区视频| 欧美性感美女热舞视频| 天天干天天色综合久久| 精品免费一区二区三区四区视频| 99久久99九九九99九| 久久av色噜噜ai换脸| 天天爱天天日天天爽| 日本少妇三级交换做爰做| 亚洲综合首页综合在线观看| 天天操天天日天天碰| 青青青免费手机视频在线观看| av天堂a亚洲va天堂va里番| 骚穴被阴茎插免费视频| 色视频免费观看网址| 欧美大鸡吧男操女啊啊啊视频| 中文字幕亚洲乱码精品无限| 黑人巨大精品一区二区在线 | 亚洲制服丝袜网站中文字幕| 自拍偷拍视频亚洲一区| 韩国毛片w妈妈的朋友7| 天天操天天搞天天操| 亚洲一区二区三区无码在线| 中文字幕日韩人妻在线三区| 人人人妻人人人妻精品少妇| 日韩在线 中文字幕| 精品视频一区二区三区◇| 韩国资源视频一区二区三区| 国产在线观看av一区| 在线免费视频999| 国产精品免费看一区二区三区| 亚洲精品国品乱码久久久久| 亚洲最强的25个城市| 久久久久夜色国产精品电影| 午夜精品久久久久久久精品乱码| 精久久久久久久久久久久| 日本少妇精品免费视频| 亚洲gay视频在线观看| 天天弄天天草天天日天天| 性色蜜桃臀x88av天美传媒| 放荡人妻极品少妇全集| 九九九九九久久久国产| 日韩一区二区在线播放观看| 七色福利视频在线观看| 国产亚洲精品啪啪视频| 91精品在线视频免费视频| xxnxx国产美女| 丰满少妇_区二区三区| 欧美日韩黄片免费在线观看| 黄色网络中文字幕日本| 日本不卡视频一二三区| 嗯~嗯~啊啊啊~高潮了软件| 不卡视频在线 欧美日韩| 久久无码高清免费视频| 亚洲综合另类欧美久久| 日韩人妻精品久久久久| 制服丝袜 中文字幕 日韩| 亚洲av综合av一去二区三区| 中文字幕一区二区人妻视频| 东京热日本一区二区三区| 亚洲一区视频中文字幕在线播放 | 人妻系列中文字幕大乳丰满人妻 | 快进来插我的逼嗯啊视频| 天天干天天弄天天日| 91中文字幕视频网站| 日韩一级视频一区二区三区| 国产亚洲精品啪啪视频| 18禁男女啪啪啪无遮挡| 国产激情视频在线观看的| 大鸡扒操大逼大片免费关看| 亚洲午夜熟女在线观看| 亚洲欧美国产人成在线| 久久亚洲国产成人精品麻豆 | 在线国产精品欧美| 亚洲成人三级黄色片| 亚洲第一区av中文字幕| 漂亮人妻口爆久久精品| 色就色综合偷拍区欧美在线| 国产精品亚洲精品亚洲| 国产精品剧情在线亚洲| 一级毛片特级毛片免费的| 国产午夜羞羞一区二区三区| 在线 激情 亚洲 视频| 婷婷色九月综合激情丁香| 操人妻人妻天天爽天天偷| 在线观看中文字幕精品av| 免费观看在线中文字幕视频| 午夜美女福利视频在线| 国产毛片特级Av片| 午夜一区二区三区视频在线观看| 亚洲高清免费在线观看视频| 国产精品网站的黄色| 国产一级一国产一级毛片| 手机看片1024精品国产| 国产精品 亚洲欧美 自拍偷拍| 日本少妇三级交换做爰做| 亚洲同性同志一二三专区| 92午夜免费福利视频www| 亚洲同性同志一二三专区 | 欧美日韩久久丝袜在线| 东京热日本一区二区三区| 最新日韩av电影在线播放| 国产高清视频www夜色资源| 欧美色区国产日韩亚洲区| 成人做爰av在线观看网站| 九九六视频,这里只有精品| 欧美黄色性视频网站| 日本小视频一区二区| 熟妇人妻丰满久久久久久久| 91亚洲最新蜜桃在线| 91美女在线观看视频| 亚洲精品一区二区gif| 免费24小时人妻视频| 青青免费观看视频| 999久久久人妻精品一区| 99色在线观看免费观看| 高清av在线婷一区二区色日韩| 精品免费一区二区三区四区视频| 黑人大巨屌操美女逼| 人妻超清中文字幕在线乱码| 黑人爆操女人免费视频| xxnxx国产美女| 亚洲日本欧美韩国另类综合| 日日躁夜夜躁狠狠操| 抽插小穴啊啊啊视频| 久久国产精品久精国产爱| 首页欧美日韩中文字幕| 蜜桃臀少妇白色紧身裤细高跟| 天天看天天爱天天日| 亚洲精品9999蜜桃| 每日更新日韩欧美在线| 午夜精品老牛av一区二区三区| 我爱搞在线观看视频| 2020国产激情视频在线观看| 一二区二区不卡视频| 天天色 天天操 天天好逼| 日本少妇三级交换做爰做| 天堂在线中文字幕av| 夏目彩春av在线看| 77亚洲视频在线观看| 亚洲精品色图1234| 韩国一级片最火爆中文字幕| 欧美久久蜜臀蜜桃资源吧| 凹凸视频一区二区在线观看| 日本清纯中文字幕版| 91精品视频在线观看视频| 操死你美女在线视频| 4438全国成人免费视频| 日本少妇人妻凌辱在线| 欧美不卡一二三区精品| 男女69视频在线观看免费| 国产极品气质外围av| 亚洲情色777中文字幕| 女同大尺度视频网站在线观看| 日本一区二区三区的资源| av 一区二区三区 熟女| 91美女在线观看视频| 国产漂亮白嫩美女在线图片| 天天在线播放日韩av| 亚洲激情噜噜噜久久久| 久久99精品热在线观看| julia人妻av一区二区三区| 91精品91久久久久| 亚洲女人自熨在线视频| jandara在线观看| 美国男的操女孩的小嫩逼| 中文字幕人妻精品精品| 国色天香一二三期区别大象| 9662av在线视频| 日本不卡 中文字幕| 69国产在线视频网站| 日本在线免费观看国产精品| 亚洲人妻系列在线视频| 亚洲一区二区三区四区入口| 一区二区三区四区 在线播放| 人妻在线中文视频视频| 亚洲熟妇在线视频观看| 欧美区日本区国产区| 亚洲免费在线不卡视频| 国产精品午夜无码AV体验区| 天天干天天日天天弄| 91中文字幕视频网站| 中文人妻av一区二区三区| 熟妇人妻丰满久久久久久久| 有码一区二区三区四区五区| 亚洲日本欧美韩国另类综合| 伦理在线观看未删减中文字幕 | 久久久精品人妻无码专区不卡| av在线中文字幕在线| 亚洲成a人片777777张柏芝| 欧美一区二区三区爽爽| 国产经典精品欧美日韩| 1区3区4区产品乱入视频| 骚穴被阴茎插免费视频| 99久久国产精品免费热| 2021国产剧情麻豆| 中文字幕熟女乱一区二区| 欧美大鸡吧男操女啊啊啊视频| 日本福利视频网站导航| 日韩av水蜜桃一区二区三区| 亚洲制服丝袜在线看| 一二区二区不卡视频| 大尺度av毛片在线网址| 老熟妇一区二区三区v∧88| 国产大桥未久一区二区| 中字幕人妻熟女人妻a62v网| 青青青国产精品视频| 亚洲唯美激情综合四射| 国产成人在线观看hd| 乌克兰美女操逼高清内射视频| 午夜国产精品免费视频| 国产白丝一区二区三区av| 久久久久国产精品二区| 欧美一区二区三区视频看| 亚洲熟女乱一区二区精品成人| 黑川堇人妻88av| 青娱乐免费最新视频| 呻吟求饶的人妻中文字幕| 欧美黑人1区2区3区| 亚洲国产精品自拍偷拍视频在线| 放荡人妻极品少妇全集| 不卡一区二区视频在线| 999久久久人妻精品一区| 五月婷婷伊人久久中文字幕| 天天干夜夜操夜夜骑| 国产,亚洲,欧美综合| 丰满少妇人妻一区二区三区蜜桃| 加勒比东京热绿帽人妻多人操| 亚洲人妻系列在线视频| 美女激情久久久久久久| 日本福利视频网站导航| 91人妻人人爽色啊啊啊| 可以免费观看日韩av| 国产精品网站亚洲发布| 亚洲国内精品久久久久久久| 国产精品内射婷婷一级| 核xp工厂精品久久亚洲| 日韩精品视频一区二区三区在线| 国产在线小视频一区二区| av网页免费在线观看| 天天早上头和脸出汗是怎么办| 欧美vs亚洲vs日韩| 亚洲女人自熨在线视频| 亚洲男人天堂最新网址大全| 在线播放 日韩 av| 最新免费在线观看污视频| 中文字幕亚洲无线乱码| 亚洲一区二区中文字幕久久| 亚洲av综合av一去二区三区 | 91麻豆精品国产在线| 女人高潮潮呻吟喷水网站| 久久久久久免费观看av| 女女抠逼白虎白丝袜| 在线观看中文字幕少妇av| 在线 激情 亚洲 视频| 最新激情中文字幕视频| 亚洲av三级电影在线观看| 亚洲最大先锋资源采集站| 婷婷一区二区三区五月丁| 国产精品久久久久久成人久| 国产精品久久久久久成人久| 视频在线+欧美十亚洲曰本| 亚洲第一页欧美第一页| 国产精品无码无卡免费观| 黄版视频在线免费观看| 丝袜美腿日韩av一区| 不卡视频在线 欧美日韩| 中国精品人妻一区二区| 青娱乐不卡视频在线| 午夜福利午夜福利影院| 亚洲成人中文无码在线| 亚洲综合熟女乱中文| 亚洲成人中文无码在线| 黑吊操欧美极品美女| 久久精品国产亚洲av热软件| 五月天男人的天堂中文字幕| 亚洲熟妇丰满多毛xxxx网站| 国产成人情侣激情视频| 国产高清视频www夜色资源| 极品风骚人妻3p视频| 天天色天天射天天日天天干| 男人和女人的逼视频| 亚洲熟妇丰满多毛xxxx网站| 中字幕人妻熟女人妻a62v网| 国产视频成人自拍蝌蚪视频| 18在线观看免费观看| av资源中文字幕在线观看| av天堂hezyo| 夫妻黄色一级性生活片| 亚洲精品色图1234| 中文字幕在线字幕乱码怎么设置 | 国产精品视频网站污污污| 色狠狠色综合久久久绯色| 国产黑色丝袜 在线日韩欧美| 天天在线播放日韩av| 日本少妇人妻中文在线| 无码人妻丰满熟妇区五路| 精品一区二区三区喷水内射高潮| 亚洲成人三级黄色片| 青青操久久综合激情| 91久久精品美女高潮喷水白浆| 区一区二区三免费观看视频| 美女网站福利在线观看| 国产极品气质外围av| 日本人妻少妇xxxxxxx| 人妻少妇的va视频| 女女抠逼白虎白丝袜| 开心五月综合激情婷婷| 青青操天堂在线观看视频| 亚洲黄色免费在线观看网站| 操操操操操操操操操网| 特级aaaaa黄色片| 国长拍拍视频免费孕妇| 亚洲全国精品女人久久久| 日本少妇丰满大bbb的小乳沟| 亚洲午夜国产末满十八岁勿进网站| 日本高清久久人人爽| 99福利一区二区视频| 亚洲妹妹我爱你在线观看| 亚洲一区亚洲二区成人福利| 呻吟求饶的人妻中文字幕| 亚洲乱熟女一区二区三区影片| 精品一区二区三区喷水内射高潮| 久久人妻诱惑我视频| 亚洲国产精品久久久久久无码| 国产成人情侣激情视频| 午夜精品久久久久久久精品乱码| 男女69视频在线观看免费| 亚洲av中文免费在线| 正在播放麻豆精品一区二区| 国产自拍偷拍视频在线免费观看| 中文字幕在线观看亚洲情色| 亚洲国产精品自拍偷拍视频在线| 五月婷婷伊人久久中文字幕| 情趣视频在线观看91| 日韩久久九九精品视频| 夜夜爽夜夜操夜夜爱| 久久久久夜色国产精品电影| 青娱乐免费视频一二三| 久久久视频在线播放| 亚洲天堂男人的天堂| 国产激情视频在线观看的| 欧美黄色性视频网站| 亚洲成人动漫av在线| 亚洲18片综合国产av| 午夜在线观看一级毛| 麻豆国产精品777777在| av大尺度一区二区三区| 在线 激情 亚洲 视频| 99久久国产精品免费热| 色欲AV蜜桃一区二区三| 91污污在线观看视频| 玖辛奈18禁同人污本子| 一区二区三区四区影片| 久久久久久高清一区| 成人人妻h在线观看| 中文字幕免费啪啪啪| 精品日本少妇久久久| 国产美女主播av在线| 国产激情一区二区视频| 中国特黄色性生活片| 亚洲av日韩久久网站| 漂亮人妻口爆久久精品| 69视频在线精品国自产拍| 一区二区三区国产在线成人av| 亚洲中文字幕无线乱码人妻精品 | 亚洲一区二区在线视频观看免费| 亚洲欧美日韩电影一区| 青青在线免费手机播放视频| 大香蕉尹人在线最新| 黑川堇人妻88av| 中文字幕综合网91| 在线观看黄页网站视频网站| av天堂新资源在线| 97精品人妻免费视频| 成熟了的熟妇毛茸茸| 午夜福利午夜福利影院| 岛国av成人午夜高清| 91九色91在线视频| 大香蕉在线欧美在线视频| 视频在线+欧美十亚洲曰本| 91 精品视频在线看| 国产熟妇色xxⅹ交白浆视频 | 啪啪啪网站免费在线看 | 岛国av成人午夜高清| 不用付费特黄特色亚洲特级黄色片| 制服丝袜中文字幕熟女人妻| 亚洲精品色图1234| 日韩精品视频一区二区三区在线 | 日本国产亚洲欧美色综合| jiee日本美女视频网站| 一区二区三区内射美女| 久久国产半精品99精品国产| 中文字幕福利视频在线一区| xxoo福利视频导航| 国产精品免费看一区二区三区| 中文字幕一区二区人妻视频| —区二区三区女厕偷拍| 制服丝袜中文字幕熟女人妻| 91精品国产人妻麻豆| 日本福利视频网站导航| www,日韩av,com| 两个奶被揉得又硬又翘怎么回事| 91激情四射婷婷综合| 精品国产av虐杀两警花| 青青青在线视频观看97| 国产大桥未久一区二区| 亚洲avav天堂av在线网毛片| 天天干天天操天天日天天日| 国产高清自拍偷拍在线| 亚洲成a人77777| 182tv精品免费在线观看| 亚洲国产中文字幕在线看| 核xp工厂精品久久亚洲| 91九色91在线视频| 视频自拍偷拍视频自拍| tushy一区二区三区视频| 琪琪日本福利伦理视频| 4438x亚洲最大的成人| 青青操天堂在线观看视频| 国产天堂av不卡网| 99re这里是国产精品首页| 天天操,天天射,天天爽| 亚洲成人自拍av在线| av激情四射五月婷婷| xxxx69在线观看视频| 波多野结衣在线一区别| 美女网站视频久久精品| 国产成人在线观看hd| 国产大桥未久一区二区| 久久国产半精品99精品国产| 天堂网成人av电影| 综合激情网,激情五月| 成年男女免费视频网站无毒| 久久精品久久久久观看99水蜜桃| 人人妻人人爽人人摸| 亚洲人人爽人人澡起碰av| 不卡一区二区视频在线| 亚洲在线观看中文字幕av| av激情四射五月婷婷| 一区二区三区四区视频精品免费| 桃色成人开心激情网| 九九热精品视频在线播放| 亚洲美女色www色| 美女把腿张开给男的捅| 熟女人妻少妇一区二区| 漂亮人妻口爆久久精品| 港台美女明星av天堂| 69久久夜色精品国产69乱电影| 黄版视频在线免费观看| 午夜情色一区二区三区| 午夜偷拍的视频久久久免费大全| 在线免费观看欧美小视频| 人妻女侠被擒受辱记| 男插女视频大全免费| 天天插天天干天天狠| 欧美一区二区三区视频看| 人妻中文字幕亚洲在线| 欧美日韩不卡视频合集| 东京热日本一区二区三区| 欧美在线观看视频欧美| av在线播放观看h| 99 re国产精品| 欧美插插插插插插| 国产探花自拍亚洲av| 青青青在线视频观看97| 97精品久久久久久无码人妻| 亚洲a级视频在线播放| 视频在线 一区二区| 午夜夫妻性生活视频| 蜜乳av一区二区三区免费观看| 男人和女人的逼视频| 中文字幕国产一区在线视频 | 日韩最近中文在线观看| 182tv精品免费在线观看| 得得爱在线视频观看| 欧美日韩亚洲tv不卡久久| 亚洲国产精品青青草| 欧美日韩高清片在线观看| 亚洲欧美精品海量播放| 亚洲av中文免费在线| 亚洲成人自拍av在线| 亚洲国产精品青青草| 91九色国产在线视频| 国内销魂老女人老泬| 快色视频在线观看免费| 天天弄天天草天天日天天| 中国精品人妻一区二区| 人人妻人人澡人人爽97| 国产人妻777人伦精品hd超碰| 91亚洲国产成人久久精品| 天天夜夜久久精品综合| 天天操天天舔天天做| av中文字幕国产精品| 国产igao激情在线视频入口 | 东京热日本一区二区三区| 成年人免费福利在线| 久久精品四虎夜夜拍拍拍| 川上优所有中文字幕在线| 国际日韩日韩日韩日韩日韩| 人妻在线中文视频视频| 国产精品igao为爱寻找激情| 三级欧美日韩一区二区三区| 国产午夜在线播放视频| 波多野结衣在线一区别| 亚洲精品一区二区gif| caopeng97在线观看视频| 亚洲美女黄色福利视频网站大全| 国产男女无套?免费网站下载| 五月天天堂视频在线| 亚洲国产电影的一区| 欧美一级特黄大片在线| 精品人妻 色中文熟女 oo| 亚洲欧美日韩电影一区| 亚洲精品色图1234| 在线观看2022av| 桃色成人开心激情网| 美利坚合众国av天堂| 欧美精品熟妇免费在线| 免费在线观看亚洲福利| av 资源在线播放| 波多野结衣在线一区别| 伊人情人成综合视频| 69精品人妻久久久久久久久久久| 欧美极品少妇高潮喷水| 日本老女人日比视频| 国产精品无码无卡免费观| 加勒比不卡在线视频| 人人妻人人爽人人爽欧美一区| 北野中文字幕一区二区| av天堂hezyo| 国产av啊啊啊啊啊啊啊| 欧美日韩一区二区三区成人影院| 久久久久九九九九九12| 男人的天堂aⅴ在线| 亚av一二三在线观看| 国产激情视频在线观看的| 日韩av熟妇在线观看| 欧美日本国产一区二区| 不卡一二三区别视频| 69精品人妻久久久久久久久久久| 亚洲欧美国产人成在线| 麻豆国产91制片厂| 91日本精产品一区二区三区| 在线能看视频你懂的| 狠狠操av一区二区三区| 欧美大胆a级视频秒播| 亚洲熟女少妇中文字幕系列| 裸日本资源在线午夜| 蜜臀久久精品久久久久久av| 自拍偷自拍亚洲精品10p| 亚洲一区二区三区四区入口| 亚洲欧洲一区二区三区在线| 日韩黄色在线观看网站上| 亚洲最大先锋资源采集站| 日本不卡 中文字幕| 日韩一区二区在线播放观看| 偷拍熟女大胆免费视频| 亚洲免费午夜污福利| 在线观看免费啪啪啪| 国产高清在线观看av| 青青在线免费手机播放视频| 91国产精品乱码久久久久久| 亚洲国产精品自产拍在线观看| 91福利高清在线播放| 人人妻人人爽人人爽欧美一区 | 女同性恋av在线播放| 干逼又爽又黄又免费的视频| 精产国品一二三77777| 黄在线看片免费人成视频| 午夜亚洲国产精品中字| 熟妇人妻丰满久久久久久久| 蜜桃臀av在线一区二区| 在线观看黄页网站视频网站| 国产精品国产三级在线高清观看| 漂亮人妻口爆久久精品| 日韩女同与成人用品电影免费看| 色老头一区二区三区四区五区| 人妻视频网站快射视频网站| 黄版视频在线免费观看| 老熟女xxxⅹhd老熟女性| 国产白丝一区二区三区av| 亚欧洲乱码视频一二三区| 黑人巨大精品一区二区在线| 中文字幕在线字幕乱码怎么设置| 欧美啪啪一区二区三区| 日韩免费黄色片在线观看| 日本东京热最新中文字幕| 日韩在线 中文字幕| 欧美最新一区二区三区| 伊人免费观看视频一| 西野翔人妻中文字幕中字在| 2020国产激情视频在线观看| 中文字幕久久久国产| 日本免费人爱做视频在线观看不卡 | 在线免费视频999| 一区二区三区婷婷中文字幕| 男插女视频大全免费| 女女抠逼白虎白丝袜| 九九六视频,这里只有精品 | 在线中文字幕人妻av| 日本黄页在线观看视频| xxxx69在线观看视频| 在线看的免费网站黄| 青青青免费手机视频在线观看| 亚洲精品9999蜜桃| 欧美日韩亚洲tv不卡久久| 好看的日本中文字幕在线观看二区| 国产精品久久久久久成人久| 精品不卡一区二区三区| 免费成人av麻豆| 精品人妻在线激情视频| 精品欧美乱码久久久| 亚洲无人区乱码中文字幕一区| 大陆中文字幕视频在线| 91精品资源在线观看| 国产三级自拍视频在线观看网站| 午夜宅男电影av网站| 亚洲天堂色综合久久| 日韩少妇免费在线播放| 亚洲国产精品一区二区第二页| 日本午夜福利免费在线播放| 亚洲天堂男人的天堂| 婷婷色综合五月天视频| 91精品91久久久久| 免费24小时人妻视频| 18岁禁一二三区免费体验| 国产天堂av不卡网| 91美女在线观看视频| 天天色天天射天天日天天干| 操人妻人妻天天爽天天偷| 熟妇高潮久久久久久久| 在线 制服 中文字幕 日韩| 亚洲同性同志一二三专区| 日本一本午夜在线播放| 91精品在线视频免费视频| 中文字幕麻绳捆绑的人妻| 久久久久性感美女偷拍视频| 女人的天堂av在线网| 2018中文字字幕人妻| 猫咪亚洲中文在线中文字幕| 亚洲字幕一区二区夜色av| 女同大尺度视频网站在线观看| av大尺度一区二区三区| av福利免费体验观看| 一区二区三区免费版在线| 91色哟哟视频在线观看| 亚洲av中文免费在线| 色哟哟亚洲乱码国产乱码精品精| 人妻熟女 亚洲 一页二页| 操死你美女在线视频| 裸日本资源在线午夜| 川上优所有中文字幕在线| 欧美日韩高清片在线观看| 两个奶被揉得又硬又翘怎么回事| 熟女一区二区三区综合| 国产91免费在线观看| 黄版视频在线免费观看| 日本不卡视频一二三区| 欧美男女一区二区三区| 在线观看中文字幕少妇av| av天堂a亚洲va天堂va里番| 92午夜免费福利视频www| av丝袜免费在线观看| 久久一级片三上悠亚| 55夜色66夜色亚洲精品| 婷婷一区二区三区五月丁| 精品欧美乱码久久久| 亚洲国产精品一区二区第二页| 91激情四射婷婷综合| 麻豆午夜激情在线观看| 青娱乐这里只有精品| 99免费观看在线视频| 在线免费视频999| 制服丝袜中文字幕熟女人妻| 亚洲av三级电影在线观看| 亚洲熟女乱一区二区精品成人| 欧美日韩高清片在线观看| 天天摸天天干夜夜操| 国产主播诱惑毛片av| 在线看日韩av不卡| 美国十次了亚洲天堂网国产| 黄色av 在线观看| 欧美日韩国产在线中文字幕| 久久久久久高清一区| 男人电影天堂在线观看| 亚洲 偷拍 自拍 欧美| 日韩女同与成人用品电影免费看| 亚洲欧美另类丝袜另类自拍| 黄色网络中文字幕日本| 91精品资源在线观看| 国产免费久久精品99re丫丫| av天堂新资源在线| 全彩漫画口工18禁| 在线国产精品欧美| 亚洲男人的天堂最新网址| 女人的天堂av在线网| 欧美日韩亚洲国产视频二区| 日本韩国欧美在线视频| 91激情四射婷婷综合| 中国精品人妻一区二区| 亚洲最大的自拍偷拍网| 亚洲无码专区中文字幕专区| 最新福利二区三区视频| 欧美vs亚洲vs日韩| 中文字幕欧美一区二区视频| 一区二区九日韩美女| 国产视频成人一区二区| 天天早上头和脸出汗是怎么办| 日韩成人精品久久久免费看| 老司机免费视频福利0| 日本黄页在线观看视频| 最近最新最好看的中文字幕| 国产激情一区二区视频| 久久av色噜噜ai换脸| 色哟哟亚洲乱码国产乱码精品精| 视频自拍偷拍视频自拍| 最近在线中文字幕免费| 亚洲成人五月婷婷久久综合| 久久精品国产亚洲av热软件| 亚洲理论在线a中文字幕97| 国产福利三级在线观看| 欧洲精品在线免费观看| 男人资源站中文字幕| 99女福利女女视频在线播放| 人妻视频网站快射视频网站| 熟女人妻aⅴ一区二区三| 性感美女极品18禁网站在线| 啪啪啪网站免费在线看| 快进来插我的逼嗯啊视频| 91精品久久久久久久久99蜜臀| 成人午夜av电影网| 青青免费观看视频| 最近日韩免费在线观看| 午夜3p福利视频合集| 全彩漫画口工18禁| 日本亚洲午夜福利一区二区三区| 亚洲综合熟女乱中文| 午夜福利在线不卡视频| 操操操操操操操操操网| 青娱乐不卡视频在线| 日韩在线 中文字幕| 亚洲人妻系列在线视频| 夜夜操夜夜爱夜夜摸| 无人区一码二码三码区别在哪| 男人和女人的逼视频| 天天干天天操天天要| 午夜福利午夜福利影院| 久久99精品热在线观看| av网页免费在线观看| 亚洲欧美日韩电影一区| 男人av一区二区三区| 亚洲欧美激情久久久| 国产男女无套?免费网站下载| 天天日天天干天天日天天干天天| 大片a免费观看在线视频观看| av天堂a亚洲va天堂va里番| 亚洲av中文无码网站| 最新日韩中文字幕啪啪啪| 人妻女侠被擒受辱记| 国产漂亮白嫩美女在线图片| 91精品麻豆91夜夜骚| 大乳丰满人妻中文字幕韩国hd| 日韩一区二区在线播放观看| 另类欧美激情校园春色| www国产亚洲精品久久久| 午夜精品小视频在线播放| 日韩国产欧美久久一区| 亚洲 偷拍 自拍 欧美| av天堂a亚洲va天堂va里番| 日本高清激情乱一区二区三区| 欧美在线视频不卡一区| 精品一区二区三区免费毛片W| 日本少妇人妻中文在线| 高潮喷水一区二区三区| 欧美情色av在线观看| 亚洲人成大片在线观看| 日本午夜福利免费在线播放| 青青草一个释放的网站| 欧美日韩高清片在线观看| av 资源在线播放| 日本少妇人妻中文在线| 亚洲图片另类综合小说| 汤姆提醒30秒中转进站口| 欧美久久一区二区伊人| 伊人久久综合国产精品| 蜜乳av中文字幕一区二区| 欧美猛少妇色ⅹⅹⅹⅹⅹ猛叫| 午夜久久人妻一级内射av网址| 无人区一码二码三码区别在哪| 黄色av 在线观看| 天天天天天天天天日日日| 玖玖资源站在线观看亚洲| 中文字幕日本一二三区| av天堂hezyo| 天天日天天亲天天操| 韩国资源视频一区二区三区| 熟妇人妻av无码中文字幕| 熟女人妻少妇一区二区| 久久久久久久岛国免费观看| 国产高清在线观看av| yy4080黄色片| 性感人妻 中文字幕| 无码人妻丰满熟妇区五路| 最新国产午夜激情视频| 亚洲欧美不卡专业视频| 三区美女视频在线观看| 精产国品一二三产品区别97| 亚洲一区在线视频观看地址| 中文字幕一区二区人妻视频| 日本一道中文字幕99| 妈妈的朋友2中文字幕在线| 婷婷六月天在线视频| 中文字幕在线免费观看成人| 亚洲AV无码久久精品国产一区老| 99热99这里免费的精品| 强乱人妻中文字幕日本| 亚洲一区二区在线激情| 五月天男人的天堂中文字幕| 欧美成人性生活视频播放| 97人妻人人揉人人躁人人夜夜爽| 人妻人妻在线视频网站| 天天插天天干天天狠| 裸日本资源在线午夜| 国产资源网站在线播放| 日韩激情文学在线视频| 人人妻人人爽人人摸| 日日夜夜免费视频精品| 丰满少妇_区二区三区| 国际日韩日韩日韩日韩日韩| 亚洲国产电影的一区| 最新日韩中文字幕免费在线观看| 美国伦理片午夜理论片| 国产av嗯嗯啊啊av| 深夜福利免费观看在线看| 亚洲综合首页综合在线观看| 亚洲欧洲无码一区2区无码| 亚洲国产精品一区51动漫| 性高潮视频在线观看日韩| 日本一本午夜在线播放| 国产,亚洲,欧美综合| 亚洲熟女乱色一区二区三区视频| 在线成人教育平台排名| 精久久久久久久久久久久 | 久久亚洲国产成人精品麻豆| 日本福利片在线播放| av里面的动作是真进去吗| av激情四射五月婷婷| 国产原创一区二区三区在线播放| 久久久人妻免费视频| 大屁股熟女一区二区视频| 久久久久久高清一区| 国产精品亚洲精品亚洲| 91精品久久久久久久99蜜月 | 99在线视频精品观看高| 精品人妻人人做人人爽| 911精产国品一二三产区区| www一区二区91| 大鸡扒操大逼大片免费关看| 麻豆午夜激情在线观看| 中文字幕在线字幕乱码怎么设置| 人妻人妻在线视频网站| lutu玩弄人妻短视频| 黄色av日韩在线观看| 欧美一级aaaaaaa片| 岛国av成人午夜高清| 妈妈的朋友中字在线免费观看| 天天综合久久无人区| 得得爱在线视频观看| 公侵犯人妻中文字幕巨| 97精品视频,全部免费| 熟妇人妻丰满久久久久久久| 999精品视频免费在线观看| 成人十欧美亚洲综合在线| 视频免费在线观看网站| 1区3区4区产品乱入视频| 亚洲综合首页综合在线观看| 裸露视频免费在线观看| 午夜一区二区三区视频在线观看| 精品久久久久久久久久久久久| 国产成人情侣激情视频| 91麻豆精品国产在线| 全球高清中文字幕av| 日本亚洲午夜福利一区二区三区| 妈妈的朋友中字在线免费观看| 亚洲欧美韩国日本一区二区| 天天天天天天天天日日日| jizzjizz国产精品传媒| 中日韩又粗又硬又大精品| 欧美亚洲另类精品第一页| 91性高湖久久久久久久久久| 成人免费视频现网站99在线观看| 亚洲 自拍 激情 另类| 国产 亚洲 欧美 自拍| 午夜野花视频在线观看| 亚洲欧美精品日韩偷拍| 亚洲成人 国产精品| 熟女一区二区三区综合| 97cao在线视频| 国产激情免费在线视频| 日韩男女视频网站在线观看| 欧洲精品在线免费观看| xxnxx国产美女| 加勒比东京热绿帽人妻多人操| 啪啪啪网站免费在线看| 9662av在线视频| 亚洲无人区乱码中文字幕一区| 亚洲欧洲一区二区三区在线| 天天看片天天摸天天操| 二十四小时日本高清在线观看| 一区二区三区国产精华液区别大吗| 台湾18禁久久久久久久激情视频| 在线观看黄页网站视频网站| 东京热日本一区二区三区| 手机看片1024精品国产| 亚洲欧美国产人成在线| 国产清纯一区二区在线观看| 天天干天天色综合久久| 性感美女极品18禁网站在线| 中文字字幕在线精品乱码| 999久久久人妻精品一区| 九九六视频,这里只有精品| 丰满少妇高潮喷水视频| 成人免费电影二区三区| 国产精品无码无卡免费观| 久久99精品热在线观看| 第一福利视频在线观看| 黑人爆操女人免费视频| 精品人妻人人做人人爽| 91精品久久久久久久久99蜜臀| 亚州av嫩草av极品在线观看| 性感美女极品18禁网站在线| 日本不卡 中文字幕| 日本人妻熟妇丰满成熟HD系列| 中文字幕欧美人妻在线.| 日本高清激情乱一区二区三区 | 青青免费观看视频| 欧美性感美女热舞视频| 亚洲欧美国产一本综合首页| 黄片操操操操操操c| 抽插小穴啊啊啊视频| 91佛爷视频在线观看| 国产成人在线观看hd| 综合久久伊人久久88| 可在线免费观看av| 中文字幕日韩首页欧美在线激情| 欧美vr专区日韩vr专区| 岳母的诱惑电影在线观看| 九色porny91国产| 中文字幕 人妻 熟女| 91精品国产91久久久久久密臀| 亚洲一区二区三区四区入口| 99久久国产精品免费消防器材| 久久sm人妻中出精品一区二区| 妈妈的朋友2中文字幕在线| 97精品视频,全部免费| 欧洲亚洲一区二区三区四区| 亚洲制服丝袜在线看| 亚洲中文字幕无线乱码人妻精品| 欧美视频亚洲视频在线| 欧美精品激情在线不卡| 在线免费视频999| 亚洲男人的天堂最新网址| 欧美精品乱码99久久蜜桃免费| 4438全国成人免费视频| 裸露视频免费在线观看| 岛国av成人午夜高清| 豆豆专区操逼性视频在线| av一区二区三区四区五区在线| 91九色91在线视频| 狠狠操深爱婷婷综合一区| 免费24小时人妻视频| 中文字幕 一区二区在线观看| 欧洲成熟女人色惰片| 丰满人妻熟女aⅴ一区| 少妇熟女天堂网av| 国产av嗯嗯啊啊av| 全球高清中文字幕av| 福利在线国产小视频| 中文字幕国产一区在线视频| 桃色成人开心激情网| 亚洲中文字幕在线视频观看二区| 伊人综合在线视频免费观看| 一区二区在线观看视频观看| 最新日韩中文字幕免费在线观看 | 日韩一区二区在线播放观看| 18福利视频在线观看| 天天干天天操天天日天天日| 美女精品久久久久久久久| 一二三四区国产在线观看| av 资源在线播放| 夜夜骚av一二三区| 亚洲gay视频在线观看| 91亚洲精品久久蜜桃| 综合激情网,激情五月| 69久久夜色精品国产69乱电影| 1级黄色片在线观看| 日本人妻熟妇丰满成熟HD系列| 亚洲 偷拍 自拍 欧美| 亚洲国产精品自拍偷拍视频在线 | 熟女人妻aⅴ一区二区三| 福利视频导航在线观看| 亚洲码av一区二区三区| 亚洲成a人片777777张柏芝| 天天看片天天摸天天操| 911美女片黄在线观看| 9662av在线视频| 5566熟女人妻人妻| 中文字幕国产一区在线视频| 亚洲乱熟女一区二区三区山| 日本成人福利电影网| 午夜精品小视频在线播放| 大鸡扒操大逼大片免费关看| 午夜3p福利视频合集| avjpm亚洲伊人久久| 国产在线小视频一区二区| 中文字幕在线免费观看人妻| 日本久久久久久黄色| 69久久夜色精品国产69乱电影| 不用付费特黄特色亚洲特级黄色片| 99国产精品久久99久久久| 亚洲黑人欧美二区三区| 91亚洲国产成人久久精品| 99精品久久精品一区二区| 懂色av之国产精品| 亚洲成人自拍图片网站| 区一区二区三免费观看视频| 午夜国产免费视频亚洲| 午夜五十路久久福利| 亚洲自拍偷拍av在线| 亚洲另类激情视频在线看| 中文字幕在线观看av观看| 两个奶被揉得又硬又翘怎么回事| 大乳人妻一区二区三区| 国产清纯一区二区在线观看| 97人妻av人人澡人人爽| 99久久国语露脸国产精品| 美女福利视频一区二区三区四区| 日韩一级视频一区二区三区| 3344永久在线观看视频下载| avtt中文字幕手机版| 天天干天天日天天弄| 亚洲国产美女主播在线观看| 亚洲成人 国产精品| 精品免费一区二区三区四区视频| 日本久久久久久黄色| 无人区一码二码三码区别在哪| 神马午夜久久电影网| xxnxx国产美女| 青青青在线视频免费播放| 久久久人妻免费视频| 亚洲欧洲无码一区2区无码| 亚洲 偷拍 自拍 欧美| 天天插天天透天天爽| 97人妻av人人澡人人爽| 2021国产在线视频| 日韩黄色在线观看网站上| 青青草原在线播放日韩| 三区美女视频在线观看| 爱搞视频在线观看视频91| 极品风骚人妻3p视频| 久操资源在线免费播放| 日本一区二区三区的资源| av无限看熟女人妻另类av| 午夜精品老牛av一区二区三区| 国产自拍偷拍在线精品| 最新福利二区三区视频| 亚洲午夜国产末满十八岁勿进网站| 最新日韩中文字幕免费在线观看| 国产美女视频带a∨黄色片| 精品视频在线观看免费99| 美女张开腿给男人桶爽的软件 | av 一区二区三区 熟女| 熟女人妻精品视频一区| 国语对白性爱三级片免费看| 不用付费特黄特色亚洲特级黄色片 | 国内精品一区二区2021在线| 亚洲欧美综合另类最新| 每日更新日韩欧美在线| 性感美女极品18禁网站在线| 欧美最新一区二区三区| 人妻少妇精品二三区| 天天操天天舔天天射天天日天天干| 国产伦理二区三区在干嘛呢| 国产美女主播av在线| 一区二区三区四区 在线播放| 亚洲乱码国产乱码精品精视频| 日本亚洲精品视频在线观看| 亚洲美女黄色福利视频网站大全| 91精品国产欧美在线| 精品视频在线观看免费99| 日本一区二区三区的资源| 日韩激情文学在线视频| 青青操久久综合激情| 在线看日韩av不卡| avtt中文字幕手机版| 国产91黑丝小视频在线观看| 午夜久久人妻一级内射av网址| 一区二区三区四区视频精品免费| 久久久久久a女人处女| av福利免费体验观看| 18禁男女啪啪啪无遮挡| 久久一级片三上悠亚| av男人站在线观看| 免费在线观看视频啪啪| 另类欧美激情校园春色| 最近最新欧美日韩精品| 360偷拍蜜桃臀69式| 女女抠逼白虎白丝袜| 美女妩媚午夜诱惑网站| 天堂网免费在线电影| 青青免费观看视频| 午夜精品一区二区三区不卡顿| 最新日韩中文字幕啪啪啪| 五月的婷婷综合视频| yy4080黄色片| 日韩免费黄色片在线观看| 天堂一区二区三区在线等| 波多野结衣在线一区别| 亚洲综合第一区二区| 中文字幕在线观看av观看| 精产国品一二三产品区别91| 午夜福利在线不卡视频| 国产白丝一区二区三区av| 好看的日本中文字幕在线观看二区| 亚洲制服丝袜美腿在线| 韩国在线播放一区二区三区| 午夜在线观看一级毛| 欧美区一区二区三视频| 午夜3p福利视频合集| 亚洲美女色www色| 国产黑色丝袜 在线日韩欧美| 性感美女极品18禁网站在线| 日韩av水蜜桃一区二区三区| 亚洲欧美成人午夜一区二区| 青青在线免费手机播放视频| 日韩加勒比精品在线看| 自拍偷拍亚洲综合第一页| 奇米网首页神马久久| aaaa级少妇高潮在线观看| 一区二区欧美 国产日韩| 日本少妇丰满大bbb的小乳沟| 高清欧美色欧美综合网站| 国产 亚洲 欧美 自拍| 麻豆国产精品777777在| 国产欧美福利在线观看| 黄片操操操操操操c| 久久视频 在线播放| 亚洲欧美国产一本综合首页| 青青青在线视频免费播放| 美国伦理片午夜理论片| 天堂一区二区三区在线等| 色哟哟亚洲乱码国产乱码精品精| 亚洲制服丝袜在线看| 免费高清av一区二区| 视频在线 一区二区| 欧美男男在线观看视频网站| 人妻少妇的va视频| 亚洲综合熟女乱中文| 污视频在线观看地址| 2019年中文字幕在线播放视频| 白白色在线免费视频发布视频 | www,日韩av,com| 欧美一区二区三区爽爽| 亚洲字幕一区二区夜色av| 最新日韩av电影在线播放 | 亚洲第一区av中文字幕| 人人妻人人澡人人爽97| 首页欧美日韩中文字幕| 亚洲国产日韩精品在线| 首页欧美日韩中文字幕| 亚洲永远av在线播放| 亚洲一区二区在线视频观看免费| 久久午夜免费鲁丝片| 成人av中文字幕在线看| 欧美黄色一区二区三区视频| 国产精品网站亚洲发布| 真人一进一出抽搐大尺度视频| 精产国品一二三产品区别97| 久草久热这里只有精品| 99久久国产精品免费消防器材|