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

    非結(jié)構(gòu)網(wǎng)格瞬態(tài)電磁場(chǎng)計(jì)算中的高效通信方法

    2022-08-09 07:16:06李岷軒江樹剛吳慶愷林中朝
    關(guān)鍵詞:時(shí)域進(jìn)程網(wǎng)格

    李岷軒,江樹剛,吳慶愷,林中朝

    (西安電子科技大學(xué) 陜西省超大規(guī)模電磁計(jì)算重點(diǎn)實(shí)驗(yàn)室,陜西 西安 710071)

    隨著無(wú)線電技術(shù)和微電子技術(shù)的飛速發(fā)展,大型設(shè)備和電子電路的電磁兼容問(wèn)題變得越來(lái)越普遍且重要。以電磁場(chǎng)仿真算法為基礎(chǔ)的電磁場(chǎng)仿真軟件作為天線設(shè)計(jì)、電磁兼容等工程問(wèn)題的基礎(chǔ)工具,具有非常重要的實(shí)際意義。在計(jì)算電磁學(xué)領(lǐng)域的眾多電磁算法中,時(shí)域有限差分(Finite Difference Time Domain,F(xiàn)DTD)方法、有限元(Finite Element Method,F(xiàn)EM)方法和矩量法(Method of Moments,MoM)是3種精度較高的經(jīng)典算法。而時(shí)域間斷伽略金(Discontinuous Galerkin Time Domain,DGTD)方法兼具時(shí)域有限差分方法的時(shí)域?qū)掝l帶特性和時(shí)域有限元(Finite Element Time Domain,F(xiàn)ETD)方法的非結(jié)構(gòu)網(wǎng)格特性,具有很高的理論精度,尤其在計(jì)算寬帶場(chǎng)路耦合問(wèn)題和復(fù)雜媒質(zhì)問(wèn)題時(shí)更具優(yōu)勢(shì)[1-7]。

    不同于時(shí)域有限元方法求解強(qiáng)耦合單元的特性,時(shí)域間斷伽略金方法使用間斷伽略金法降低了網(wǎng)格間的耦合性,求解過(guò)程中計(jì)算各單元的獨(dú)立矩陣,避免了計(jì)算大型稀疏矩陣,計(jì)算復(fù)雜度遠(yuǎn)小于時(shí)域有限元方法[8]。時(shí)域間斷伽略金方法中常用的顯式迭代方案有l(wèi)eapfrog和ADI等[9-10],在時(shí)域間斷伽略金方法的顯式迭代并行計(jì)算過(guò)程中,各個(gè)單元在迭代更新時(shí)所需的相鄰單元信息為前一時(shí)刻信息,因此需要與該單元的相鄰單元進(jìn)行信息交換。在并行迭代計(jì)算過(guò)程中,這一信息交換過(guò)程是時(shí)域間斷伽略金方法的主要通信熱點(diǎn)。相鄰單元的數(shù)量會(huì)明顯影響并行計(jì)算過(guò)程中的通信量,即通信負(fù)載。將單元?jiǎng)澐纸o各并行進(jìn)程的過(guò)程決定了各個(gè)進(jìn)程的計(jì)算量,即各進(jìn)程的計(jì)算負(fù)載。通信負(fù)載與計(jì)算負(fù)載共同決定了算法的并行效率。受制于區(qū)域劃分方法的局限性[11-12],與采用結(jié)構(gòu)性網(wǎng)格進(jìn)行離散的時(shí)域有限差分[13-14]方法相比,在時(shí)域間斷伽略金方法并行計(jì)算的過(guò)程中需要在計(jì)算負(fù)載和通信負(fù)載間做取舍[15-16]。同時(shí)在并行迭代過(guò)程中,進(jìn)程間的數(shù)據(jù)收發(fā)過(guò)程非常頻繁,因此在大規(guī)模并行狀態(tài)下難以取得很高的并行效率[15]。筆者將給出一種適用于大規(guī)模并行的時(shí)域間斷伽略金通信策略。

    1 并行DGTD原理與通信策略

    無(wú)源麥克斯韋方程組的微分形式中,電場(chǎng)E和磁場(chǎng)H表示為

    (1)

    (2)

    (3)

    其中,Z為單元i的阻抗,Y為單元i的導(dǎo)納;Z+為單元i相鄰單元的阻抗,Y+為單元i相鄰單元的導(dǎo)納;E+、H+為單元i相鄰單元的電場(chǎng)和磁場(chǎng)。由于數(shù)值通量的存在,求解式(2)的過(guò)程中需要相鄰單元的電磁場(chǎng)值。對(duì)式(2)進(jìn)行伽略金檢驗(yàn)后得到[9]:

    (4)

    (5)

    其中,Γlp為與單元i屬于相同進(jìn)程的相鄰單元面,Γcp為與單元i屬于不同進(jìn)程的相鄰單元面。以式(5)為基礎(chǔ),采用不同時(shí)間差分近似方法可以得到關(guān)于ei和hi的并行迭代方程。同時(shí),在不同進(jìn)程的交界面單元處,各進(jìn)程的邊界單元需要與相鄰進(jìn)程單元交換信息。

    (6)

    除邊緣進(jìn)程外(被分配區(qū)域處于模型邊緣的進(jìn)程),各進(jìn)程的通信量Pmessage與計(jì)算量Pamount之比為

    (7)

    其中,αcpn為相鄰進(jìn)程數(shù),對(duì)于立方體均勻分割的情況,αcpn可視為6。γ為每個(gè)相鄰網(wǎng)格所需的通信量,β為各個(gè)網(wǎng)格每次迭代時(shí)的計(jì)算量。式(7)表明,對(duì)每個(gè)進(jìn)程來(lái)說(shuō),其通信量與進(jìn)程數(shù)正相關(guān),計(jì)算量與進(jìn)程數(shù)負(fù)相關(guān),當(dāng)計(jì)算規(guī)模(進(jìn)程數(shù))擴(kuò)大時(shí),程序的計(jì)算效率會(huì)因?yàn)橥ㄐ帕康奶嵘档汀K惴看蔚牡臅r(shí)間滿足

    (8)

    其中,T為每次迭代的總時(shí)間,Tm為每次迭代的通信時(shí)間,Tc為每次迭代的計(jì)算時(shí)間,Sm為計(jì)算平臺(tái)每個(gè)通信周期的通信速度,γp為各進(jìn)程在每個(gè)通信周期中的平均通信量,Sc為計(jì)算平臺(tái)的運(yùn)算速度,N為整個(gè)通信過(guò)程的非自鎖點(diǎn)對(duì)點(diǎn)通信次數(shù)(通信周期)。顯然,要減少算法的迭代時(shí)間,需要提高算法的計(jì)算速度和通信速度。使用圖形處理器(Graphic Processor Unit,GPU)和多核處理器(Many Integerated Core,MIC)卡等異構(gòu)加速技術(shù)和中央處理器(Central Processing Unit,CPU)多線程技術(shù)能夠直接而有效地提升算法的計(jì)算速度[15-17];筆者著重討論提升算法通信速度的方式。根據(jù)文獻(xiàn)[10]中并行效率的定義,n個(gè)進(jìn)程計(jì)算時(shí)的并行效率Rn可以表示為

    (9)

    其中,Tc為基準(zhǔn)核數(shù)的計(jì)算時(shí)間,c為基準(zhǔn)核數(shù)。視Sm、γp、Sc、β為常數(shù),可以得到:

    (10)

    (11)

    (12)

    式(10)表明,在均勻劃分的理想情況下,點(diǎn)對(duì)點(diǎn)通信的DGTD方法迭代過(guò)程的并行效率與并行規(guī)模和通信周期呈負(fù)相關(guān)。其中,Cm是一次迭代中方形結(jié)構(gòu)一個(gè)面上所有網(wǎng)格的通信時(shí)間,Cc是單進(jìn)程情況下整個(gè)模型的總迭代計(jì)算時(shí)間。顯然,在大規(guī)模并行時(shí)想要取得較好的并行效率,就需要減小每次通信的通信周期。

    基于消息傳遞接口(Message Passing Interface,MPI)中的MPI_ISEND和MPI_IRECV兩個(gè)接口,可以方便地完成相鄰進(jìn)程的信息交換[12]。為了降低算法在通信時(shí)的通信周期,需要在通信過(guò)程開始之前規(guī)定好各進(jìn)程點(diǎn)對(duì)點(diǎn)通信的順序。最簡(jiǎn)單的通信順序?yàn)檫M(jìn)程號(hào)由小到大依次通信的順序通信方式。以圖1(a)所示的進(jìn)程關(guān)聯(lián)情況為例,其進(jìn)程間的相鄰關(guān)系可以用圖1(b)所示的關(guān)聯(lián)矩陣表示。圖1(b)中列號(hào)為本地進(jìn)程號(hào),行號(hào)為相鄰進(jìn)程號(hào),進(jìn)程間有相鄰關(guān)系的元素為1,無(wú)相鄰關(guān)系的元素為0。采用順序通信方式后,各個(gè)進(jìn)程的通信順序即通信拓?fù)淙鐖D1(c)所示,各個(gè)矩陣元素值代表對(duì)應(yīng)進(jìn)程之間所屬的通信周期數(shù);該策略通信過(guò)程中的實(shí)際通信情況如圖2所示,圖中由虛線箭頭連接的方框表示一組收發(fā)過(guò)程,橫軸為通信周期數(shù)。

    實(shí)際應(yīng)用過(guò)程中,一個(gè)進(jìn)程往往與多個(gè)進(jìn)程相關(guān)聯(lián),順序通信方式下總通信周期數(shù)近似為進(jìn)程總數(shù),每個(gè)通信周期內(nèi)的通信量遠(yuǎn)小于通信緩沖區(qū)大小,造成了通信資源的浪費(fèi)。

    當(dāng)一對(duì)相鄰進(jìn)程進(jìn)行通信時(shí),與這兩個(gè)進(jìn)程沒(méi)有相鄰關(guān)系的另一對(duì)相鄰進(jìn)程稱為這兩個(gè)進(jìn)程的無(wú)關(guān)通信組。在同一個(gè)通信周期內(nèi),幾個(gè)無(wú)關(guān)通信組能夠同時(shí)進(jìn)行通信操作,從而避免通信帶寬的浪費(fèi)。如圖3(a)所示,按進(jìn)程號(hào)從小到大的順序,依次分析每個(gè)通信周期內(nèi)的無(wú)關(guān)通信組,能夠得到通信周期最小的通信拓?fù)?,稱為最小周期通信拓?fù)?。以下給出獲取最小周期通信拓?fù)涞膫未a。

    最小周期通信拓?fù)鋫未a:

    DO WHILE NotAllCommunicationFinished

    NowPeriod++;

    DO send_rank=0,mpi_size

    IF comm_check(send_rank).equal..false. THEN

    DO recv_rank=ConnectProcessor0,ConnectProcessorN;

    IF comm_check(recv_rank).equal..false. THEN

    comm_check(send_rank)=.true.;CommMatrix(send_rank,recv_rank)=NowPeriod;

    comm_check(recv_rank)=.true.;CommMatrix(recv_rank,send_rank)=NowPeriod;exit;

    ENDIF

    ENDDO

    ENDIF

    ENDDO

    ENDDO。

    最小周期通信的實(shí)際通信過(guò)程如圖3(b)所示。相對(duì)于點(diǎn)對(duì)點(diǎn)通信中簡(jiǎn)單的順序通信方式,最小周期通信方式的總通信周期更少;根據(jù)式(10),其理論并行效率也更高。

    2 數(shù)值算例

    以圖4(a)所示的矩形波導(dǎo)縫隙天線輻射特性計(jì)算為例,驗(yàn)證方法的有效性。矩形波導(dǎo)尺寸為710 mm×13.2 mm×6.5 mm,使用同軸線在波導(dǎo)中心處饋電,采用頻帶為14~20 GHz的調(diào)制高斯脈沖激勵(lì)源激勵(lì),計(jì)算時(shí)使用由Intel Xeon Gold 6140 CPU @ 2.30 GHz處理器組成的計(jì)算集群;該集群共16節(jié)點(diǎn),每節(jié)點(diǎn)含兩塊CPU處理器及72 GB內(nèi)存;每個(gè)處理器核心數(shù)為18,最大線程數(shù)為36。計(jì)算時(shí),給每個(gè)處理器核心分配兩個(gè)進(jìn)程。其中,縫隙天線16 GHz的xOz面和yOz面增益如圖4(b)和圖4(c)所示。表1給出了使用順序通信和最小周期通信的計(jì)算情況。可以看出,在不同計(jì)算核心數(shù)并行規(guī)模的情況下,順序通信的通信周期遠(yuǎn)大于最小周期通信,導(dǎo)致并行效率大大降低,這一點(diǎn)在大規(guī)模計(jì)算時(shí)更為突出。

    表1 矩形波導(dǎo)縫隙天線計(jì)算情況(迭代663 217步)

    文獻(xiàn)[15]采用非阻塞的順序通信策略,在計(jì)算1.6×106網(wǎng)格量的模型時(shí),由6進(jìn)程擴(kuò)展至30進(jìn)程(5倍)時(shí)效率降至90%以下。文中算例網(wǎng)格量為1.1×106,與文獻(xiàn)[15]中算例有一定相似性。文中算例在核心數(shù)為936的計(jì)算過(guò)程中,每個(gè)進(jìn)程平均計(jì)算的網(wǎng)格量為1 176,每對(duì)進(jìn)程的單次通信的網(wǎng)格量約為121個(gè)單元的信息。以120進(jìn)程為基準(zhǔn),采用最小周期通信策略擴(kuò)展8倍的并行效率達(dá)到了約89%,同時(shí)其通信周期約為順序通信的3%,這表明,采用最小周期通信策略能有效降低程序的通信周期,從而提升算法的并行效率。

    以圖5(a)所示的汽車模型為例,驗(yàn)證并行規(guī)模較大時(shí),最小周期通信策略的并行效率。汽車尺寸為4.8 m×1.8 m×1.45 m,計(jì)算頻率為1 GHz,計(jì)算平臺(tái)為神威·太湖之光超級(jí)計(jì)算機(jī),激勵(lì)源為水平極化平面波,入射方向沿z軸負(fù)向,極化方向沿x軸正向,計(jì)算θ=90°,φ由0°到360°的xOy面雙站散射截面(Radar Cross Section,RCS)。需要注意的是,神威·太湖之光超級(jí)計(jì)算機(jī)采用由申威26010異構(gòu)眾核處理器組成的計(jì)算節(jié)點(diǎn),該處理器由1個(gè)主核和64個(gè)從核組成,稱為一個(gè)核組,每個(gè)節(jié)點(diǎn)包含4個(gè)核組[18]。計(jì)算過(guò)程中,每個(gè)核組作為一個(gè)進(jìn)程共同參與計(jì)算。DGTD方法在申威26010異構(gòu)眾核處理器上運(yùn)行時(shí)需要對(duì)主核和從核分別進(jìn)行編程[13],而筆者主要討論核組間的通信過(guò)程,對(duì)此不多敘述。

    圖5(b)展示了DGTD與MoM的RCS計(jì)算結(jié)果。該算例迭代10 000步的計(jì)算時(shí)間和并行效率如圖6所示。以800進(jìn)程為基準(zhǔn),當(dāng)并行規(guī)模擴(kuò)展到8 000進(jìn)程(并行規(guī)模擴(kuò)展10倍)時(shí),并行效率約為70.38%。根據(jù)式(8),圖6中的通信時(shí)間受通信周期數(shù)和進(jìn)程數(shù)增加的共同影響而逐漸上升,這是影響其并行效率的主要因素。

    3 結(jié)束語(yǔ)

    在并行規(guī)模較大時(shí),點(diǎn)對(duì)點(diǎn)通信的DGTD方法的并行效率主要取決于每次迭代時(shí)通信過(guò)程的通信周期數(shù)。筆者采用最小周期通信策略能夠有效降低通信時(shí)的通信周期,使算法在進(jìn)程間進(jìn)行點(diǎn)對(duì)點(diǎn)通信時(shí)取得更好的并行效率,從而提升DGTD方法在大規(guī)模并行計(jì)算時(shí)的計(jì)算效率,進(jìn)而使其更高效地應(yīng)用于工程問(wèn)題的計(jì)算中。

    猜你喜歡
    時(shí)域進(jìn)程網(wǎng)格
    用全等三角形破解網(wǎng)格題
    債券市場(chǎng)對(duì)外開放的進(jìn)程與展望
    反射的橢圓隨機(jī)偏微分方程的網(wǎng)格逼近
    基于時(shí)域信號(hào)的三電平逆變器復(fù)合故障診斷
    重疊網(wǎng)格裝配中的一種改進(jìn)ADT搜索方法
    基于極大似然準(zhǔn)則與滾動(dòng)時(shí)域估計(jì)的自適應(yīng)UKF算法
    基于曲面展開的自由曲面網(wǎng)格劃分
    基于時(shí)域逆濾波的寬帶脈沖聲生成技術(shù)
    基于時(shí)域波形特征的輸電線雷擊識(shí)別
    社會(huì)進(jìn)程中的新聞學(xué)探尋
    精品久久久久久久久久久久久| 性插视频无遮挡在线免费观看| 91麻豆精品激情在线观看国产| 亚洲自拍偷在线| 欧美又色又爽又黄视频| 村上凉子中文字幕在线| 国产av在哪里看| 久久精品夜夜夜夜夜久久蜜豆| 男女做爰动态图高潮gif福利片| 一区福利在线观看| 亚洲av美国av| a级毛片免费高清观看在线播放| 亚洲精品粉嫩美女一区| 久久久久久国产a免费观看| 99在线人妻在线中文字幕| 亚洲无线在线观看| 一夜夜www| 美女cb高潮喷水在线观看| 久久久国产成人精品二区| 中文字幕av在线有码专区| 三级国产精品欧美在线观看| 女生性感内裤真人,穿戴方法视频| 欧美性猛交╳xxx乱大交人| 国产伦人伦偷精品视频| 俺也久久电影网| 99精品久久久久人妻精品| 日韩 亚洲 欧美在线| 99久国产av精品| 一本综合久久免费| 免费看a级黄色片| a级毛片a级免费在线| 亚洲成人久久爱视频| 免费人成在线观看视频色| 亚洲精品在线美女| 欧美黄色片欧美黄色片| 国产精华一区二区三区| 欧美最黄视频在线播放免费| 国产午夜福利久久久久久| 丰满人妻一区二区三区视频av| 偷拍熟女少妇极品色| 在线免费观看不下载黄p国产 | 国产欧美日韩一区二区精品| 国产精品电影一区二区三区| 大型黄色视频在线免费观看| 国产精品精品国产色婷婷| 亚洲人成伊人成综合网2020| 一个人看视频在线观看www免费| 成人国产一区最新在线观看| 亚洲美女黄片视频| 在线观看免费视频日本深夜| 精品欧美国产一区二区三| 国产v大片淫在线免费观看| 美女被艹到高潮喷水动态| 国内精品一区二区在线观看| 在线观看av片永久免费下载| eeuss影院久久| 亚洲,欧美精品.| 一个人免费在线观看电影| 最近在线观看免费完整版| 欧美极品一区二区三区四区| 午夜免费男女啪啪视频观看 | 日韩有码中文字幕| 国产大屁股一区二区在线视频| 国产伦一二天堂av在线观看| 又粗又爽又猛毛片免费看| 久久精品国产清高在天天线| 丝袜美腿在线中文| 亚洲精品456在线播放app | 搡老熟女国产l中国老女人| 国产精品综合久久久久久久免费| 久久久久久久久大av| 18禁在线播放成人免费| 国产亚洲精品久久久com| 丰满人妻熟妇乱又伦精品不卡| 特级一级黄色大片| 每晚都被弄得嗷嗷叫到高潮| 好男人电影高清在线观看| 日韩成人在线观看一区二区三区| 丁香欧美五月| 国产精品一区二区三区四区久久| 久久久久久大精品| 老鸭窝网址在线观看| 别揉我奶头~嗯~啊~动态视频| 美女高潮的动态| 老女人水多毛片| 成人欧美大片| 免费无遮挡裸体视频| 成人av一区二区三区在线看| av欧美777| www.999成人在线观看| av视频在线观看入口| 亚洲欧美激情综合另类| 亚洲,欧美,日韩| 18禁黄网站禁片免费观看直播| 日本a在线网址| 国产精品亚洲av一区麻豆| 国产探花极品一区二区| 一进一出好大好爽视频| 一本一本综合久久| 99精品在免费线老司机午夜| 日韩欧美精品v在线| 一区二区三区四区激情视频 | 国产成年人精品一区二区| 最近在线观看免费完整版| 真人做人爱边吃奶动态| 国产三级在线视频| 一边摸一边抽搐一进一小说| 又黄又爽又刺激的免费视频.| 不卡一级毛片| 色综合站精品国产| 欧美色视频一区免费| 非洲黑人性xxxx精品又粗又长| 亚洲国产日韩欧美精品在线观看| 精品久久久久久,| 久久6这里有精品| 直男gayav资源| 国产午夜福利久久久久久| 国产伦在线观看视频一区| 亚洲精品一卡2卡三卡4卡5卡| 一a级毛片在线观看| 欧美精品国产亚洲| 亚洲国产色片| 欧美激情国产日韩精品一区| 黄片小视频在线播放| 日韩国内少妇激情av| 18禁黄网站禁片免费观看直播| 天天躁日日操中文字幕| 丁香六月欧美| 男女下面进入的视频免费午夜| 俺也久久电影网| 日日夜夜操网爽| 日韩人妻高清精品专区| 搡老妇女老女人老熟妇| 精品久久久久久成人av| 国产爱豆传媒在线观看| 一夜夜www| 精品人妻一区二区三区麻豆 | 日韩大尺度精品在线看网址| 亚洲乱码一区二区免费版| 久久6这里有精品| 国产 一区 欧美 日韩| 悠悠久久av| 深爱激情五月婷婷| 国产精品永久免费网站| 欧美日韩黄片免| 欧美性猛交黑人性爽| 国产69精品久久久久777片| 国产精品国产高清国产av| 一进一出好大好爽视频| 99久久精品国产亚洲精品| 黄色配什么色好看| 国产精品日韩av在线免费观看| 97超级碰碰碰精品色视频在线观看| 又黄又爽又免费观看的视频| 波多野结衣高清作品| 欧美xxxx性猛交bbbb| 伊人久久精品亚洲午夜| 亚洲熟妇熟女久久| 成熟少妇高潮喷水视频| 成人毛片a级毛片在线播放| av在线天堂中文字幕| 国产精品久久电影中文字幕| 99久久久亚洲精品蜜臀av| 美女xxoo啪啪120秒动态图 | a级毛片a级免费在线| 老司机福利观看| 熟女电影av网| 久久久精品欧美日韩精品| 欧美日韩亚洲国产一区二区在线观看| 国产精品综合久久久久久久免费| 中亚洲国语对白在线视频| 国产色婷婷99| 2021天堂中文幕一二区在线观| 欧美黑人巨大hd| 高潮久久久久久久久久久不卡| 亚洲精品在线观看二区| 午夜老司机福利剧场| 国产欧美日韩精品一区二区| 欧美在线黄色| 18禁黄网站禁片午夜丰满| 亚洲国产精品合色在线| 免费高清视频大片| 中文字幕免费在线视频6| 国产午夜福利久久久久久| 日韩欧美免费精品| 啦啦啦观看免费观看视频高清| 亚洲精品一区av在线观看| 欧美黑人巨大hd| 最新在线观看一区二区三区| 级片在线观看| 国产精品人妻久久久久久| 在线观看免费视频日本深夜| 人人妻,人人澡人人爽秒播| 亚洲av美国av| 国产亚洲欧美98| 三级男女做爰猛烈吃奶摸视频| 亚洲国产欧美人成| 亚洲不卡免费看| av在线蜜桃| 狠狠狠狠99中文字幕| 精品一区二区三区av网在线观看| 欧美在线黄色| 成人性生交大片免费视频hd| av福利片在线观看| 观看美女的网站| 乱码一卡2卡4卡精品| 日韩精品青青久久久久久| 色播亚洲综合网| 窝窝影院91人妻| 男女下面进入的视频免费午夜| 麻豆成人午夜福利视频| av天堂中文字幕网| 少妇裸体淫交视频免费看高清| 一本精品99久久精品77| 精品久久久久久久久亚洲 | 午夜精品一区二区三区免费看| 日韩欧美精品v在线| 久久九九热精品免费| 一本综合久久免费| 国产高清视频在线观看网站| 欧美高清性xxxxhd video| 久久精品人妻少妇| 久久性视频一级片| 欧美乱色亚洲激情| 国产亚洲欧美在线一区二区| 欧洲精品卡2卡3卡4卡5卡区| 免费看a级黄色片| 久久性视频一级片| 国产精品1区2区在线观看.| 九九热线精品视视频播放| 直男gayav资源| 精品人妻1区二区| 18禁裸乳无遮挡免费网站照片| 在线播放无遮挡| 91麻豆精品激情在线观看国产| 夜夜夜夜夜久久久久| 亚洲中文字幕日韩| 69av精品久久久久久| av专区在线播放| 搡女人真爽免费视频火全软件 | 亚洲最大成人av| 中文在线观看免费www的网站| 欧美日本亚洲视频在线播放| av女优亚洲男人天堂| 国内精品一区二区在线观看| 黄色视频,在线免费观看| 国产 一区 欧美 日韩| 搡老岳熟女国产| 在线看三级毛片| 亚洲经典国产精华液单 | 伊人久久精品亚洲午夜| 露出奶头的视频| 国产色爽女视频免费观看| 我的老师免费观看完整版| 亚洲av第一区精品v没综合| 看免费av毛片| 国产精品av视频在线免费观看| 国产成人av教育| 又紧又爽又黄一区二区| 一区二区三区四区激情视频 | 国产高清有码在线观看视频| 亚洲国产日韩欧美精品在线观看| 久久久久久久午夜电影| 亚洲人成网站高清观看| 国产免费av片在线观看野外av| 国产亚洲欧美98| 永久网站在线| 亚洲精品影视一区二区三区av| 中文字幕av在线有码专区| 色吧在线观看| 日本在线视频免费播放| 3wmmmm亚洲av在线观看| 97超级碰碰碰精品色视频在线观看| 色哟哟哟哟哟哟| 欧美成人a在线观看| 国产亚洲精品av在线| 亚州av有码| 亚洲人成网站在线播| 婷婷精品国产亚洲av在线| 禁无遮挡网站| 久久香蕉精品热| 亚洲国产精品成人综合色| 日韩欧美在线乱码| 成人特级av手机在线观看| 国产大屁股一区二区在线视频| www日本黄色视频网| 欧美xxxx性猛交bbbb| 动漫黄色视频在线观看| 成人无遮挡网站| 日韩国内少妇激情av| 露出奶头的视频| 啦啦啦观看免费观看视频高清| 亚洲综合色惰| 99久久成人亚洲精品观看| 婷婷色综合大香蕉| 国产单亲对白刺激| 色在线成人网| 国产69精品久久久久777片| 午夜福利高清视频| 国产淫片久久久久久久久 | 少妇裸体淫交视频免费看高清| a级一级毛片免费在线观看| 丰满人妻一区二区三区视频av| 免费黄网站久久成人精品 | av在线老鸭窝| 宅男免费午夜| av在线老鸭窝| 婷婷六月久久综合丁香| 精品久久久久久成人av| 国内揄拍国产精品人妻在线| 欧美另类亚洲清纯唯美| 一区福利在线观看| 亚洲国产精品久久男人天堂| 男女做爰动态图高潮gif福利片| 蜜桃亚洲精品一区二区三区| av在线观看视频网站免费| 日韩有码中文字幕| 日韩欧美 国产精品| 日韩欧美精品免费久久 | 国产精品嫩草影院av在线观看 | 成人国产一区最新在线观看| 国产成人av教育| 成人av一区二区三区在线看| 亚洲av免费高清在线观看| 亚洲人成网站在线播放欧美日韩| 十八禁人妻一区二区| 精品福利观看| 制服丝袜大香蕉在线| 看黄色毛片网站| 免费高清视频大片| 美女被艹到高潮喷水动态| www.999成人在线观看| 国产色爽女视频免费观看| 他把我摸到了高潮在线观看| 亚洲天堂国产精品一区在线| 国产精品免费一区二区三区在线| 国产综合懂色| 色综合亚洲欧美另类图片| 亚洲不卡免费看| 久久精品国产自在天天线| 成年女人毛片免费观看观看9| 欧美日本亚洲视频在线播放| 久久精品国产99精品国产亚洲性色| 人人妻,人人澡人人爽秒播| 国产成人av教育| 日韩人妻高清精品专区| 99热6这里只有精品| 国产高清视频在线播放一区| 一进一出抽搐动态| 日韩欧美国产一区二区入口| 日本撒尿小便嘘嘘汇集6| 亚洲精品亚洲一区二区| 青草久久国产| 色综合婷婷激情| 久久精品久久久久久噜噜老黄 | 欧美午夜高清在线| 禁无遮挡网站| 人妻制服诱惑在线中文字幕| 欧美精品国产亚洲| 如何舔出高潮| 国内少妇人妻偷人精品xxx网站| 欧美成人一区二区免费高清观看| 日韩欧美在线乱码| 美女黄网站色视频| 首页视频小说图片口味搜索| 国内少妇人妻偷人精品xxx网站| 欧美成人一区二区免费高清观看| 亚洲av中文字字幕乱码综合| 久久精品国产清高在天天线| 很黄的视频免费| 国产一区二区三区视频了| 国产精品爽爽va在线观看网站| 国产精品永久免费网站| 欧美又色又爽又黄视频| 中国美女看黄片| 久久精品国产亚洲av涩爱 | 欧美性感艳星| 乱码一卡2卡4卡精品| 亚洲乱码一区二区免费版| 午夜激情福利司机影院| 村上凉子中文字幕在线| 老司机福利观看| 少妇丰满av| 久久伊人香网站| 亚洲综合色惰| 啦啦啦观看免费观看视频高清| АⅤ资源中文在线天堂| 中文字幕高清在线视频| av在线天堂中文字幕| 国产伦一二天堂av在线观看| 久久久久久久精品吃奶| av福利片在线观看| 久久午夜亚洲精品久久| 大型黄色视频在线免费观看| 国产69精品久久久久777片| 亚洲精品日韩av片在线观看| 蜜桃久久精品国产亚洲av| 99久久精品热视频| 久久久久精品国产欧美久久久| 亚洲av.av天堂| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 欧美色欧美亚洲另类二区| 最好的美女福利视频网| 亚洲av不卡在线观看| 亚洲美女视频黄频| 国产精品一区二区免费欧美| 九色成人免费人妻av| 少妇的逼水好多| 中文字幕熟女人妻在线| 又黄又爽又免费观看的视频| www.www免费av| 日韩中文字幕欧美一区二区| 成人性生交大片免费视频hd| 国产三级在线视频| 亚洲人成网站在线播| 欧美午夜高清在线| 国产精品电影一区二区三区| 18禁在线播放成人免费| 国产乱人伦免费视频| 高清毛片免费观看视频网站| 国产av不卡久久| 亚洲中文日韩欧美视频| 亚洲七黄色美女视频| 国语自产精品视频在线第100页| 亚洲av熟女| av福利片在线观看| 国产欧美日韩一区二区精品| 午夜免费成人在线视频| 精品99又大又爽又粗少妇毛片 | 久久久久九九精品影院| 每晚都被弄得嗷嗷叫到高潮| 两人在一起打扑克的视频| 亚洲国产精品sss在线观看| 亚洲久久久久久中文字幕| 我的女老师完整版在线观看| 国产av麻豆久久久久久久| 嫩草影院精品99| 亚洲欧美日韩东京热| 久久6这里有精品| 在线a可以看的网站| 日韩亚洲欧美综合| 亚洲精品一区av在线观看| 久久中文看片网| 国产大屁股一区二区在线视频| 免费高清视频大片| 欧美日韩国产亚洲二区| 亚洲精华国产精华精| 不卡一级毛片| 桃色一区二区三区在线观看| 国产精品av视频在线免费观看| av在线老鸭窝| 国产一区二区三区在线臀色熟女| 久久久久国产精品人妻aⅴ院| 亚洲精品久久国产高清桃花| 亚洲成人久久性| 成年免费大片在线观看| av欧美777| 又黄又爽又刺激的免费视频.| 色综合站精品国产| 成人国产综合亚洲| 精品欧美国产一区二区三| 久久久久亚洲av毛片大全| 国产单亲对白刺激| 宅男免费午夜| 国产精品人妻久久久久久| 亚洲欧美日韩卡通动漫| 999久久久精品免费观看国产| 国产精品美女特级片免费视频播放器| 成人三级黄色视频| 观看免费一级毛片| 色综合站精品国产| 两性午夜刺激爽爽歪歪视频在线观看| 97热精品久久久久久| 人妻夜夜爽99麻豆av| 色视频www国产| 久久精品夜夜夜夜夜久久蜜豆| 亚洲av成人精品一区久久| 国产一区二区三区在线臀色熟女| 国产乱人伦免费视频| 1024手机看黄色片| 亚洲国产精品999在线| 精品人妻一区二区三区麻豆 | 国产精品综合久久久久久久免费| 久久精品国产99精品国产亚洲性色| 精品国产三级普通话版| 成人欧美大片| 亚洲欧美日韩卡通动漫| av在线蜜桃| 欧美3d第一页| 在线免费观看的www视频| 精品久久久久久久人妻蜜臀av| 在线观看美女被高潮喷水网站 | 一个人免费在线观看电影| 女同久久另类99精品国产91| 51国产日韩欧美| 国产一区二区在线观看日韩| 亚洲av成人精品一区久久| 国产精品爽爽va在线观看网站| 国产毛片a区久久久久| 美女cb高潮喷水在线观看| 男人舔奶头视频| 午夜福利成人在线免费观看| 国产高清视频在线观看网站| 91狼人影院| 亚洲精品久久国产高清桃花| 神马国产精品三级电影在线观看| 亚洲久久久久久中文字幕| 99在线人妻在线中文字幕| 99国产精品一区二区三区| 亚洲精品影视一区二区三区av| 成人毛片a级毛片在线播放| 特大巨黑吊av在线直播| 日本精品一区二区三区蜜桃| 亚洲,欧美,日韩| 男插女下体视频免费在线播放| 婷婷精品国产亚洲av在线| 99久久成人亚洲精品观看| 久久久久久大精品| 久久国产精品影院| 老司机午夜十八禁免费视频| 99久久精品热视频| 18+在线观看网站| 一区二区三区免费毛片| 神马国产精品三级电影在线观看| www.熟女人妻精品国产| 久久久精品欧美日韩精品| 久久国产精品影院| 欧美色欧美亚洲另类二区| 一进一出好大好爽视频| 国产亚洲精品av在线| 午夜老司机福利剧场| 精品人妻1区二区| 日韩欧美国产在线观看| 美女xxoo啪啪120秒动态图 | 很黄的视频免费| 色综合亚洲欧美另类图片| 午夜日韩欧美国产| 在线观看午夜福利视频| 久久久久亚洲av毛片大全| av黄色大香蕉| 偷拍熟女少妇极品色| 精品国产亚洲在线| 国产v大片淫在线免费观看| 欧美精品国产亚洲| 久久久久国产精品人妻aⅴ院| 亚洲国产色片| 国产高清激情床上av| 麻豆一二三区av精品| 日韩av在线大香蕉| 日韩欧美 国产精品| 精品久久久久久成人av| 亚洲无线在线观看| 在线观看免费视频日本深夜| 国产69精品久久久久777片| 国产成年人精品一区二区| 伊人久久精品亚洲午夜| 成人性生交大片免费视频hd| 亚洲片人在线观看| 国产亚洲精品综合一区在线观看| 男插女下体视频免费在线播放| 国产真实伦视频高清在线观看 | 两个人的视频大全免费| 校园春色视频在线观看| a级毛片免费高清观看在线播放| 一级黄色大片毛片| 亚洲一区高清亚洲精品| 长腿黑丝高跟| h日本视频在线播放| 免费在线观看成人毛片| 88av欧美| 99久久九九国产精品国产免费| 能在线免费观看的黄片| 丰满的人妻完整版| 国产精品99久久久久久久久| 真人一进一出gif抽搐免费| 在线观看一区二区三区| 国产在线男女| 日韩高清综合在线| 亚洲精品一区av在线观看| 精品一区二区三区av网在线观看| 欧美激情在线99| 日韩精品青青久久久久久| 国产 一区 欧美 日韩| 日本a在线网址| 国产乱人视频| 最近最新中文字幕大全电影3| 日韩欧美免费精品| 国产高潮美女av| 国产av一区在线观看免费| 我要搜黄色片| 老熟妇仑乱视频hdxx| 国产精品永久免费网站| 亚洲av一区综合| 国产高潮美女av| 尤物成人国产欧美一区二区三区| 深爱激情五月婷婷| 99久久久亚洲精品蜜臀av| 国内久久婷婷六月综合欲色啪| 免费大片18禁| 99久久久亚洲精品蜜臀av| 国内久久婷婷六月综合欲色啪| 日韩欧美免费精品| 日本免费一区二区三区高清不卡| 中文字幕人妻熟人妻熟丝袜美| 嫩草影院新地址| 国内少妇人妻偷人精品xxx网站| 亚洲天堂国产精品一区在线| а√天堂www在线а√下载| 国产爱豆传媒在线观看| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 日韩大尺度精品在线看网址| 日本 欧美在线| 亚洲国产精品sss在线观看| 麻豆一二三区av精品| 国产亚洲精品久久久com|