arxiv: v2 [math.pr] 9 Dec 2016
|
|
- Damon Townsend
- 6 years ago
- Views:
Transcription
1 Omnithermal Perfet Simulation for Multi-server Queues Stephen B. Connor 3th Deember 206 arxiv: v2 [math.pr] 9 De 206 Abstrat A number of perfet simulation algorithms for multi-server First Come First Served queues have reently been developed. hose of Connor and Kendall (205) and Blanhet, Pei, and Sigman (205) use dominated Coupling from the Past (domcfp) to sample from the equilibrium distribution of the Kiefer-Wolfowitz workload vetor for stable M/G/ and GI/GI/ queues respetively, using random assignment queues as dominating proesses. In this note we answer a question posed by Connor and Kendall (205), by demonstrating how these algorithms may be modified in order to arry out domcfp simultaneously for a range of values of (the number of servers), at minimal extra ost. Keywords and phrases: Dominated Coupling from the Past; First Come First Served disipline [FCFS]; Kiefer-Wolfowitz workload vetor; perfet simulation; M/G/ queue; Random Assignment disipline [RA]; sandwihing; stohasti ordering Mathematis Subjet Classifiation: Primary 65C05; Seondary 60K25; 60J05; 68U20 Introdution During the past five years there have been a number of signifiant advanes in perfet simulation methods for multi-server queues. Sigman (20) pioneered the use of dominated Coupling from the Past (domcfp) for super-stable M/G/ queues with First Come First Served (FCFS) disipline. ( Super-stable means that the queue would remain stable even if all but one of the servers were removed.) he limitation to super-stable queues is neessitated by Sigman s use of a stable M/G/ queue as dominating proess in the domcfp algorithm. Connor and Kendall (205) subsequently showed how to generalise this idea to work for stable M/G/ queues, by using as dominating proess an M/G/ queue with random assignment (RA) disipline (under whih the servers are independent). hey desribe two algorithms (outlined in Setion 2 below) and ompare their effiieny; they show that their Algorithm, whih requires waiting for the dominating proess to empty, is signifiantly less effiient than Algorithm 2, whih makes use of sandwihing proesses (in ommon with many other domcfp algorithms). Blanhet, Dong, and Pei (205) were the first authors to show how to perform perfet simulation for multi-server queues with general inter-arrival time and servie time distributions (i.e. relaxing the assumption of Exponential inter-arrival times). Rather than use a random assignment queue as dominating proess, they make use of a so-alled vaation system.
2 (his idea is also employed by Blanhet and Chen (206) to sample from the equilibrium of a generalized Jakson Network of single-server queues.) However, Blanhet et al. (205) have sine demonstrated how to make the random assignment dominating proess work in this setting. he hard part here is working out how to simulate the dominating proess in reverse-time; with renewal, as opposed to Poisson arrivals, the servers in the RA model are no longer independent. hese piees of work all serve to demonstrate that perfet simulation is a pratial and effiient method for simulating from a wide lass of multi-server queueing systems. Connor and Kendall (205) ask a very natural question: is it possible to arry out dominated CFP simultaneously for M/G/ queues with a range of, the number of servers? he authors refer to this as omnithermal dominated CFP, borrowing a term used to desribe Grimmett (995) s oupling of random-luster proesses for all values of a speifi parameter, and applied to CFP in Propp and Wilson (996). he potentially diffiult issue in the queueing ontext is that of deteting a time at whih we an be sure that the appropriate sandwihing proesses will oalese for all in the range being onsidered. In this paper we show how suh oalesene may be deteted with the aid of a simple riterion that uses information about the sandwihing proesses only for the queue with the least number of servers. he outline of the paper is as follows. In Setion 2 we reall the definition of the Kiefer- Wolfowitz workload proess assoiated to a multi-server FCFS queue, and then sketh the two perfet simulation algorithms of Connor and Kendall (205). In Setion 3 we present a natural partial order between Kiefer-Wolfowitz vetors of potentially different lengths; we subsequently use this to determine a ondition whih ensures that the termination time of Connor and Kendall (205) s Algorithm 2 is monotoni in the number of servers. In Setion 4 we use this ondition to produe an Omnithermal Algorithm, and briefly report on the results of applying this to some M/M/ queues. Finally, in Setion 5, we indiate how our results may be used to perform omnithermal perfet simulation for queues with general renewal input, or in the situation where we are interested in saling the distribution of servie durations, rather than hanging the number of servers. 2 Dominated CFP for M/G/ queues Consider a general / / FCFS queue. We denote the Kiefer-Wolfowitz workload vetor (Kiefer and Wolfowitz, 955) at time t 0 by V(t) = (V (, t), V (2, t),..., V (, t)), where V (, t) V (2, t)... represent the ordered amounts of residual work in the system for the servers at time t. Customer n arrives at time t n (for 0 t t 2...), with inter-arrival times denoted by n = t n+ t n (where we set t 0 = 0). Customer n brings with it a servie duration S n. Observing V just before arrival of the n th ustomer (but definitely after the arrival of the (n ) st ustomer) generates the proess W n : in the ase t n < t n we have W n = V(t n ). his satisfies the well-known reursion W n+ = R(W n + S n e n f) +, for n 0, () where e = (, 0, 0,..., 0), f = (,,..., ), R plaes the oordinates of a vetor in inreasing order, and + replaes negative oordinates of a vetor by zeros. For simpliity of exposition we shall primarily disuss M/G/ queues in what follows (i.e. inter-arrival times are Exponential). However, our method for performing omnithermal 2
3 perfet simulation for these queues applies equally well to GI/GI/ queues using an algorithm of Blanhet et al. (205), as will be observed in Setion 5. Let the arrival rate be λ > 0, and let servie durations S n be i.i.d. with mean /µ and E [ S 2] <. (As explained in Connor and Kendall (205), this seond moment ondition is required in order to guarantee finite mean run-time of their perfet simulation algorithms.) he queue is stable if and only if λ/(µ) <, and so we restrit attention to this senario. Connor and Kendall (205) propose two domcfp algorithms for sampling from the equilibrium distribution of the Kiefer-Wolfowitz workload vetor for a stable M/G/ queue X. Both of these algorithms use as dominating proess an M/G/ queue Y with Random Assignment (RA) servie disipline. hat is, ustomers in Y are alloated upon arrival to a uniformly hosen server; this renders the servers independent, whih allows us to easily simulate a stationary version of the dominating proess in reverse-time, as required by domcfp. It is possible to arrange for X to be path-wise dominated by Y as long as the two queues are oupled by assigning servie durations in order of initiation of servie. (Under FCFS ustomers initiate servie in the same order in whih they arrive, but this is typially not the ase for other servie disiplines.) he preise statement of this domination an be found in Connor and Kendall (205), an abridged version of whih is reprodued here for onveniene. heorem (heorem 3.3 of Connor and Kendall (205)). Consider a -server queueing system viewed as a funtion of (a) the sequene of arrival times 0 t t 2 t 3... and (b) the sequene of servie durations S, S 2, S 3,... assigned in order of initiation of servie. Consider the following different alloation rules, in some ases varying over time:. / / [RA]; 2. / / [RA] until a speified non-random time, then swithing to / / [F CF S]; 3. / / [RA] until a speified non-random time, 0, then swithing to / / [F CF S]; 4. / / [F CF S]; hen ase k dominates ase k +, in the sense that the m th initiation of servie in ase k + ours no later than the m th initiation of servie in ase k, and the m th departure in ase k + ours no later than the m th departure in ase k. Moreover, for all times t the Kiefer-Wolfowitz workload vetor for ase 3 dominates (oordinate-by-oordinate) that of ase 4, with similar domination holding for ases 2 and 3 for all t. We an now summarise the two domcfp algorithms of Connor and Kendall (205). Algorithm. Construt a stationary M/G/ [RA] proess bakwards in time until it empties at some time < 0; 2. Use this to reate a forwards in time trajetory of an M/G/ [RA] queue Y started from empty at time ; 3. Use the sequenes of arrival times and servie durations in Y to onstrut an M/G/ [F CF S] queue X that is dominated by Y over [, 0]; 4. Return X 0. 3
4 Algorithm 2. Fix a bakoff (or inspetion) time < 0, and onstrut a path of the stationary M/G/ [RA] queue Y over the time period [, 0]; 2. Construt sandwihing proesses L and U over [, 0] as follows: (a) L and U both evolve as the Kiefer-Wolfowitz vetors of M/G/ [F CF S] queues, using the same sequenes of arrival times and servie durations as Y over (, 0]; (b) L is empty, while U is instantiated using the same residual workloads present in Y ; 3. Chek for oalesene: if L 0 = U 0. return this value; else set 2 and return to Step Connor and Kendall (205) provide more details for eah of the steps outlined above, and demonstrate that Algorithm 2, although more ompliated to desribe, is in general signifiantly faster than Algorithm. 3 Omnithermal perfet simulation In this setion we onsider the following question: is it possible to adapt the domcfp algorithms outlined in Setion 2 in order to simultaneously sample from the equilibrium of M/G/( + m) queues, where m ranges over some subset of N? As pointed out by Connor and Kendall (205), it is straightforward to aomplish this in a relatively effiient manner using Algorithm : one an emptying time has been established for the M/G/ queue, then all M/G/( + m) queues of interest an be started from empty at this time and run over [, 0] using the same arrival times and servie durations; a simple workload domination argument shows that their values at time 0 will form a single perfet sample from the required set of equilibrium distributions. However, given the signifiantly faster run-time of Algorithm 2, a far more interesting question is whether or not one an produe a omparably effiient omnithermal domcfp algorithm using sandwihing proesses. Suppose that we have implemented Algorithm 2, and have obtained one equilibrium sample for the M/G/ queue. hat is, we have established some bakoff time < 0, along with sequenes of arrival times and servie durations, suh that L 0 = U 0. Our first observation is the following: suppose that we use these sequenes to produe new FCFS proesses L +m and U +m over [, 0] in the manner desribed in step 2 of Algorithm 2. (hat is, L +m is empty, and U +m ontains the residual workloads present in Y, the M/G/ [RA] proess.) hen the workload vetor for Ut +m dominates (oordinate-by-oordinate) that of L +m t for all t [, 0], and if L +m 0 = U0 +m then this value will be a perfet draw from the equilibrium of the M/G/( + m) queue, as required. his follows from heorem : due to the way in whih it is instantiated, U +m is a queueing system that hanges from M/G/ [RA] to M/G/(+m) [F CF S] at time. But the former of these an be thought of as an M/G/(+m) system with a random alloation rule that uniformly distributes jobs amongst only of the ( + m) servers; sine this is less effiient than the FCFS disipline, the proof of heorem 3.3 in Connor and Kendall (205) holds with this slightly modified setup. his observation implies that, given the arrival times and servie durations used in Algorithm 2 with servers, we ould just onstrut sandwihing proesses L +m and U +m over [, 0] and see whether they oalese. If they do, then we have obtained a sample from the 4
5 required distribution; if not, then we need to extend the dominating proess Y for this sample further into the past (setting 2 ), and then hek again for oalesene. But this method is not as lean as we would like: in partiular, it is likely that the extent to whih any single sample path of Y needs to be extended will vary with the value of m. (As will be shown in the next setion, oalesene of L +m and U +m over [, 0] does not imply oalesene of L +n and U +n over the same interval for all n > m.) Assuming that we want to obtain samples for a range of values of m N, this method is therefore rather ineffiient. Ideally we would like to use Algorithm 2 to produe a sample for the M/G/ queue, and then re-use the path of Y for this sample in order to draw from the equilibrium of M/G/( + m) for any desired m N. 3. Comparing queues with different numbers of servers Suppose that we have two FCFS queues, eah seeing the same set of arrival times and assoiated servie durations. We first of all need to show that the workload vetor with fewer servers dominates that of the other, with respet to a ertain natural partial order. Definition 2. For V R and V +m R +m, we write V +m V if and only if V +m (k + m) V (k), k =,...,. hus if V and V +m are workload vetors, V +m V if and only if eah of the busiest servers in V +m has no more work remaining than the orresponding server in V. Proposition 3. Let V and V +m be Kiefer-Wolfowitz workload vetors for an M/G/ and an M/G/( + m) FCFS queue respetively. Suppose that V0 +m V0 and that eah queue sees the same set of arrival times and assoiated servie durations. hen Vt +m Vt for all t 0. Proof. It is lear that the ordering between V and V +m will hold until the first arrival time, t. Furthermore, one we show that Vt +m Vt the result will follow simply by indution. Let S denote the servie duration attahed to the arrival at time t. Reall that the effet of this arrival is that S is added to any outstanding work at the first (least busy) oordinate in V and V +m, and then the resulting vetors are eah reordered in inreasing order. Suppose that after this reordering has taken plae, the oordinate with value S + V t () (the amount of work now at the server to whih the arrival at t was alloated) is loated in position i of Vt, et. If i +m m then the result is trivial (sine the last oordinates of V +m are unhanged by the arrival at time t ). So suppose that i +m > m. hen for k < min{i, i +m m} we have V +m t (k + m) = V +m (k + m + ) Vt (k + ) = Vt (k). t Similarly, for k > max{i, i +m m} we have V +m t (k + m) = Vt + (k + m) Vt (k) = Vt (k). here are now two ases to onsider, depending on whih of i and i +m m is larger. Case : i i +m m. hen for k = i,..., i +m m: V +m t (k + m) V +m (i +m ) = V +m () + S V t () + S = Vt (i ) Vt (k). t t Here the first and last inequalities hold sine the oordinates of the workload vetors at time t are arranged in inreasing order; the middle inequality follows from Vt +m V t. 5
6 Case 2: i +m m < i. Note that the following relations hold between V t and V t : { Vt Vt (k) = (k + ) k = i +m m,..., i and V +m t (k + m) = For k = i +m m,..., i it follows that V +m t Vt () + S { k = i Vt +m () + S k = i+m m Vt +m (k + m) k = i+m m +,..., i. (k + m) V +m (k + m + ) = V +m (k + m + ) Vt (k + ) = Vt (k). t t he proof is ompleted by observing that when k = i, V +m t 3.2 Coalesene (k + m) = V +m (i + m) Vt (i ) = Vt (i ) Vt (i ). t Suppose one again that we have used Algorithm 2 to obtain a single perfet sample from the M/G/ queue: this yields a bakoff time < 0 and a sequene of arrival times and assoiated servie durations suh that the sandwihing proesses U and L oalese over [, 0]. Define D to be the non-negative vetor-valued proess given by the oordinate-wise differene between U and L : D t = U t L t, t 0. Let be the oalesene time for this realisation: = inf{t > : D t = 0} < 0. We shall write L t for the number of ustomers in L t, and A t for the set of oordinates where U t and L t agree: A t = {k : D t (k) = 0, k }. We are interested in the question of whether oalesene of U and L implies oalesene of U +m and L +m (instantiated at time as desribed in Setion 2) over the same period. he following example shows that this is not guaranteed. Example 4. Consider sandwihing proesses for two and three server systems (i.e. = 2 and m = ), as desribed above. Suppose that U 2 and U 3 are both initiated at time 0 with a single servie duration of length, and that these queues proeed to see arrival times/servies (t, S) as follows: (0.,.2), (0.3,.8), (0.8, 5). he evolution of these proesses viewed at arrival times is as follows: t 0 = 0 t = 0. t 2 = 0.3 t 3 = 0.8 U 2 (0.0,.0) (0.9,.2) (.0, 2.5) (2.0, 5.5) L 2 (0.0, 0.0) (0.0,.2) (.0,.8) (.3, 5.5) If there are no further arrivals within the next two units of time, we see that U 2 and L 2 will oalese at time 2 = 2.8 (sine it will take two more units of time for their first oordinates to agree, and their seond oordinates are already mathed). However, feeding the same sequene of arrival times/servies to U 3 and L 3, we see that they will not oalese before time 2 : 6
7 t 0 = 0 t = 0. t 2 = 0.3 t 3 = 0.8 U 3 (0.0, 0.0,.0) (0.0, 0.9,.2) (0.7,.0,.8) (0.5,.3, 5.2) L 3 (0.0, 0.0, 0.0) (0.0, 0.0,.2) (0.0,.0,.8) (0.5,.3, 5.0) Furthermore, if we were to onsider sandwihing proesses for a four-server system, these would oalese by time 2 using the above sequene of arrivals: t 0 = 0 t = 0. t 2 = 0.3 t 3 = 0.8 U 4 (0.0, 0.0, 0.0,.0) (0.0, 0.0, 0.9,.2) (0.0, 0.7,.0,.8) (0.2, 0.5,.3, 5.0) L 4 (0.0, 0.0, 0.0, 0.0) (0.0, 0.0, 0.0,.2) (0.0, 0.0,.0,.8) (0.0, 0.5,.3, 5.0) A simple, and intuitively obvious, ondition whih guarantees that the sandwihing proesses L +m and U +m will oalese by time is that no ustomer arriving at the lower proess L during the period [, ] has to wait to ommene servie: Proposition 5. If L t for all t [, ] then D +m m N. = 0 (and so +m ) for any Proof. Sine L +m is started from empty, the ondition implies that its first m oordinates are identially zero over [, ]. Moreover, oalesene of U and L implies that there must exist an empty server in both of these proesses at time (see Connor and Kendall (205)); i.e. U () = L +m () = 0. Sine U U +m, Proposition 3 ensures that U U, and so the first m oordinates of U +m (and of L +m ) must also be equal to zero. It remains to show that the last oordinates of U +m and L +m agree. Sine no server in L ever has more than one ustomer to deal with at any moment, the same is true for L +m, and so L +m t (k + m) = L t(k) for all k and for all t [, ]. hen by the domination established in Proposition 3, and the fat that U = L, for all k, as required. U +m (k + m) U (k) = L (k) = L+m (k + m), he ondition of Proposition 5 is rather strong, and an in fat be weakened, as we now show. heorem 6. Suppose that, for all arrival times τ [, ], if D τ () = 0 (equivalently, A τ ) then U τ () = 0. hen +m for any m N. Remark 7. he ondition of Proposition 5 is stronger than that of heorem 6. o see this, suppose that L t for all t [, ]. If at some arrival time τ [, ] we have D τ () = 0 but U τ () > 0, then there must be at least ustomers in L τ (sine L τ () = U τ () > 0). But then the ustomer arriving at time τ would fore L τ = +, whih would break our initial assumption. herefore if D τ () = 0 it must be the ase that U τ () = 0. Note that in Example 4 the ondition of Proposition 5 learly fails, and so therefore does the ondition of heorem 6. (In fat U 2 t 3 () = L 2 t 3 () > 0.) he key to proving heorem 6 is to onsider the time until oalesene of the sandwihing proesses U and L when viewed at time t, under the assumption of there being no further 7
8 arrivals in (t, ). Note that this is just the time taken for U to lear all work in oordinates whih disagree with those in L at time t. Let us write C t for this quantity: Ct = max U k / A t (k) = Ut (n t), (2) t where we define n t = max{ k : k / A t}. It is lear that C t dereases deterministially at unit rate until it either hits zero (at whih point L and U oalese) or a new ustomer arrives. Consider then what happens to C if there is an arrival at time τ with assoiated servie duration S. Let k U and k L be the oordinates satisfying U τ (k U ) = U τ () + S and L τ (k L ) = L τ () + S. (hat is, the arriving job gets alloated to the server with the least work in eah of U τ and L τ, and then when the workload vetors are reordered this job finds itself in position k U in U τ and k L in L τ.) o be expliit k U = min{k : U τ () + S U τ (k + ), k < }, with k U = if the minimum above is taken over the empty set. Note that this onvention that the new job is plaed at the lowest oordinate possible, after reordering, in U τ allows us to deal with the possibility that U τ () + S = U τ (k + ) for some k <, i.e. of U τ ontaining two equal but non-zero oordinates; when arrivals are Poisson this possibility ours with probability zero, of ourse. In partiular, this implies that U τ (k U ) > U τ (k U ). (3) here are two ases to onsider when assessing the impat of an arrival on C, depending on whether or not the servers with least workload in U τ and L τ are in agreement. Case : A τ (i) Suppose first that k U n τ. Sine U τ (k) = L τ (k) for all k > n τ, it must be the ase that k L = k U A τ. So n τ = n τ and C τ = U τ (n τ ) = U τ (n τ ) = C τ. (ii) Alternatively, if k U < n τ then n τ = n τ, and so C τ = C τ one again. hus there is no hange to C if the arriving ustomer finds U τ () = L τ (). Case 2: / A τ (i) Suppose that k U n τ. Sine Uτ (k) = L τ (k) for all k > n τ, it must be the ase that k L k U. We laim that k U / A τ, and so n τ = k U ; it then follows that Cτ = Uτ (n τ ) = Uτ (k U ) = Uτ () + S > Cτ. o see that n τ = k U, we need to show that L τ (k U ) < Uτ (k U ). Notie that L τ (k U ) = max{l τ (k U ), L τ () + S} and U τ (k U ) = U τ () + S. Clearly L τ () + S < Uτ () + S (sine / A τ ). Uτ (k U ) < Uτ (k U ) thanks to (3). Furthermore, L τ (k U ) (ii) Alternatively, if k U < n τ then k L n τ also, and so n τ = n τ. hus C τ = C τ. hus when / A τ, C τ = max{u τ (n τ ), U τ () + S}. 8
9 In summary, we see that the time to oalesene of L and U inreases exatly when there is an arrival at time τ, with / A τ and k U n τ. hat is, { Cτ Cτ if A τ = max{cτ, Uτ () + S} if / A (4) τ. he next result is key to the proof of heorem 6: it shows that, under the same assumption as the theorem, the time to oalesene with + m servers is dominated by the time to oalesene with servers. Lemma 8. Suppose that at arrival time τ, if A τ then Uτ () = 0. Suppose also that C τ for some m N. hen Cτ +m Cτ. C +m τ Proof. We onsider the two possible senarios seen by the ustomer arriving at time τ.. A +m τ. 2. / A +m τ and / A τ. (Note that the third possibility, that / A +m τ and A τ, is exluded by our assumption. Indeed, if A τ then our assumption would fore L τ () = Uτ () = 0. So the arrival at time τ would find a server empty in U, whih would mean that it must also find a server empty in U +m. But that would imply that Aτ +m.) Let s treat these two possibilities in order.. Sine A +m t, we know from (4) that the oalesene time for L+m and U +m is unhanged by the new arrival. In addition, the oalesene time for L and U annot derease due to this arrival. So C +m τ = C +m τ C τ C τ. 2. Here the arrival potentially affets the time-to-oalesene of both pairs of sandwihing proesses. However, C +m τ = max{cτ +m, U τ +m () + S} max{c τ, Uτ () + S} = Cτ, where the inequality follows from the seond assumption of the Lemma, and the previously established fat that Uτ +m U τ. We an now omplete the proof of heorem 6. Reall that the sandwihing proesses L +m and U +m are started at time < 0 with L +m empty and U +m instantiated using the same set of residual workloads that are present in U. Now, it is lear that departures in U +m our no later than in U, and sine C is simply the time taken for all ustomers present in U to depart it follows that C +m C. Given that the assumption of heorem 6 holds, Lemma 8 tells us that this ordering is preserved for all t [, ]: Ct +m Ct, t [, ]. But sine L and U oalese at time < 0, we see that C +m = C = 0, and so +m, as laimed. 9
10 4 Simulations he result of heorem 6 provides us with a reipe for performing omnithermal perfet simulation for M/G/( + m) queues, for m = 0,, 2,.... Omnithermal Algorithm. Use Algorithm 2 of Connor and Kendall (205) to establish a bakoff time < 0 and upper and lower sandwihing proesses U and L over [, 0] suh that U 0 = L 0 : (i) Calulate the oalesene time [, 0] of U and L, and hek to see whether the ondition of heorem 6 is satisfied for all arrival times in [, ]; (ii) While the ondition is not satisfied, set 2, and use Algorithm 2 of Connor and Kendall (205) to extend the simulation of the sandwihing proesses over the new window [, 0]. 2. For eah required m N, onstrut L +m over [, 0], using the same sequene of arrival times and servies as in the onstrution of L. Return L0 +m as a perfet equilibrium draw of the Kiefer-Wolfowitz vetor for the M/G/( + m) queue of interest. (Note that there is no need to also onstrut U +m : heorem 6 guarantees that U0 +m = L +m 0.) Remark 9. If we are alled upon to use part (ii) of the algorithm and extend the simulations of U and L further into the past, we are guaranteed that these new sandwihing proesses (Ũ and L, say) will still oalese by time : this follows from heorem 5. of Connor and Kendall (205), whih implies that L t L t Ũ t U t, t [, 0]. We applied the Omnithermal Algorithm to the three sets of simulations presented in Figure 4 of Connor and Kendall (205); these all onern M/M/ queues with λ = and µ = 2. 5,000 runs were performed for eah set of parameters (λ = 0, 30, 50): for eah run we reorded the original bakoff time (required for oalesene of U and L ), and the number of times that this was doubled in part (ii) in order to satisfy the ondition of heorem 6. When λ = 30 (respetively 50) we found that only 4 (respetively 0) of the original bakoff times needed to be extended; when λ = 0 this number rose to 08 (2.% of the sample). Furthermore, those runs whih needed extending did not require to be doubled a signifiant number of times: no run required us to double more than three times, with most of the runs requiring only a single additional bakoff in order to satisfy the ondition of heorem 6 see Figure. In addition, we used the Omnithermal Algorithm to investigate the effet of hanging server number for an M/M/ queue with arrival rate λ =.2, servie rate µ = and number of servers {2, 3, 4}. We ran the algorithm 5,000 times, for whih 94 runs failed the ondition of heorem 6 and needed extending further into the past. Figure 2 shows the mean value of eah oordinate of the Kiefer-Wolfowitz workload vetor from this sample. he effet on the value of the first oordinate (whih represents the mean waiting time of a ustomer arriving in equilibrium) of inreasing the number of servers from two to three an be seen to be substantial: the mean value drops from 0.57 to just Further detail is provided in Figure 3, where we show the effet on the distribution funtion of the remaining workload at the first two oordinates of the workload vetors for the same set of simulations. 0
11 70 Frequeny Original value of Number of additional bakoffs required Figure : Simulation results for 5,000 runs of the Omnithermal Algorithm for an M/M/ queue with λ = = 0 and µ = runs failed the ondition of heorem 6 and needed extending further into the past; the figure shows the number of times eah of these 08 runs needed to be extended (using the usual binary bakoff sheme). No run needed to be extended more than three times. 5 Conluding omments We have shown how the effiient Algorithm 2 of Connor and Kendall (205) for M/G/ queues may be modified to allow for omnithermal perfet simulation; our new algorithm uses a simple test to determine whether or not the dominating proess used for the -server algorithm needs to be extended further into the past in order to allow for simultaneous sampling from M/G/( + m) queues for any desired m N. Furthermore, we have provided numerial evidene whih suggests that, at least in the ase of moderate traffi intensity, the Omnithermal Algorithm requires minimal additional omputational expense. We onlude by noting two extensions of our algorithm. Our first observation is that the oalesene arguments underpinning Setion 3 do not rely in any way on the distribution of inter-arrival times. As noted in the introdution, Blanhet et al. (205) have reently shown how to implement domcfp for GI/GI/ queues using a random assignment dominating.4.2 Number of servers Mean workload Workload vetor oordinate Figure 2: Mean of eah oordinate of the workload vetor for an M/M/ queue with λ =.2, µ = and {2, 3, 4}. (Results from 5,000 runs of the Omnithermal Algorithm.)
12 .0.0 Distribution funtion Number of servers Distribution funtion Number of servers Workload (first oordinate) Workload (seond oordinate) Figure 3: Distribution funtions for workload at first (left) and seond (right) oordinates of the workload vetor, for the set of simulations presented in Figure 2. proess with upper and lower sandwihing proesses in the style of Algorithm 2 above. It is therefore possible to perform omnithermal simulation for these queues, by using their algorithm in plae of Algorithm 2 in step of the Omnithermal Algorithm. he seond natural extension is to omnithermal simulation of queues whih are more stable than the M/G/ on aount of having shorter servie durations, but the same number of servers. (We ould equivalently onsider queues with longer inter-arrival times of ourse; however, for the domination arguments of Setion 3 to hold it is essential that the two systems being ompared have the same set of arrival times. It is therefore more onvenient to adjust the servie durations instead.) Suppose that servie times in the more stable system are distributed as βs for some β (0, ]: this is equivalent to the servie times being distributed as S, but with eah server now ompleting work at rate β. So we an ompare the two systems as in Setion 3, feeding both the same sets of arrival times and servie durations, but with the time-to-oalesene in (2) replaed by C β t = βu β t (nβ t ) (where we have one again used the supersript to indiate the parameter that varies between the queues under onsideration). In a similar manner to Example 4, it is easy to onjure up a set of arrival times and servie durations suh that the system ompleting work at rate β has C β t > Ct for some values of t. However, we note that in this new setting equation (4) beomes { Cτ β C β τ if A β τ = max{c β τ, β(u β τ () + S)} if / Aβ τ. Using this, it is a simple exerise to hek that if the ondition of heorem 6 is satisfied, then the sandwihing proesses for the faster-working system will oalese no later than do U and L for the original M/G/ queue. In other words, we an perform omnithermal simulation in this setting by simply replaing step 2 of the Omnithermal Algorithm with the following variant: 2. For eah required β (0, ), onstrut L β over [, 0], using the same set of arrival times and servies as in the onstrution of L. Return L β 0 as a perfet equilibrium draw of the Kiefer-Wolfowitz vetor for the M/G/ queue, in whih work is ompleted at rate β. 2
13 (Step in whih we possibly extend some simulations further into the past does not hange at all.) In onlusion, omnithermal perfet simulation for multi-server queues under FCFS alloation may be performed in a pratial and effiient manner; moreover, one we have done the work of extending any runs of the dominating proess further into the past, the resulting simulations may be used to produe perfet samples when either the number of servers or the rate at whih work is performed is inreased. Referenes Blanhet, J. and X. Chen (206). arxiv Perfet Sampling of Generalized Jakson Network. Blanhet, J., J. Dong, and Y. Pei (205). Perfet Sampling of GI/GI/ Queues. arxiv Blanhet, J., Y. Pei, and K. Sigman (205). Exat sampling for some multi-dimensional queueing models with renewal input. arxiv Connor, S. B. and W. S. Kendall (205). Perfet simulation of M/G/ queues. Advanes in Applied Probability 47 (4), Grimmett, G. (995). he stohasti random-luster proess and the uniqueness of randomluster measures. he Annals of Probability 23 (4), Kiefer, J. and J. Wolfowitz (955). On the theory of queues with many servers. ransations of the Amerian Mathematial Soiety 8 (), 8. Propp, J. G. and D. B. Wilson (996). Exat sampling with oupled Markov hains and appliations to statistial mehanis. Random Strutures and Algorithms 9, Sigman, K. (20). Exat Simulation of the Stationary Distribution of the FIFO M/G/ Queue. Journal of Applied Probability 48A,
The Hanging Chain. John McCuan. January 19, 2006
The Hanging Chain John MCuan January 19, 2006 1 Introdution We onsider a hain of length L attahed to two points (a, u a and (b, u b in the plane. It is assumed that the hain hangs in the plane under a
More informationA Queueing Model for Call Blending in Call Centers
A Queueing Model for Call Blending in Call Centers Sandjai Bhulai and Ger Koole Vrije Universiteit Amsterdam Faulty of Sienes De Boelelaan 1081a 1081 HV Amsterdam The Netherlands E-mail: {sbhulai, koole}@s.vu.nl
More informationOmnithermal perfect simulation for multi-server queues
Omnithermal perfect simulation for multi-server queues Stephen Connor stephen.connor@york.ac.uk LMS-EPSRC Durham Symposium July-August 2017 Dominated CFTP in a nutshell Suppose that we re interested in
More informationMillennium Relativity Acceleration Composition. The Relativistic Relationship between Acceleration and Uniform Motion
Millennium Relativity Aeleration Composition he Relativisti Relationship between Aeleration and niform Motion Copyright 003 Joseph A. Rybzyk Abstrat he relativisti priniples developed throughout the six
More informationA Characterization of Wavelet Convergence in Sobolev Spaces
A Charaterization of Wavelet Convergene in Sobolev Spaes Mark A. Kon 1 oston University Louise Arakelian Raphael Howard University Dediated to Prof. Robert Carroll on the oasion of his 70th birthday. Abstrat
More informationComplexity of Regularization RBF Networks
Complexity of Regularization RBF Networks Mark A Kon Department of Mathematis and Statistis Boston University Boston, MA 02215 mkon@buedu Leszek Plaskota Institute of Applied Mathematis University of Warsaw
More informationThe Effectiveness of the Linear Hull Effect
The Effetiveness of the Linear Hull Effet S. Murphy Tehnial Report RHUL MA 009 9 6 Otober 009 Department of Mathematis Royal Holloway, University of London Egham, Surrey TW0 0EX, England http://www.rhul.a.uk/mathematis/tehreports
More informationMaximum Entropy and Exponential Families
Maximum Entropy and Exponential Families April 9, 209 Abstrat The goal of this note is to derive the exponential form of probability distribution from more basi onsiderations, in partiular Entropy. It
More informationRemark 4.1 Unlike Lyapunov theorems, LaSalle s theorem does not require the function V ( x ) to be positive definite.
Leture Remark 4.1 Unlike Lyapunov theorems, LaSalle s theorem does not require the funtion V ( x ) to be positive definite. ost often, our interest will be to show that x( t) as t. For that we will need
More informationComputer Science 786S - Statistical Methods in Natural Language Processing and Data Analysis Page 1
Computer Siene 786S - Statistial Methods in Natural Language Proessing and Data Analysis Page 1 Hypothesis Testing A statistial hypothesis is a statement about the nature of the distribution of a random
More informationMethods of evaluating tests
Methods of evaluating tests Let X,, 1 Xn be i.i.d. Bernoulli( p ). Then 5 j= 1 j ( 5, ) T = X Binomial p. We test 1 H : p vs. 1 1 H : p>. We saw that a LRT is 1 if t k* φ ( x ) =. otherwise (t is the observed
More informationOptimization of Statistical Decisions for Age Replacement Problems via a New Pivotal Quantity Averaging Approach
Amerian Journal of heoretial and Applied tatistis 6; 5(-): -8 Published online January 7, 6 (http://www.sienepublishinggroup.om/j/ajtas) doi:.648/j.ajtas.s.65.4 IN: 36-8999 (Print); IN: 36-96 (Online)
More informationDIGITAL DISTANCE RELAYING SCHEME FOR PARALLEL TRANSMISSION LINES DURING INTER-CIRCUIT FAULTS
CHAPTER 4 DIGITAL DISTANCE RELAYING SCHEME FOR PARALLEL TRANSMISSION LINES DURING INTER-CIRCUIT FAULTS 4.1 INTRODUCTION Around the world, environmental and ost onsiousness are foring utilities to install
More informationNonreversibility of Multiple Unicast Networks
Nonreversibility of Multiple Uniast Networks Randall Dougherty and Kenneth Zeger September 27, 2005 Abstrat We prove that for any finite direted ayli network, there exists a orresponding multiple uniast
More informationTaste for variety and optimum product diversity in an open economy
Taste for variety and optimum produt diversity in an open eonomy Javier Coto-Martínez City University Paul Levine University of Surrey Otober 0, 005 María D.C. Garía-Alonso University of Kent Abstrat We
More informationJournal of Inequalities in Pure and Applied Mathematics
Journal of Inequalities in Pure and Applied Mathematis A NEW ARRANGEMENT INEQUALITY MOHAMMAD JAVAHERI University of Oregon Department of Mathematis Fenton Hall, Eugene, OR 97403. EMail: javaheri@uoregon.edu
More informationEstimating the probability law of the codelength as a function of the approximation error in image compression
Estimating the probability law of the odelength as a funtion of the approximation error in image ompression François Malgouyres Marh 7, 2007 Abstrat After some reolletions on ompression of images using
More informationDiscrete Bessel functions and partial difference equations
Disrete Bessel funtions and partial differene equations Antonín Slavík Charles University, Faulty of Mathematis and Physis, Sokolovská 83, 186 75 Praha 8, Czeh Republi E-mail: slavik@karlin.mff.uni.z Abstrat
More informationOn Component Order Edge Reliability and the Existence of Uniformly Most Reliable Unicycles
Daniel Gross, Lakshmi Iswara, L. William Kazmierzak, Kristi Luttrell, John T. Saoman, Charles Suffel On Component Order Edge Reliability and the Existene of Uniformly Most Reliable Uniyles DANIEL GROSS
More informationWhere as discussed previously we interpret solutions to this partial differential equation in the weak sense: b
Consider the pure initial value problem for a homogeneous system of onservation laws with no soure terms in one spae dimension: Where as disussed previously we interpret solutions to this partial differential
More informationSingular Event Detection
Singular Event Detetion Rafael S. Garía Eletrial Engineering University of Puerto Rio at Mayagüez Rafael.Garia@ee.uprm.edu Faulty Mentor: S. Shankar Sastry Researh Supervisor: Jonathan Sprinkle Graduate
More informationHankel Optimal Model Order Reduction 1
Massahusetts Institute of Tehnology Department of Eletrial Engineering and Computer Siene 6.245: MULTIVARIABLE CONTROL SYSTEMS by A. Megretski Hankel Optimal Model Order Redution 1 This leture overs both
More informationmax min z i i=1 x j k s.t. j=1 x j j:i T j
AM 221: Advaned Optimization Spring 2016 Prof. Yaron Singer Leture 22 April 18th 1 Overview In this leture, we will study the pipage rounding tehnique whih is a deterministi rounding proedure that an be
More informationControl Theory association of mathematics and engineering
Control Theory assoiation of mathematis and engineering Wojieh Mitkowski Krzysztof Oprzedkiewiz Department of Automatis AGH Univ. of Siene & Tehnology, Craow, Poland, Abstrat In this paper a methodology
More informationLikelihood-confidence intervals for quantiles in Extreme Value Distributions
Likelihood-onfidene intervals for quantiles in Extreme Value Distributions A. Bolívar, E. Díaz-Franés, J. Ortega, and E. Vilhis. Centro de Investigaión en Matemátias; A.P. 42, Guanajuato, Gto. 36; Méxio
More informationProduct Policy in Markets with Word-of-Mouth Communication. Technical Appendix
rodut oliy in Markets with Word-of-Mouth Communiation Tehnial Appendix August 05 Miro-Model for Inreasing Awareness In the paper, we make the assumption that awareness is inreasing in ustomer type. I.e.,
More informationAverage Rate Speed Scaling
Average Rate Speed Saling Nikhil Bansal David P. Bunde Ho-Leung Chan Kirk Pruhs May 2, 2008 Abstrat Speed saling is a power management tehnique that involves dynamially hanging the speed of a proessor.
More informationChapter 8 Hypothesis Testing
Leture 5 for BST 63: Statistial Theory II Kui Zhang, Spring Chapter 8 Hypothesis Testing Setion 8 Introdution Definition 8 A hypothesis is a statement about a population parameter Definition 8 The two
More informationLecture 7: Sampling/Projections for Least-squares Approximation, Cont. 7 Sampling/Projections for Least-squares Approximation, Cont.
Stat60/CS94: Randomized Algorithms for Matries and Data Leture 7-09/5/013 Leture 7: Sampling/Projetions for Least-squares Approximation, Cont. Leturer: Mihael Mahoney Sribe: Mihael Mahoney Warning: these
More informationRobust Recovery of Signals From a Structured Union of Subspaces
Robust Reovery of Signals From a Strutured Union of Subspaes 1 Yonina C. Eldar, Senior Member, IEEE and Moshe Mishali, Student Member, IEEE arxiv:87.4581v2 [nlin.cg] 3 Mar 29 Abstrat Traditional sampling
More informationOn the Licensing of Innovations under Strategic Delegation
On the Liensing of Innovations under Strategi Delegation Judy Hsu Institute of Finanial Management Nanhua University Taiwan and X. Henry Wang Department of Eonomis University of Missouri USA Abstrat This
More informationIndian Institute of Technology Bombay. Department of Electrical Engineering. EE 325 Probability and Random Processes Lecture Notes 3 July 28, 2014
Indian Institute of Tehnology Bombay Department of Eletrial Engineering Handout 5 EE 325 Probability and Random Proesses Leture Notes 3 July 28, 2014 1 Axiomati Probability We have learned some paradoxes
More informationMath 151 Introduction to Eigenvectors
Math 151 Introdution to Eigenvetors The motivating example we used to desrie matrixes was landsape hange and vegetation suession. We hose the simple example of Bare Soil (B), eing replaed y Grasses (G)
More informationA model for measurement of the states in a coupled-dot qubit
A model for measurement of the states in a oupled-dot qubit H B Sun and H M Wiseman Centre for Quantum Computer Tehnology Centre for Quantum Dynamis Griffith University Brisbane 4 QLD Australia E-mail:
More informationOrdered fields and the ultrafilter theorem
F U N D A M E N T A MATHEMATICAE 59 (999) Ordered fields and the ultrafilter theorem by R. B e r r (Dortmund), F. D e l o n (Paris) and J. S h m i d (Dortmund) Abstrat. We prove that on the basis of ZF
More informationModel-based mixture discriminant analysis an experimental study
Model-based mixture disriminant analysis an experimental study Zohar Halbe and Mayer Aladjem Department of Eletrial and Computer Engineering, Ben-Gurion University of the Negev P.O.Box 653, Beer-Sheva,
More informationThe shape of a hanging chain. a project in the calculus of variations
The shape of a hanging hain a projet in the alulus of variations April 15, 218 2 Contents 1 Introdution 5 2 Analysis 7 2.1 Model............................... 7 2.2 Extremal graphs.........................
More informationA RUIN MODEL WITH DEPENDENCE BETWEEN CLAIM SIZES AND CLAIM INTERVALS
A RUIN MODEL WITH DEPENDENCE BETWEEN CLAIM SIZES AND CLAIM INTERVALS Hansjörg Albreher a, Onno J. Boxma b a Graz University of Tehnology, Steyrergasse 3, A-8 Graz, Austria b Eindhoven University of Tehnology
More informationOn the Quantum Theory of Radiation.
Physikalishe Zeitshrift, Band 18, Seite 121-128 1917) On the Quantum Theory of Radiation. Albert Einstein The formal similarity between the hromati distribution urve for thermal radiation and the Maxwell
More informationOn the density of languages representing finite set partitions
On the density of languages representing finite set partitions Nelma Moreira Rogério Reis Tehnial Report Series: DCC-04-07 Departamento de Ciênia de Computadores Fauldade de Ciênias & Laboratório de Inteligênia
More informationAharonov-Bohm effect. Dan Solomon.
Aharonov-Bohm effet. Dan Solomon. In the figure the magneti field is onfined to a solenoid of radius r 0 and is direted in the z- diretion, out of the paper. The solenoid is surrounded by a barrier that
More information23.1 Tuning controllers, in the large view Quoting from Section 16.7:
Lesson 23. Tuning a real ontroller - modeling, proess identifiation, fine tuning 23.0 Context We have learned to view proesses as dynami systems, taking are to identify their input, intermediate, and output
More informationSimplified Buckling Analysis of Skeletal Structures
Simplified Bukling Analysis of Skeletal Strutures B.A. Izzuddin 1 ABSRAC A simplified approah is proposed for bukling analysis of skeletal strutures, whih employs a rotational spring analogy for the formulation
More information3 Tidal systems modelling: ASMITA model
3 Tidal systems modelling: ASMITA model 3.1 Introdution For many pratial appliations, simulation and predition of oastal behaviour (morphologial development of shorefae, beahes and dunes) at a ertain level
More informationCapacity Pooling and Cost Sharing among Independent Firms in the Presence of Congestion
Capaity Pooling and Cost Sharing among Independent Firms in the Presene of Congestion Yimin Yu Saif Benjaafar Graduate Program in Industrial and Systems Engineering Department of Mehanial Engineering University
More informationSimplification of Network Dynamics in Large Systems
Simplifiation of Network Dynamis in Large Systems Xiaojun Lin and Ness B. Shroff Shool of Eletrial and Computer Engineering Purdue University, West Lafayette, IN 47906, U.S.A. Email: {linx, shroff}@en.purdue.edu
More information22.54 Neutron Interactions and Applications (Spring 2004) Chapter 6 (2/24/04) Energy Transfer Kernel F(E E')
22.54 Neutron Interations and Appliations (Spring 2004) Chapter 6 (2/24/04) Energy Transfer Kernel F(E E') Referenes -- J. R. Lamarsh, Introdution to Nulear Reator Theory (Addison-Wesley, Reading, 1966),
More informationBäcklund Transformations: Some Old and New Perspectives
Bäklund Transformations: Some Old and New Perspetives C. J. Papahristou *, A. N. Magoulas ** * Department of Physial Sienes, Helleni Naval Aademy, Piraeus 18539, Greee E-mail: papahristou@snd.edu.gr **
More informationLecture 3 - Lorentz Transformations
Leture - Lorentz Transformations A Puzzle... Example A ruler is positioned perpendiular to a wall. A stik of length L flies by at speed v. It travels in front of the ruler, so that it obsures part of the
More informationPhysical Laws, Absolutes, Relative Absolutes and Relativistic Time Phenomena
Page 1 of 10 Physial Laws, Absolutes, Relative Absolutes and Relativisti Time Phenomena Antonio Ruggeri modexp@iafria.om Sine in the field of knowledge we deal with absolutes, there are absolute laws that
More informationMeasuring & Inducing Neural Activity Using Extracellular Fields I: Inverse systems approach
Measuring & Induing Neural Ativity Using Extraellular Fields I: Inverse systems approah Keith Dillon Department of Eletrial and Computer Engineering University of California San Diego 9500 Gilman Dr. La
More informationWavetech, LLC. Ultrafast Pulses and GVD. John O Hara Created: Dec. 6, 2013
Ultrafast Pulses and GVD John O Hara Created: De. 6, 3 Introdution This doument overs the basi onepts of group veloity dispersion (GVD) and ultrafast pulse propagation in an optial fiber. Neessarily, it
More informationIMPEDANCE EFFECTS OF LEFT TURNERS FROM THE MAJOR STREET AT A TWSC INTERSECTION
09-1289 Citation: Brilon, W. (2009): Impedane Effets of Left Turners from the Major Street at A TWSC Intersetion. Transportation Researh Reord Nr. 2130, pp. 2-8 IMPEDANCE EFFECTS OF LEFT TURNERS FROM THE
More informationFair Integrated Scheduling of Soft Real-time Tardiness Classes on Multiprocessors
Fair Integrated Sheduling of Soft Real-time Tardiness Classes on Multiproessors UmaMaheswari C. Devi and James H. Anderson Department of Computer Siene, The University of North Carolina, Chapel Hill, NC
More informationMOLECULAR ORBITAL THEORY- PART I
5.6 Physial Chemistry Leture #24-25 MOLECULAR ORBITAL THEORY- PART I At this point, we have nearly ompleted our rash-ourse introdution to quantum mehanis and we re finally ready to deal with moleules.
More informationarxiv:gr-qc/ v2 6 Feb 2004
Hubble Red Shift and the Anomalous Aeleration of Pioneer 0 and arxiv:gr-q/0402024v2 6 Feb 2004 Kostadin Trenčevski Faulty of Natural Sienes and Mathematis, P.O.Box 62, 000 Skopje, Maedonia Abstrat It this
More information12.1 Events at the same proper distance from some event
Chapter 1 Uniform Aeleration 1.1 Events at the same proper distane from some event Consider the set of events that are at a fixed proper distane from some event. Loating the origin of spae-time at this
More informationSensor management for PRF selection in the track-before-detect context
Sensor management for PRF seletion in the tra-before-detet ontext Fotios Katsilieris, Yvo Boers, and Hans Driessen Thales Nederland B.V. Haasbergerstraat 49, 7554 PA Hengelo, the Netherlands Email: {Fotios.Katsilieris,
More informationCOMBINED PROBE FOR MACH NUMBER, TEMPERATURE AND INCIDENCE INDICATION
4 TH INTERNATIONAL CONGRESS OF THE AERONAUTICAL SCIENCES COMBINED PROBE FOR MACH NUMBER, TEMPERATURE AND INCIDENCE INDICATION Jiri Nozika*, Josef Adame*, Daniel Hanus** *Department of Fluid Dynamis and
More informationTight bounds for selfish and greedy load balancing
Tight bounds for selfish and greedy load balaning Ioannis Caragiannis Mihele Flammini Christos Kaklamanis Panagiotis Kanellopoulos Lua Mosardelli Deember, 009 Abstrat We study the load balaning problem
More informationMetric of Universe The Causes of Red Shift.
Metri of Universe The Causes of Red Shift. ELKIN IGOR. ielkin@yande.ru Annotation Poinare and Einstein supposed that it is pratially impossible to determine one-way speed of light, that s why speed of
More informationGeneralized Dimensional Analysis
#HUTP-92/A036 7/92 Generalized Dimensional Analysis arxiv:hep-ph/9207278v1 31 Jul 1992 Howard Georgi Lyman Laboratory of Physis Harvard University Cambridge, MA 02138 Abstrat I desribe a version of so-alled
More informationEvaluation of effect of blade internal modes on sensitivity of Advanced LIGO
Evaluation of effet of blade internal modes on sensitivity of Advaned LIGO T0074-00-R Norna A Robertson 5 th Otober 00. Introdution The urrent model used to estimate the isolation ahieved by the quadruple
More informationThe Impact of Information on the Performance of an M/M/1 Queueing System
The Impat of Information on the Performane of an M/M/1 Queueing System by Mojgan Nasiri A thesis presented to the University of Waterloo in fulfillment of the thesis requirement for the degree of Master
More informationRelativity in Classical Physics
Relativity in Classial Physis Main Points Introdution Galilean (Newtonian) Relativity Relativity & Eletromagnetism Mihelson-Morley Experiment Introdution The theory of relativity deals with the study of
More informationLECTURE NOTES FOR , FALL 2004
LECTURE NOTES FOR 18.155, FALL 2004 83 12. Cone support and wavefront set In disussing the singular support of a tempered distibution above, notie that singsupp(u) = only implies that u C (R n ), not as
More informationA NETWORK SIMPLEX ALGORITHM FOR THE MINIMUM COST-BENEFIT NETWORK FLOW PROBLEM
NETWORK SIMPLEX LGORITHM FOR THE MINIMUM COST-BENEFIT NETWORK FLOW PROBLEM Cen Çalışan, Utah Valley University, 800 W. University Parway, Orem, UT 84058, 801-863-6487, en.alisan@uvu.edu BSTRCT The minimum
More informationV. Interacting Particles
V. Interating Partiles V.A The Cumulant Expansion The examples studied in the previous setion involve non-interating partiles. It is preisely the lak of interations that renders these problems exatly solvable.
More informationModal Horn Logics Have Interpolation
Modal Horn Logis Have Interpolation Marus Kraht Department of Linguistis, UCLA PO Box 951543 405 Hilgard Avenue Los Angeles, CA 90095-1543 USA kraht@humnet.ula.de Abstrat We shall show that the polymodal
More informationEXACT TRAVELLING WAVE SOLUTIONS FOR THE GENERALIZED KURAMOTO-SIVASHINSKY EQUATION
Journal of Mathematial Sienes: Advanes and Appliations Volume 3, 05, Pages -3 EXACT TRAVELLING WAVE SOLUTIONS FOR THE GENERALIZED KURAMOTO-SIVASHINSKY EQUATION JIAN YANG, XIAOJUAN LU and SHENGQIANG TANG
More informationChapter 2 Linear Elastic Fracture Mechanics
Chapter 2 Linear Elasti Frature Mehanis 2.1 Introdution Beginning with the fabriation of stone-age axes, instint and experiene about the strength of various materials (as well as appearane, ost, availability
More informationREFINED UPPER BOUNDS FOR THE LINEAR DIOPHANTINE PROBLEM OF FROBENIUS. 1. Introduction
Version of 5/2/2003 To appear in Advanes in Applied Mathematis REFINED UPPER BOUNDS FOR THE LINEAR DIOPHANTINE PROBLEM OF FROBENIUS MATTHIAS BECK AND SHELEMYAHU ZACKS Abstrat We study the Frobenius problem:
More informationMultiPhysics Analysis of Trapped Field in Multi-Layer YBCO Plates
Exerpt from the Proeedings of the COMSOL Conferene 9 Boston MultiPhysis Analysis of Trapped Field in Multi-Layer YBCO Plates Philippe. Masson Advaned Magnet Lab *7 Main Street, Bldg. #4, Palm Bay, Fl-95,
More informationCan Learning Cause Shorter Delays in Reaching Agreements?
Can Learning Cause Shorter Delays in Reahing Agreements? Xi Weng 1 Room 34, Bldg 2, Guanghua Shool of Management, Peking University, Beijing 1871, China, 86-162767267 Abstrat This paper uses a ontinuous-time
More informationMost results in this section are stated without proof.
Leture 8 Level 4 v2 he Expliit formula. Most results in this setion are stated without proof. Reall that we have shown that ζ (s has only one pole, a simple one at s =. It has trivial zeros at the negative
More informationLightpath routing for maximum reliability in optical mesh networks
Vol. 7, No. 5 / May 2008 / JOURNAL OF OPTICAL NETWORKING 449 Lightpath routing for maximum reliability in optial mesh networks Shengli Yuan, 1, * Saket Varma, 2 and Jason P. Jue 2 1 Department of Computer
More informationarxiv:physics/ v4 [physics.gen-ph] 9 Oct 2006
The simplest derivation of the Lorentz transformation J.-M. Lévy Laboratoire de Physique Nuléaire et de Hautes Energies, CNRS - IN2P3 - Universités Paris VI et Paris VII, Paris. Email: jmlevy@in2p3.fr
More informationON A PROCESS DERIVED FROM A FILTERED POISSON PROCESS
ON A PROCESS DERIVED FROM A FILTERED POISSON PROCESS MARIO LEFEBVRE and JEAN-LUC GUILBAULT A ontinuous-time and ontinuous-state stohasti proess, denoted by {Xt), t }, is defined from a proess known as
More informationParallel disrete-event simulation is an attempt to speed-up the simulation proess through the use of multiple proessors. In some sense parallel disret
Exploiting intra-objet dependenies in parallel simulation Franeso Quaglia a;1 Roberto Baldoni a;2 a Dipartimento di Informatia e Sistemistia Universita \La Sapienza" Via Salaria 113, 198 Roma, Italy Abstrat
More informationBoundary value problems for the one-dimensional Willmore equation
Boundary value problems for the one-dimensional Willmore equation Klaus Dekelnik and Hans-Christoph Grunau Fakultät für Mathematik Otto-von-Guerike-Universität Postfah D-396 Magdeburg November 7, 6 Abstrat
More informationTHE METHOD OF SECTIONING WITH APPLICATION TO SIMULATION, by Danie 1 Brent ~~uffman'i
THE METHOD OF SECTIONING '\ WITH APPLICATION TO SIMULATION, I by Danie 1 Brent ~~uffman'i Thesis submitted to the Graduate Faulty of the Virginia Polytehni Institute and State University in partial fulfillment
More informationON THE MOVING BOUNDARY HITTING PROBABILITY FOR THE BROWNIAN MOTION. Dobromir P. Kralchev
Pliska Stud. Math. Bulgar. 8 2007, 83 94 STUDIA MATHEMATICA BULGARICA ON THE MOVING BOUNDARY HITTING PROBABILITY FOR THE BROWNIAN MOTION Dobromir P. Kralhev Consider the probability that the Brownian motion
More informationCase I: 2 users In case of 2 users, the probability of error for user 1 was earlier derived to be 2 A1
MUTLIUSER DETECTION (Letures 9 and 0) 6:33:546 Wireless Communiations Tehnologies Instrutor: Dr. Narayan Mandayam Summary By Shweta Shrivastava (shwetash@winlab.rutgers.edu) bstrat This artile ontinues
More informationPerturbation Analyses for the Cholesky Factorization with Backward Rounding Errors
Perturbation Analyses for the holesky Fatorization with Bakward Rounding Errors Xiao-Wen hang Shool of omputer Siene, MGill University, Montreal, Quebe, anada, H3A A7 Abstrat. This paper gives perturbation
More informationSimple FIR Digital Filters. Simple FIR Digital Filters. Simple Digital Filters. Simple FIR Digital Filters. Simple FIR Digital Filters
Simple Digital Filters Later in the ourse we shall review various methods of designing frequeny-seletive filters satisfying presribed speifiations We now desribe several low-order FIR and IIR digital filters
More informationA simple expression for radial distribution functions of pure fluids and mixtures
A simple expression for radial distribution funtions of pure fluids and mixtures Enrio Matteoli a) Istituto di Chimia Quantistia ed Energetia Moleolare, CNR, Via Risorgimento, 35, 56126 Pisa, Italy G.
More informationCoding for Random Projections and Approximate Near Neighbor Search
Coding for Random Projetions and Approximate Near Neighbor Searh Ping Li Department of Statistis & Biostatistis Department of Computer Siene Rutgers University Pisataay, NJ 8854, USA pingli@stat.rutgers.edu
More informationCollinear Equilibrium Points in the Relativistic R3BP when the Bigger Primary is a Triaxial Rigid Body Nakone Bello 1,a and Aminu Abubakar Hussain 2,b
International Frontier Siene Letters Submitted: 6-- ISSN: 9-8, Vol., pp -6 Aepted: -- doi:.8/www.sipress.om/ifsl.. Online: --8 SiPress Ltd., Switzerland Collinear Equilibrium Points in the Relativisti
More informationA Functional Representation of Fuzzy Preferences
Theoretial Eonomis Letters, 017, 7, 13- http://wwwsirporg/journal/tel ISSN Online: 16-086 ISSN Print: 16-078 A Funtional Representation of Fuzzy Preferenes Susheng Wang Department of Eonomis, Hong Kong
More informationSURFACE WAVES OF NON-RAYLEIGH TYPE
SURFACE WAVES OF NON-RAYLEIGH TYPE by SERGEY V. KUZNETSOV Institute for Problems in Mehanis Prosp. Vernadskogo, 0, Mosow, 75 Russia e-mail: sv@kuznetsov.msk.ru Abstrat. Existene of surfae waves of non-rayleigh
More informationSensitivity Analysis in Markov Networks
Sensitivity Analysis in Markov Networks Hei Chan and Adnan Darwihe Computer Siene Department University of California, Los Angeles Los Angeles, CA 90095 {hei,darwihe}@s.ula.edu Abstrat This paper explores
More informationModeling of Threading Dislocation Density Reduction in Heteroepitaxial Layers
A. E. Romanov et al.: Threading Disloation Density Redution in Layers (II) 33 phys. stat. sol. (b) 99, 33 (997) Subjet lassifiation: 6.72.C; 68.55.Ln; S5.; S5.2; S7.; S7.2 Modeling of Threading Disloation
More informationTail-robust Scheduling via Limited Processor Sharing
Performane Evaluation Performane Evaluation 00 200) 22 Tail-robust Sheduling via Limited Proessor Sharing Jayakrishnan Nair a, Adam Wierman b, Bert Zwart a Department of Eletrial Engineering, California
More informationReliability Guaranteed Energy-Aware Frame-Based Task Set Execution Strategy for Hard Real-Time Systems
Reliability Guaranteed Energy-Aware Frame-Based ask Set Exeution Strategy for Hard Real-ime Systems Zheng Li a, Li Wang a, Shuhui Li a, Shangping Ren a, Gang Quan b a Illinois Institute of ehnology, Chiago,
More informationRESEARCH ON RANDOM FOURIER WAVE-NUMBER SPECTRUM OF FLUCTUATING WIND SPEED
The Seventh Asia-Paifi Conferene on Wind Engineering, November 8-1, 9, Taipei, Taiwan RESEARCH ON RANDOM FORIER WAVE-NMBER SPECTRM OF FLCTATING WIND SPEED Qi Yan 1, Jie Li 1 Ph D. andidate, Department
More informationNEW MEANS OF CYBERNETICS, INFORMATICS, COMPUTER ENGINEERING, AND SYSTEMS ANALYSIS
Cybernetis and Systems Analysis, Vol. 43, No. 5, 007 NEW MEANS OF CYBERNETICS, INFORMATICS, COMPUTER ENGINEERING, AND SYSTEMS ANALYSIS ARCHITECTURAL OPTIMIZATION OF A DIGITAL OPTICAL MULTIPLIER A. V. Anisimov
More informationRelativistic Dynamics
Chapter 7 Relativisti Dynamis 7.1 General Priniples of Dynamis 7.2 Relativisti Ation As stated in Setion A.2, all of dynamis is derived from the priniple of least ation. Thus it is our hore to find a suitable
More informationarxiv:gr-qc/ v7 14 Dec 2003
Propagation of light in non-inertial referene frames Vesselin Petkov Siene College, Conordia University 1455 De Maisonneuve Boulevard West Montreal, Quebe, Canada H3G 1M8 vpetkov@alor.onordia.a arxiv:gr-q/9909081v7
More informationFINITE WORD LENGTH EFFECTS IN DSP
FINITE WORD LENGTH EFFECTS IN DSP PREPARED BY GUIDED BY Snehal Gor Dr. Srianth T. ABSTRACT We now that omputers store numbers not with infinite preision but rather in some approximation that an be paed
More informationExtending LMR for anisotropic unconventional reservoirs
Extending LMR for anisotropi unonventional reservoirs Maro A. Perez Apahe Canada Ltd Summary It has beome inreasingly advantageous to haraterize rok in unonventional reservoirs within an anisotropi framework.
More information