• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    多GPU異構(gòu)模型實(shí)現(xiàn)放射治療中卷積/積分算法的快速計(jì)算

    2021-12-22 08:20:42賴佳路
    核技術(shù) 2021年12期
    關(guān)鍵詞:模體線程個(gè)數(shù)

    賴佳路 宋 瑩 周 莉 白 雪 侯 氫

    1(四川大學(xué)華西醫(yī)院放療科 成都610041)2(四川大學(xué)原子核科學(xué)技術(shù)研究所輻射物理及技術(shù)教育部重點(diǎn)實(shí)驗(yàn)室 成都610064)

    如今調(diào)強(qiáng)技術(shù)(Intensity Modulated Radiation Therapy,IMRT)已經(jīng)廣泛應(yīng)用于臨床,調(diào)強(qiáng)計(jì)劃制定的質(zhì)量和速度受限于劑量計(jì)算的精度與速度。目前,放療中光子線劑量計(jì)算算法主要有蒙特卡羅算法[1](Monte Carlo,MC)、卷 積/積 分 算 法[2-5](Convolution/Superposition,CS)以及有限筆形束算法[6](Finite Size Pencil Beam,F(xiàn)SPB)。MC算法被稱為劑量計(jì)算的“金標(biāo)準(zhǔn)”,它能計(jì)算各種復(fù)雜條件下的劑量分布[7]。由于MC計(jì)算時(shí)間太長,該方法極少應(yīng)用于臨床。FSPB是目前應(yīng)用于臨床劑量計(jì)算的方法,其計(jì)算速度快,但不能對(duì)三維非均勻介質(zhì)進(jìn)行修正,當(dāng)介質(zhì)非均勻梯度很大時(shí)候帶來誤差較大[8]。CS算 法 模 型 來 源 于Boyer[9]、Mohan[10]和Mackie[11]等在各自研究中提出計(jì)算模型。CS算法能進(jìn)行三維非均勻性修正,其計(jì)算精度僅次于MC。盡管CS算法計(jì)算速度遠(yuǎn)遠(yuǎn)快于MC,但該算法計(jì)算速度要完全應(yīng)用于臨床較困難。1989年Ahnesj?[4]提出了筒串卷積疊加算法(Collapsed Cone Convolution,CCC),該算法對(duì)CS進(jìn)行了簡化,縮短了計(jì)算時(shí)間。目前CCC算法已經(jīng)應(yīng)用于Pinnacle、RayStation治療計(jì)劃系統(tǒng)。由于簡化后的CCC算法計(jì)算精度稍遜于CS算法,本文針對(duì)的是原始的CS算法。目前,圖形處理器(Graphics Processing Unit,GPU)已經(jīng)廣泛應(yīng)用于并行計(jì)算及醫(yī)學(xué)影像處理[12]。王先良等[13]借助單顆GPU對(duì)CS算法進(jìn)行了加速,與僅用CPU計(jì)算相比速度提升了60倍,單野計(jì)算時(shí)間達(dá)到1 min左右。雖然加速效果不錯(cuò),但該速度應(yīng)用于調(diào)強(qiáng)計(jì)劃顯然不足。基于此,本文詳細(xì)研究了采用1~7塊GPU對(duì)CS算法的加速情況,探討通過CPU+多GPU異構(gòu)計(jì)算模型讓CS算法應(yīng)用于臨床調(diào)強(qiáng)計(jì)劃的可行性。

    1 CS算法原理和CPU+多GPU環(huán)境搭建

    1.1 CS算法原理

    CS算法以點(diǎn)核為能量沉積核(Energy Deposit Kernel,EDK),主要涉及兩個(gè)步驟[14]。首先,計(jì)算原射線在單位質(zhì)量內(nèi)釋放的總能量(Total Energy Released Per Unit Mass,Terma)。再 將Terma和EDK進(jìn)行卷積/積分得到最終沉積的能量。其原理如圖1所示,均勻介質(zhì)中任意劑量沉積點(diǎn)r處沉積的能量可通過式(1)計(jì)算。

    圖1 卷積/積分算法原理圖Fig.1 Schematic diagram of convolution/superposition algorithm

    式中:r′為光子源到光子散射點(diǎn)的矢量;r0為光子源到模體表面的距離;r為光子源到能量沉積點(diǎn)的位置矢量;D(r)為r處的劑量;T(r′)為原射線在位移處r′處產(chǎn)生的Terma;keff(r-r′)為有效能量沉積核,它代表原始射線在r′處與介質(zhì)發(fā)生作用,產(chǎn)生的次級(jí)電子在位移r處單位體積產(chǎn)生的能量沉積分?jǐn)?shù)。Terma可通過式(2)和式(3)計(jì)算[4]得到:

    式中:Ψ(Ei)為能量Ei對(duì)應(yīng)的能量注量;uˉ(r)為點(diǎn)r處的平均線性衰減系數(shù);u(Ei,r)為點(diǎn)r對(duì)應(yīng)能量為Ei的線性衰減系數(shù)。

    1.2 CUDA架構(gòu)

    計(jì)算機(jī)統(tǒng)一設(shè)備架構(gòu)(Compute Unified Device Architecture,CUDA)為GPU提供了一整套通用并行計(jì)算解決方案[15]。該架構(gòu)由“中央”處理器(Central Processing Unit,CPU)和“協(xié)”處理器GPU組成。CPU稱為主機(jī)端,GPU稱為設(shè)備端。在CUDA架構(gòu)中,GPU上運(yùn)行的函數(shù)叫內(nèi)核函數(shù)(Kernel)。Kernel函數(shù)以線程網(wǎng)格(Grid)的形式組織,Grid為多個(gè)可以被并行執(zhí)行的線程塊(Block)集合。Block包含大小相同的線程(Thread),Thread是函數(shù)執(zhí)行的最基本單元。GPU的計(jì)算核心是流多處理器(SM),每個(gè)SM以線程束(warp)為單位執(zhí)行,同一個(gè)warp的線程執(zhí)行不同數(shù)據(jù)的相同指令,這種執(zhí)行方式能有效隱藏線程的延遲和等待問題[16]。在CUDA編程模型中,編程者可以通過三個(gè)預(yù)定義變量[17]blockDim、blockIdx和threadIdx實(shí)現(xiàn)數(shù)組元素與線程的對(duì)應(yīng),從而實(shí)現(xiàn)線程間并行計(jì)算。在多GPU架構(gòu)中,主機(jī)只有一個(gè),設(shè)備可以有多個(gè)。借助CUDA 4.0工具包,可以通過cudaSetDevice()函數(shù)在各個(gè)GPU之間切換。

    2 CPU+多GPU實(shí)現(xiàn)CS算法加速

    CS算法的程序?qū)崿F(xiàn)過程主要涉及6個(gè)步驟:

    1)相關(guān)計(jì)算參數(shù)讀??;

    2)創(chuàng)建虛擬模體;

    3)計(jì)算Terma;

    4)計(jì)算卷積/積分;

    5)劑量插值;

    6)計(jì)算所有射束產(chǎn)生的總劑量;。其中:步驟2)~5)這4個(gè)模塊(下文簡稱為M2、M3、M4、M5)計(jì)算數(shù)據(jù)量大、數(shù)據(jù)相關(guān)性低、數(shù)據(jù)有相同的執(zhí)行程序。該4個(gè)模塊適合用GPU進(jìn)行加速。下面將對(duì)4個(gè)模塊進(jìn)行詳細(xì)介紹。

    2.1 創(chuàng)建虛擬模體

    CS算法中的EDK來自MC方法在水中模擬得到的數(shù)據(jù)。模擬過程中,射束垂直入射到模體。因此計(jì)算過程中我們需要?jiǎng)?chuàng)建一個(gè)與射束中心軸方向垂直的虛擬體模,如圖2所示。

    圖2 創(chuàng)建虛擬模體原理圖Fig.2 Schematic diagram of creating virtual phantom

    式中:θG、θT、θC分別代表機(jī)架、治療床以及準(zhǔn)直器角度,Matrixv→b(θG,θT,θC)為射束坐標(biāo)系到人體坐標(biāo)系的轉(zhuǎn)換矩陣。為了得到更精確的CT密度,需要細(xì)化虛擬模體的體元。細(xì)化后的體元多達(dá)數(shù)百萬以上。與CPU循環(huán)計(jì)算方式不同,借助GPU可以讓一個(gè)線程負(fù)責(zé)一個(gè)體元的計(jì)算。

    2.2 計(jì)算Terma

    計(jì)算Terma采用射線追蹤的方法[18]。射線從點(diǎn)源出發(fā)沿著固定的方向依次穿過多個(gè)體元,在射野平面加入射束的強(qiáng)度信息后計(jì)算每個(gè)體元的Terma。為了保留射束發(fā)散的物理性質(zhì),把射束看成連續(xù)的在入射表面等面積的錐形束,這樣可以得到更精確的計(jì)算結(jié)果。為了保證追蹤過程每個(gè)體元都至少有一條射線穿過,需要將射束細(xì)化,細(xì)化后的射束數(shù)目較多,且相互獨(dú)立,該過程非常適合GPU并行計(jì)算。

    2.3 計(jì)算卷積/積分

    卷積計(jì)算中將要計(jì)算的體元作為中心,計(jì)算所有位于積分核半徑Radcut之內(nèi)對(duì)其有影響的Terma產(chǎn)生的劑量總和。卷積/積分中計(jì)算時(shí)間和精度受Radcut大小影響較大。參考陳炳周的研究結(jié)果[19],本文積分核半徑選擇3 cm。計(jì)算過程中將需要計(jì)算的體元r作為追蹤的起點(diǎn),沿著r′-r方向進(jìn)行射線追蹤,然后得到對(duì)應(yīng)等效長度的積分核,對(duì)于密度不均勻的介質(zhì),需要對(duì)其做等效厚度修正,計(jì)算等效厚度使用射線追蹤的方法,三維空間中點(diǎn)(x,y,z)的等效深度d(x,y,z)可以通過下面公式進(jìn)行計(jì)算:

    式中:ρ(x′,y′,z′)是體元(x′,y′,z′)相對(duì)于水的CT值;l(x′,y′,z′)表示射線穿過體元(x′,y′,z′)的幾何長度。CS算法需要將點(diǎn)核離散化[20],為了實(shí)現(xiàn)精度與速度的平衡,本文中,頂點(diǎn)角φ和方位角θ都取48個(gè)。在計(jì)算積分之前,用光線追蹤方法[21]建立一個(gè)體元密度和坐標(biāo)無關(guān)的查找表,在等效厚度計(jì)算的時(shí)候,只需要按角度和步數(shù)Nsteps就可以簡單地得到在x、y、z、r方向的增量,將r方向的增量Δr與體元的電子密度ρ相乘便是射線穿過體元的等效長度。該過程計(jì)算量大,非常適合GPU進(jìn)行加速。

    2.4 劑量插值

    由于機(jī)架角度變化會(huì)改變射束坐標(biāo)系,為了處理不同射束坐標(biāo)系計(jì)算得到的劑量,計(jì)算得到的虛擬模體中劑量還應(yīng)插值到初始模體中去,便于劑量的分析、顯示等。該過程類似于模塊2,同樣適合并行計(jì)算。

    3 測(cè)試結(jié)果與分析

    本文實(shí)現(xiàn)了CS算法基于CUDA架構(gòu)的多GPU并行計(jì)算模型。測(cè)試平臺(tái)參數(shù)如表1所示。測(cè)試模體是一個(gè)大小為30 cm×30 cm×30 cm密度非均勻模體,體元密度值是相對(duì)于水為0~1.2之間的隨機(jī)數(shù),體元尺寸大小設(shè)定為0.2 cm×0.2 cm×0.2 cm,總的體元數(shù)目約為337萬,該體元數(shù)目遠(yuǎn)遠(yuǎn)超過我們臨床中常見部位的體元數(shù)目,如頭部約160萬、胸部約239萬、腹部約250萬。虛擬模體體元的大小和初始模體體元大小設(shè)置相同。射野的大小為10 cm×10 cm,野的分辨率為0.2 cm×0.2 cm。測(cè)試機(jī)架角度為0°和45°。光子能量采用6 MV,GPU采用1~7個(gè)。本文定義了加速比TiG/C,表示第i個(gè)模塊用CPU計(jì)算花費(fèi)時(shí)間與使用GPU花費(fèi)的時(shí)間之比,式如(6)所示,為總的加速因子。

    本文GPU計(jì)算時(shí)間包括在各個(gè)GPU上分配內(nèi)存以及將數(shù)據(jù)拷貝到各個(gè)GPU上進(jìn)行并行計(jì)算所用時(shí)間,由于本實(shí)驗(yàn)中內(nèi)核啟動(dòng)極其快,所以本文忽略掉了內(nèi)核啟動(dòng)時(shí)間。

    表2和表3分別給出了機(jī)架角度為0°和45°時(shí),CS算法4個(gè)模塊單用CPU和采用1~7個(gè)GPU的計(jì)算時(shí)間及加速比。

    多GPU并行計(jì)算的環(huán)境比傳統(tǒng)的單GPU復(fù)雜,本文測(cè)試環(huán)境中使用的是共享系統(tǒng)多GPU,即一個(gè)單獨(dú)的系統(tǒng)中包含多個(gè)GPU,這些GPU通過共享一個(gè)CPU的RAM進(jìn)行連接,各個(gè)GPU數(shù)據(jù)之間不能直接傳輸,必須通過CPU內(nèi)存進(jìn)行中轉(zhuǎn)。隨著GPU使用個(gè)數(shù)的增加,需要在更多的GPU上為核函數(shù)所需數(shù)據(jù)分配內(nèi)存,同時(shí)需要將CPU數(shù)據(jù)分別傳向更多的GPU。當(dāng)多個(gè)GPU掛載到總線上后,總線傳輸壓力增大[22],由于本文研究內(nèi)容需要數(shù)據(jù)在CPU和GPU之間頻繁進(jìn)行交換,數(shù)據(jù)傳輸時(shí)間延遲便成了本文研究過程中的一個(gè)關(guān)鍵問題。特別是對(duì)于在CPU上計(jì)算時(shí)間本身較短的情況。表2和表3中的M2、M3、M5三個(gè)模塊采用多個(gè)GPU計(jì)算時(shí)間說明了這種情況。該三個(gè)模塊在CPU計(jì)算時(shí)間均在1 s之內(nèi),隨著GPU使用個(gè)數(shù)的增加一方面數(shù)據(jù)傳輸總時(shí)間增加,另一方面內(nèi)核啟動(dòng)總時(shí)間也會(huì)相應(yīng)增加,Dana Schaa與David Kaeli對(duì)共享總線結(jié)構(gòu)的單GPU和多GPU的傳輸速率做過測(cè)試,結(jié)果表明單GPU的吞吐率甚至?xí)哂诙郍PU[23]。針對(duì)本文研究情況,建議M2、M3、M5三塊模塊使用1~2個(gè)GPU進(jìn)行計(jì)算。對(duì)于M4模塊,隨著GPU使用個(gè)數(shù)的增加加速倍數(shù)越來越高,分析原因是數(shù)據(jù)傳輸延遲和內(nèi)核啟動(dòng)時(shí)間與模塊整個(gè)計(jì)算時(shí)間相比完全可以忽略不計(jì),當(dāng)采用7個(gè)GPU時(shí)候加速倍數(shù)已經(jīng)高達(dá)180倍以上。

    表2 機(jī)架角度為0°時(shí),2~5四個(gè)模塊計(jì)算時(shí)間和加速比Table 2 Execution time and speedup factors for 2~5 module when gantry angle was 0°

    表3 機(jī)架角度為45°時(shí),2~5四個(gè)模塊計(jì)算時(shí)間和加速比Table 3 Execution time and speedup factors for 2~5 module when gantry angle was 45°

    對(duì)于模塊M4經(jīng)過代碼優(yōu)化后,加速比達(dá)到了198倍以上,主要原因是研究過程中取消了M4中大量的邏輯判斷語句。該模塊采用CPU計(jì)算時(shí),為了減少計(jì)算時(shí)間會(huì)加入一些邏輯判斷語句來減少循環(huán)次數(shù)。GPU與CPU相比并不擅長邏輯運(yùn)算,考慮到GPU中擁有大量線程,因此沒必要通過邏輯判斷來減少循環(huán)次數(shù)。通過測(cè)試不難發(fā)現(xiàn),提高CS算法計(jì)算加速比可以通過兩種方式:一是對(duì)程序代碼進(jìn)行優(yōu)化,盡量減少邏輯判斷語句;二是對(duì)程序中各個(gè)模塊的GPU使用個(gè)數(shù)進(jìn)行優(yōu)化。表4是優(yōu)化代碼和各個(gè)模塊GPU使用個(gè)數(shù)后計(jì)算結(jié)果,對(duì)模塊M2和M5采用1個(gè)GPU,對(duì)模塊M3采用兩個(gè)GPU,對(duì)模塊M4采用7個(gè)GPU且取消M4中邏輯判斷語句。通過優(yōu)化組合,最終加速倍數(shù)達(dá)到了207倍,計(jì)算時(shí)間縮短至9 s。上文提到臨床中頭部、胸部和腹部擁有更少的體元數(shù)目,這些部位若采用本文的方法將花費(fèi)更少的計(jì)算時(shí)間,因此加速后的CS算法能夠滿足臨床要求。

    表4 通過優(yōu)化代碼和GPU使用個(gè)數(shù)優(yōu)化后四個(gè)模塊總的加速比Table 4 The speedup factors for 4 modules after optimizing code and number of GPU used

    由于不是所有的CUDA函數(shù)都滿足IEEE 754標(biāo)準(zhǔn),因此有必要對(duì)結(jié)果準(zhǔn)確性進(jìn)行驗(yàn)證。驗(yàn)證過程中將CPU計(jì)算的劑量結(jié)果作為參考劑量,將采用不同GPU個(gè)數(shù)計(jì)算后的劑量結(jié)果作為評(píng)估劑量,對(duì)二者進(jìn)行γ因子分析[24]。分析過程中距離標(biāo)準(zhǔn)選擇3 mm,百分劑量差標(biāo)準(zhǔn)選擇3%。經(jīng)過評(píng)估,無論采用多少個(gè)GPU進(jìn)行計(jì)算,GPU計(jì)算結(jié)果的γ通過率都是100%。說明使用不同個(gè)數(shù)GPU計(jì)算后的結(jié)果滿足精度要求。值得一提的是勾成俊等[25]借助單個(gè)GPU讓電子束劑量計(jì)算加速了98倍,在后續(xù)工作中,我們將研究多GPU在電子束劑量計(jì)算中加速情況。

    4 結(jié)語

    本文研究采用不同個(gè)數(shù)GPU對(duì)CS算法的加速結(jié)果表明:并行計(jì)算的加速倍數(shù)并非完全與GPU使用個(gè)數(shù)呈線性關(guān)系,還應(yīng)充分考慮數(shù)據(jù)傳輸問題。對(duì)于在CPU上計(jì)算時(shí)間較短的程序代碼,僅用一兩個(gè)GPU進(jìn)行加速即可。此外將CPU執(zhí)行代碼改為GPU執(zhí)行代碼過程中還應(yīng)該考慮到邏輯運(yùn)算問題,在編寫GPU計(jì)算代碼過程中應(yīng)該盡量減少邏輯判斷語句,而不是簡單地對(duì)CPU代碼進(jìn)行移植??傊?,本文研究結(jié)果表明:通過程序代碼和GPU使用個(gè)數(shù)的優(yōu)化,CS算法能完全應(yīng)用于臨床調(diào)強(qiáng)計(jì)劃。

    猜你喜歡
    模體線程個(gè)數(shù)
    怎樣數(shù)出小正方體的個(gè)數(shù)
    基于Matrix Profile的時(shí)間序列變長模體挖掘
    等腰三角形個(gè)數(shù)探索
    怎樣數(shù)出小木塊的個(gè)數(shù)
    植入(l, d)模體發(fā)現(xiàn)若干算法的實(shí)現(xiàn)與比較
    怎樣數(shù)出小正方體的個(gè)數(shù)
    基于網(wǎng)絡(luò)模體特征攻擊的網(wǎng)絡(luò)抗毀性研究
    淺談linux多線程協(xié)作
    基于模體演化的時(shí)序鏈路預(yù)測(cè)方法
    Linux線程實(shí)現(xiàn)技術(shù)研究
    国产精品免费视频内射| 亚洲成国产人片在线观看| 亚洲成人国产一区在线观看| 欧美日本中文国产一区发布| 后天国语完整版免费观看| 男女高潮啪啪啪动态图| 在线观看免费视频日本深夜| 亚洲少妇的诱惑av| 欧美日韩亚洲综合一区二区三区_| 成人黄色视频免费在线看| 校园春色视频在线观看| 亚洲av五月六月丁香网| 一级a爱片免费观看的视频| 久久久久久久午夜电影 | 可以在线观看毛片的网站| 一边摸一边做爽爽视频免费| 亚洲视频免费观看视频| av网站免费在线观看视频| 国产97色在线日韩免费| av电影中文网址| 老汉色∧v一级毛片| a级片在线免费高清观看视频| 国产精品久久久久久人妻精品电影| 精品无人区乱码1区二区| 日日爽夜夜爽网站| 亚洲国产欧美日韩在线播放| 成年人免费黄色播放视频| 久久九九热精品免费| 国产高清videossex| 手机成人av网站| 一级a爱片免费观看的视频| 国产精品自产拍在线观看55亚洲| 免费在线观看亚洲国产| 可以在线观看毛片的网站| 妹子高潮喷水视频| 88av欧美| 国产av又大| 久久影院123| 757午夜福利合集在线观看| 久久人人精品亚洲av| 成人亚洲精品一区在线观看| 视频区图区小说| 日日爽夜夜爽网站| 精品电影一区二区在线| 国产精品爽爽va在线观看网站 | 亚洲国产精品999在线| 欧美老熟妇乱子伦牲交| 亚洲成a人片在线一区二区| 免费在线观看黄色视频的| 最新在线观看一区二区三区| 日韩大尺度精品在线看网址 | 亚洲av片天天在线观看| 欧美不卡视频在线免费观看 | 黄片小视频在线播放| 成在线人永久免费视频| 真人一进一出gif抽搐免费| 色综合欧美亚洲国产小说| 国产精品 欧美亚洲| 色尼玛亚洲综合影院| 色婷婷av一区二区三区视频| 在线观看免费日韩欧美大片| 桃红色精品国产亚洲av| 国产麻豆69| 国产亚洲欧美精品永久| 人妻丰满熟妇av一区二区三区| 亚洲自偷自拍图片 自拍| 精品国产亚洲在线| 欧美日本亚洲视频在线播放| 国产欧美日韩一区二区三| 国产免费现黄频在线看| 亚洲精华国产精华精| 老司机午夜十八禁免费视频| 国产精品一区二区在线不卡| 欧美乱码精品一区二区三区| 亚洲精品在线美女| 欧美日本中文国产一区发布| 在线av久久热| www日本在线高清视频| 黄色 视频免费看| 不卡一级毛片| 波多野结衣高清无吗| 男女下面进入的视频免费午夜 | 欧美日韩亚洲国产一区二区在线观看| 精品国产乱子伦一区二区三区| 国产免费男女视频| 亚洲成人久久性| 亚洲国产看品久久| 亚洲精品粉嫩美女一区| 久久久久久久精品吃奶| 99热国产这里只有精品6| 国产成人啪精品午夜网站| 久久 成人 亚洲| 激情视频va一区二区三区| 天天躁夜夜躁狠狠躁躁| 一级a爱片免费观看的视频| 性色av乱码一区二区三区2| 99久久综合精品五月天人人| 亚洲三区欧美一区| 午夜福利免费观看在线| 精品日产1卡2卡| 一边摸一边抽搐一进一出视频| 亚洲精品久久成人aⅴ小说| 久久久久国内视频| 91精品国产国语对白视频| 嫩草影视91久久| 久久久精品欧美日韩精品| 亚洲成a人片在线一区二区| 欧美人与性动交α欧美软件| 亚洲国产看品久久| 亚洲一卡2卡3卡4卡5卡精品中文| 久久精品亚洲精品国产色婷小说| 精品一品国产午夜福利视频| 在线观看免费午夜福利视频| 亚洲va日本ⅴa欧美va伊人久久| 99re在线观看精品视频| 亚洲一卡2卡3卡4卡5卡精品中文| 午夜91福利影院| 国产在线观看jvid| 亚洲一区二区三区欧美精品| 真人做人爱边吃奶动态| 亚洲国产精品999在线| 免费人成视频x8x8入口观看| 免费在线观看日本一区| 亚洲国产精品999在线| av超薄肉色丝袜交足视频| 级片在线观看| 99久久人妻综合| 欧美成人性av电影在线观看| 50天的宝宝边吃奶边哭怎么回事| 亚洲成人免费av在线播放| 黄频高清免费视频| 多毛熟女@视频| 夜夜看夜夜爽夜夜摸 | 亚洲欧洲精品一区二区精品久久久| 身体一侧抽搐| 男女之事视频高清在线观看| 怎么达到女性高潮| 无遮挡黄片免费观看| 亚洲aⅴ乱码一区二区在线播放 | 露出奶头的视频| 欧美中文日本在线观看视频| 亚洲精品成人av观看孕妇| 99国产精品免费福利视频| 人人妻人人澡人人看| 一本大道久久a久久精品| 男女下面进入的视频免费午夜 | 亚洲成人免费电影在线观看| 日本 av在线| 夜夜爽天天搞| 黄色a级毛片大全视频| 伊人久久大香线蕉亚洲五| 丝袜美腿诱惑在线| 久久人人爽av亚洲精品天堂| 天天躁夜夜躁狠狠躁躁| av网站免费在线观看视频| 亚洲欧美日韩高清在线视频| 国产精品久久久av美女十八| 久久精品亚洲av国产电影网| 国产伦一二天堂av在线观看| 美女高潮喷水抽搐中文字幕| 久久久久国产一级毛片高清牌| 九色亚洲精品在线播放| 在线天堂中文资源库| 桃色一区二区三区在线观看| 激情视频va一区二区三区| 嫩草影院精品99| 欧美日韩亚洲国产一区二区在线观看| 女性被躁到高潮视频| 亚洲国产欧美日韩在线播放| 亚洲伊人色综图| 99精品欧美一区二区三区四区| av超薄肉色丝袜交足视频| aaaaa片日本免费| 国产1区2区3区精品| 亚洲精品一卡2卡三卡4卡5卡| 天天躁狠狠躁夜夜躁狠狠躁| 久久久精品欧美日韩精品| 午夜老司机福利片| 日韩大尺度精品在线看网址 | 天天躁夜夜躁狠狠躁躁| 成年人免费黄色播放视频| 欧美精品一区二区免费开放| 丰满人妻熟妇乱又伦精品不卡| 亚洲国产精品合色在线| 国产精品久久久av美女十八| 一夜夜www| 婷婷丁香在线五月| 久久久久久免费高清国产稀缺| 嫩草影院精品99| www.熟女人妻精品国产| 涩涩av久久男人的天堂| 黄片大片在线免费观看| 一本大道久久a久久精品| 亚洲第一青青草原| 国产单亲对白刺激| ponron亚洲| 亚洲国产精品一区二区三区在线| 午夜日韩欧美国产| 91精品三级在线观看| 少妇的丰满在线观看| 久久久久久免费高清国产稀缺| 国产精品av久久久久免费| 久久久久精品国产欧美久久久| 国产欧美日韩一区二区三区在线| 1024视频免费在线观看| 无限看片的www在线观看| 免费在线观看完整版高清| 亚洲欧美激情综合另类| 午夜两性在线视频| 亚洲黑人精品在线| 久久精品亚洲精品国产色婷小说| 亚洲精品国产精品久久久不卡| 国产精品亚洲av一区麻豆| 一边摸一边做爽爽视频免费| 亚洲av熟女| 脱女人内裤的视频| 午夜福利欧美成人| 日本vs欧美在线观看视频| 亚洲人成网站在线播放欧美日韩| 满18在线观看网站| 精品卡一卡二卡四卡免费| 国产精品影院久久| 久久中文看片网| 一区二区三区激情视频| 国产免费av片在线观看野外av| 欧美在线一区亚洲| 国产一区二区三区综合在线观看| 午夜福利一区二区在线看| 亚洲aⅴ乱码一区二区在线播放 | 在线观看www视频免费| 在线永久观看黄色视频| 男人操女人黄网站| netflix在线观看网站| 一a级毛片在线观看| 日韩大码丰满熟妇| 人人妻人人添人人爽欧美一区卜| 欧美一区二区精品小视频在线| 淫秽高清视频在线观看| 三级毛片av免费| 亚洲一区二区三区不卡视频| 久久亚洲精品不卡| 国产精华一区二区三区| 亚洲欧美日韩高清在线视频| 大型黄色视频在线免费观看| 亚洲片人在线观看| 19禁男女啪啪无遮挡网站| 性色av乱码一区二区三区2| 午夜免费激情av| 国产成人欧美在线观看| 91九色精品人成在线观看| 91字幕亚洲| 久久久久久久午夜电影 | 欧美一区二区精品小视频在线| 国产免费现黄频在线看| 中文字幕色久视频| 亚洲一区中文字幕在线| avwww免费| 1024视频免费在线观看| 淫妇啪啪啪对白视频| 国产精品综合久久久久久久免费 | 80岁老熟妇乱子伦牲交| 亚洲五月婷婷丁香| 丁香六月欧美| 午夜影院日韩av| 日韩av在线大香蕉| 正在播放国产对白刺激| 少妇粗大呻吟视频| 一边摸一边抽搐一进一小说| 老司机在亚洲福利影院| 51午夜福利影视在线观看| 亚洲色图综合在线观看| 亚洲激情在线av| 精品欧美一区二区三区在线| 久久国产乱子伦精品免费另类| 欧美色视频一区免费| 手机成人av网站| 叶爱在线成人免费视频播放| 色播在线永久视频| 成人亚洲精品一区在线观看| 国产1区2区3区精品| 亚洲,欧美精品.| 欧美日韩亚洲国产一区二区在线观看| 涩涩av久久男人的天堂| 日本vs欧美在线观看视频| 99久久久亚洲精品蜜臀av| 久久天躁狠狠躁夜夜2o2o| 一本综合久久免费| 嫩草影院精品99| av国产精品久久久久影院| 午夜福利一区二区在线看| 亚洲第一欧美日韩一区二区三区| 国产三级在线视频| 色婷婷av一区二区三区视频| 婷婷精品国产亚洲av在线| 后天国语完整版免费观看| 又紧又爽又黄一区二区| 欧美在线一区亚洲| 国产精品国产高清国产av| 亚洲av电影在线进入| 50天的宝宝边吃奶边哭怎么回事| 美女高潮喷水抽搐中文字幕| 两个人看的免费小视频| 欧美激情 高清一区二区三区| 精品午夜福利视频在线观看一区| 新久久久久国产一级毛片| 久久香蕉激情| 国产精品综合久久久久久久免费 | 日本一区二区免费在线视频| 午夜激情av网站| 波多野结衣高清无吗| 12—13女人毛片做爰片一| 国产av一区二区精品久久| 69av精品久久久久久| 国产日韩一区二区三区精品不卡| 亚洲国产欧美日韩在线播放| 欧美激情极品国产一区二区三区| 自拍欧美九色日韩亚洲蝌蚪91| 久久久久国产精品人妻aⅴ院| 深夜精品福利| 日本撒尿小便嘘嘘汇集6| 激情在线观看视频在线高清| 免费在线观看完整版高清| 精品久久久久久电影网| 美女国产高潮福利片在线看| 欧美激情 高清一区二区三区| av片东京热男人的天堂| 亚洲五月色婷婷综合| 国产精品一区二区三区四区久久 | 欧美国产精品va在线观看不卡| 亚洲激情在线av| 国产欧美日韩一区二区精品| 亚洲av片天天在线观看| 国产主播在线观看一区二区| 99riav亚洲国产免费| 午夜精品久久久久久毛片777| 日本一区二区免费在线视频| 不卡av一区二区三区| 亚洲男人天堂网一区| 美女 人体艺术 gogo| 精品国产国语对白av| 日本五十路高清| 日日干狠狠操夜夜爽| 黄色a级毛片大全视频| 国产在线观看jvid| 亚洲色图 男人天堂 中文字幕| 色播在线永久视频| 久久天堂一区二区三区四区| 97超级碰碰碰精品色视频在线观看| 一级毛片精品| 亚洲成人免费电影在线观看| 国产精品爽爽va在线观看网站 | 看黄色毛片网站| 精品一品国产午夜福利视频| 欧美日韩视频精品一区| 99国产精品一区二区三区| 亚洲精品美女久久av网站| av网站在线播放免费| 丝袜在线中文字幕| 香蕉丝袜av| 亚洲在线自拍视频| 亚洲成人精品中文字幕电影 | 国产99久久九九免费精品| 久久国产精品男人的天堂亚洲| 18美女黄网站色大片免费观看| 久久久久久久久免费视频了| 亚洲美女黄片视频| 成年人黄色毛片网站| 美女国产高潮福利片在线看| 欧美在线一区亚洲| 国产精品98久久久久久宅男小说| 18禁裸乳无遮挡免费网站照片 | 亚洲九九香蕉| 两性午夜刺激爽爽歪歪视频在线观看 | 美女国产高潮福利片在线看| 午夜福利,免费看| 十分钟在线观看高清视频www| 亚洲色图av天堂| 午夜91福利影院| 窝窝影院91人妻| 成人三级黄色视频| 桃红色精品国产亚洲av| 女人爽到高潮嗷嗷叫在线视频| 欧美日本中文国产一区发布| av天堂久久9| 男女下面进入的视频免费午夜 | 黑人操中国人逼视频| 色在线成人网| 99香蕉大伊视频| 少妇 在线观看| 精品福利永久在线观看| 亚洲av成人一区二区三| 亚洲自拍偷在线| 如日韩欧美国产精品一区二区三区| 日日夜夜操网爽| 欧美日韩亚洲高清精品| 亚洲av成人av| 69精品国产乱码久久久| 一区二区三区国产精品乱码| 最好的美女福利视频网| 久久精品亚洲av国产电影网| 国产精品久久久久成人av| 成年人免费黄色播放视频| 91麻豆av在线| 少妇粗大呻吟视频| 日日干狠狠操夜夜爽| 丝袜美足系列| 久久天堂一区二区三区四区| 欧美中文日本在线观看视频| 亚洲精品国产区一区二| 国产精品 国内视频| 精品久久久精品久久久| 99精品欧美一区二区三区四区| 亚洲国产中文字幕在线视频| 伊人久久大香线蕉亚洲五| 国产精品成人在线| 美女国产高潮福利片在线看| 日韩视频一区二区在线观看| 精品一区二区三区av网在线观看| 国产真人三级小视频在线观看| 亚洲一卡2卡3卡4卡5卡精品中文| 欧美日韩亚洲高清精品| 人人妻人人澡人人看| 日本五十路高清| 亚洲第一欧美日韩一区二区三区| 久久狼人影院| 可以免费在线观看a视频的电影网站| 老汉色∧v一级毛片| 一区二区三区精品91| 又紧又爽又黄一区二区| 免费搜索国产男女视频| 亚洲av五月六月丁香网| 两性夫妻黄色片| 国产一区二区三区视频了| 女警被强在线播放| 一级a爱片免费观看的视频| 午夜成年电影在线免费观看| 男女床上黄色一级片免费看| 色播在线永久视频| 欧美在线一区亚洲| 国产99白浆流出| bbb黄色大片| 天天躁狠狠躁夜夜躁狠狠躁| 深夜精品福利| 国产真人三级小视频在线观看| 三级毛片av免费| 欧美日韩精品网址| 国内毛片毛片毛片毛片毛片| 在线观看免费高清a一片| 深夜精品福利| 久久午夜综合久久蜜桃| 久久精品aⅴ一区二区三区四区| 精品无人区乱码1区二区| 免费在线观看影片大全网站| 午夜免费激情av| 在线观看66精品国产| 男女下面插进去视频免费观看| 欧美日韩国产mv在线观看视频| av有码第一页| 黄片播放在线免费| 欧美午夜高清在线| 人妻久久中文字幕网| 亚洲成人国产一区在线观看| 成人精品一区二区免费| 免费女性裸体啪啪无遮挡网站| 国产精品1区2区在线观看.| 日韩高清综合在线| 嫩草影视91久久| 免费在线观看黄色视频的| 亚洲自拍偷在线| 国产精品永久免费网站| 久久精品影院6| 美国免费a级毛片| 狠狠狠狠99中文字幕| 久热爱精品视频在线9| 制服诱惑二区| 天天躁狠狠躁夜夜躁狠狠躁| 最新美女视频免费是黄的| 国产精品日韩av在线免费观看 | 亚洲人成网站在线播放欧美日韩| 视频区图区小说| 啦啦啦免费观看视频1| 国产日韩一区二区三区精品不卡| 亚洲人成77777在线视频| 欧美精品啪啪一区二区三区| 欧美精品一区二区免费开放| 国产三级在线视频| 国产人伦9x9x在线观看| 国产成人精品在线电影| 国产成+人综合+亚洲专区| 电影成人av| 亚洲成人久久性| 国产成人欧美在线观看| 免费在线观看完整版高清| 国产亚洲精品综合一区在线观看 | 999精品在线视频| 国产亚洲欧美98| 日韩欧美一区视频在线观看| 日日摸夜夜添夜夜添小说| 午夜亚洲福利在线播放| 激情视频va一区二区三区| 日本wwww免费看| 亚洲成人免费av在线播放| www.精华液| 黄色视频不卡| 在线十欧美十亚洲十日本专区| 桃红色精品国产亚洲av| 涩涩av久久男人的天堂| 最新美女视频免费是黄的| 一个人免费在线观看的高清视频| 亚洲欧美精品综合一区二区三区| 欧美日韩一级在线毛片| 日本欧美视频一区| 精品一区二区三区av网在线观看| 免费观看人在逋| 成人特级黄色片久久久久久久| 日本撒尿小便嘘嘘汇集6| 另类亚洲欧美激情| 黑人操中国人逼视频| 午夜福利在线观看吧| 纯流量卡能插随身wifi吗| 好看av亚洲va欧美ⅴa在| 最近最新免费中文字幕在线| 不卡一级毛片| 婷婷丁香在线五月| 成年人免费黄色播放视频| 757午夜福利合集在线观看| 真人一进一出gif抽搐免费| 国产在线观看jvid| 免费观看精品视频网站| 精品免费久久久久久久清纯| 麻豆久久精品国产亚洲av | 国产精品二区激情视频| 老司机深夜福利视频在线观看| 日韩高清综合在线| 91成人精品电影| 欧美久久黑人一区二区| 黄频高清免费视频| 黄色 视频免费看| 成人三级做爰电影| 国产三级黄色录像| 亚洲专区国产一区二区| 99久久国产精品久久久| 国产精品国产av在线观看| 99久久99久久久精品蜜桃| 午夜福利在线观看吧| 亚洲伊人色综图| 国产av一区在线观看免费| 久久久国产精品麻豆| 久久国产精品影院| 国产亚洲欧美98| 美女高潮喷水抽搐中文字幕| 国产精品乱码一区二三区的特点 | 夜夜躁狠狠躁天天躁| 久久精品国产亚洲av高清一级| 12—13女人毛片做爰片一| 欧美另类亚洲清纯唯美| 90打野战视频偷拍视频| 亚洲人成网站在线播放欧美日韩| 成人国语在线视频| 天堂√8在线中文| 女性被躁到高潮视频| 99在线视频只有这里精品首页| 在线观看一区二区三区激情| 久久久久久久久免费视频了| 久久中文字幕人妻熟女| 777久久人妻少妇嫩草av网站| 国产精华一区二区三区| 国产精品综合久久久久久久免费 | 成年人免费黄色播放视频| 看黄色毛片网站| 久久人妻福利社区极品人妻图片| 一二三四在线观看免费中文在| 如日韩欧美国产精品一区二区三区| 天天躁夜夜躁狠狠躁躁| 美女扒开内裤让男人捅视频| 99热只有精品国产| 国产熟女xx| 久久精品国产亚洲av高清一级| 亚洲五月色婷婷综合| 嫁个100分男人电影在线观看| 亚洲成人久久性| 99久久综合精品五月天人人| 男女床上黄色一级片免费看| 9191精品国产免费久久| 在线看a的网站| 老汉色∧v一级毛片| 桃色一区二区三区在线观看| 欧美黑人精品巨大| 黄色怎么调成土黄色| 黄色视频,在线免费观看| 亚洲va日本ⅴa欧美va伊人久久| 久久精品成人免费网站| 每晚都被弄得嗷嗷叫到高潮| 首页视频小说图片口味搜索| 少妇被粗大的猛进出69影院| 可以在线观看毛片的网站| 久久久久久久精品吃奶| 在线观看www视频免费| 热99re8久久精品国产| avwww免费| 色综合婷婷激情| 成年人黄色毛片网站| 欧美黑人精品巨大| 国产伦一二天堂av在线观看| 精品熟女少妇八av免费久了| 91麻豆精品激情在线观看国产 | 久久久久九九精品影院| 亚洲国产精品999在线| 91精品国产国语对白视频| 大香蕉久久成人网| 男女下面进入的视频免费午夜 | 一边摸一边抽搐一进一小说| 亚洲全国av大片| 亚洲情色 制服丝袜| 精品福利永久在线观看|