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

    High Throughput Scheduling Algorithms for Input Queued Packet Switches

    2022-11-09 08:16:28ChithraDeviJemiFlorinabelandNarayananPrasanth
    Computers Materials&Continua 2022年1期

    R.Chithra Devi,D.Jemi Florinabel and Narayanan Prasanth

    1Department of IT,Dr.Sivanthi Aditanar College of Engineering,Tiruchendur,628205,India

    2Department of CSE,Dr.Sivanthi Aditanar College of Engineering,Tiruchendur,628205,India

    3School of CSE,Vellore Institute of Technology,Vellore,632014,India

    Abstract: The high-performance computing paradigm needs high-speed switching fabrics to meet the heavy traffic generated by their applications.These switching fabrics are efficiently driven by the deployed scheduling algorithms.In this paper,we proposed two scheduling algorithms for input queued switches whose operations are based on ranking procedures.At first,we proposed a Simple 2-Bit(S2B)scheme which uses binary ranking procedure and queue size for scheduling the packets.Here,the Virtual Output Queue(VOQ)set with maximum number of empty queues receives higher rank than other VOQ’s.Through simulation,we showed S2B has better throughput performance than Highest Ranking First (HRF) arbitration under uniform,and non-uniform traffic patterns.To further improve the throughput-delay performance,an Enhanced 2-Bit(E2B)approach is proposed.This approach adopts an integer representation for rank,which is the number of empty queues in a VOQ set.The simulation result shows E2B outperforms S2B and HRF scheduling algorithms with maximum throughput-delay performance.Furthermore,the algorithms are simulated under hotspot traffic and E2B proves to be more efficient.

    Keywords: Crossbar switch;input queued switch;virtual output queue;scheduling algorithm;high performance computing

    1 Introduction

    In recent years,a lot of the commercial and scientific application requires high-performance computing (HPC).Modeling the environmental issues,product designing,materials,and biological research,managing Internet traffic are such example applications [1].HPC is a paradigm that solves problems that demand a high degree of computations in a shorter time with maximum accuracy [2].HPC aggregates resources from multiple machines or computes nodes to build a powerful infrastructure such as a high-end cluster,supercomputer,grid,or cloud to resolve complex problems [3].An HPC cluster consist of a large number of worker nodes and servers interconnected through a high-speed network using a Gigabit Ethernet preferably InfiniBand or Myrinet.The performance of the cluster largely depends on how fast the data communication was established between their components.As data are mostly residing in data centers,which are accessed by various resources in the cluster might create huge IP traffic around them.In some centers,Storage Area Networks (SAN) are used for carrying data between these cluster resources despite its implementation complexity.However,a widely preferred option for these data centers is to use high-speed switches [4].These switches can provide high-speed data transmission with maximum throughput and minimum latency on any kind of traffic [5].As the number of compute node get increases,the traffic around the data center too increase.As a result,the number of ports in the switch fabric has been increased however it is not a viable solution because it might minimize the data throughput [6].

    In order to design a fast as well as cheap switch,crossbar fabric is preferred.Here,between input and output ports all possible permutation of connections can be designed.A crossbar switch fabric can be classified as Input Queued (IQ),Output Queued (OQ),Combined Input and Output Queued (CIOQ) and Crosspoint Queued (CQ) switches [7].Among all the types of crossbars,input queued switch remains simple to implement and more effective because the memory bandwidth required for packet transmission at each timeslot is very less.This makes it highly preferred and most suitable for high-speed communications in Internet routers or data centers [7-9].The structure of the IQ switch is shown in Fig.1 and is discussed detailed in Section 3.A switch is driven by its arbitration scheme,which makes it more powerful.The performance of a switch is purely based on its employed scheduling algorithm.In IQ switch,iterative scheduling is preferred by various scheduling schemes such as iSLIP,RRM,etc.where scheduling is done in three phases: request,grant and accept.A matching request is sent by each input to the output during the request phase.One of the received requests is selected by the output during the grant phase.Respective Input will accept one of the available grants received from the output in the last phase.The same set of handshake operations will be performed at each iteration to establish a maximum matching.The drawback of the scheme,it is not possible to achieve maximum matching at every timeslot and hence suffers from average throughput performance.

    Some classes of algorithms such as DRRM,Grant Aware (GA) designed to work with two phases i.e.,request and grant [10].But still,a considerable amount of time is spent in finalizing the maximal matching,and hence the switching is not possible within a timeslot.In the case of HPC,switches are designed to handle huge traffic.Therefore the deployed arbitration scheme should be capable to find maximum matching at each timeslot in order to achieve maximum throughput-delay performance.For instance,on a 10 GPS line carrying 32-byte packets,5 ns is about to be one timeslot.A timeslot is allotted to make scheduling decisions as well as for data transmission.Completing all the iterations and making scheduling decisions within the minimum duration of 5 ns is a highly challenging task.This time-critical task offer severe headache to scheduling schemes used in high-speed switches.In some instances,the severity of task might lead to additional timeslot or in need of speedup.This will further reduce the switch performance or increase the implementation complexity,for the latter case.Therefore a high-performance scheduling algorithm is required to attain maximum matching at every timeslot to meet the demand of compute-intensive applications.

    Figure 1:Input queued switch with virtual output queue

    The rest of the paper is discussed as follows: Several iterative algorithms are discussed in the next section.In Section 3,we introduce the S2B and E2B algorithms with necessary examples.In Section 4,the performance of the proposed algorithms is evaluated and compared with existing schemes.In Section 5,we concluded the paper.

    2 Related Works

    One of the remarkable scheduling scheme in the history of IQ is iSLIP proposed by Nick McKeown [11],achieves 100% through under uniform traffic but offers poor performance under non-uniform traffic [11,12].It is an iterative scheduler which took three phases for every queue matching.A matching is computed in O (log n) iterations but largely suffers from longer queuing delay and matching overhead.Iterative Least Recently Used (iLRU) uses a similar procedure of iSLIP except for the choice of next scheduled cell.Here,the least recently used cell will be given maximum priority for the next schedule compared to RR fashion used in iSLIP.Parallel Iterative Matching (PIM) is an iterative maximal matching scheduler where the input-output matching is done in two phases i.e.,accept and grant by using random selection arbiters.It provides less than 65% of throughput under heavy load which is the limiting factor along with its high implementation cost and unfairness among the queues during switch schedule [13].Performance of iSLIP,iLRU and PIM are compared with 16 × 16 switch with single iteration as well as N iterations under Bernoulli arrivals.For single iteration,iSLIP provides maximum throughput performance because of its ability to desynchronize the output arbiters however for N iterations,all three schedulers can achieve maximum throughput irrespective of their poor delay performance.Providing N iterations for each schedule will introduce sever overhead issues.For irregular as well as bursty arrivals,their overall performance is not,as similar to earlier cases.

    Modified Round Robin Algorithm (MRRA) proposed by (Shanmugam) made an attempt to reduce the number of iterations at each timeslot to achieve maximum matching.The throughput performance of MRRA is similar to iSLIP under Bernoulli traffic however under bursty and their self-similar traffic model,they offer high delay.Authors preferred to use speedup in the future to improve the throughput and delay performance.On real-time applications,Dynamic Round Robin Matching (DRRM) is a widely used scheduler among high-speed switches because of its high arbitration rate.It is an iterative scheduler which took only two phases for every queue matching.It is commonly preferred because it is simple to implement,offer high efficiency and fairness.However,its throughput under heavy load is less than 50% which limits the scheduler usage in compute-intensive applications [14].

    To further reduce the scheduling overhead and complexity,few schemes are proposed with single iterations.Synchronous Round Robin (SRR) is the first of its kind,a distributed scheduler for Input Queued switches [15].It is a single bit single request algorithm which means only one request will be sent by an input during the request phase.Either global RR or longest queue strategy will be used to generate the 1-request.SRR performance under non-uniform traffic is average,hence not suited for high-speed networks.In [16],the author proposed three algorithms to support multiple server architecture.They are designed in such a way they can switch more than one fixed-sized cell in a timeslot.However,the hardware implementation and complexity seems to be stumbling block for practical.π-RGA is another single iteration algorithm that uses multiple bits during the request phase.The multi-bit is used to convey the history of information about the VOQ and is used while making scheduling decisions.It provides better performance under bursty traffic,however for other traffic,the performance is poor and its implementation is very difficult for high-speed networks.

    Another interesting single bit algorithm is the Round Robin with Longest Queue First(RR.LQF) [17].It uses a global RR scheduler and LQF scheduler to switch the cell from the longest queue.At every slot,a single bit is sent to the output indicating whether the VOQ receives a new cell or not.A counter is required at the output port to record the bit information however it increases the implementation as well as maintenance complexity.Some class of algorithms made an approach to introduce bonus schemes in order to improve the priority of waiting packets [18,19],thereby they avoided starvation and achieves better throughput delay performance.However,the scheme is most suited to buffered crossbar switches rather than IQ switch.

    As of now,most of the algorithm adopts a queue-based approach for scheduling the packets.But Highest Ranking First (HRF) uses a rank-based approach instead of using queue details for scheduling.It is one of the most successful single bit algorithm proposed by [20].In the basic model,every VOQ of the respective input is ranked based on the length of the queue.The rank is sent to the output during the request phase.The queue holding the highest rank will be granted permission by the output.On the accept phase,input accepts the received grant.If input receives more than one grant,the length of the queue will be considered to select a grant.In the next model,the RR scheduler is used to avoid packet delay,especially when working with a heavy load.These classes of the algorithm can offer better performance however the implementation complexity is relatively high.So in the third model,a three rank approach is used to rank all the VOQ.They made an attempt to reduce the complexity however it still looks very complicated for practical implementation in high-speed networks.HRF throughput-delay-performance under non-uniform traffic or hotspot traffic is not good.Moreover,their ranking procedure still based on queue length,and in most cases (under non-uniform traffic),it hugely suffers from packet delay which may lead to starvation.

    For example,consider a 4 × 4 switch in Fig.2,rank is calculated for all the VOQ in each input.(1,3,4,2) is the rank of the VOQ’s in the input 1 and is sent to the respective output.Similarly (2,1,3,4),(2,1,0,0) and (2,1,0,0) are the ranks of VOQ in inputs 2,-4 respectively.As per the HRF algorithm,output 1 receives the rank (1,2,2,2) from all the inputs and it grants the queue with the highest rank i.e.,VOQ1 is granted from input 1.This process continues for all the input matching,however for the given example,on a single iteration (in one timeslot) only 50% matching is possible.This is because inputs 1 and 2 receives two grants and they can execute only one.This situation continues in the further timeslots too,with a maximum matching ratio in single iteration per timeslot is less than 60%.This severely increases the delay ratio of switched packets and hence more attention is required.We propose a class of algorithms that enhances the ranking procedure to increase the input-output matching ratio as well as reduces the packet delay and thereby improves the schedule performance.

    Figure 2:A VOQ based input queued switch with ranking procedure

    3 Proposed Works

    We consider an N × N crossbar switch,a complete bipartite graph with N input ports and N output ports.A switch is said to be Input Queued if a queue is deployed on the input side as shown in Fig.1.Packets arriving from different sources are stored in an input queue,in specific,in their respective Virtual Output Queue (VOQ) [21,22].IQ fabric is integrated with VOQ’s on the input side.Each input i consist of a set of VOQ called as VOQ_Set (VS),where the packets are stored in a specific order.The number of queues in the VOQ_Set is the same as the number of Output ports and packets stored on the first queue of VOQ1 is destined only to a first output port and so on.Throughout our work,it is assumed that we use only fixed-sized packets (or cells)for scheduling and are reassembled on the output side.Packets are selected from each VOQ_Set based on the employed scheduling scheme.At each timeslot,the matched packets are switched from the input port to the output port of the crossbar.We made sure that the transmission speed of the switch fabric (line rate) as well as on the input and output ports are the same i.e.,1 cell per slot.The input to the switch can be generated through various arrivals.Whenever a packet arrives,its destination information can be identified from its header.Then the packet is directed to the respective input queue for scheduling purposes.

    3.1 A Simple 2-Bit(S2B)Approach

    As HRF failed to consider the empty VOQ,switch suffers from latency issues.We proposed a Simple Two Bit scheduling approach for IQ switches.This simple scheme considers empty VOQ as a vital factor for the ranking process.Like other iterative algorithms,this scheme works in three phases.During the request phase,VOQ in each input sends two bits (R,P) to output.R is the rank of the VOQ_Set,which is a binary value,states whether the respective VOQ_Set has empty VOQ or not (say 0 or 1,0-non-empty and 1-empty queue) and P is an integer which specifies the number of packets in the Queue.In the grant phase,R-bit is compared to ensure the input with empty queues got the preference.That is,input-R bit with value 1 will be given preference for the grant,in case,if R-bit is the same for more than one input,then P value will be considered i.e.,input with maximum P value will be granted.If still not able to grant,the RR scheme will be used to grant an input.During the accept phase,input accepts the grant and gets ready for switching.If an input receives more than one grant from the output,it accepts the grant which has a larger P-bit value.

    Figure 3:A simple 2-bit arbitration scheme

    For example,consider the 4 × 4 switch in Fig.3,inputs is equipped with VOQ which has fixed-sized packets for switching.In the request phase,VOQ1 of the input 1 send the value pair(0,5),similarly,output 1 receives 4 pair of VOQ values [(0,5),(0,4),(1,2),(1,3)] from VOQ1 of all the input.In the grant phase,at step 1,output considers the VOQ 1 from inputs 3 and 4[(1,2),(1,3) alone because of R-bit.For example,consider the 4 × 4 switch in Fig.3,inputs is equipped with VOQ which has fixed-sized packets for switching.In the request phase,VOQ1 of the input 1 send the value pair (0,5),similarly,output 1 receives 4 pair of VOQ values [(0,5),(0,4),(1,2),(1,3)] from VOQ1 of all the input.In the grant phase,at step 1,output considers the VOQ 1 from inputs 3 and 4 [(1,2),(1,3) alone because of R-bit.In step 2,P-bit is compared,and therefore input 4 i.e.,(1,3) is granted.In the accept phase,input 4 accepts the grant of output 1,and matching is confirmed (I4->O1).Similarly,other matchings in the same slot are[(I1->O4),(I3->O2)] at iteration1,and (I2->O3) at iteration 2.This shows a maximum matching is highly possible at every timeslot through this simple 2-bit approach.

    Algorithm

    Algorithm 1: Simple Two Bit Approach 1: procedure Request (input)2: R=VOQ set status of input i 3: Compute(R)4: Calculate P[j]=length of VOQ[j] at input i 5: if P[j] >0 then 6: send (R,P[j]) to output i 7: end if 8: end procedure 9: procedure Grant (rank,length)10: Gi[r,l]=request from input i 11: if Gi[r]==1 then 12: Check the rank of other inputs with same status 13: If true,send grant to input with max (Gi[l])14: else send grant to input i 15: else if Gi[r]==0 then 16: Check the rank of other inputs with same status 17: If true,send grant to input with max (Gi[l])18: else send grant to input i 19: end if 20: end procedure 21: procedure Accept (out)22: A[i]=No.of grant from output i 23: if A >1 then 24: Max(P[j]) is accepted 25: Reset other input VOQ [j] to 0 26: else if A==1 then 27: accept the grant from output i 28: Reset other input VOQ [j] to 0 29: else 30: call next iteration 31: Compute (R)32: Calculate P[j]=length of VOQ [j] at input i(Continued)

    33: end if 34: e nd procedure 35: procedure Compute(R)36: Check for empty queue-if true return 1 else 0 37: end procedure

    3.2 Enhanced 2-Bit(E2B)Approach

    A simple 2-bit approach can improve the latency performance compared to the HRF algorithm.S2B uses ranking procedures as well as queue lengths to make scheduling decisions.Ranks used are binary values which imply whether a VOQ_Set has empty queues or not.The ranking procedure can be further enhanced by considering the number of empty queues in the VOQ_SET by using an integer representation instead of binary.The rank defined by this integer indicates the number of empty queues in the VOQ_SET.Consider Fig.4,VOQ1 of input 1 send value pairs(1,5) to output in the request phase.Here,VOQ1 rank is 1 and length is 5.Similarly,VOQ1 of inputs 2-4 send the value pairs [(1,4),(2,4),(3,3)] respectively to output1.In the grant phase,the rank of the VOQ (R-bit) is compared,input having the highest rank will be preferred.If more than one VOQ holds a similar rank then P-bit is compared and the input having a larger P-bit value will be granted.

    Figure 4:Enhanced 2-bit arbitration scheme

    In our example,output1 grant input 4 (I4->O1) and is accepted in the last phase.However,if an input receives more than one grant,input checks for larger P-bit and accepts the respective output.In our example,already output 1 is granted for input 4,output 2 grants the input 3,output 3 grants Input 2 and output 4 grants input 1.For the given example,the maximum matching [(I1->O4),(I2->O3),(I3->O2),(I4->01)] is achieved with a single iteration of a timeslot.However,in some cases maximum matching is possible with two iterations.The time complexity of E2B is O(N).The proposed arbitration schemes are better suited to IQ switches than the existing schemes because of the following:

    · Probability to achieve maximum matching is very high because it prioritizes the VOQ_Set with maximum empty queues for scheduling

    · Matching is possible in the first iteration itself (most of the cases) of a timeslot,hence reduces the overhead issue

    · It is designed in such a way to achieve maximum throughput performance under nonuniform traffic where other algorithms failed.

    · Most suited to real time high-performance applications.

    Algorithm

    Algorithm 2: Enhanced Two Bit Approach 1: procedure Request (input)2: R=VOQ_Set status of input i 3: Compute (R)4: Calculate P[j]=length of VOQ[j] at input i 5: if P[j] >0 then 6: send (R,P[j]) to output i 7: end if 8: end procedure 9: procedure Grant (rank,length)10: Gi[r,l]=request from input i of 11: Compute m=max (Gi[r])12: c=number of occurrence of m in Gi[r]13: If (c==1) then 14: grant corresponding input i 15: else 16: Compute max(Gi [l]) from m and grant corresponding input i 17: end if 18: end procedure 19: procedure Accept (out)20: A[i]=No.of grant from output i 21: if A >1 then 22: accept the grant with Max (P[j])23: Reset other input VOQ [i] to 0 24: else if A==1 then 25: accept the grant from output i 26: Reset other input VOQ [i] to 0 27: else 28: call next iteration 29: Compute(R)30: Calculate P[j]=length of VOQ [j] at input i 31: end if 32: end procedure(Continued)

    33: procedure Compute(R)34: R=No.of empty VOQ in input i 35: return R 36: end procedure?

    4 Performance Analysis

    4.1 Simulation Setup

    A simulation environment is developed to measure the throughput and delay performance of S2B and E2B scheduling schemes.Both the algorithms executed with a load 10,000 fixedsized packets under various traffic patterns and their performances are noted.The performance of S2B and E2B are compared with the HRF scheme with respect to throughput and delay metrics.Throughout our simulation,a 4 × 4 VOQ based IQ crossbar switch is used with 4 input and 4 output ports.In order to avoid packet drop on the input side,VOQ is designed with infinite buffer size.The scheduler is tested with various load sizes and the traffic patterns (Bernoulli) used are independent of each other.Uniform,non-uniform,and hotspot are three types of workload generated for the simulations.Uniform traffic ensures the workload remains the same for all the output ports at each simulation however those packets can be distributed through any VOQ.In the case of non-uniform traffic,for every simulation,the load size remains different from every output ports.

    4.2 Simulation Results

    At first,we investigate the throughput performance of HRF,S2B,and E2B scheduling schemes under uniform traffic distribution.The throughput is analysed under various load strengths (ranging from 10% to 100%) and the result is depicted in Fig.5.At 40% of the load,HRF managed to deliver 85% of throughput,however,it drops below 80% when maximum load is offered.The effect of prioritizing the empty VOQ during the scheduling process provides a considerable impact on S2B throughput performance,which manages a 5% improvement than HRF at all load levels.

    Further enhancement in S2B leads to E2B which achieves a maximum throughput of 100%for all the load levels.E2B is also tested with various uniform workloads and it achieves 100%in all the cases.This is achieved because it provides 100% matching in most of the switching cases with one or more iterations in a slot.Fig.6 shows the throughput performance of HRF,S2B,and E2B under non-uniform traffic.Until 40% of the load,HRF and S2B provide similar performance,however,when the load increases,the difference in throughput reaches at a maximum of 10%.This is the point where most of the schedulers struggle to deliver an acceptable throughput when the load get increases (>50%) under non-uniform traffic.However,it is understood that the decrease in throughput performance is due to minimal matching at each slot.If the amount of matching can be improved then the overall performance will increase.E2B is designed in such a way to increase the number of matchings at each slot.Fig.6 also shows the throughput performance of E2B which is approximately 95% when maximum load is offered.E2B is simulated with various non-uniform load structures (up to 100,000 packets) and in all the cases its throughput performance is not less than 94% however for few load structures it achieves 100% throughput too.To our knowledge,this is the first IQ scheduler which can attain 100%throughput under certain non-uniform load patterns.

    Figure 5:Throughput performance of HRF,S2B and E2B under uniform traffic

    Figure 6:Throughput performance of HRF,S2B,and E2B under non-uniform traffic

    In a switch fabric,a packet can arrive and depart within the same slot in the best possible case for delay analysis.However,it is not possible for all the packets in real-time scenarios.Delay measure is defined as the total number of timeslots spent by a packet in the switch or,it is the difference between the output port times to the input port time of a packet.Fig.7 shows the delay performance of HRF,S2B,and E2B scheduling schemes under uniform traffic.All three schemes provide similar delay performance under any workload.This is because all these schemes adopted ranking based approaches and hence they attain similar performances.Moreover,there is no much difference to separate these schemes under non-uniform traffic too,and is depicted in Fig.8.Therefore it is understood the change of approach in the ranking process by S2B and E2B has no delay improvement over HRF.

    4.3 Hotspot Traffic

    In the real world high-performance computing environment,the study shows there are numerous situations where traffic is targeted at a specific destination.In an interconnection network,numerous application faces hotspot traffic,for example,multiple clients access a server is one such scenario that will generate a lot of traffic in and around the ports.For summer vacation,a number of students from a particular lab is attempting to reserve an airline ticket at a discounted price at a specific period of time is a perfect scenario for hotspot traffic.When n number of inputs targeting a single output port at a given period of time,then the traffic occurred in the crossbar interconnection is said to be hotspot traffic.From the study,we understood that most of the scheduling schemes are not simulated with hotspot traffic however we analysed the throughput and delay performance of HRF,S2B,and E2B algorithms.For our simulation,uniform traffic is generated with a hotspot fraction f=0.8 and their throughput performance is depicted in Fig.9.

    Figure 7:Delay performance of HRF,S2B and E2B under uniform traffic

    Figure 8:Delay performance of HRF,S2B and E2B under non-uniform traffic

    The distribution of load has less impact on the overall throughput performance of S2B and HRF schemes.With maximum load,the throughput performance of HRF is 70% and S2Bc is 65% however the method of prioritizing the queue with empty slots for switching gives significant improvement.Therefore E2B achieves 95% throughput when maximum load is offered.We conducted numerous simulations by changing the load structures and hotspot fraction however the minimum achieved throughput at any scenario is 95%.To our knowledge,this is the maximum achieved throughput by any IQ scheduler under hotspot traffic.Fig.10 shows the delay performance of HRF,S2B,and E2B under hotspot traffic.Up to 50% of the load,all three algorithms have similar delay however E2B performs better when the load exceeds.At maximum load,E2B is comparatively better than S2B and HRF by a minimum of 5%.

    In summary,E2B outperforms HRF and S2B with maximum throughput under uniform,nonuniform and hotspot traffic.E2B delay performance under uniform and non-uniform is almost similar to S2B and HRF however under hotspot traffic it offers better delay.

    Figure 9:Throughput performance of HRF,S2B and E2B under hotspot traffic

    Figure 10:Delay performance of HRF,S2B and E2B under hotspot traffic

    5 Conclusion

    In this paper,we first proposed a simple 2-bit scheduling algorithm for input queued switches.It uses a simple ranking procedure and VOQ length for scheduling the packets.The simulation result shows S2B offers better throughput performance for all the tested traffic patterns.To further improve the performance,an extended 2-bit scheme is proposed and simulated under various traffic patterns.We showed that E2B throughput-delay performance is far better than S2B and HRF under heavy loads and hence it is most suited for high-performance computing environments.As a future work,scheduling schemes can be designed in a non-iterative fashion to reduce the overhead involved.

    Funding Statement: The authors received no specific funding for this study.

    Conflicts of Interest: The authors declare that they have no conflicts of interest to report regarding the present study.

    欧美日韩黄片免| 日韩视频一区二区在线观看| 国产精品成人在线| 777米奇影视久久| 亚洲av成人一区二区三| 欧美日韩黄片免| 亚洲视频免费观看视频| 永久免费av网站大全| 日韩一卡2卡3卡4卡2021年| 91老司机精品| 999久久久精品免费观看国产| 日韩欧美一区视频在线观看| 色老头精品视频在线观看| 亚洲精品成人av观看孕妇| 午夜福利视频在线观看免费| 99精国产麻豆久久婷婷| 国产熟女午夜一区二区三区| 久久热在线av| 999久久久精品免费观看国产| 国产不卡av网站在线观看| 国产福利在线免费观看视频| 久久久久视频综合| av电影中文网址| 亚洲成人免费电影在线观看| 亚洲三区欧美一区| 亚洲天堂av无毛| 首页视频小说图片口味搜索| 2018国产大陆天天弄谢| 天堂8中文在线网| 国产免费一区二区三区四区乱码| 麻豆乱淫一区二区| 美女国产高潮福利片在线看| 午夜免费成人在线视频| 蜜桃在线观看..| 欧美日韩中文字幕国产精品一区二区三区 | 亚洲 国产 在线| 久久久久久久大尺度免费视频| 午夜福利视频在线观看免费| 亚洲国产精品999| 免费高清在线观看日韩| 美女大奶头黄色视频| 国产一区二区三区av在线| 桃花免费在线播放| 久久午夜综合久久蜜桃| 王馨瑶露胸无遮挡在线观看| 亚洲欧美精品综合一区二区三区| 在线亚洲精品国产二区图片欧美| 在线观看一区二区三区激情| 欧美黄色片欧美黄色片| 国产av精品麻豆| 亚洲一区二区三区欧美精品| 菩萨蛮人人尽说江南好唐韦庄| 多毛熟女@视频| 真人做人爱边吃奶动态| 亚洲国产精品成人久久小说| 成年av动漫网址| 成年动漫av网址| 黄色视频不卡| 交换朋友夫妻互换小说| 日韩有码中文字幕| 欧美日韩亚洲国产一区二区在线观看 | 亚洲精品在线美女| videos熟女内射| a级片在线免费高清观看视频| 99精国产麻豆久久婷婷| 这个男人来自地球电影免费观看| 国产麻豆69| 成人三级做爰电影| 啦啦啦中文免费视频观看日本| 欧美精品啪啪一区二区三区 | 一个人免费在线观看的高清视频 | 看免费av毛片| 国产精品免费视频内射| 99久久人妻综合| 99精品欧美一区二区三区四区| 国产精品av久久久久免费| 久久精品亚洲av国产电影网| 亚洲免费av在线视频| 久久影院123| 欧美激情久久久久久爽电影 | 亚洲精品久久午夜乱码| 岛国在线观看网站| 亚洲精品一区蜜桃| 亚洲国产精品999| 黄片大片在线免费观看| 亚洲国产欧美在线一区| 精品久久久精品久久久| 好男人电影高清在线观看| 精品乱码久久久久久99久播| av网站免费在线观看视频| av片东京热男人的天堂| 91老司机精品| 亚洲精品成人av观看孕妇| 亚洲欧美色中文字幕在线| 欧美激情高清一区二区三区| 亚洲人成电影免费在线| 亚洲精品国产色婷婷电影| 9191精品国产免费久久| 飞空精品影院首页| 在线精品无人区一区二区三| 99国产精品一区二区蜜桃av | 久久女婷五月综合色啪小说| 黄片大片在线免费观看| 大码成人一级视频| 精品少妇久久久久久888优播| 秋霞在线观看毛片| videos熟女内射| 亚洲中文av在线| 午夜福利在线免费观看网站| 久久国产精品大桥未久av| 国产成人精品在线电影| 久久人妻熟女aⅴ| 人成视频在线观看免费观看| 日韩有码中文字幕| 免费av中文字幕在线| 一本综合久久免费| 999久久久精品免费观看国产| 精品国产一区二区三区四区第35| 国产深夜福利视频在线观看| 欧美精品亚洲一区二区| 亚洲情色 制服丝袜| 最新在线观看一区二区三区| av在线播放精品| 首页视频小说图片口味搜索| 中亚洲国语对白在线视频| 亚洲精品美女久久av网站| videos熟女内射| 午夜免费成人在线视频| 香蕉丝袜av| av福利片在线| 国产av国产精品国产| 黑人欧美特级aaaaaa片| 夜夜骑夜夜射夜夜干| 一区二区三区激情视频| av在线老鸭窝| 国产野战对白在线观看| 性高湖久久久久久久久免费观看| 欧美日韩一级在线毛片| 欧美成狂野欧美在线观看| 黄色视频不卡| 美女高潮喷水抽搐中文字幕| 成在线人永久免费视频| 69av精品久久久久久 | 黄色视频,在线免费观看| 热99久久久久精品小说推荐| 90打野战视频偷拍视频| 亚洲av欧美aⅴ国产| 成人三级做爰电影| 国产在视频线精品| 欧美日韩av久久| 人妻 亚洲 视频| 精品国产乱码久久久久久男人| 国产一区二区三区在线臀色熟女 | 激情视频va一区二区三区| 日日爽夜夜爽网站| 免费不卡黄色视频| 美女大奶头黄色视频| 99re6热这里在线精品视频| 亚洲专区中文字幕在线| 国产在视频线精品| 国产激情久久老熟女| 免费高清在线观看视频在线观看| 国产一区二区激情短视频 | 嫁个100分男人电影在线观看| 日本黄色日本黄色录像| 日本91视频免费播放| 水蜜桃什么品种好| 精品第一国产精品| 久久久精品国产亚洲av高清涩受| 免费女性裸体啪啪无遮挡网站| 亚洲成人国产一区在线观看| 韩国精品一区二区三区| 老熟妇乱子伦视频在线观看 | 一区二区三区乱码不卡18| 国产激情久久老熟女| 午夜精品久久久久久毛片777| 国产伦人伦偷精品视频| 高清视频免费观看一区二区| 老熟女久久久| 亚洲精品第二区| 亚洲av成人一区二区三| 精品一区在线观看国产| 热99re8久久精品国产| 国产成人精品久久二区二区91| 99国产精品免费福利视频| 国产一区有黄有色的免费视频| 丰满饥渴人妻一区二区三| 国产免费一区二区三区四区乱码| 一级毛片电影观看| 国产国语露脸激情在线看| 男女免费视频国产| 美女大奶头黄色视频| 精品国产乱码久久久久久男人| 欧美日本中文国产一区发布| tube8黄色片| 每晚都被弄得嗷嗷叫到高潮| 色老头精品视频在线观看| 中文字幕av电影在线播放| 欧美日韩亚洲国产一区二区在线观看 | 少妇被粗大的猛进出69影院| 大香蕉久久网| 亚洲国产日韩一区二区| videos熟女内射| 建设人人有责人人尽责人人享有的| 大陆偷拍与自拍| 巨乳人妻的诱惑在线观看| 久久久久久人人人人人| 亚洲一区中文字幕在线| 亚洲色图综合在线观看| 亚洲avbb在线观看| 操出白浆在线播放| a级毛片在线看网站| 色94色欧美一区二区| 亚洲第一青青草原| 无遮挡黄片免费观看| 天堂俺去俺来也www色官网| 亚洲欧美清纯卡通| 啦啦啦视频在线资源免费观看| svipshipincom国产片| 成人国产av品久久久| 久久久欧美国产精品| 精品人妻一区二区三区麻豆| 建设人人有责人人尽责人人享有的| 不卡av一区二区三区| 国精品久久久久久国模美| 欧美精品av麻豆av| 日韩制服骚丝袜av| 又大又爽又粗| 国产男女超爽视频在线观看| 午夜福利,免费看| 亚洲国产看品久久| 自线自在国产av| 国产精品久久久久久精品古装| 久久久久精品人妻al黑| 女人久久www免费人成看片| 国产免费av片在线观看野外av| 一个人免费看片子| 日韩三级视频一区二区三区| 99热国产这里只有精品6| 国产高清国产精品国产三级| 中文字幕高清在线视频| 久久久精品区二区三区| 欧美黄色片欧美黄色片| 亚洲国产av新网站| 亚洲欧美清纯卡通| 国产精品久久久久久精品古装| 搡老岳熟女国产| 啪啪无遮挡十八禁网站| 777米奇影视久久| 免费日韩欧美在线观看| 亚洲精品一区蜜桃| 亚洲av欧美aⅴ国产| 两个人看的免费小视频| 男女午夜视频在线观看| 亚洲欧美一区二区三区黑人| 欧美另类亚洲清纯唯美| 亚洲欧美激情在线| 韩国精品一区二区三区| 天天影视国产精品| 国产欧美亚洲国产| 日韩视频一区二区在线观看| 美女大奶头黄色视频| 亚洲国产av影院在线观看| 久久99热这里只频精品6学生| 欧美人与性动交α欧美软件| 精品人妻一区二区三区麻豆| 99热国产这里只有精品6| 亚洲五月色婷婷综合| 国产av精品麻豆| 久久青草综合色| 午夜福利乱码中文字幕| 可以免费在线观看a视频的电影网站| 亚洲精品国产区一区二| 男女下面插进去视频免费观看| 国产成人精品久久二区二区免费| 国产在线免费精品| 免费在线观看日本一区| 亚洲国产精品999| 亚洲精品国产一区二区精华液| 亚洲精品一区蜜桃| 国产精品秋霞免费鲁丝片| 999久久久精品免费观看国产| 久久久久久久久久久久大奶| 欧美+亚洲+日韩+国产| 免费观看av网站的网址| av电影中文网址| 国产精品九九99| 正在播放国产对白刺激| 在线永久观看黄色视频| 亚洲精品国产一区二区精华液| 99久久人妻综合| 麻豆国产av国片精品| 亚洲五月色婷婷综合| 久久99一区二区三区| 男男h啪啪无遮挡| 99国产精品免费福利视频| 男女高潮啪啪啪动态图| 一区福利在线观看| 国产男人的电影天堂91| 一个人免费看片子| 国产精品亚洲av一区麻豆| 婷婷色av中文字幕| 久久九九热精品免费| 老司机在亚洲福利影院| av网站免费在线观看视频| 亚洲专区字幕在线| 久久女婷五月综合色啪小说| 亚洲成人手机| 性高湖久久久久久久久免费观看| 日日爽夜夜爽网站| av有码第一页| 真人做人爱边吃奶动态| videosex国产| 脱女人内裤的视频| 亚洲男人天堂网一区| 久久久久精品国产欧美久久久 | 久久性视频一级片| 激情视频va一区二区三区| 亚洲激情五月婷婷啪啪| 欧美精品亚洲一区二区| 色精品久久人妻99蜜桃| 在线 av 中文字幕| 久久精品国产亚洲av高清一级| 91九色精品人成在线观看| 女性生殖器流出的白浆| 热re99久久国产66热| 久久久久精品人妻al黑| 亚洲色图综合在线观看| av有码第一页| 大片免费播放器 马上看| 一本一本久久a久久精品综合妖精| 一进一出抽搐动态| 在线 av 中文字幕| 在线观看www视频免费| 精品一品国产午夜福利视频| 电影成人av| 国产不卡av网站在线观看| 大型av网站在线播放| 亚洲精品国产色婷婷电影| 国产高清视频在线播放一区 | 真人做人爱边吃奶动态| 久久狼人影院| 在线观看免费午夜福利视频| 欧美国产精品va在线观看不卡| 欧美黑人欧美精品刺激| 国产成人欧美在线观看 | 高清av免费在线| 国产无遮挡羞羞视频在线观看| 成在线人永久免费视频| √禁漫天堂资源中文www| 黄色 视频免费看| 最近最新中文字幕大全免费视频| 国产高清videossex| 好男人电影高清在线观看| 伊人亚洲综合成人网| 国产片内射在线| 亚洲国产精品成人久久小说| 国产精品香港三级国产av潘金莲| 视频区图区小说| 建设人人有责人人尽责人人享有的| 亚洲av国产av综合av卡| 动漫黄色视频在线观看| 欧美人与性动交α欧美精品济南到| 国产一级毛片在线| 99精品久久久久人妻精品| 老司机靠b影院| 久久精品亚洲av国产电影网| 老司机靠b影院| 久久久国产成人免费| 999久久久精品免费观看国产| 久久九九热精品免费| 久久久久久免费高清国产稀缺| 91国产中文字幕| 另类精品久久| 国产极品粉嫩免费观看在线| 久久人人97超碰香蕉20202| 欧美日韩福利视频一区二区| 好男人电影高清在线观看| 国产精品秋霞免费鲁丝片| 99久久人妻综合| 99久久精品国产亚洲精品| 十八禁高潮呻吟视频| 亚洲国产av新网站| 波多野结衣av一区二区av| 成年女人毛片免费观看观看9 | 91九色精品人成在线观看| 久久中文字幕一级| a级毛片黄视频| 电影成人av| 秋霞在线观看毛片| 国产欧美日韩一区二区精品| 日韩制服丝袜自拍偷拍| 爱豆传媒免费全集在线观看| 日韩精品免费视频一区二区三区| 国产人伦9x9x在线观看| 高清视频免费观看一区二区| 久久国产亚洲av麻豆专区| 一个人免费在线观看的高清视频 | 啦啦啦 在线观看视频| 少妇猛男粗大的猛烈进出视频| 国产精品.久久久| 9191精品国产免费久久| 狠狠狠狠99中文字幕| 老司机在亚洲福利影院| 91国产中文字幕| 美女视频免费永久观看网站| 一进一出抽搐动态| 一个人免费看片子| svipshipincom国产片| 1024香蕉在线观看| 国产av国产精品国产| 久久女婷五月综合色啪小说| 不卡av一区二区三区| 亚洲黑人精品在线| 在线观看www视频免费| 久久香蕉激情| 精品亚洲成国产av| 久久久国产一区二区| 国产精品久久久人人做人人爽| 亚洲精品国产色婷婷电影| 高清视频免费观看一区二区| 中文字幕高清在线视频| 国产免费福利视频在线观看| 99国产精品一区二区三区| 亚洲精品中文字幕一二三四区 | 他把我摸到了高潮在线观看 | 蜜桃国产av成人99| 少妇人妻久久综合中文| 天天影视国产精品| 中文字幕最新亚洲高清| 国产三级黄色录像| 99国产精品一区二区三区| av超薄肉色丝袜交足视频| 一区二区三区激情视频| 一区二区av电影网| 国产精品久久久av美女十八| 黑人猛操日本美女一级片| 人妻一区二区av| 少妇猛男粗大的猛烈进出视频| 亚洲九九香蕉| 欧美日韩av久久| 国产黄频视频在线观看| 黑人巨大精品欧美一区二区蜜桃| 亚洲国产欧美在线一区| 亚洲,欧美精品.| 黄色毛片三级朝国网站| 香蕉丝袜av| 午夜激情av网站| 国产精品成人在线| 人人妻,人人澡人人爽秒播| 91麻豆av在线| 一本一本久久a久久精品综合妖精| 日韩视频一区二区在线观看| 我要看黄色一级片免费的| 久久久久久久久久久久大奶| 美女视频免费永久观看网站| 久久久久久亚洲精品国产蜜桃av| 人妻 亚洲 视频| 黄色怎么调成土黄色| 免费高清在线观看日韩| a 毛片基地| 日本av免费视频播放| av线在线观看网站| 午夜久久久在线观看| 欧美黑人欧美精品刺激| 天天躁日日躁夜夜躁夜夜| 中文字幕另类日韩欧美亚洲嫩草| 亚洲精品成人av观看孕妇| av福利片在线| 啦啦啦免费观看视频1| 国产片内射在线| 欧美激情极品国产一区二区三区| 国产成+人综合+亚洲专区| 男女边摸边吃奶| 女人被躁到高潮嗷嗷叫费观| 狂野欧美激情性xxxx| 久久天堂一区二区三区四区| 精品高清国产在线一区| 交换朋友夫妻互换小说| 黄色 视频免费看| 十分钟在线观看高清视频www| 午夜福利乱码中文字幕| 窝窝影院91人妻| 免费女性裸体啪啪无遮挡网站| 亚洲五月色婷婷综合| 欧美日韩亚洲综合一区二区三区_| 岛国在线观看网站| 亚洲天堂av无毛| 亚洲伊人色综图| 精品一区二区三区四区五区乱码| 国产深夜福利视频在线观看| 亚洲七黄色美女视频| bbb黄色大片| 亚洲专区国产一区二区| 老鸭窝网址在线观看| 亚洲专区国产一区二区| 桃红色精品国产亚洲av| 亚洲成人免费电影在线观看| 日韩大码丰满熟妇| 秋霞在线观看毛片| 一级,二级,三级黄色视频| 777久久人妻少妇嫩草av网站| 欧美在线一区亚洲| 欧美激情久久久久久爽电影 | 国产在线一区二区三区精| 国产成人一区二区三区免费视频网站| 男人舔女人的私密视频| 国产熟女午夜一区二区三区| 日韩一区二区三区影片| 视频在线观看一区二区三区| a 毛片基地| a在线观看视频网站| 悠悠久久av| 久久精品久久久久久噜噜老黄| 妹子高潮喷水视频| 免费高清在线观看视频在线观看| 在线十欧美十亚洲十日本专区| 午夜福利在线免费观看网站| 热re99久久国产66热| 亚洲性夜色夜夜综合| 午夜福利视频精品| 国产精品一区二区在线观看99| 亚洲欧美一区二区三区久久| 大陆偷拍与自拍| av不卡在线播放| 国产成人精品无人区| 国产精品99久久99久久久不卡| 免费看十八禁软件| 国产亚洲一区二区精品| 99热国产这里只有精品6| 午夜福利视频在线观看免费| 91麻豆av在线| av欧美777| svipshipincom国产片| 另类亚洲欧美激情| 国产一区二区 视频在线| 人人妻人人爽人人添夜夜欢视频| 免费在线观看日本一区| 女人久久www免费人成看片| 免费一级毛片在线播放高清视频 | 国产真人三级小视频在线观看| 法律面前人人平等表现在哪些方面 | 黄频高清免费视频| 欧美 亚洲 国产 日韩一| 狠狠狠狠99中文字幕| 欧美日韩精品网址| 最新的欧美精品一区二区| 丰满饥渴人妻一区二区三| 免费观看a级毛片全部| 好男人电影高清在线观看| 99热网站在线观看| 午夜影院在线不卡| 亚洲国产中文字幕在线视频| 久久天堂一区二区三区四区| 淫妇啪啪啪对白视频 | 菩萨蛮人人尽说江南好唐韦庄| 亚洲成人免费av在线播放| 成人国产av品久久久| 亚洲精品一区蜜桃| 一本—道久久a久久精品蜜桃钙片| 精品一区二区三卡| av一本久久久久| 99re6热这里在线精品视频| 真人做人爱边吃奶动态| 精品国产一区二区三区久久久樱花| 久久国产精品影院| 欧美黄色淫秽网站| 一本一本久久a久久精品综合妖精| 超碰成人久久| 亚洲天堂av无毛| 亚洲精品久久成人aⅴ小说| 黄色毛片三级朝国网站| 亚洲一区二区三区欧美精品| 国产亚洲av片在线观看秒播厂| 久久国产精品人妻蜜桃| 肉色欧美久久久久久久蜜桃| 妹子高潮喷水视频| 欧美另类亚洲清纯唯美| 国产精品一区二区在线观看99| 亚洲中文字幕日韩| 婷婷丁香在线五月| 午夜福利免费观看在线| 天堂俺去俺来也www色官网| 精品欧美一区二区三区在线| 国产亚洲欧美在线一区二区| 秋霞在线观看毛片| 日韩中文字幕视频在线看片| 亚洲中文av在线| 国产日韩欧美视频二区| 精品久久蜜臀av无| 亚洲av日韩精品久久久久久密| 麻豆乱淫一区二区| 777米奇影视久久| 久久精品人人爽人人爽视色| 黄网站色视频无遮挡免费观看| 777米奇影视久久| videosex国产| 亚洲精品在线美女| 免费久久久久久久精品成人欧美视频| 精品久久久久久久毛片微露脸 | 亚洲精品中文字幕一二三四区 | 黄色视频,在线免费观看| 超碰97精品在线观看| 欧美中文综合在线视频| 美女脱内裤让男人舔精品视频| 悠悠久久av| 久久久久精品国产欧美久久久 | 亚洲av成人一区二区三| 久久九九热精品免费| 久久精品aⅴ一区二区三区四区| 91九色精品人成在线观看| 国产熟女午夜一区二区三区| 久久毛片免费看一区二区三区| 亚洲自偷自拍图片 自拍| 视频区图区小说|