On the Control of Fork-Join Networks

Size: px
Start display at page:

Download "On the Control of Fork-Join Networks"

Transcription

1 On the Contol of Fok-Join Netwoks Ehun Özkan1 and Amy R. Wad 2 Abstact In a fok-join pocessing netwok, jobs aive, and then fok into tasks, some of which can be pocessed sequentially and some of which can be done in paallel, accoding to a set of deteministic pecedence constaints. Befoe a job can depat the netwok, the tasks must be joined togethe, which gives ise to synchonization constaints in the netwok. When thee ae multiple job types that shae esouces, the contol decisions occu at any seve that can pocess moe than one job type. We focus on a pototypical fok-join netwok with two job classes, two fok opeations, one shaed seve, two join opeations, and holding costs. The fok opeations ae fist, followed by the simultaneous pocessing of type a b jobs by a dedicated seve and a shaed seve, and, finally, a join opeation. We solve the contol poblem fo this netwok in an asymptotic egime in which the shaed seve is in heavy taffic. We show that a cµ-type static pioity policy is asymptotically optimal when the shaed seve is in some sense slow at pocessing the moe expensive type a jobs. Othewise, an asymptotically optimal contol is a state-dependent slow depatue pacing contol unde which the shaed seve pioitizes the type b jobs wheneve the numbes of type a jobs waiting to be joined in the downsteam buffes ae balanced. Finally, by consideing a boade class of fok-join netwoks, we see that the depatue pacing idea seems to be obust. Keywods: Fok-join pocessing netwok; Scheduling Contol; Diffusion Appoximation; Asymptotic Optimality. AMS Classification: Pimay 60K25, 90B22, 90B36, 93E20; Seconday 60F17, 60J70. 1 Intoduction Netwoks in which pocessing of jobs occus both sequentially and in paallel ae pevalent in many application domains, such as compute systems Xia et al. [2007], healthcae Amony et al. [2015], manufactuing Dalley and Geshwin [1992], poject management Adle et al. [1995], and the justice system Lason et al. [1993]. The paallel pocessing of jobs gives ise to synchonization constaints that can be a main eason fo job delay. Although delays in fok-join netwoks can be appoximated unde the common fist-come-fist-seved FCFS scheduling discipline Nguyen [1993, 1994], thee is no eason to believe FCFS scheduling minimizes delay. Ou objective in this pape is to devise contol policies that minimize delay o, moe geneally, holding costs in fok-join netwoks with multiple custome classes that shae pocessing esouces. To do this, we caefully study the pototype netwok in Figue 1. Fom that, we lean how to develop good contol policies fo a boade class of fok-join netwoks. The fok-join netwok in Figue 1 seves to illustate why fok-join netwok contol is difficult. In that netwok, thee ae two aiving job types a and b, seven seves numbeed 1 to 7, and ten buffes numbeed 1 to 10. We assume h a is the cost pe unit time to hold a type a job, and h b to hold a type b job. Type a b jobs ae fist pocessed at seve 1 2, then fok into two jobs, one that must be pocessed at seve 3 5 and the othe at seve 4, and finally join togethe to complete thei pocessing at seve 6 7. Thee is synchonization because the pocessing at seve 6 7 cannot begin until thee is at least one job in both buffes 7 and 8 9 and 10. Seve 4 pocesses both job types, but can only seve one job at a time. The contol decision is to decide which job type seve 4 should pioitize. This could be done by always pioitizing the moe expensive job type, in accodance with the well-known cµ-ule. Then, if h a µ A h b µ B whee µ A µ B is the ate at which 1 Mashall School of Business, Univesity of Southen Califonia, Ehun.Ozkan.2018@mashall.usc.edu. 2 Mashall School of Business, Univesity of Southen Califonia, amywad@mashall.usc.edu. 1

2 Type jobs Type jobs Figue 1: An example of a fok-join pocessing netwok. seve 4 pocesses type a b jobs, seve 4 always pefes to wok on a type a job ove a type b job. Howeve, when thee ae multiple jobs waiting at buffes 8 and 10, and no jobs waiting at buffe 9, it may be pefeable to have seve 4 wok on a type b job instead of a type a job and especially if also no jobs ae waiting at buffe 7. This is because seve 4 can pevent the join seve 7 fom idling without being the cause of seve 6 s foced idling. Seve 3 is the cause. The fok-join netwok contol poblem is too difficult to solve exactly, and so we seach fo an asymptotic solution. We do this unde the assumption that seve 4 is in heavy taffic. Othewise, the seve 4 scheduling contol has negligible impact on the delay of type a and type b jobs. We futhe assume the seves 6 and 7 ae in light taffic, which focuses attention on when the equied simultaneous pocessing of jobs at those seves foces idling. The seves 1, 2, 3, and 5 can all be in eithe light o heavy taffic. In the afoementioned heavy taffic egime, we fomulate and solve an appoximating diffusion contol poblem DCP. The DCP solution matches the numbe of jobs in buffe 4 to that in buffe 3, except when the total numbe of jobs waiting fo pocessing by seve 4 is too small fo that to be possible. The implication is that when seve 3 is in light taffic, so that buffe 3 is empty, buffe 4 is empty and all jobs waiting to be pocessed by seve 4 ae type b jobs. Othewise, when seve 3 is in heavy taffic, the contol at seve 4 must caefully pace its pocessing of type a jobs to pevent getting ahead of seve 3. Ou poposed policy is motivated by the obsevation fom the DCP solution that thee is no eason to have fewe type a jobs in buffe 4 than in buffe 3. If seve 3 can pocess jobs at least as quickly as seve 4 can pocess type a jobs, then the contol unde which seve 4 gives static pioity to type a jobs pefoms well. Othewise, we intoduce a slow depatue pacing SDP contol in which seve 4 slows its pocessing of type a jobs to match the depatue pocess of type a jobs fom the buffe 4 to the one fom the buffe 3. SDP is a obust idea that looks to be applicable to moe geneal fok-join netwoks. To see this, 2

3 we fomulate and solve appoximating DCPs fo a fok-join netwok with an abitay numbe of foks cf. Figue 5 in Section 9.1 and a fok-join netwok with thee job types cf. Figue 6 in Section 9.2. In each case, the DCP solutions suggest that, depending on the pocessing capacities of the seves, the seves in the netwok that pocess moe than one job type should eithe give static pioity to the moe expensive job types o slow the depatue pocess of these moe expensive jobs in ode to sometimes pioitize the less expensive jobs. This pevents unnecessay foced idling of the downsteam join seves that pocess the less expensive job types, without sacificing the speed at which the moe expensive job types depat the netwok. We pove that ou poposed policy fo the fok-join netwok in Figue 1 is asymptotically optimal in heavy taffic. To do this, we fist show that the DCP solution povides a stochastic lowe bound on the holding cost unde any policy at evey time instant. This is a stong objective, in line with that in Ata and Kuma [2005]. Then, we pove a weak convegence esult that implies the afoementioned lowe bound is achieved unde ou poposed policy. The weak convegence esult when the netwok opeates unde the SDP contol is a majo technical challenge fo the pape. This is because the SDP contol is a dynamic contol that depends on the netwok state. In ode to obtain the weak convegence, we must caefully constuct andom intevals on which we know the job type seve 4 is pioitizing. Although this idea is simila in spiit to the andom inteval constuction in Bell and Williams [2001], Ghamami and Wad [2013], the poof to show convegence on the intevals is much diffeent, due to the desied matching of the type a job depatue pocesses fom the seves 3 and 4. When we piece those intevals togethe, we see the DCP solution aise. That equies a esult on the ate of convegence of a light taffic GI/GI/1 queue that is diffeent fom that in Dai and Weiss [2002]. The emainde of this pape is oganized as follows. We conclude this section with a liteatue eview and a summay of ou mathematical notation. Section 2 specifies ou model and Section 3 povides ou asymptotic famewok. We constuct and solve an appoximating DCP in Section 4. We intoduce the SDP contol in Section 5, and specify when the poposed policy is SDP and when it is static pioity. Section 6 poves that the DCP solution povides a lowe bound on the pefomance of any contol, and Section 7 poves weak convegence unde the poposed policy. Section 8 povides extensive simulation esults. In Section 9 we constuct and solve appoximating DCPs fo a boade class of fok-join netwoks. Section 10 makes concluding emaks and poposes a futue eseach diection. We sepaate out ou ate of convegence esult fo a light taffic GI/GI/1 queue in the appendix, as that is a esult of inteest in its own ight. We also povide in the appendix poofs of esults that use moe standad methodology as well as moe detailed simulation esults. 1.1 Liteatue Review The inspiation fo this wok came fom the papes Nguyen [1993, 1994]. Nguyen [1993] establishes that a feedfowad FCFS fok-join netwok with one job type and single-seve stations in heavy taffic can be appoximated by a eflected Bownian motion RBM, and Nguyen [1994] extends this esult to include multiple job types. The dimension of the eflected Bownian motion equals the numbe of stations, and its state space is a polyhedal egion. In contast to the eflected Bownian motion appoximation fo feedfowad queueing netwoks Haison [1996], Haison and Van Mieghem [1997], Haison [1998, 2006], the effect of the synchonization constaints in fok-join netwoks is to incease the numbe of faces defining the state space. Also in contast to feedfowad queueing netwoks, that numbe is inceased futhe when moving fom the single job type to multiple job type scenaio. Although delay estimates fo fok-join netwoks follow fom the esults of Nguyen [1993, 1994], they leave open the question of whethe and how much delays can be shotened by scheduling jobs in a non-fcfs ode. To solve the scheduling poblem, we follow the standad Bownian machiney poposed in 3

4 Haison [1996]. This is typically done by fist intoducing a heavy-taffic asymptotic egime in which esouces ae almost fully utilized and the buffe content pocesses can be appoximated by a function of a Bownian motion, and second fomulating an appoximating Bownian contol poblem. Often, the dimension of the appoximating Bownian contol poblem can be educed by showing its equivalence with a so-called wokload fomulation Haison and Van Mieghem [1997], Haison [2006]. The intiguing diffeence when the undelying netwok is a fok-join netwok is that the join seves must be in light taffic to aive at an equivalent wokload fomulation. The issue is that othewise the appoximating poblem is non-linea. This light taffic assumption is asymptotically equivalent to the assumption that pocessing times ae instantaneous. Ou simulation esults suggest that ou poposed contol that is asymptotically optimal when the join seves ae in light taffic also pefoms vey well when the join seves ae in heavy taffic. The papes Pedasani et al. [2014a,b] ae some of the few studies we find that conside the contol of fok-join pocessing netwoks. In both papes, thee ae multiple job classes, but in Pedasani et al. [2014a] the seves can coopeate on the pocessing of jobs and in Pedasani et al. [2014b] they cannot. Thei focus is on finding obust polices in the discete-time setting that do not depend on system paametes and ae ate stable. They do not detemine whethe o not thei poposed policies minimize delay, which is ou focus. The pape Guvich and Wad [2014] seeks to minimize delay, but in the context of a matching queue netwok that has only joins and no foks. An essential question to answe when thinking about contols fo multiclass fok-join netwoks, as can be seen fom the papes Lu and Pang [2015a,b], Ata et al. [2012], is whethe o not the jobs being joined ae exchangeable; that is, whethe o not a task foked fom one job can be late joined with a task foked fom a diffeent job. Exchangeability is geneally tue in the manufactuing setting, because thee is no diffeence between pats with the same specifications, and geneally false in healthcae settings, because all papewok and test esults associated with one patient cannot be joined with anothe patient. The papes Lu and Pang [2015a,b] develop heavy taffic appoximations fo a non-exchangeable fok-join netwok with one aival steam that foks into aival steams to multiple many-seve queues, and then must be joined togethe to poduce one depatue steam. The heavy-taffic appoximation fo the non-exchangeable netwok is diffeent than fo the exchangeable netwok, and the non-exchangeability assumption inceases the poblem difficulty. Thei focus, diffeent than ous, is on the effect of coelation in the sevice times, and thee is no contol. The pape Ata et al. [2012] looks at a fok-join netwok in which thee is no contol decision if jobs ae exchangeable, and shows that the pefomance of the exchangeable netwok lowe bounds the pefomance of the non-exchangeable netwok. Then, they popose a contol fo the non-exchangeable netwok that achieves pefomance vey close to the exchangeable netwok. In compaison to the afoementioned papes, the exchangeability assumption is ielevant in ou case. This is because we assume head-of-line pocessing fo each job type, so that the exact same type a b jobs foked fom seve 1 2 ae the ones joined at seve Notation and Teminology The set of nonnegative integes is denoted by N, the set of stictly positive integes ae denoted by N +, and N := N {+ }. The k dimensional k N + Euclidean space is denoted by R k, R + denotes [0, +, and 0 k is the zeo vecto in R k. Fo x, y R, x y := max{x, y}, x y := min{x, y}, and x + := x 0. Fo any x R, x x denotes the geatest smallest intege which is smalle geate than o equal to x. The escipt denotes the tanspose of a matix o vecto. Fo each k N +, D k denotes the the space of all ω : R + R k which ae ight continuous with left limits on R +. Let 0 D be such that 0t = 0 fo all t R +. Fo ω D and T R +, we let ω T := 0 t T ωt. We conside D k endowed with the usual Skookhod J 1 topology cf. Chapte 3 of Billingsley [1999]. Let BD k denote the Boel σ-algeba on D k associated with Skookhod J 1 4

5 topology. By Theoem of Whitt [2002], BD k coincides with the Kolmogoov σ-algeba geneated by the coodinate pojections. Fo stochastic pocesses W, R +, and W whose sample paths ae in D k fo some k N +, W = W means that the pobability measues induced by W on D k, BD k weakly convege to the one induced by W on D k, BD k as. Fo x, y D, x y, x y, and x + ae pocesses in D such that x yt := xt yt, x yt := xt yt, and x + t := xt + fo all t R +. Fo x D, we define the mappings Ψ, Φ : D D such that fo all t 0, Ψxt := xs +, Φxt := xt + Ψxt, 0 s t whee Φ is the one-sided, one-dimensional eflection map, which is intoduced by Skookhod [1961]. Let Z = {1, 2,..., n} and X i be a pocess in D fo each i Z. Then X i, i Z denotes the pocess X 1, X 2,..., X n in D n. We denote e as the deteministic identity pocess in D such that et = t fo all t 0. We abbeviate the phase unifomly in compact intevals by u.o.c., almost suely by a.s.. We let a.s. denote almost sue convegence and = d denote equal in distibution. I denotes the indicato function and BM q θ, Σ denotes a Bownian motion with dift vecto θ and covaiance matix Σ, which stats at point q. The big-o notation is denoted by O, i.e., if x : R + R and y : R + R ae two functions, then xt = Oyt as t if and only if thee exist constants C and t 0 such that xt C yt fo all t t 0. Lastly, o p is the little-o in pobability notation, i.e., if {X n, n N} and {Y n, n N} ae sequences of andom vaiables, then X n = o p Y n if and only if X n / Y n conveges in pobability to 0. 2 Model Desciption We conside the contol of the fok-join pocessing netwok depicted in Figue 1. In this netwok, thee ae 2 job types, 7 seves, 10 buffes, and 8 activities. The set of job types is denoted by J, whee J = {a, b} and a and b denote the type a and type b jobs, espectively. The set of seves is denoted by S, whee S = {1, 2,..., 7}. The set of buffes is denoted by K, whee K = {1, 2,..., 10}, and the set of activities is denoted by A whee A = {1, 2, 3, A, B, 5, 6, 7}. Except fo seve 4, each seve is associated with a single activity. Seve 4 is associated with two activities, denoted by A and B, which ae pocessing type a jobs fom buffe 4 and type b jobs fom buffe 5, espectively. Both of the seves 6 and 7 deplete jobs fom 2 diffeent buffes. Note that these two seves ae join seves and pocess jobs wheneve both of the coesponding buffes ae nonempty. Hence, both of the seves 6 and 7 ae associated with a single activity, namely activities 6 and 7, espectively. Let s : A S be a function such that sj denotes the seve in which activity j, j A takes place. Let f : K\{1, 2} A be a function such that fk denotes the activity which feeds buffe k, k K\{1, 2}. Lastly, let d : K A be a function such that dk denotes the activity which depletes buffe k, k K. Fo example, sa = sb = 4, f4 = 1 and d4 = A. 2.1 Stochastic Pimitives We assume that all the andom vaiables and stochastic pocesses ae defined in the same complete pobability space Ω, F, P, E denotes the expectation unde P, and PA, B := PA B. We associate the extenal aival time of each job and the pocess time of each job in the coesponding activities with the sequence of andom vaiables { v j i, j J A} i=1 and the stictly positive constants {λ j, j J } and {µ j, j A}. We assume that fo each j J A, { v j i} i=1 is a stictly positive, independent and identically distibuted i.i.d. sequence of andom vaiables mutually independent of { v k i} i=1 fo all k J A\{j}, E[ v j1] = 1, and the vaiance of v j 1, denoted by Va v j 1, is σj 2. Fo j J, let v ji := v j i/λ j be the inteaival time between the i 1st and ith type j job. Then, λ j and σ j ae the aival ate and the coefficient of vaiation of the inteaival times of the type j jobs, whee j J. Fo j {1, 3, A, 6} j {2, B, 5, 7}, let 5

6 v j i := v j i/µ j be the sevice time of the ith incoming type a b job associated with the activity j. Then, µ j and σ j ae the sevice ate and the coefficient of vaiation of the sevice times elated to activity j, j A. Let fo each j J A, V j 0 := 0 and V j n := n v j i n N +, S j t := {n N : V j n t}. 1 i=1 Then, S j is a enewal pocess fo each j J A. If j J, S j t counts the numbe of extenal type j aivals until time t; if j A, S j t counts the numbe of sevice completions associated with the activity j until time t given that the coesponding seve woks continuously on this activity duing [0, t]. 2.2 Scheduling Contol and Netwok Dynamics Let T j t, j A, denote the cumulative amount of time seve sj devotes fo activity j duing [0, t]. Then, a scheduling contol is defined by the two dimensional sevice time allocation pocess T A, T B. Although a scheduling contol indiectly affects T 6, T 7, since we do not have any diect contol on seves 6 and 7, we exclude T 6, T 7 fom the definition of the scheduling contol. Let, I sj t := t T j t, j A\{A, B}, 2a I 4 t := t T A t T B t, denote the cumulative idle time of the seves duing the inteval [0, t]. Fo any j A, S j T j t denotes the total numbe of sevice completions elated to activity j in seve sj up to time t. Fo any k K, let Q k t be the numbe of jobs waiting in buffe k at time t, t 0, including the jobs that ae being seved. Then, fo all t 0, Q 1 t := S a t S 1 T 1 t 0, Q 2 t := S b t S 2 T 2 t 0, 3a 2b Q k t := S fk T fk t S dk T dk t 0, k K\{1, 2}. 3b Fo simplicity, we assume that initially all of the buffes ae empty, i.e., Q k 0 = 0 fo each k K. Late, we elax this assumption in Remak 11. We have V j S j T j t T j t < V j S j T j t + 1, fo all j A and t 0, 4 which implies that we conside only the head-of-the-line HL policies, whee jobs ae pocessed in the FCFS ode within each buffe. In this netwok, a task associated with a specific job cannot join a task oiginating in anothe job unde the HL policies; that is, ecalling ou liteatue eview cf. Section 1.1, the notion of exchangeability is not pesent. It is staightfowad to see that the wok conseving policies ae moe efficient than the othes in this netwok. Hence, we ensue that all of the seves wok in a wok-conseving fashion by the following constaints: Fo all t 0, I j is nondeceasing and I j 0 = 0 fo all j S, I sdk t inceases if and only if Q k t = 0, fo all k {1, 2, 3, 6}, I 4 t inceases if and only if Q 4 t Q 5 t = 0, I 6 t inceases if and only if Q 7 t Q 8 t = 0, I 7 t inceases if and only if Q 9 t Q 10 t = 0. 5a 5b 5c 5d 5e 6

7 A scheduling policy T := T A, T B is admissible if it satisfies the following conditions: Fo any T i, i A, I j, j S, and Q k, k K satisfying 2, 3, 4, and 5, T j t F, t 0 and j {A, B}, 6a T j is continuous and nondeceasing with T j 0 = 0, j {A, B}, I 4 is continuous and nondeceasing with I 4 0 = 0. Conditions 6a 6c imply that we conside a wide ange of scheduling policies including the ones which can anticipate the futue. 2.3 The Objective A natual objective is to minimize the discounted expected total holding cost. Let h a and h b denote the holding cost ate pe job pe unit time fo a type a and b job, espectively; and δ > 0 be the discount paamete. Moeove, let 6b 6c Z a t := Q 3 t + Q 4 t + Q 7 t + Q 8 t, Z b t := Q 5 t + Q 6 t + Q 9 t + Q 10 t. Then Z a t and Z b t denote the total numbe of type a and b jobs in the system except jobs waiting in buffes 1 and 2. Since Q 1 t and Q 2 t ae independent of the scheduling policy, we exclude these pocesses fom the definitions of Z a t and Z b t. Then, the discounted expected total holding cost unde admissible policy T is [ ] J T = E e δt h a Z a t + h b Z b t dt, 7 0 and ou objective is to find an admissible policy which minimizes 7. Anothe natual objective is to minimize the expected total cost up to time t, t R +, which is [ t ] J T = E h a Z a s + h b Z b sds. 8 0 Yet anothe possible objective is to minimize the long-un aveage cost pe unit time, [ 1 t ] J T = lim t t E h a Z a s + h b Z b s ds. 9 We focus on a moe challenging objective which is minimizing 0 P h a Z a t + h b Z b t > x, fo all t R + and x > It is possible to see that any policy that minimizes 10 also minimizes the objectives 7, 8, and 9. In this specific netwok, fo all t 0 Q 3 t + Q 7 t = Q 4 t + Q 8 t, Q 5 t + Q 9 t = Q 6 t + Q 10 t. 11 By 11, a policy is optimal unde the objective 10 if and only if it is optimal unde the objective of minimizing P h a Q 3 t + Q 7 t + h b Q 6 t + Q 10 t > x, fo all t R + and x > We will focus on the objective 12 fom this point fowad. 7

8 3 Asymptotic Famewok It is vey difficult to analyze the system descibed in Section 2 exactly. Even if we can accomplish this vey challenging task, it is even less likely that the optimal contol policy will be simple enough to be expessed by a few paametes. Theefoe, we focus on finding an asymptotically optimal contol policy unde diffusion scaling and the assumption that seve 4 is in heavy taffic. We fist intoduce a sequence of fok-join systems and pesent the main assumptions done fo this study in Section 3.1. Then we fomally define the fluid and diffusion scaled pocesses and pesent convegence esults fo the diffusion scaled wokload facing seve 4 and the diffusion scaled queue length pocesses associated with the seves 1, 2, 3, and 5 in Section 3.2. The question left open is to detemine what should be the elationship between the contol and the numbe of type a and b jobs in the wokload facing seve A Sequence of Fok-Join Systems We conside a sequence of fok-join systems indexed by whee though a sequence of values in R +. Each queueing system has the same stuctue defined in Section 2 except that the aival and sevice ates depend on. To be moe pecise, in the th system, we associate the extenal aival time of each job and the pocess time of each job in the coesponding activities with the sequence of andom vaiables { v j i, j J A} i=1, which we have defined in Section 2.1, and the stictly positive constants {λ j, j J } and {µ j, j A} such that v j i := v ji/λ j, j J is the inteaival time between the i 1st and ith type j job and vj i := v ji/µ j, j {1, 3, A, 6} j {2, B, 5, 7} is the sevice time of the ith incoming type a b job associated with the activity j in the th system. Theefoe, λ j, j J and µ j, j A ae the aival ates and sevice ates in the th system wheeas the coefficient of vaiations ae the same with the oiginal system defined in Section 2. Fom this point fowad, we will use the escipt to show the dependence of the stochastic pocesses to the th queueing system. Next, we pesent ou assumptions elated to the system paametes. The fist one is elated to cost paametes. Assumption 1 Without loss of geneality, we assume that h a µ A h b µ B. This assumption implies that it is moe expensive to keep the type a jobs than the type b jobs in seve 4. Second, we make the following assumptions elated to the stochastic pimitives of the netwok. Assumption 2 Thee exists a non-empty open neighbohood, O, aound 0 such that fo all α O, E[e αv j1 ] <, fo all j J A. Assumption 2 is the exponential moment assumption fo the inteaival and sevice time pocesses. This assumption is common in the queueing liteatue, cf. Haison [1998], Bell and Williams [2001], Maglaas [2003], Meyn [2003]. Remak 1 We will not need exponential moment assumption fo all inteaival and sevice time pocesses. We state the weakest moment conditions fo each inteaival and sevice time pocess equied in the poofs in Remak 10. Ou final assumption concens the convegence of the aival and sevice ates, and sets up heavy taffic asymptotic egime. Assumption 3 8

9 1. λ j λ j > 0 fo all j J as. 2. µ j µ j fo all j A as. 3. λ a /µ A + λ b /µ B = λ a/µ A + λ b /µ B 1 θ 4/µ A as, whee θ 4 is a constant in R. 5. As, 6. λ a < µ 6 and λ b < µ 7. λ a µ 1 θ 1 R { }, λ b µ 2 θ 2 R { }, 13a λ a µ 3 θ 3 R { }, λ b µ 5 θ 5 R { }. 13b Pats 3 and 4 of Assumption 3 ae the heavy taffic assumptions fo seve 4. Note that, if seve 4 was in light taffic, the contol of the netwok would not be inteesting. By Pat 3, we have µ A > λ a. If µ A = λ a, the fok-join netwok is stable in the diffusion limit only if λ b = 0, which is an uninteesting case to conside. Similaly, µ B > λ b. Pat 5 states that each of the seves 1, 2, 3, and 5 can be eithe in light o heavy taffic. Specifically fo each i {1, 2, 3, 5}, if θ i =, then seve i is in light taffic. On the othe hand, if θ i R, then seve i is in heavy taffic. Pat 6 states that the two join seves ae in light taffic. Note that, Ata et al. [2012], Guvich and Wad [2014], Lu and Pang [2015a,b] assume that the sevice pocesses in the join seves ae instantaneous. Hence, Pat 6 of Assumption 3 extends the assumptions on the join seves made in the liteatue. Fo simplicity, we assume that Q k 0 = 0, fo all k K and R +. Late, we elax this assumption in Remak 11. Assumptions 1, 2, and 3 ae assumed thoughout the pape. 3.2 Fluid and Diffusion Scaled Pocesses In this section, we pesent the fluid and diffusion scaled pocesses. Fo all t 0, the fluid scaled pocesses ae defined as S jt := 2 S j 2 t, j J A, T jt := 2 T j 2 t, j A, 14a Q jt := 2 Q j 2 t, j K, I jt := 2 I j 2 t, j S, 14b and the diffusion scaled pocesses ae defined as Ŝ j t := 1 S j 2 t λ j 2 t, j J, Ŝ j t := 1 S j 2 t µ j 2 t, j A, 15a T j t := 1 T j 2 t, j A, Î j t := 1 I j 2 t, j S, 15b Q jt := 1 Q j 2 t, j K. 15c By the functional cental limit theoem FCLT, cf. Theoem of Whitt [2002], we have the following weak convegence esult. Ŝ j, j J A = Sj, j J A, 16 whee S j is a one-dimensional Bownian Motion fo each j J A such that S d j = BM0 0, λ j σj 2 fo j J, S d j = BM0 0, µ j σj 2 fo j A and each S j is mutually independent of S i, i J A\{j}. 9

10 Fo t 0, let us define the wokload pocess up to a constant scale facto W4 t := Q 4t + µ A µ Q 5t. 17 B Then, W4 t/µ A is the expected time to deplete buffes 4 and 5 at time t, if no moe jobs aive afte time t in the th system. Let W 4t := 2 W4 2 t and Ŵ 4 t := 1 W4 2 t denote the fluid and diffusion scaled wokload pocesses, espectively. Next, we pesent the convegence of the fluid scaled pocesses unde any wok-conseving policy. Poposition 1 Unde any wok-conseving policy cf. 5, as Q k, k K, W 4, T j, j A a.s. Q k, k K, W 4, T j, j A u.o.c., whee Q k = 0 fo all k K, W 4 = 0, and T j t = λ a /µ j t fo j {1, 3, A, 6} and T j t = λ b /µ j t fo j {2, B, 5, 7} fo all t 0. The poof of Poposition 1 is pesented in Appendix B.1. We will use Poposition 1 to pove convegence esults fo the diffusion scaled pocesses. Consideing Assumption 3 Pat 5, let H := {i {1, 2, 3, 5} : θ i R}, i.e., H is the set of seves which ae in heavy taffic among the seves 1, 2, 3, and 5. Let H be the cadinality of the set H, and fo each i {1, 2, 3, 5}, let { 1, if seve i is in heavy taffic, i.e., i H, χ i := 0, if seve i is in light taffic, i.e., i / H. Fo all t 0, let S 1t := χ 1 Ŝ 1 T 1t λ a µ 1t + 1 χ 1 Ŝ a t Q 1t, 18a S 2t := χ 2 Ŝ 2 T 2t λ b µ 2t + 1 χ 2 Ŝ b t Q 2t. 18b Then, by using 18, we define the so-called netput pocess fo each buffe as X k t := Ŝ l t dk Ŝ T dkt + λ l µ dk t, fo k, l {1, a, 2, b}, 19a X k t := S l t dk Ŝ T dkt + λ i µ dk t, fo k, l, i {3, 1, a, 6, 2, b}, 19b X k t := S l t dk Ŝ T dkt + µ λ i dk λ i t, fo k, l, i {4, 1, a, 5, 2, b}, µ dk µ dk X k t := Ŝ l t Q i t Q jt Ŝ dk T dkt + λ l µ dk t, fo k, l, i, j {7, a, 1, 3, 8, a, 1, 4, 9, b, 2, 5, 10, b, 2, 6}. 19c 19d 10

11 Let θ := Q := Q 1 Q 2 Q 3 Q 6 Ŵ 4 θ 1 θ 2 θ 3 χ 1 θ 1 θ 5 χ 2 θ 2 θ 4 χ 1 θ 1 χ 2 µ A µb θ 2 X 1 X 2, X := X 3, Î := X 6 X 4 + µ A µ X B 5 Î 1 Î 2 Î 3 Î 5 Î 4, 20 µ µ , R := χ 1 µ 1 0 µ 3 0 0, 21 0 χ 2 µ 2 0 µ 5 0 χ 1 µ µ 1 χ A 2 µ µ B µ A and let R be a 5 5 matix which is the component-wise limit of R. Then, we have Q = X + R Î, 22 whee Let us define Σ := Q k = X k + µ dkî sdk, k {7, 8, 9, 10} λ a σa 2 + σ1 2 0 λ aσ1 2 0 λ a σ λ b σb 2 + σ2 2 0 λ bσ2 2 µ A µb λ b σ λ a σ1 2 0 Cov3, 3 0 Cov3, λ b σ2 2 0 Cov6, 6 Cov4, 6 4 λ a σ1 2 µ A µb λ b σ2 2 Cov3, 4 Cov4, 6 Cov4, 4 Cov3, 3 := λ a χ1 σ χ 1 σa 2 + σ3 2, Cov3, 4 := λa χ1 σ χ 1 σa 2, Cov4, 6 := µ A λ b χ2 σ χ 2 σ 2 b, Cov6, 6 := λb χ2 σ χ 2 σb 2 µ + σ2 5, B Cov4, 4 := λ a χ1 σ χ 1 σa 2 + σa 2 2 µa + λ b χ2 σ χ 2 σb 2 + B σ2. 25 Then, we have the following weak convegence esult. Poposition 2 Unde any wok-conseving policy cf. 5, Q 1, Q 2, Q 3, Q 6, Ŵ 4 = Q1, Q 2, Q 3, Q 6, W 4, 26 whee Q i = 0 fo each i / H and Qi, i H, W4 is a semimatingale eflected Bownian motion SRBM associated with the data P H, θ H, Σ H, R H, 0 H. P H is the nonnegative othant in R H ; θ H is an H -dimensional vecto deived fom the vecto θ cf. 21 by deleting the ows coesponding to each i, i / H; Σ H and R H ae H H -dimensional matices deived fom Σ cf. 24 and R cf. 21 by deleting the ows and columns coesponding to each i, i / H, espectively; and 0 H is the oigin in P H. The state space of the SRBM is P H ; θ H and Σ H ae the dift vecto and the covaiance matix of the undelying Bownian motion of the SRBM, espectively; R H is the eflection matix; and 0 H is the stating point of the SRBM. The fomal definition of a SRBM can be found in Definition 3.1 of Williams [1998]. The poof of Poposition 2 is pesented in Appendix B µ B 24

12 4 The Appoximating Diffusion Contol Poblem In this section, we constuct an appoximating diffusion contol poblem DCP and it will help us to guess policies which ae good candidates fo asymptotic optimality. Paallel with the objective 12, conside the objective of minimizing P h a Q 3 t + Q 7t + h b Q 6 t + Q 10t > x, t 0, x > 0, 27 fo some. Motivated by Assumption 3 Pat 6, let us petend that the sevice pocesses at seves 6 and 7 happen instantaneously. Since the diffusion scaled queue length pocess weakly conveges to 0 in a light taffic queue, we believe that consideing instantaneous sevice ates in the downsteam seves will not change the behavio of the system in the limit. In this case, jobs can accumulate in buffe 7 8 only at the times buffe 8 7 is empty. Similaly, jobs can accumulate in buffe 10 9 only at the times buffe 9 10 is empty. By this fact and 11, Q 7 = Q 4 Q 3 +, Q 8 = Q 3 Q 4 +, Q 9 = Q 6 Q 5 +, Q 10 = Q 5 Q By 28, the objective 27 is equivalent to minimizing P h a Q 3t + Q 4 Q h b Q 6t + Q 5 Q + 6 > x, t 0, x > Fom the objective 29, we need appoximations fo Q 3, Q 4, Q 5, and Q 6 and we will achieve this goal by letting. By 15c, 17, and Poposition 2, we know that Q 3, Q 6, weakly conveges to Q3, Q 6 and Q 4 + µ A /µ B Q 5 weakly conveges to W 4. At this point, let us assume that Q 4, Q 5 = Q4, Q 5. Then we constuct the following DCP: Fo any x > 0 and t 0, min P h a Q 3 t + Q4 t Q + 3 t + h b Q 6 t + Q5 t Q + 6 t > x, s.t. Q4 t + µ A µ B Q5 t = W 4 t, 30 Q k t 0, fo all k {4, 5}. Intuitively, we want to minimize the total cost by splitting the total wokload in the seve 4 to buffes 4 and 5 in the DCP 30. Now, we will conside DCP 30 path-wise. Let ω ω Ω denote a sample path of the pocesses in DCP 30 and fo any F : Ω D, F ω t denote the value of the pocess F at time t in the sample path ω. Then, conside the following optimization poblem fo each ω Ω and t 0. min h a Q4 ω t Q 3 ω t + + hb Q5 ω t Q 6 ω t +, 31a s.t. Q4 ω t + µ A µ B Q5 ω t = W 4 ω t, 31b Q k ω t 0, fo all k {4, 5}. 31c Note that, we exclude the tem h a Q3 ω t +h b Q6 ω t fom the objective function 31a because Q 3 ω t and Q 6 ω t ae independent of the decision vaiables Q 4 ω t and Q 5 ω t. Although the objective function 31a is nonlinea, the optimization poblem 31 is easy to solve because it has linea constaints. Moeove, Lemma 1, whose poof is pesented in Appendix B.3, povides a closed-fom solution fo

13 Lemma 1 Conside the optimization poblem min h a q 4 q h b q 5 q 6 +, s.t. q 4 + µ A µ B q 5 = w 4, q 4 0, q 5 0, whee q 4 and q 5 ae the decision vaiables, all the paametes ae nonnegative, and h a µ A h b µ B. Then q 4 = q 3 w 4 and q 5 = µ B /µ A w 4 q 3 + is an optimal solution of this poblem. Theefoe, by Lemma 1, we see that an optimal solution of the optimization poblem 31 is Q 4 ω t = W 4 ω t Q 3 ω t and Q 5 ω t = µ B /µ A W4 ω t Q 3 ω t + fo all ω and t 0. This esult and 28 imply the following poposition. Poposition 3 Q4, Q 5, Q 7, Q 8, Q 9, Q 10 = Q 3 W 4, µ B µ A W4 Q 3 +, 0, Q3 W 4 +, Q 6 µ B W4 µ Q + + 3, A µb W4 µ Q Q6 A 32 is an optimal solution of the DCP 30. Note that the ight hand side of 32 is independent of the scheduling policies by Poposition 2. Theefoe, a contol policy in which the coesponding pocesses weakly convege to the ight hand side of 32 is a good candidate fo an asymptotically optimal policy. The DCP solution in Poposition 3 matches the content level of buffe 4 to that of buffe 3, except when the buffe 3 content level exceeds the total wok facing seve 4 that is, the combined contents of buffes 4 and 5. This ensues that seve 4 neve causes seve 6 to idle because of the join opeation, as is evidenced by the fact that buffe 7 is always empty wheeas buffe 8 sometimes has a positive content level. At the same time, seve 4 pevents any unnecessay idling of seve 7 by devoting its emaining effot to pocessing the contents of buffe 5. Sometimes, that effot is sufficient to keep up with seve 5 and pevent the contents of buffe 10 fom building and sometimes it is not. That is why sometimes buffe 9 has a positive content in the DCP solution and sometimes buffe 10 does but neve both simultaneously. In the next section, we fomally intoduce the poposed policy. 5 Poposed Policy Ou objective is to popose a policy unde which the diffusion scaled queue-length pocesses tack the DCP solution given in Poposition 3. This is because the DCP solution povides a lowe bound on the asymptotic pefomance of any admissible policy, as we will pove in Section 6 see Theoem 2. The key obsevation fom the DCP solution in Poposition 3 is that thee is no eason fo the depatue pocess of the moe expensive type a jobs fom seve 4 to exceed that of seve 3. Instead of eve letting seve 4 get ahead, it is pefeable to have seve 4 wok on type b jobs, so as to pevent as much foced idling at seve 7 due to the join opeation as possible. The only time thee 13

14 should have been moe cumulative type a job depatues fom seve 4 than fom seve 3 is when the total numbe of jobs facing seve 4 is less than that facing seve 3. In that case, seve 4 can outpace seve 3 without focing additional idling at seve 7. The intuition in the peceding paagaph motivates the following depatue pacing policy, in which seve 4 gives pioity to type a jobs when the numbe of type a jobs in buffe 4 exceeds that in buffe 3 and gives pioity to the type b jobs in buffe 5 othewise. Definition 1 Slow Depatue Pacing SDP Policy. The allocation pocess T A, T B satisfies I Q 3 t < Q 4 t d t T A t = 0, I Q 3 t Q 4 t, Q 5 t > 0 d t T B t = 0, I Q 3 t Q 4 t, Q 4 t + Q 5 t > 0 di 4 t = 0, 33a 33b 33c togethe with 2, 3, 4, and 5. It is possible to see that T A, T B that satisfies 33 also satisfies 6, and so is admissible. If Q 3 t < Q 4 t, 33a ensues that seve 4 gives static pioity to buffe 4. If Q 3 t Q 4 t and Q 5 t > 0, 33b ensues that seve 4 gives static pioity to buffe 5. 33c ensues a wok-conseving contol policy in seve 4 when Q 3 t Q 4 t. When µ 3 < µ A, so that seve 3 is the slowe seve, we use the slow depatue pacing policy to detemine when seve 4 can allocate effot to pocessing type b jobs without inceasing type a job delay. Othewise, when µ 3 µ A, thee is almost neve exta pocessing powe to allocate to type b jobs, and so a static pioity policy will pefom similaly to the slow depatue pacing policy see Remak 12 egading ou numeical esults. Definition 2 Static Pioity Policy. The allocation pocess T A, T B satisfies 0 0 I Q 4 t > 0 d t T A t = 0, I Q 4 t + Q 5 t > 0 di 4 t = 0, 34a 34b togethe with 2, 3, 4, and 5. It is possible to see that T A, T B that satisfies 34 also satisfies 6, and so is admissible. 34a ensues that seve 4 gives static pioity to buffe 4 and 34b ensues a wok-conseving contol policy in seve 4. The poposed policy is the SDP policy in Definition 1 when µ 3 < µ A and is the static pioity policy in Definition 2 when µ 3 µ A. We have the following weak convegence esult associated with the poposed policy. Theoem 1 Unde the poposed policy, Q k, k K, Ŵ 4 = Qk, k K, W 4, whee Q1, Q 2, Q 3, Q 6, W 4 is defined in Poposition 2 and Q4, Q 5, Q 7, Q 8, Q 9, Q 10 Poposition 3. is defined in 14

15 The poof of Theoem 1 is pesented in Section 7. Theoem 1, the continuous mapping theoem see, fo example Theoem of Whitt [2002], and Theoem of Whitt [2002] establish the asymptotic behavio of the objective function 12 unde the poposed policy. Coollay 1 Unde the poposed policy, fo all t 0 and x > 0, we have lim P h a Q 3 t + Q 7t = P + h b Q 6 t + Q 10t h a Q3 t + h b > x µb Q 6 t + W4 t µ Q t Q6 t > x A. 35 Remak 2 The poposed policy is a peemptive policy. Howeve, it is often pefeed to use a nonpeemptive policy. To specify a non-peemptive policy, we must specify which type of job seve 4 chooses to pocess each time seve 4 becomes fee, and thee ae both type a and type b jobs waiting in buffes 4 and 5. The non-peemptive vesion of the SDP policy has seve 4 choose to seve a type a job when the numbe of jobs in buffe 4 exceeds that in buffe 3, and to seve a type b job othewise. The non-peemptive vesion of the static pioity policy has seve 4 always choose a type a job. We expect the pefomance of the non-peemptive vesion of ou poposed policy to be indistinguishable fom ou poposed policy in ou asymptotic egime, and we veify that the fome policy pefoms vey well by ou numeic esults in Section 8. Remak 3 In the classical open pocessing netwoks, if a seve is in light taffic, then the coesponding diffusion scaled buffe length pocess conveges to 0 see Theoem 6.1 of Chen and Mandelbaum [1991]. Howeve, we see in Theoem 1 that although the seves 6 and 7 ae in light taffic, Q 9 and Q 10 ae non-zeo pocesses moeove Q 8 is a non-zeo pocess when seve 3 is in heavy taffic. Theefoe, even though a join seve has moe than enough pocessing capacity, significant amount of jobs can accumulate in the coesponding buffes due to the synchonization equiements between the jobs in diffeent buffes. This makes the contol of fok-join netwoks moe challenging than the one of classical open-pocessing netwoks. Remak 4 We will pove that the SDP and static pioity policies ae asymptotically optimal when λ a µ 3 < µ A and µ 3 > λ a, espectively cf. Sections 7.1 and 7.2 and Remak 6. Hence, both SDP and static pioity policies ae asymtotically optimal when λ a < µ 3 < µ A and this implies that thee ae many othe asymptotically optimal contol policies in addition to the one that we popose. On the othe hand, the simulation expeiments show that the static pioity policy pefoms pooly when λ a µ 3 but the SDP policy pefoms vey well even when µ 3 > µ A cf. Section 8.2. Theefoe, the pefomance of the poposed policy that we constuct is obust with espect to the paametes λ a, µ 3, and µ A, which has moe pactical appeal. In the following section, we fomally define asymptotic optimality and pove that the poposed policy is asymptotically optimal given Theoem 1 and Coollay 1. 6 Asymptotic Optimality In this section, we pove that the DCP solution cf. Poposition 3 is a lowe bound fo all admissible policies with espect to the objective function 12, i.e., we have the following esult. 15

16 Theoem 2 Let {T, 0} be an abitay sequence of admissible policies. Then fo all t 0 and x > 0, we have lim inf P h a Q 3 t + Q 7t P + h b Q 6 t + Q 10t h a Q3 t + h b > x µb Q 6 t + W4 t µ Q t Q6 t > x A Theoem 2 togethe with Coollay 1 state that the poposed policy is asymptotically optimal.. 36 Remak 5 In Section 2.2, we state that the objective 12 implies the objectives 7, 8, and 9. Although 35 and 36 imply asymptotic optimality of the sequences of contol policies with espect to the objectives 7 and 8, they do not necessaily imply the same esult elated to the objective 9. That is because we need to change the ode of the limits with espect to and t to get the desied esult. Howeve, fo that, we need additional esults such as unifom convegence of the elated pocesses. 6.1 Poof of Theoem 2 Let us conside the tem in the left hand side of 36. By 11 and the fact that Q k 0 fo each k K, fo all t 0, Q 7t Q 4t Q 3t +, Q 10t Q 5t Q 6t +. Theefoe, it is enough to pove fo all t 0 and x > 0, lim inf h P a Q 3t + Q 4 t Q + 3t + h b Q 6t + Q 5 t Q + 6t h a Q3 t + h b P > x µb Q 6 t + W4 t µ Q t Q6 t > x A. 37 By 16, Poposition 1, and Theoem of Whitt [2002] joint convegence when one limit is deteministic, we have Ŝ j, j J A, T i, i A = Sj, j J A, T i, i A. 38 Now, we use Skookhod s epesentation theoem cf. Theoem of Whitt [2002] to obtain the equivalent distibutional epesentations of the pocesses in 38 fo which we use the same symbols and call Skookhod epesented vesions such that all Skookhod epesented vesions of the pocesses ae defined in the same pobability space and the weak convegence in 38 is eplaced by almost sue convegence. Then we have Ŝ j, j J A, T a.s. i, i A Sj, j J A, T i, i A, u.o.c. as. 39 We will conside the Skookhod epesented vesions of these pocesses fom this point fowad and pove 37 with espect to these pocesses. By 3, 15, 17, 19, 22, 39, and Poposition 2, we have Q 3, Ŵ 4, Q 6 a.s. Q 3, W 4, Q 6, u.o.c., 40 16

17 whee Ŵ4 = Q 4 + µ A Q µ 5, a.s., 41 B and all of the pocesses in 40 and 41 have the same distibution with the oiginal ones. Then by Fatou s lemma, the tem in the left hand side of 37 is geate than o equal to P lim inf [h a Q 3t + Q 4 t Q + 3t + h b Q 6t + Q 5 t Q ] + 6t > x. 42 Fo each sufficiently lage such that h a µ A h bµ B, we will find a path-wise lowe bound fo the tem [h a Q 3t + Q 4 t Q + 3t + h b Q 6t + Q 5 t Q ] + 6t. 43 Fom this point fowad, we will conside the sample paths in which Q k t, k {3, 4, 5, 6} ae finite fo all and t. By 40 and 41, these sample paths occu with pobability one. Let ω be a sample path and ω t be defined as in Section 4. By 41, 43, and the fact that Q 3 t and Q 6 t ae independent of the contol policy, we constuct the following optimization poblem. Fo each ω in Ω except a null set and t 0 min h a Q 4 ω t Q 3ω t + + hb Q 5 ω t Q 6ω t +, 44a s.t. Q 4ω t + µ A Q µ 5ω t = Ŵ 4 ω t, 44b B Q 4ω t 0, Q 5 ω t 0, 44c whee Q 4 ω t and Q 5 ω t ae the decision vaiables. The optimization poblem 44 has the same stuctue with the one pesented in Lemma 1. Theefoe, we can use Lemma 1 to solve 44 and in the optimal solution Q 4ω t = Q 3ω t Ŵ 4 ω t, Q 5 ω t = µ B Ŵ µ 4 ω t Q + 3ω t. 45 A Theefoe, by 45, a path-wise lowe bound on 43 unde the admissible policy T is µ h a Q 3 t + h b Q 6t + B Ŵ µ 4 t Q + + 3t Q 6 t. 46 A When we take the lim inf of the tem in 46, by 40 and the continuous mapping theoem specifically we use the continuity of the mapping + and Theoem 4.1 of Whitt [1980], which shows the continuity of addition, 42 is geate than o equal to µb P h a Q3 t + h b Q 6 t + W4 t µ Q t Q6 t > x. 47 A Note that the lowe bound in 47 is independent of contol by Poposition 2 and 40. Theefoe, 47 poves 36 fo the Skookhod epesented vesions of the pocesses. Since these pocesses have the same joint distibution with the oiginal ones, 47 also poves Theoem 2. 17

18 7 Weak Convegence Poof In this section, we pove Theoem 1. We conside the cases λ a µ 3 < µ A and µ 3 µ A sepaately. Note that the poposed policy is the SDP policy in 33 in the fist case and the static pioity policy in 34 in the second case. The poof of the second case is staightfowad, but the poof of the fist case is complicated because the poposed policy is a continuous-eview and state-dependent policy in this case. 7.1 Case I: λ a µ 3 < µ A Slow Depatue Pacing Policy The following esult plays a cucial ole in the weak convegence of Q 4 and Q 5 unde the poposed policy. Poposition 4 Unde the poposed policy, fo all ɛ > 0 and T > 0, lim P Q 4 Q 3 Ŵ 4 T > ɛ = The poof of Poposition 4 is pesented in Sections 7.3 and 7.4. By 26, 48, and Theoem of Whitt [2002] convegence-togethe theoem, we have the following joint convegence esult. Q 1, Q 2, Q 3, Q 4, Q 6, Ŵ 4 = Q1, Q 2, Q 3, Q 3 W 4, Q 6, W By 17 and 49, we have Q 5 = µ B µ A W 4 Q At this point, we invoke the Skookhod epesentation theoem again fo all the pocesses in 38, 49, and 50, and we will use the same symbols again. Then, we can eplace the weak convegence in 38, 49, and 50 with almost sue convegence in u.o.c. fo the Skookhod epesented vesions of the pocesses. Next, we will conside seve 6. Let Ẑ 6 := Q 7 Q 8. Then, by 19d and 23, Ẑ 6t = Û 6 t + µ 6Î 6t, Û 6 t := Ŝ at Q 1t Q 3t Q 4t Ŝ 6T 6t + λ a µ 6t, 51 whee Û 6 is defined in 51. Since Î 6 t inceases only if Ẑ 6 t = 0, we have a Skookhod poblem with espect to Û 6. By the uniqueness of the solution of the Skookhod poblem cf. Poposition B.1 of Bell and Williams [2001] µ 6Î 6 = ΨÛ 6 and Ẑ 6 = ΦÛ 6 fo each. By 39 and the fact that Q 1, Q 3, Q a.s. 4 Q1, Q 3, Q 3 W 4 µ 6Î 6 + λ a µ 6e, Ẑ6 a.s. S a + Q 1 + Q 3 + S 6 T 6, 0, u.o.c., 52 by Lemma 6.4 ii of Chen and Yao [2001]. By 19d, 23, and 52, Q 7 = Ŝ a Q 1 Q 3 Ŝ 6T 6t + λ a µ 6e + µ 6Î a.s. 6 S a Q 1 Q 3 S 6 T 6 S a + Q 1 + Q 3 + S 6 T 6 = 0 Q 8 = Ŝ a Q 1 Q 4 Ŝ 6T 6t + λ a µ 6e + µ 6Î 6 a.s. u.o.c., S a Q 1 Q 3 W 4 S 6 T 6 S a + Q 1 + Q 3 + S 6 T 6 = Q3 W 4 + u.o.c. By the same way, we can deive the following esult fo seve 7. Q 9, Q a.s. 10 Q6 µ B W4 µ Q + + µb 3, W4 A µ Q Q6, u.o.c. A 18

19 Since the Skookhod epesented vesion of the pocesses have the same joint distibution with the oiginal ones, when the Skookhod epesented vesions of the pocesses convege almost suely u.o.c., then the oiginal pocesses weakly convege and we get the desied esult. 7.2 Case II: µ 3 µ A Static Pioity Policy When µ 3 µ A, seve 3 is in light taffic because Assumption 3 Pat 3 implies µ A > λ a. Then Q 3 = Q 3 = 0 by Poposition 2. Since seve 4 gives static pioity to buffe 4 ove buffe 5 when µ 3 µ A in the poposed policy, then buffe 4 acts like a light taffic queue and Q 4 = 0 = Q 3 W 4. This implies that all of the wokload in seve 4 accumulates in buffe 5 and Q 5 = µ B/µ A W 4 by 17 and Poposition 2. The convegence poof of all othe pocesses is vey simila to the one pesented in Section 7.1. Remak 6 It is staightfowad to see that the poof pesented above holds when seve 3 is in light taffic λ a < µ 3. Hence, Theoem 1 and the Coollay 1 holds unde the static pioity policy wheneve seve 3 is in light taffic. Theefoe, as stated in Remak 4, the static pioity policy is asymptotically optimal wheneve seve 3 is in light taffic. Remak 7 When µ 3 µ A, we do not need the Poposition 4 to get the desied weak convegence esult. Howeve, the poof of Poposition 4 is staightfowad in this case. By Theoem 1 and the continuous mapping theoem, Q 4 Q 3 Ŵ 4 = 0, and this clealy gives us the desied esult. Fom this point fowad, we will only conside the case λ a µ 3 < µ A and we pove Poposition 4 unde this case in the next two sections. The key is the constuction of andom intevals such that in any given inteval, we know if seve 4 is giving pioity to the buffe 4 o to the buffe 5. We define the intevals in Section 7.3 and then use them to pove Poposition 4 in Section Inteval Constuction The intevals ae defined as shown in Figue 2. In the th system, let τn : Ω R + {+ } be such that τ0 := 0 and τ 2n 1 := inf{t > τ 2n 2 : Q 3t = Q 4t 1}, n N +, 53a τ 2n := inf{t > τ 2n 1 : Q 3t = Q 4t}, n N +. 53b Fo completeness, if τn 0 = + fo some n 0 N +, we define τn := + fo all n n 0. We also define fo all n N, inf { t [τ τ, 2n := 2n, τ 2n+1 : Q 5 t = 0}, if Q 5 t = 0 fo some t [τ 2n, τ 2n+1 whee τ2n <, 54 +, othewise. Note that, it is possible that the numbe of jobs in buffe 5 neve becomes 0 duing the inteval [τ2n, τ 2n+1,. In this case, we define τ2n := in 54 fo completeness. Fo each n N +, we call [τ2n 1, τ 2n as an up inteval because Q 4 t > Q 3 t fo all t [τ2n 1, τ 2n by constuction. On the othe hand, fo each n N, we call [τ 2n, τ 2n+1 as a down inteval because Q 4 t Q 3 t fo all t [τ 2n, τ 2n+1 by constuction. Moeove, fo each n N, if τ, 2n <, we call [τ 2n, τ,, 2n as a down 1 inteval and [τ2n, τ 2n+1 as a down 2 inteval. Thus, we sepaate a down inteval into two disjoint intevals. The poposed policy gives static peemptive pioity to buffe 4 ove buffe 5 duing the up intevals and gives static peemptive pioity to buffe 5 duing the down intevals. 19

arxiv: v2 [math.pr] 14 Jun 2016

arxiv: v2 [math.pr] 14 Jun 2016 On the Contol of Fok-Join Netwoks Ehun Özkan1 and Amy R. Wad 2 axiv:1505.04470v2 [math.pr] 14 Jun 2016 Abstact Netwoks in which the pocessing of jobs occus both sequentially and in paallel ae pevalent

More information

4/18/2005. Statistical Learning Theory

4/18/2005. Statistical Learning Theory Statistical Leaning Theoy Statistical Leaning Theoy A model of supevised leaning consists of: a Envionment - Supplying a vecto x with a fixed but unknown pdf F x (x b Teache. It povides a desied esponse

More information

Dynamic Control of N-Systems with Many Servers: Asymptotic Optimality of a Static Priority Policy in Heavy Traffic

Dynamic Control of N-Systems with Many Servers: Asymptotic Optimality of a Static Priority Policy in Heavy Traffic Dynamic Contol of N-Systems with Many Seves: Asymptotic Optimality of a Static Pioity Policy in Heavy Taffic Tolga Tezcan Industial and Entepise Systems Engineeing Univesity of Illinois at Ubana-Champaign

More information

ON INDEPENDENT SETS IN PURELY ATOMIC PROBABILITY SPACES WITH GEOMETRIC DISTRIBUTION. 1. Introduction. 1 r r. r k for every set E A, E \ {0},

ON INDEPENDENT SETS IN PURELY ATOMIC PROBABILITY SPACES WITH GEOMETRIC DISTRIBUTION. 1. Introduction. 1 r r. r k for every set E A, E \ {0}, ON INDEPENDENT SETS IN PURELY ATOMIC PROBABILITY SPACES WITH GEOMETRIC DISTRIBUTION E. J. IONASCU and A. A. STANCU Abstact. We ae inteested in constucting concete independent events in puely atomic pobability

More information

Validity of heavy-traffic steady-state approximations in multiclass queueing networks: The case of queue-ratio disciplines

Validity of heavy-traffic steady-state approximations in multiclass queueing networks: The case of queue-ratio disciplines Validity of heavy-taffic steady-state appoximations in multiclass queueing netwos: The case of queue-atio disciplines Itai Guvich Kellogg School of Management, Nothwesten Univesity, Evanston, IL 628 email:

More information

A Bijective Approach to the Permutational Power of a Priority Queue

A Bijective Approach to the Permutational Power of a Priority Queue A Bijective Appoach to the Pemutational Powe of a Pioity Queue Ia M. Gessel Kuang-Yeh Wang Depatment of Mathematics Bandeis Univesity Waltham, MA 02254-9110 Abstact A pioity queue tansfoms an input pemutation

More information

PAPER 39 STOCHASTIC NETWORKS

PAPER 39 STOCHASTIC NETWORKS MATHEMATICAL TRIPOS Pat III Tuesday, 2 June, 2015 1:30 pm to 4:30 pm PAPER 39 STOCHASTIC NETWORKS Attempt no moe than FOUR questions. Thee ae FIVE questions in total. The questions cay equal weight. STATIONERY

More information

Fractional Zero Forcing via Three-color Forcing Games

Fractional Zero Forcing via Three-color Forcing Games Factional Zeo Focing via Thee-colo Focing Games Leslie Hogben Kevin F. Palmowski David E. Robeson Michael Young May 13, 2015 Abstact An -fold analogue of the positive semidefinite zeo focing pocess that

More information

Lecture 28: Convergence of Random Variables and Related Theorems

Lecture 28: Convergence of Random Variables and Related Theorems EE50: Pobability Foundations fo Electical Enginees July-Novembe 205 Lectue 28: Convegence of Random Vaiables and Related Theoems Lectue:. Kishna Jagannathan Scibe: Gopal, Sudhasan, Ajay, Swamy, Kolla An

More information

Control of Patient Flow in Emergency Departments, or Multiclass Queues with Deadlines and Feedback

Control of Patient Flow in Emergency Departments, or Multiclass Queues with Deadlines and Feedback Submitted to Opeations Reseach manuscipt Please, povide the mansucipt numbe! Contol of Patient Flow in Emegency Depatments, o Multiclass Queues with Deadlines and Feedback Junfei Huang National Univesity

More information

Absorption Rate into a Small Sphere for a Diffusing Particle Confined in a Large Sphere

Absorption Rate into a Small Sphere for a Diffusing Particle Confined in a Large Sphere Applied Mathematics, 06, 7, 709-70 Published Online Apil 06 in SciRes. http://www.scip.og/jounal/am http://dx.doi.og/0.46/am.06.77065 Absoption Rate into a Small Sphee fo a Diffusing Paticle Confined in

More information

arxiv:math/ v1 [math.pr] 24 Aug 2005

arxiv:math/ v1 [math.pr] 24 Aug 2005 The Annals of Applied Pobability 5, Vol 15, No 3, 1887 1935 DOI: 1114/1551655 c Institute of Mathematical Statistics, 5 axiv:math/5846v1 [mathpr] 4 Aug 5 A LARGE DEVIATIONS APPROACH TO ASYMPTOTICALLY OPTIMAL

More information

New problems in universal algebraic geometry illustrated by boolean equations

New problems in universal algebraic geometry illustrated by boolean equations New poblems in univesal algebaic geomety illustated by boolean equations axiv:1611.00152v2 [math.ra] 25 Nov 2016 Atem N. Shevlyakov Novembe 28, 2016 Abstact We discuss new poblems in univesal algebaic

More information

arxiv: v1 [math.co] 4 May 2017

arxiv: v1 [math.co] 4 May 2017 On The Numbe Of Unlabeled Bipatite Gaphs Abdullah Atmaca and A Yavuz Ouç axiv:7050800v [mathco] 4 May 207 Abstact This pape solves a poblem that was stated by M A Haison in 973 [] This poblem, that has

More information

Internet Appendix for A Bayesian Approach to Real Options: The Case of Distinguishing Between Temporary and Permanent Shocks

Internet Appendix for A Bayesian Approach to Real Options: The Case of Distinguishing Between Temporary and Permanent Shocks Intenet Appendix fo A Bayesian Appoach to Real Options: The Case of Distinguishing Between Tempoay and Pemanent Shocks Steven R. Genadie Gaduate School of Business, Stanfod Univesity Andey Malenko Gaduate

More information

Do Managers Do Good With Other People s Money? Online Appendix

Do Managers Do Good With Other People s Money? Online Appendix Do Manages Do Good With Othe People s Money? Online Appendix Ing-Haw Cheng Haison Hong Kelly Shue Abstact This is the Online Appendix fo Cheng, Hong and Shue 2013) containing details of the model. Datmouth

More information

DIFFUSION LIMITS FOR SHORTEST REMAINING PROCESSING TIME QUEUES

DIFFUSION LIMITS FOR SHORTEST REMAINING PROCESSING TIME QUEUES Stochastic Systems 2011, Vol. 1, No. 1, 1 16 DOI: 10.1214/10-SSY016 DIFFUSION LIMITS FOR SHORTEST REMAINING PROCESSING TIME QUEUES By H. Chistian Gomoll, Lukasz Kuk and Ambe L. Puha Univesity of Viginia,

More information

Web-based Supplementary Materials for. Controlling False Discoveries in Multidimensional Directional Decisions, with

Web-based Supplementary Materials for. Controlling False Discoveries in Multidimensional Directional Decisions, with Web-based Supplementay Mateials fo Contolling False Discoveies in Multidimensional Diectional Decisions, with Applications to Gene Expession Data on Odeed Categoies Wenge Guo Biostatistics Banch, National

More information

The Substring Search Problem

The Substring Search Problem The Substing Seach Poblem One algoithm which is used in a vaiety of applications is the family of substing seach algoithms. These algoithms allow a use to detemine if, given two chaacte stings, one is

More information

A Multivariate Normal Law for Turing s Formulae

A Multivariate Normal Law for Turing s Formulae A Multivaiate Nomal Law fo Tuing s Fomulae Zhiyi Zhang Depatment of Mathematics and Statistics Univesity of Noth Caolina at Chalotte Chalotte, NC 28223 Abstact This pape establishes a sufficient condition

More information

3.1 Random variables

3.1 Random variables 3 Chapte III Random Vaiables 3 Random vaiables A sample space S may be difficult to descibe if the elements of S ae not numbes discuss how we can use a ule by which an element s of S may be associated

More information

arxiv: v1 [math.pr] 6 May 2010

arxiv: v1 [math.pr] 6 May 2010 axiv:1005.1035v1 [math.pr] 6 May 2010 DIFFUSION LIMITS FOR SHORTEST REMAINING PROCESSING TIME QUEUES H. Chistian Gomoll, Lukasz Kuk, Ambe L. Puha Univesity of Viginia, Maia Cuie-Sk lodowska Univesity,

More information

Pearson s Chi-Square Test Modifications for Comparison of Unweighted and Weighted Histograms and Two Weighted Histograms

Pearson s Chi-Square Test Modifications for Comparison of Unweighted and Weighted Histograms and Two Weighted Histograms Peason s Chi-Squae Test Modifications fo Compaison of Unweighted and Weighted Histogams and Two Weighted Histogams Univesity of Akueyi, Bogi, v/noduslód, IS-6 Akueyi, Iceland E-mail: nikolai@unak.is Two

More information

On the Quasi-inverse of a Non-square Matrix: An Infinite Solution

On the Quasi-inverse of a Non-square Matrix: An Infinite Solution Applied Mathematical Sciences, Vol 11, 2017, no 27, 1337-1351 HIKARI Ltd, wwwm-hikaicom https://doiog/1012988/ams20177273 On the Quasi-invese of a Non-squae Matix: An Infinite Solution Ruben D Codeo J

More information

Bayesian Congestion Control over a Markovian Network Bandwidth Process

Bayesian Congestion Control over a Markovian Network Bandwidth Process Bayesian Congestion Contol ove a Makovian Netwok Bandwidth Pocess Paisa Mansouifad,, Bhaska Kishnamachai, Taa Javidi Ming Hsieh Depatment of Electical Engineeing, Univesity of Southen Califonia, Los Angeles,

More information

EM Boundary Value Problems

EM Boundary Value Problems EM Bounday Value Poblems 10/ 9 11/ By Ilekta chistidi & Lee, Seung-Hyun A. Geneal Desciption : Maxwell Equations & Loentz Foce We want to find the equations of motion of chaged paticles. The way to do

More information

Efficiency Loss in a Network Resource Allocation Game

Efficiency Loss in a Network Resource Allocation Game Efficiency Loss in a Netwok Resouce Allocation Game Ramesh Johai johai@mit.edu) John N. Tsitsiklis jnt@mit.edu) June 11, 2004 Abstact We exploe the popeties of a congestion game whee uses of a congested

More information

Temporal-Difference Learning

Temporal-Difference Learning .997 Decision-Making in Lage-Scale Systems Mach 17 MIT, Sping 004 Handout #17 Lectue Note 13 1 Tempoal-Diffeence Leaning We now conside the poblem of computing an appopiate paamete, so that, given an appoximation

More information

DIFFUSION LIMITS FOR SHORTEST REMAINING PROCESSING TIME QUEUES UNDER NONSTANDARD SPATIAL SCALING

DIFFUSION LIMITS FOR SHORTEST REMAINING PROCESSING TIME QUEUES UNDER NONSTANDARD SPATIAL SCALING Submitted to the Annals of Applied Pobability axiv: axiv:1407.3837 DIFFUSION LIMITS FOR SHORTEST REMAINING PROCESSING TIME QUEUES UNDER NONSTANDARD SPATIAL SCALING By Ambe L. Puha Califonia State Univesity

More information

Unobserved Correlation in Ascending Auctions: Example And Extensions

Unobserved Correlation in Ascending Auctions: Example And Extensions Unobseved Coelation in Ascending Auctions: Example And Extensions Daniel Quint Univesity of Wisconsin Novembe 2009 Intoduction In pivate-value ascending auctions, the winning bidde s willingness to pay

More information

Central Coverage Bayes Prediction Intervals for the Generalized Pareto Distribution

Central Coverage Bayes Prediction Intervals for the Generalized Pareto Distribution Statistics Reseach Lettes Vol. Iss., Novembe Cental Coveage Bayes Pediction Intevals fo the Genealized Paeto Distibution Gyan Pakash Depatment of Community Medicine S. N. Medical College, Aga, U. P., India

More information

10/04/18. P [P(x)] 1 negl(n).

10/04/18. P [P(x)] 1 negl(n). Mastemath, Sping 208 Into to Lattice lgs & Cypto Lectue 0 0/04/8 Lectues: D. Dadush, L. Ducas Scibe: K. de Boe Intoduction In this lectue, we will teat two main pats. Duing the fist pat we continue the

More information

6 PROBABILITY GENERATING FUNCTIONS

6 PROBABILITY GENERATING FUNCTIONS 6 PROBABILITY GENERATING FUNCTIONS Cetain deivations pesented in this couse have been somewhat heavy on algeba. Fo example, detemining the expectation of the Binomial distibution (page 5.1 tuned out to

More information

Optimal Multiplexing on a Single Link: Delay and Buer. Requirements. Abhay Parekh y

Optimal Multiplexing on a Single Link: Delay and Buer. Requirements. Abhay Parekh y Optimal Multiplexing on a Single Link: Delay and Bue Requiements Leonidas Geogiadis Aistotle Univesity of Thessaloniki Dept. of Electical and Compute Eng. P.O. Box 435 Thessaloniki, 54006 - GREECE Roch

More information

Chapter 5 Linear Equations: Basic Theory and Practice

Chapter 5 Linear Equations: Basic Theory and Practice Chapte 5 inea Equations: Basic Theoy and actice In this chapte and the next, we ae inteested in the linea algebaic equation AX = b, (5-1) whee A is an m n matix, X is an n 1 vecto to be solved fo, and

More information

Stanford University CS259Q: Quantum Computing Handout 8 Luca Trevisan October 18, 2012

Stanford University CS259Q: Quantum Computing Handout 8 Luca Trevisan October 18, 2012 Stanfod Univesity CS59Q: Quantum Computing Handout 8 Luca Tevisan Octobe 8, 0 Lectue 8 In which we use the quantum Fouie tansfom to solve the peiod-finding poblem. The Peiod Finding Poblem Let f : {0,...,

More information

On the global uniform asymptotic stability of time-varying dynamical systems

On the global uniform asymptotic stability of time-varying dynamical systems Stud. Univ. Babeş-Bolyai Math. 59014), No. 1, 57 67 On the global unifom asymptotic stability of time-vaying dynamical systems Zaineb HajSalem, Mohamed Ali Hammami and Mohamed Mabouk Abstact. The objective

More information

Surveillance Points in High Dimensional Spaces

Surveillance Points in High Dimensional Spaces Société de Calcul Mathématique SA Tools fo decision help since 995 Suveillance Points in High Dimensional Spaces by Benad Beauzamy Januay 06 Abstact Let us conside any compute softwae, elying upon a lage

More information

An Application of Fuzzy Linear System of Equations in Economic Sciences

An Application of Fuzzy Linear System of Equations in Economic Sciences Austalian Jounal of Basic and Applied Sciences, 5(7): 7-14, 2011 ISSN 1991-8178 An Application of Fuzzy Linea System of Equations in Economic Sciences 1 S.H. Nassei, 2 M. Abdi and 3 B. Khabii 1 Depatment

More information

Efficiency Loss in a Network Resource Allocation Game: The Case of Elastic Supply

Efficiency Loss in a Network Resource Allocation Game: The Case of Elastic Supply Efficiency Loss in a Netwok Resouce Allocation Game: The Case of Elastic Supply axiv:cs/0506054v1 [cs.gt] 14 Jun 2005 Ramesh Johai (johai@stanfod.edu) Shie Manno (shie@mit.edu) John N. Tsitsiklis (jnt@mit.edu)

More information

Hypothesis Test and Confidence Interval for the Negative Binomial Distribution via Coincidence: A Case for Rare Events

Hypothesis Test and Confidence Interval for the Negative Binomial Distribution via Coincidence: A Case for Rare Events Intenational Jounal of Contempoay Mathematical Sciences Vol. 12, 2017, no. 5, 243-253 HIKARI Ltd, www.m-hikai.com https://doi.og/10.12988/ijcms.2017.7728 Hypothesis Test and Confidence Inteval fo the Negative

More information

Duality between Statical and Kinematical Engineering Systems

Duality between Statical and Kinematical Engineering Systems Pape 00, Civil-Comp Ltd., Stiling, Scotland Poceedings of the Sixth Intenational Confeence on Computational Stuctues Technology, B.H.V. Topping and Z. Bittna (Editos), Civil-Comp Pess, Stiling, Scotland.

More information

ANA BERRIZBEITIA, LUIS A. MEDINA, ALEXANDER C. MOLL, VICTOR H. MOLL, AND LAINE NOBLE

ANA BERRIZBEITIA, LUIS A. MEDINA, ALEXANDER C. MOLL, VICTOR H. MOLL, AND LAINE NOBLE THE p-adic VALUATION OF STIRLING NUMBERS ANA BERRIZBEITIA, LUIS A. MEDINA, ALEXANDER C. MOLL, VICTOR H. MOLL, AND LAINE NOBLE Abstact. Let p > 2 be a pime. The p-adic valuation of Stiling numbes of the

More information

Multihop MIMO Relay Networks with ARQ

Multihop MIMO Relay Networks with ARQ Multihop MIMO Relay Netwoks with ARQ Yao Xie Deniz Gündüz Andea Goldsmith Depatment of Electical Engineeing Stanfod Univesity Stanfod CA Depatment of Electical Engineeing Pinceton Univesity Pinceton NJ

More information

Functions Defined on Fuzzy Real Numbers According to Zadeh s Extension

Functions Defined on Fuzzy Real Numbers According to Zadeh s Extension Intenational Mathematical Foum, 3, 2008, no. 16, 763-776 Functions Defined on Fuzzy Real Numbes Accoding to Zadeh s Extension Oma A. AbuAaqob, Nabil T. Shawagfeh and Oma A. AbuGhneim 1 Mathematics Depatment,

More information

On the integration of the equations of hydrodynamics

On the integration of the equations of hydrodynamics Uebe die Integation de hydodynamischen Gleichungen J f eine u angew Math 56 (859) -0 On the integation of the equations of hydodynamics (By A Clebsch at Calsuhe) Tanslated by D H Delphenich In a pevious

More information

On a quantity that is analogous to potential and a theorem that relates to it

On a quantity that is analogous to potential and a theorem that relates to it Su une quantité analogue au potential et su un théoème y elatif C R Acad Sci 7 (87) 34-39 On a quantity that is analogous to potential and a theoem that elates to it By R CLAUSIUS Tanslated by D H Delphenich

More information

Journal of Inequalities in Pure and Applied Mathematics

Journal of Inequalities in Pure and Applied Mathematics Jounal of Inequalities in Pue and Applied Mathematics COEFFICIENT INEQUALITY FOR A FUNCTION WHOSE DERIVATIVE HAS A POSITIVE REAL PART S. ABRAMOVICH, M. KLARIČIĆ BAKULA AND S. BANIĆ Depatment of Mathematics

More information

6 Matrix Concentration Bounds

6 Matrix Concentration Bounds 6 Matix Concentation Bounds Concentation bounds ae inequalities that bound pobabilities of deviations by a andom vaiable fom some value, often its mean. Infomally, they show the pobability that a andom

More information

Chapter 3: Theory of Modular Arithmetic 38

Chapter 3: Theory of Modular Arithmetic 38 Chapte 3: Theoy of Modula Aithmetic 38 Section D Chinese Remainde Theoem By the end of this section you will be able to pove the Chinese Remainde Theoem apply this theoem to solve simultaneous linea conguences

More information

Lecture 8 - Gauss s Law

Lecture 8 - Gauss s Law Lectue 8 - Gauss s Law A Puzzle... Example Calculate the potential enegy, pe ion, fo an infinite 1D ionic cystal with sepaation a; that is, a ow of equally spaced chages of magnitude e and altenating sign.

More information

PROBLEM SET #1 SOLUTIONS by Robert A. DiStasio Jr.

PROBLEM SET #1 SOLUTIONS by Robert A. DiStasio Jr. POBLM S # SOLUIONS by obet A. DiStasio J. Q. he Bon-Oppenheime appoximation is the standad way of appoximating the gound state of a molecula system. Wite down the conditions that detemine the tonic and

More information

Asymptotically Lacunary Statistical Equivalent Sequence Spaces Defined by Ideal Convergence and an Orlicz Function

Asymptotically Lacunary Statistical Equivalent Sequence Spaces Defined by Ideal Convergence and an Orlicz Function "Science Stays Tue Hee" Jounal of Mathematics and Statistical Science, 335-35 Science Signpost Publishing Asymptotically Lacunay Statistical Equivalent Sequence Spaces Defined by Ideal Convegence and an

More information

Chem 453/544 Fall /08/03. Exam #1 Solutions

Chem 453/544 Fall /08/03. Exam #1 Solutions Chem 453/544 Fall 3 /8/3 Exam # Solutions. ( points) Use the genealized compessibility diagam povided on the last page to estimate ove what ange of pessues A at oom tempeatue confoms to the ideal gas law

More information

Encapsulation theory: radial encapsulation. Edmund Kirwan *

Encapsulation theory: radial encapsulation. Edmund Kirwan * Encapsulation theoy: adial encapsulation. Edmund Kiwan * www.edmundkiwan.com Abstact This pape intoduces the concept of adial encapsulation, wheeby dependencies ae constained to act fom subsets towads

More information

Goodness-of-fit for composite hypotheses.

Goodness-of-fit for composite hypotheses. Section 11 Goodness-of-fit fo composite hypotheses. Example. Let us conside a Matlab example. Let us geneate 50 obsevations fom N(1, 2): X=nomnd(1,2,50,1); Then, unning a chi-squaed goodness-of-fit test

More information

A Converse to Low-Rank Matrix Completion

A Converse to Low-Rank Matrix Completion A Convese to Low-Rank Matix Completion Daniel L. Pimentel-Alacón, Robet D. Nowak Univesity of Wisconsin-Madison Abstact In many pactical applications, one is given a subset Ω of the enties in a d N data

More information

Multiple Criteria Secretary Problem: A New Approach

Multiple Criteria Secretary Problem: A New Approach J. Stat. Appl. Po. 3, o., 9-38 (04 9 Jounal of Statistics Applications & Pobability An Intenational Jounal http://dx.doi.og/0.785/jsap/0303 Multiple Citeia Secetay Poblem: A ew Appoach Alaka Padhye, and

More information

Physics 211: Newton s Second Law

Physics 211: Newton s Second Law Physics 211: Newton s Second Law Reading Assignment: Chapte 5, Sections 5-9 Chapte 6, Section 2-3 Si Isaac Newton Bon: Januay 4, 1643 Died: Mach 31, 1727 Intoduction: Kinematics is the study of how objects

More information

Fresnel Diffraction. monchromatic light source

Fresnel Diffraction. monchromatic light source Fesnel Diffaction Equipment Helium-Neon lase (632.8 nm) on 2 axis tanslation stage, Concave lens (focal length 3.80 cm) mounted on slide holde, iis mounted on slide holde, m optical bench, micoscope slide

More information

Chapter 3 Optical Systems with Annular Pupils

Chapter 3 Optical Systems with Annular Pupils Chapte 3 Optical Systems with Annula Pupils 3 INTRODUCTION In this chapte, we discuss the imaging popeties of a system with an annula pupil in a manne simila to those fo a system with a cicula pupil The

More information

Lab #4: Newton s Second Law

Lab #4: Newton s Second Law Lab #4: Newton s Second Law Si Isaac Newton Reading Assignment: bon: Januay 4, 1643 Chapte 5 died: Mach 31, 1727 Chapte 9, Section 9-7 Intoduction: Potait of Isaac Newton by Si Godfey Knelle http://www.newton.cam.ac.uk/at/potait.html

More information

2 Governing Equations

2 Governing Equations 2 Govening Equations This chapte develops the govening equations of motion fo a homogeneous isotopic elastic solid, using the linea thee-dimensional theoy of elasticity in cylindical coodinates. At fist,

More information

CENTER FOR MULTIMODAL SOLUTIONS FOR CONGESTION MITIGATION (CMS)

CENTER FOR MULTIMODAL SOLUTIONS FOR CONGESTION MITIGATION (CMS) Final Repot to the CENTER FOR MULTIMODAL SOLUTIONS FOR CONGESTION MITIGATION (CMS) CMS Poect Numbe: _8-4_ Title: Chaacteizing the Tadeoffs and Costs Associated with Tanspotation Congestion in Supply Chains

More information

1 Explicit Explore or Exploit (E 3 ) Algorithm

1 Explicit Explore or Exploit (E 3 ) Algorithm 2.997 Decision-Making in Lage-Scale Systems Mach 3 MIT, Sping 2004 Handout #2 Lectue Note 9 Explicit Exploe o Exploit (E 3 ) Algoithm Last lectue, we studied the Q-leaning algoithm: [ ] Q t+ (x t, a t

More information

Vanishing lines in generalized Adams spectral sequences are generic

Vanishing lines in generalized Adams spectral sequences are generic ISSN 364-0380 (on line) 465-3060 (pinted) 55 Geomety & Topology Volume 3 (999) 55 65 Published: 2 July 999 G G G G T T T G T T T G T G T GG TT G G G G GG T T T TT Vanishing lines in genealized Adams spectal

More information

ON THE INVERSE SIGNED TOTAL DOMINATION NUMBER IN GRAPHS. D.A. Mojdeh and B. Samadi

ON THE INVERSE SIGNED TOTAL DOMINATION NUMBER IN GRAPHS. D.A. Mojdeh and B. Samadi Opuscula Math. 37, no. 3 (017), 447 456 http://dx.doi.og/10.7494/opmath.017.37.3.447 Opuscula Mathematica ON THE INVERSE SIGNED TOTAL DOMINATION NUMBER IN GRAPHS D.A. Mojdeh and B. Samadi Communicated

More information

Macro Theory B. The Permanent Income Hypothesis

Macro Theory B. The Permanent Income Hypothesis Maco Theoy B The Pemanent Income Hypothesis Ofe Setty The Eitan Beglas School of Economics - Tel Aviv Univesity May 15, 2015 1 1 Motivation 1.1 An econometic check We want to build an empiical model with

More information

Relating Branching Program Size and. Formula Size over the Full Binary Basis. FB Informatik, LS II, Univ. Dortmund, Dortmund, Germany

Relating Branching Program Size and. Formula Size over the Full Binary Basis. FB Informatik, LS II, Univ. Dortmund, Dortmund, Germany Relating Banching Pogam Size and omula Size ove the ull Binay Basis Matin Saueho y Ingo Wegene y Ralph Wechne z y B Infomatik, LS II, Univ. Dotmund, 44 Dotmund, Gemany z ankfut, Gemany sauehof/wegene@ls.cs.uni-dotmund.de

More information

Method for Approximating Irrational Numbers

Method for Approximating Irrational Numbers Method fo Appoximating Iational Numbes Eic Reichwein Depatment of Physics Univesity of Califonia, Santa Cuz June 6, 0 Abstact I will put foth an algoithm fo poducing inceasingly accuate ational appoximations

More information

On the ratio of maximum and minimum degree in maximal intersecting families

On the ratio of maximum and minimum degree in maximal intersecting families On the atio of maximum and minimum degee in maximal intesecting families Zoltán Lóánt Nagy Lale Özkahya Balázs Patkós Máté Vize Mach 6, 013 Abstact To study how balanced o unbalanced a maximal intesecting

More information

To Feel a Force Chapter 7 Static equilibrium - torque and friction

To Feel a Force Chapter 7 Static equilibrium - torque and friction To eel a oce Chapte 7 Chapte 7: Static fiction, toque and static equilibium A. Review of foce vectos Between the eath and a small mass, gavitational foces of equal magnitude and opposite diection act on

More information

On the Poisson Approximation to the Negative Hypergeometric Distribution

On the Poisson Approximation to the Negative Hypergeometric Distribution BULLETIN of the Malaysian Mathematical Sciences Society http://mathusmmy/bulletin Bull Malays Math Sci Soc (2) 34(2) (2011), 331 336 On the Poisson Appoximation to the Negative Hypegeometic Distibution

More information

Localization of Eigenvalues in Small Specified Regions of Complex Plane by State Feedback Matrix

Localization of Eigenvalues in Small Specified Regions of Complex Plane by State Feedback Matrix Jounal of Sciences, Islamic Republic of Ian (): - () Univesity of Tehan, ISSN - http://sciencesutaci Localization of Eigenvalues in Small Specified Regions of Complex Plane by State Feedback Matix H Ahsani

More information

Alternative Tests for the Poisson Distribution

Alternative Tests for the Poisson Distribution Chiang Mai J Sci 015; 4() : 774-78 http://epgsciencecmuacth/ejounal/ Contibuted Pape Altenative Tests fo the Poisson Distibution Manad Khamkong*[a] and Pachitjianut Siipanich [b] [a] Depatment of Statistics,

More information

working pages for Paul Richards class notes; do not copy or circulate without permission from PGR 2004/11/3 10:50

working pages for Paul Richards class notes; do not copy or circulate without permission from PGR 2004/11/3 10:50 woking pages fo Paul Richads class notes; do not copy o ciculate without pemission fom PGR 2004/11/3 10:50 CHAPTER7 Solid angle, 3D integals, Gauss s Theoem, and a Delta Function We define the solid angle,

More information

DonnishJournals

DonnishJournals DonnishJounals 041-1189 Donnish Jounal of Educational Reseach and Reviews. Vol 1(1) pp. 01-017 Novembe, 014. http:///dje Copyight 014 Donnish Jounals Oiginal Reseach Pape Vecto Analysis Using MAXIMA Savaş

More information

Math 2263 Solutions for Spring 2003 Final Exam

Math 2263 Solutions for Spring 2003 Final Exam Math 6 Solutions fo Sping Final Exam ) A staightfowad appoach to finding the tangent plane to a suface at a point ( x, y, z ) would be to expess the cuve as an explicit function z = f ( x, y ), calculate

More information

Information Retrieval Advanced IR models. Luca Bondi

Information Retrieval Advanced IR models. Luca Bondi Advanced IR models Luca Bondi Advanced IR models 2 (LSI) Pobabilistic Latent Semantic Analysis (plsa) Vecto Space Model 3 Stating point: Vecto Space Model Documents and queies epesented as vectos in the

More information

Lifting Private Information Retrieval from Two to any Number of Messages

Lifting Private Information Retrieval from Two to any Number of Messages Lifting Pivate Infomation Retieval fom Two to any umbe of Messages Rafael G.L. D Oliveia, Salim El Rouayheb ECE, Rutges Univesity, Piscataway, J Emails: d746@scaletmail.utges.edu, salim.elouayheb@utges.edu

More information

A NEW VARIABLE STIFFNESS SPRING USING A PRESTRESSED MECHANISM

A NEW VARIABLE STIFFNESS SPRING USING A PRESTRESSED MECHANISM Poceedings of the ASME 2010 Intenational Design Engineeing Technical Confeences & Computes and Infomation in Engineeing Confeence IDETC/CIE 2010 August 15-18, 2010, Monteal, Quebec, Canada DETC2010-28496

More information

Stochastic Analysis of Periodic Real-Time Systems *

Stochastic Analysis of Periodic Real-Time Systems * Stochastic Analysis of Peiodic Real-Time Systems * José Luis Díaz Daniel F Gacía Kanghee Kim Chang-Gun Lee Lucia Lo Bello José Maía López Sang Lyul Min Oazio Miabella Abstact This pape descibes a stochastic

More information

An Exact Solution of Navier Stokes Equation

An Exact Solution of Navier Stokes Equation An Exact Solution of Navie Stokes Equation A. Salih Depatment of Aeospace Engineeing Indian Institute of Space Science and Technology, Thiuvananthapuam, Keala, India. July 20 The pincipal difficulty in

More information

q i i=1 p i ln p i Another measure, which proves a useful benchmark in our analysis, is the chi squared divergence of p, q, which is defined by

q i i=1 p i ln p i Another measure, which proves a useful benchmark in our analysis, is the chi squared divergence of p, q, which is defined by CSISZÁR f DIVERGENCE, OSTROWSKI S INEQUALITY AND MUTUAL INFORMATION S. S. DRAGOMIR, V. GLUŠČEVIĆ, AND C. E. M. PEARCE Abstact. The Ostowski integal inequality fo an absolutely continuous function is used

More information

Numerical Integration

Numerical Integration MCEN 473/573 Chapte 0 Numeical Integation Fall, 2006 Textbook, 0.4 and 0.5 Isopaametic Fomula Numeical Integation [] e [ ] T k = h B [ D][ B] e B Jdsdt In pactice, the element stiffness is calculated numeically.

More information

d 2 x 0a d d =0. Relative to an arbitrary (accelerating frame) specified by x a = x a (x 0b ), the latter becomes: d 2 x a d 2 + a dx b dx c

d 2 x 0a d d =0. Relative to an arbitrary (accelerating frame) specified by x a = x a (x 0b ), the latter becomes: d 2 x a d 2 + a dx b dx c Chapte 6 Geneal Relativity 6.1 Towads the Einstein equations Thee ae seveal ways of motivating the Einstein equations. The most natual is pehaps though consideations involving the Equivalence Pinciple.

More information

MODULE 5a and 5b (Stewart, Sections 12.2, 12.3) INTRO: In MATH 1114 vectors were written either as rows (a1, a2,..., an) or as columns a 1 a. ...

MODULE 5a and 5b (Stewart, Sections 12.2, 12.3) INTRO: In MATH 1114 vectors were written either as rows (a1, a2,..., an) or as columns a 1 a. ... MODULE 5a and 5b (Stewat, Sections 2.2, 2.3) INTRO: In MATH 4 vectos wee witten eithe as ows (a, a2,..., an) o as columns a a 2... a n and the set of all such vectos of fixed length n was called the vecto

More information

( ) [ ] [ ] [ ] δf φ = F φ+δφ F. xdx.

( ) [ ] [ ] [ ] δf φ = F φ+δφ F. xdx. 9. LAGRANGIAN OF THE ELECTROMAGNETIC FIELD In the pevious section the Lagangian and Hamiltonian of an ensemble of point paticles was developed. This appoach is based on a qt. This discete fomulation can

More information

arxiv: v1 [physics.pop-ph] 3 Jun 2013

arxiv: v1 [physics.pop-ph] 3 Jun 2013 A note on the electostatic enegy of two point chages axiv:1306.0401v1 [physics.pop-ph] 3 Jun 013 A C Tot Instituto de Física Univesidade Fedeal do io de Janeio Caixa Postal 68.58; CEP 1941-97 io de Janeio,

More information

On the ratio of maximum and minimum degree in maximal intersecting families

On the ratio of maximum and minimum degree in maximal intersecting families On the atio of maximum and minimum degee in maximal intesecting families Zoltán Lóánt Nagy Lale Özkahya Balázs Patkós Máté Vize Septembe 5, 011 Abstact To study how balanced o unbalanced a maximal intesecting

More information

As is natural, our Aerospace Structures will be described in a Euclidean three-dimensional space R 3.

As is natural, our Aerospace Structures will be described in a Euclidean three-dimensional space R 3. Appendix A Vecto Algeba As is natual, ou Aeospace Stuctues will be descibed in a Euclidean thee-dimensional space R 3. A.1 Vectos A vecto is used to epesent quantities that have both magnitude and diection.

More information

Physics 2A Chapter 10 - Moment of Inertia Fall 2018

Physics 2A Chapter 10 - Moment of Inertia Fall 2018 Physics Chapte 0 - oment of netia Fall 08 The moment of inetia of a otating object is a measue of its otational inetia in the same way that the mass of an object is a measue of its inetia fo linea motion.

More information

Constrained Langevin Equations and (Bio)Chemical Reaction Networks

Constrained Langevin Equations and (Bio)Chemical Reaction Networks Constained Langevin Equations and BioChemical Reaction Netwoks David F. Andeson 1, Desmond J. Higham 2, Saul C. Leite 3, and Ruth J. Williams 4 Abstact Stochastic effects play an impotant ole in modeling

More information

Bounds on the performance of back-to-front airplane boarding policies

Bounds on the performance of back-to-front airplane boarding policies Bounds on the pefomance of bac-to-font aiplane boading policies Eitan Bachmat Michael Elin Abstact We povide bounds on the pefomance of bac-to-font aiplane boading policies. In paticula, we show that no

More information

Review: Electrostatics and Magnetostatics

Review: Electrostatics and Magnetostatics Review: Electostatics and Magnetostatics In the static egime, electomagnetic quantities do not vay as a function of time. We have two main cases: ELECTROSTATICS The electic chages do not change postion

More information

B. Spherical Wave Propagation

B. Spherical Wave Propagation 11/8/007 Spheical Wave Popagation notes 1/1 B. Spheical Wave Popagation Evey antenna launches a spheical wave, thus its powe density educes as a function of 1, whee is the distance fom the antenna. We

More information

Compactly Supported Radial Basis Functions

Compactly Supported Radial Basis Functions Chapte 4 Compactly Suppoted Radial Basis Functions As we saw ealie, compactly suppoted functions Φ that ae tuly stictly conditionally positive definite of ode m > do not exist The compact suppot automatically

More information

COLLAPSING WALLS THEOREM

COLLAPSING WALLS THEOREM COLLAPSING WALLS THEOREM IGOR PAK AND ROM PINCHASI Abstact. Let P R 3 be a pyamid with the base a convex polygon Q. We show that when othe faces ae collapsed (otated aound the edges onto the plane spanned

More information

Encapsulation theory: the transformation equations of absolute information hiding.

Encapsulation theory: the transformation equations of absolute information hiding. 1 Encapsulation theoy: the tansfomation equations of absolute infomation hiding. Edmund Kiwan * www.edmundkiwan.com Abstact This pape descibes how the potential coupling of a set vaies as the set is tansfomed,

More information

9.1 The multiplicative group of a finite field. Theorem 9.1. The multiplicative group F of a finite field is cyclic.

9.1 The multiplicative group of a finite field. Theorem 9.1. The multiplicative group F of a finite field is cyclic. Chapte 9 Pimitive Roots 9.1 The multiplicative goup of a finite fld Theoem 9.1. The multiplicative goup F of a finite fld is cyclic. Remak: In paticula, if p is a pime then (Z/p) is cyclic. In fact, this

More information