Green s function Monte Carlo algorithms for elliptic problems
|
|
- Suzan Miles
- 5 years ago
- Views:
Transcription
1 Mathematics and Computes in Simulation 63 (2003) Geen s function Monte Calo algoithms fo elliptic poblems I.T. Dimov, R.Y. Papancheva Cental Laboatoy fo Paallel Pocessing, Depatment of Paallel Algoithms, Bulgaian Academy of Sciences, Acad. G. Bonchev St., bl. 25A, 1113 Sofia, Bulgaia Received 2 Septembe 2002; eceived in evised fom 2 July 2003; accepted 28 July 2003 Abstact In many lage-scale poblems, one is inteested to obtain diectly an appoximate value of a functional of the solution. Hee, we conside a special class of gid-fee Monte Calo algoithms fo diect computing of linea functionals of the solution of an elliptic bounday-value poblem. Such kind of poblems appea in envionmental sciences, computational physics and financial mathematics. To ceate the algoithms, we use the Geen s function analysis and define the conditions unde which the integal tansfomation kenel of the integal epesentation fo the bounday-value poblem unde consideation is non-negative. This analysis is done fo a possible set of densities, and it is used to geneate thee diffeent gid-fee Monte Calo algoithms based on diffeent choices of the density of the adius of the balls used in Monte Calo simulation. Only one of the geneated algoithms was known befoe. We shall call it Sipin s algoithm. It was poposed and studied by Sipin. The aim of this wok is to study the two new algoithms poposed hee and based on two othe (than in Sipin s algoithm) possible choices of the densities. The algoithms ae descibed and analyzed. The pefomed numeical tests show that the efficiency of one of the new algoithms, which is based on a constant density is highe than the efficiency of Sipin s algoithm IMACS. Published by Elsevie B.V. All ights eseved. Keywods: Monte Calo method; Selection methods; Density function 1. Intoduction In this pape, new Monte Calo algoithms fo computing functionals of solutions of elliptic boundayvalue poblems ae poposed and studied. Thee ae many lage-scale poblems in computational physics, biology and envionmental sciences, whee it is impotant to be able to compute functionals of the solution without computing the solution itself. It is that because some impotant paametes which can be measued in physics, and some effects on the live matte ae in fact functionals. Fo instance, the damage effect of Coesponding autho. Tel.: ; fax: addesses: ivdimov@bas.bg (I.T. Dimov), umi@canto.bas.bg (R.Y. Papancheva). URL: /$ IMACS. Published by Elsevie B.V. All ights eseved. doi: /s (03)
2 588 I.T. Dimov, R.Y. Papancheva / Mathematics and Computes in Simulation 63 (2003) dange pollution is measued by: f(x)u(x) dx, Ω whee x (x 1,...,x d ) Ω R d, u(x) is the concentation of the pollutant and f(x) is given sensitivity function. Anothe example could be given fom the statistical physics, whee paametes chaacteized a multi-paticle system like the mean-velocity V m and the enegy E ae defined by: V m = f(x)v(x) dx, E = Ω Ω f(x)v 2 (x) dx, whee f(x)(x Ω R d ) is the distibution function. The algoithms poposed in this pape can be efficiently implemented in calculating functionals of the solution. Conside the functional J(u): J(u) (g, u) = g(x)u(x) dx, (1) Ω whee Ω R 3 and x = (x 1,x 2,x 3 ) Ω is a point in the Euclidean space R 3. The functions u(x) and g(x) belong to the Banach space X and to the adjoint space X, espectively, and u(x) is the unique solution of the following Fedholm integal equation in an opeato fom: u = K(u) + f. (2) The algoithms unde consideation use a local integal epesentation of the solution by Geen s function. The Monte Calo appoach with the use of Geen s function is well known (see [8]). The basic poblem is the estimation of the computational complexity and the constuction of Monte Calo algoithms with lowe computational cost. We also conside poblems, in which a linea functional of the solution is seached fo and a ough solution is acceptable. Fo pactical computations, it means that the elative eo is about 5 10%. Fo example, such kind of poblems ae the pollution poblems in envionmental sciences. It is known that the most efficient computational method fo diect computing of functionals is the Monte Calo method (see [2,6]). This pape deals with an analysis of the choice of the density function p() in Lévy s function L p (y, x). We use the Geen s function analysis to ceate the algoithms and define the conditions unde which the integal tansfomation kenel of the integal epesentation fo the consideed bounday-value poblem is non-negative. This analysis is done fo a possible set of densities. The esult is used to geneate diffeent gid-fee Monte Calo algoithms based on diffeent choices of the density of the adius of the balls used in Monte Calo simulation. Only one of the geneated algoithms was known befoe. We call it Sipin s algoithms. It was poposed by Sipin and studied in [9]. Since in the Monte Calo simulation a selection technique is used, it is impotant to obtain an estimation fo the efficiency of the algoithm. The selection algoithm in geneal case woks as follows. Let us assume that we have to calculate a andom vaiable ξ with density function f 1 (x)/f 1 and an easily computed function f 2 (x) exists, with
3 Ω I.T. Dimov, R.Y. Papancheva / Mathematics and Computes in Simulation 63 (2003) f 1 (x) f 2 (x) fo all x Ω. Hee, F 1 = f 1 (x) dx and F 2 = f 2 (x) dx. Ω Let ξ be a andom vaiable with density function f 2 (x)/f 2, and γ be an unifomly distibuted andom vaiable in the inteval [0, 1]. It has been shown [8] that if: γf 2 (ξ) f 1 (ξ), (3) then ξ is the needed andom vaiable with density function f 1 (x)/f 1. This algoithm allows us to use the function f 2 (x) instead of the oiginal function f 1 (x). Howeve, some vaiables will be ejected, namely when condition (3) will not be fulfilled. We use the following definition fo the efficiency. Definition 1.1. The efficiency of the selection algoithm is defined as the atio: E = F 1 F 2. In 1984, Sipin poved that the efficiency of the selection gid-fee Monte Calo algoithm is geate than 1/2 [9]. An impoved esult fo the efficiency of the Sipin s algoithm was obtained by Dimov in 1988 [3]: E = 1 + α 2 + α, fo some explicit constant α. The last esult was impoved again by Dimov and Guov in 1998 [4]: E = 1 + α 2 + α ε R, fo some ε R, with 0 <ε R < 1. It seems that this impovement is exact, i.e. it could not be impoved again (but the last statement is still not poved). In all these woks [3,4,9], the density function is fixed p() = e k. Hee, we should answe the question if a diffeent function to be taken, how this affects the algoithm efficiency. The aim of this wok is to study the two new algoithms poposed hee and based on two othe (than in Sipin s algoithm) possible choices of the densities. The algoithms ae descibed and analyzed. The pefomed numeical tests show that the efficiency of one of the new algoithms which is based on a constant density is highe than the efficiency of Sipin s algoithm. The pape is oganized as follows. In Section 2, we descibe the fomulation of the poblem and give the conditions fo existing an unique solution u(x) C 2 (Ω) C( Ω).InSection 3, the set of possible densities poviding the conditions, unde which the integal tansfomation kenel is non-negative is defined. Two impotant cases ae consideed: a case of a constant density and a case of an exponential density. We find a set of paametes of the poblem fo which both constant and exponential densities can be applied. In Section 4, we descibe and analyze the algoithms geneated by the chosen densities. Numeical tests ae pesented in Section 5. The esults of the numeical simulation illustate the efficiency of the poposed algoithms. In Section 6, we give some concluding emaks.
4 590 I.T. Dimov, R.Y. Papancheva / Mathematics and Computes in Simulation 63 (2003) Fomulation of the poblem Let X = L 1 (Ω). Then, X = L (Ω) ([8], p. 148). Conside the case when K is an odinay linea integal tansfom: K(u)(x) = k(x, y)u(y) dy. Ω In this case, (2) obtains the fom: u(x) = k(x, y)u(y) dy + f(x). Ω The coesponding Neumann seies conveges when the condition: K(u) L(L1 ) = sup k(x, y) dx q<1, y Ω Ω holds. We have to calculate the functional J(u), whee u(x) is the solution of the following bounday-value poblem: M(u)(x) = Φ(x), x Ω, Ω R 3, (4) u(x) = ψ(x), x Ω, whee the opeato M is defined by: 3 ( 2 M = + b i (x) ) + c(x). x i i=1 x 2 i Hee follows the definition fo the class A (k,λ) of domains Ω. Definition 2.1. The domain Ω belongs to the class A (k,λ) if fo any point x Ω the bounday Ω can be epesented as a function z 3 = σ(z 1,z 2 ) in the neighbohood of x with k σ C (0,λ) (Ω), i.e. k σ(y) k σ(y ) Q y y λ, (5) whee the vectos y (z 1,z 2 ) and y λ (0, 1]. (z 1,z 2 ) ae two-dimensional vectos, Q is a constant, and The definition of the class C (0,λ) (Ω) can be found in [12]. If in the closed domain Ω A (1,λ) the coefficients of the opeato M satisfy the conditions: b i (i = 1, 2, 3), c C (0,λ) ( Ω), c 0 and Φ C (0,λ) (Ω) C( Ω), ψ C( Ω), the poblem (4) and (5) has an unique solution u(x) C 2 (Ω) C( Ω). The conditions fo uniqueness of the solution can be found in [1,11]. The fist step in studying the Monte Calo appoach is obtaining an integal epesentation of the solution in the fom: u(x) = k(x, y)u(y) dy + f(x). B(x)
5 I.T. Dimov, R.Y. Papancheva / Mathematics and Computes in Simulation 63 (2003) Fom an algoithmic point of view, the domain B(x) must be chosen in such a way that the coodinates of the bounday points y B(x) could be easily calculated. So we denote by B(x) the ball: B(x) = B R (x) ={y : = x y R(x)}, whee R(x) is the adius of the ball. We seek a epesentation of the integal kenel k(x, y), using Lévy s function and the adjoint opeato M fo the initial diffeential opeato M. Hee follows a shot definition fo the Lévy s function [12]. Definition 2.2. Let, H(x, y) = 1 4π x y. Evey function L(x, y) continuous in vaiables x and y, fo x and y in Ω, and x y, togethe with its fist and second deivatives with espect to the x i is called a Lévy s function if it satisfies bounds of the following type fo some λ>0 and = x y : L H = O( λ 1 ), [L H] x i = O( λ 2 ), 2 [L H] x i x j = O( λ 3 ), such bounds holding unifomly in evey bounded domain in Ω. The explicit fom of the Lévy s function fo the poblem (4) and (5) can be found in [12]: R ( 1 L p (y, x) = µ p (R) 1 ) p(ρ) dρ, R, (6) ρ whee the following notation is used: p() is a density function: µ p (R) = [4πq p (R)] 1, q p (R) = R 0 p(ρ) dρ. It is known [4] that if the vecto-function b(x) satisfies the conditions b j (x) C(Ω) (j = 1, 2, 3), and div b(x) = 0, then the adjoint opeato M, applied on functions v(x) C 2 (Ω), such that: v(x) x i = v(x) = 0, fo any x Ω, i = 1, 2, 3, has the following fom: M = 3 ( 2 x 2 i=1 i b i (x) ) + c(x). x i
6 592 I.T. Dimov, R.Y. Papancheva / Mathematics and Computes in Simulation 63 (2003) Accoding the esults of Mianda [12], when the Lévy s function is used, the following integal epesentation fo the solution u(x) of the bounday-value poblem is tue: u(x) = (u(y)my L p(y, x) + L p (y, x)φ(y)) dy T(x) + 3 T(x) i=1 [( Lp (y, x) u(y) n i u(y) L ) ] p(y, x) b i (y)u(y)l p (y, x) dy, (7) y i y i whee n = (n 1,n 2,n 3 ) is the exteio nomal to the bounday T(x). Fomula (7) holds fo any domain T(x) A (1,λ) and theefoe it is emaining tue fo evey ball B(x), lying inside the domain Ω. It is shown that it is possible to constuct the function L p (y, x) in such a way that My L p(y, x) is non-negative in B(x), and such that L p (y, x) and its deivatives vanish on B(x). Finally, the integal epesentation of the solution has the following fom: u(x) = My L p(y, x)u(y) dy + L p (y, x)φ(y) dy, (8) whee B(x) M y L p(y, x) = µ p (R) p() + µ p(r) 2 2 B(x) R p(ρ) µ p (R)c(y) ρ [ 3 c(y) + b i (y) y i x i i=1 dρ ] R p(ρ) dρ. Fo solving ou poblem, we use a Monte Calo pocedue that is called Walk on the Balls (WOB). To ensue the convegence of the pocess, we intoduce the ε-stip of the bounday. This pocess is simila to the well-known spheical pocess [7,10,13,14]. The pocess stats at point ξ 0 = x Ω, which is chosen coespondingly with some initial density function π(x). In the paticula case, when π(x) = δ(x x 0 ) the Makov s chain stats at the point x 0 Ω. The next andom point is detemined with tansition density function p(x, y), which is popotional to the kenel of the integal equation (2) inside maximal ball B(x) Ω. This pocess teminates when the point falls into the ε-stip of the bounday. In spheical coodinates, p(x, y) is equal to: p(, w) = p() p (w). We aim at analyzing the oppotunities fo a choice of density function p(). Ou goal is to conside such density functions, in which the choice of bigge adii of the next ball in the Makov chain is moe likely. One may expect that in this case the length of the chains will be shote and the pocess moe effective. 3. Geen s function analysis Hee, we intoduce the following notations: b = sup b(x), x Ω c = sup c(x), x Ω and denote by R the adius of the maximal ball lying inside Ω. b(x) = (b 1 (x), b 2 (x), b 3 (x)),
7 and I.T. Dimov, R.Y. Papancheva / Mathematics and Computes in Simulation 63 (2003) In [5] was poved that the conditions: My L p(y, x) 0, fo any y B(x), (9) L p (y, x) = L p(y, x) = 0, y i fo any y B(x), i = 1, 2, 3, (10) ae satisfied fo: p() = e k, (11) whee k b + Rc. The density function defined in (11) coesponds to the Sipin s algoithm. Ou aim is to conside the oppotunities and the conditions fo a diffeent choice of the density function p(). The following poblem can be fomulated: find the class of possible densities fo which the integal tansfomation kenel M y L p(y, x) is non-negative and the conditions (10) ae satisfied. The condition: L p (y, x) = 0, fo any y B(x), is obvious and follows diectly fom (6). Since, L p (y, x) y i we obtain: = L p = x i y i y i 3 µ p L p (y, x) = 0, fo any y B(x). y i We shall pove the following theoem. R p(ρ) dρ, Theoem 3.1. If the density function p() satisfies the inequality: p() (b + R4 ) R c p(ρ) dρ, (12) then the integal tansfomation kenel M y L p(y, x) is non-negative fo any y B(x). Poof. My L p(y, x) can be witten as follows [4]: My L p(y, x) = µ p(r) Γ 2 p (y, x), whee ( R R ) p(ρ) 3 Γ p (y, x) = p() + c(y) p(ρ) dρ dρ + b i (y) y i x i ρ i=1 R p(ρ) dρ.
8 594 I.T. Dimov, R.Y. Papancheva / Mathematics and Computes in Simulation 63 (2003) Since ρ 0, it follows that 1 /ρ 0, and theefoe, R R p(ρ) p(ρ) dρ dρ 0. (13) ρ Fo any y: 0 c(y) c. As x y =, we obtain that: 3 i=1 b i (y) y i x i b. Taking in account (13) (15), we obtain: ( R Γ p (y, x) p() c p(ρ) dρ Since 0 ρ R, it follows: 1 ρ 1 R, R p(ρ) ρ and theefoe, ( R R ) p(ρ) ( p(ρ) dρ dρ 1 ) R ρ R Fo 0 R, one has: ( 1 ) R R 4. (14) (15) ) R dρ b p(ρ) dρ. (16) p(ρ) dρ. (17) Finally, using the obtained extemum and the inequalities (16) and (17) we can wite: Γ p (y, x) p() c R R R ( p(ρ) dρ b p(ρ) dρ = p() b + c R ) R p(ρ) dρ. 4 4 Hence the esult. The inequality (12) defines a class of possible densities p(), which can be used in a Monte Calo simulation. We shall conside two possible densities: p() = p = const fo a given set of paametes of the poblem (this case is descibed in Section 3.1) and p() = e m,m>0 (descibed in Section 3.2) Case A: a constant density In this subsection, we descibe the conditions which have to be fulfilled in ode to use a constant density p() = p = const.
9 I.T. Dimov, R.Y. Papancheva / Mathematics and Computes in Simulation 63 (2003) Fo this choice of the function p() is tue that: p() (b + R4 ) R [ c p(ρ) dρ p 1 (b + R4 ) ] c R. Consequently, if (b + (R/4)c )R 1, then M y L p(y, x) 0. Lemma 3.1. If the paametes of the poblem satisfy the inequality: 1 (b Rc )R, and the density function p() is chosen to be: p() = p = const, then the condition (9) is tue. Remak. If the adius of the maximal ball, lying inside the domain Ω satisfies the inequality: R 2 c ( b 2 + c b ), then the condition (9) is tue Case B: an exponential density While in Section 3.1, the adius is unifomly distibuted in [0,R], hee, since the function p() = e m,m>0 is stictly inceasing, the choice of bigge adii will be moe pobable. Let us intoduce a paamete k, defined by the expession: k = b Rc, and conside the density function p() = e k. Fom the inequality (12) one can obtain that: k(r ) ln 2. Consequently if the condition kr ln 2 is fulfilled, then the inequality (12) is tue fo any. Lemma 3.2. If the paametes of the poblem satisfy the inequality: kr ln 2, and the density function p() is chosen to be: p() = e k, k = b Rc, then the condition (9) is fulfilled.
10 596 I.T. Dimov, R.Y. Papancheva / Mathematics and Computes in Simulation 63 (2003) Monte Calo algoithms Conside the poblem of evaluating the functional: J(u) (g, u) = g(x)u(x) dx, Ω whee u(x) is the solution of the integal equation (8). The Monte Calo algoithm fo solving this poblem can be defined as a ball pocess. To ensue the convegence of the pocess, we intoduce the ε-stip of the bounday, i.e. whee Ω ε ={x Ω : B(x) = B ε (x)}, B ε (x) ={y Ω : x y ε}. The andom vaiable, whose mathematical expectation coincides with J(u) is: Θ[g] = g(ξ 0) Q j f(ξ i ), π(ξ 0 ) whee j=0 k(ξ j 1,ξ j ) Q 0 = 1, Q j = Q j 1, j = 1, 2, 3,..., p(ξ j 1,ξ j ) and ξ 0,ξ 1,... is a Makov chain in Ω with initial density function π(x) and tansition densities p(x, y), which ae toleant to g(x) and k(x, y), espectively [2,8,15]. We use the following definition of the toleant density function. Definition 4.1. It is said that the density p(x) is toleant to the function f(x) if p(x) > 0 in all points x Ω, with f(x) 0. In [4] was poved that the integal tansfomation kenel in local integal epesentation can be used as a tansition density function in the Makov pocess. In ou case: { M y L p (y, x), when x Ω\ Ω, p(x, y) = k(x, y) = 0, when x Ω, and { B(x) f(x) = L p(y, x)φ(y) dy, ψ(y), when x Ω\ Ω, when x Ω, The function p(x, y) can be expessed in spheical coodinates as [4]: p(, w) = sin θ 4π p() q p (R) p (w),
11 I.T. Dimov, R.Y. Papancheva / Mathematics and Computes in Simulation 63 (2003) whee [ ] b(x + w) cos(b, w) + c(x + w) R c(x + R w)2 p(ρ) p (w) = 1 + p(ρ) dρ p() p() ρ dρ. Fo simulating andom vaiable with density function p (w) we use the selection algoithm [4]. Since, b R p (w) 1 + p(ρ) dρ = h(), p() the function h() can be used as a majoant. Hee follows the algoithm fo one andom walk [4]: 1. Calculate the adius R(x) of the maximal ball of cente x lying inside Ω. 2. Compute a ealization of the andom vaiable τ with the density: p() q p (R). 3. Calculate the function: b R h() = 1 + p(ρ) dρ. p() 4. Compute independent ealizations w j of an unit isotopic vecto in R Constuct independent ealizations γ j of an unifomly distibuted andom vaiable in the inteval [0, 1]. 6. Repeat steps 4 and 5 until find the paamete j 0 given by: j 0 = min{j : h()γ j p (w j )}. The andom vecto w j0 has the density p (w). 7. Calculate the andom point y, using the following fomula: y = x + w j0. The value = y x is the adius of the ball lying inside Ω and having a cente at x. 8. Stop the algoithm when the andom pocess eaches the ε-stip Ω ε, i.e. y Ω ε.ify/ Ω ε then the algoithm has to be epeated fo x = y Case A: p() = p = const Hee, p() q p (R) = 1 R, and consequently fo sampling the jump in WOB pocess one can use the fomula: = Rγ, whee γ is an unifomly distibuted andom vaiable in the inteval [0, 1].
12 598 I.T. Dimov, R.Y. Papancheva / Mathematics and Computes in Simulation 63 (2003) Fo majoant function h(),wehave: h() = 1 + b (R ), and fo p (w) we obtain: [ 3 ] ( ) R p (w) = 1 + b i (x + w)w i + c(x + w) (R ) c(x + w) 2 ln. i= Case B: p() = e k,k = b + (R/4)c In this case: p() q p (R) = and consequently, kek e kr 1, = 1 k ln(1 + (ekr 1)γ), whee γ is an unifomly distibuted andom vaiable in the inteval [0, 1]. Fo majoant function h(),wehave: h() = 1 + b k (ek(r ) 1) and fo p (w) we obtain: [ 3 p (w) = 1 + b i (x + w)w i + c(x + w) i=1 ] 1 c(x + R w)2 e kρ k (ek(r ) 1) e k ρ dρ. 5. Numeical tests Numeical tests ae pefomed in ode to compae the poposed algoithms on a poblem, which is simila to some eal-life poblems. The examples we conside hee deal with the following bounday-value poblem: 3 ( 2 u(x) i=1 x 2 i + b i (x) u(x) ) + c(x)u(x) = 0 in Ω = [0, 1] 3, x i u(x 1,x 2,x 3 ) = e a 1x 1 +a 2 x 2 +a 3 x 3, (x 1,x 2,x 3 ) Ω. In ou tests: b 1 (x) = a 2 a 3 (x 2 x 3 ), b 2 (x) = a 3 a 1 (x 3 x 1 ), b 3 (x) = a 1 a 2 (x 1 x 2 ), (18)
13 I.T. Dimov, R.Y. Papancheva / Mathematics and Computes in Simulation 63 (2003) Table 1 CPU time in seconds fo the case when ε = 0.1 and a i = 0.25 (i = 1, 2, 3) fo diffeent values of the numbes of the Makov chains N N p() = const p() = e k p() = e k and c(x) = (a a2 2 + a2 3 ), whee a 1,a 2,a 3 ae paametes of the poblem. We conside thee cases fo the coefficients: the fist case, when a i = 0.25 (i = 1, 2, 3); the second case, when a i = 0.5 (i = 1, 2, 3); the thid case, when a i = 1 (i = 1, 2, 3). Such kind of poblems appea in envionmental mathematics and descibe the pollution tanspot due to advection and diffusion and take into account the deposition of pollution in aeas fee of emission souces (fo moe details, see [16]). We also use diffeent values to define the ε-stips: ε = 0.01; ε = 0.05; ε = 0.1. The numeical expeiments estimate the solution at the point: x 0 = (0.5, 0.5, 0.5), i.e. g(x) = δ(x x 0 ). Some numeical esults ae pesented in Table 1. We compae the computational time in cases when p() = const (algoithm A), p() = e k (algoithm B) and p() = e k. The last case is used as a basic algoithm. One can see, that the algoithm A is much faste than the basic algoithm. It can be also seen, that the algoithm B is even a little bit faste than the basic one. The eason fo such a esult could be that the numbe of opeations of the poposed algoithms is less, o the aveage length of the Makov chains is less. Table 2 shows the aveage length of the Makov chain fo all thee algoithms. One can see, that the length of the Makov chains is pactically the same (fo the both of poposed algoithms the length of the Makov chains is a little bit less). It means that the eason fo deceasing the computational time fo the Table 2 Aveage lengths of the Makov chains in case when ε = 0.1 and a i = 0.25 (i = 1, 2, 3) fo diffeent values of N N p() = const p() = e k p() = e k
14 600 I.T. Dimov, R.Y. Papancheva / Mathematics and Computes in Simulation 63 (2003) Table 3 Aveage lengths of the Makov chains in case when N = 200,000 ε a i (i = 1, 2, 3) p() = const p() = e k p() = e k Table 4 Maximal lengths of the Makov chains in case when N = 200,000 ε a i (i = 1, 2, 3) p() = const p() = e k p() = e k algoithm A is mostly the fact that this algoithm is less time consuming. Ou numeical tests show that the length of the Makov chains (the numbe of Monte Calo iteations) does not depend on the numbe N of the Makov chains. Such a esult can be expected, because the numbe of the Monte Calo iteations depends on the adius of convegence of the iteative opeato and the value of the ε-stip. As smalle is the adius of convegence of the iteative opeato, smalle is the numbe of iteations (that is the length of the Makov chain fo the algoithms unde consideation). As smalle is the ε-stip of the domain, lage p()=const p()=e -k p()=e k CPU N Fig. 1. Dependence of the CPU time fom the numbe of andom tajectoies in case when ε = 0.1, a = 0.25.
15 I.T. Dimov, R.Y. Papancheva / Mathematics and Computes in Simulation 63 (2003) is the numbe of iteation (the aveage length of the Makov chain). Some of ou numeical tests pesented in Table 3 illustate these statements. The fist line pesents the esults of the aveage length of the Makov chain (which coesponds to the aveage numbe of the Monte Calo iteations) in the case when ε = 0.1 and a i = 0.25 (i = 1, 2, 3). The length of the chain is aveaged on a lage numbe of Makov chains (N = 200,000). Tests ae pefomed fo diffeent initial value of the pseudo-andom numbe geneato. The aveage length of the Makov chains pactically does not depend on the initial value of the pseudo-andom numbe geneato. Table 3 contains the same paamete obtained fo the smalle ε-stip (ε = 0.01). In this case, one can obseve a significant inceasing of the aveage numbe of Monte Calo iteations. Some numeical tests wee pefomed in ode to study how the aveage length of the Makov chain depends of the paametes of the poblem. In ou tests we vay the value of the coefficients a i (i = 1, 2, 3). It coesponds to vaiation of the advection pat of the poblem, because these coefficients define the advection coefficients b i (i = 1, 2, 3) (see (18)). The numeical esults show that such kind of changes in the opeato (espectively, in the Monte Calo iteative opeato) does not change a lot the numbe of the Monte Calo iteations. Fom the esults pesented in Table 3, one can see that the aveage length of the Makov chain does not change significantly when the values of a i ae changed fom 0.25 to 0.5. Fo these algoithms it is also impotant to know the maximal length of the Makov chain. Table 4 pesents the esults of length of Makov chains with the maximum states, which coesponds to the maximal numbe of Monte Calo iteations. The Makov chain with the maximal length is chosen among a lage numbe of chains (N = 200,000). The esults ae obtained fo ε = 0.1 and a i = 0.25 (i = 1, 2, 3) Eo p()=const p()=e^{k} p()=e^{-k} N Fig. 2. Convegence of the algoithms fo small values of N, ε = 0.1, a = 0.25, u(x 0 ) =
16 602 I.T. Dimov, R.Y. Papancheva / Mathematics and Computes in Simulation 63 (2003) Simila esults ae obtained fo ε = One can see, that thee is significant dependence of the maximal length of the Makov chain fom the value of the ε-stip of the domain. The dependence of the CPU time fom the numbe of andom tajectoies (numbe of the Makov chains) is shown in Fig. 1. One can see, that the dependence is linea. It is because the computational wok inceases linealy with the inceasing of the numbe of andom tajectoies N. The diffeence between the exact and the appoximate solutions is pesented in Figs. 2 and 3 as a paamete, called Eo. The convegence of all thee algoithms is pactically the same. One can expect such a esult, because the convegence depends on the adius of convegence of the iteative opeato of the Monte Calo method, which is the same fo all thee algoithms. Some esults of the convegence of the algoithms unde consideation ae pesented in Fig. 3. One can see, that fo values of N lage than 1000 all algoithms convege. The speed of convegence is pactically the same fo all thee algoithms. As the exact solution u(x) is geate than 1 fom Figs. 2 and 3, one can also see that the elative eo is always less than 5%, which satisfies the equiements fo the needed accuacy of the algoithms fo solving this kind of poblems. Unde elative eo, we assume the atio among absolute eo and the exact solution, i.e. Eo = u(x 0) ũ(x 0 ), u(x 0 ) whee ũ(x 0 ) is the appoximate value of the solution u(x) at the point x 0 Ω p()=const p()=e^{k} p()=e^{-k} 0.01 Eo N Fig. 3. Vaiation of the eo fo lage values of N, ε = 0.1, a = 0.25, u(x 0 ) =
17 6. Concluding emaks I.T. Dimov, R.Y. Papancheva / Mathematics and Computes in Simulation 63 (2003) We study two new gid-fee Monte Calo algoithms poposed in this wok. The algoithms unde consideation ae based on two diffeent choices of the density of the adius of the balls used in Monte Calo simulations. The algoithms ae descibed and analyzed. We obtain integal epesentations fo constant and exponential densities and pove that the integal tansfomation kenels of the coesponding epesentations can be used as density functions of the gid-fee Monte Calo algoithms. The pefomed numeical tests show that: The numbe of the Monte Calo iteations (the length of the Makov chain) depends mainly of the value of the ε-stip of the domain. The aveage length of the Makov chains pactically does not depend of the numbe of chains. Fo the both of poposed algoithms this vaiation is smalle than the vaiation fo the Sipin s algoithm consideed as a basic algoithm. The algoithms unde consideation convege; the speed of convegence is pactically the same as the speed of convegence of the Sipin s algoithm. One can also see, that the elative eo is always less than 5%, which satisfies the equiements fo the needed accuacy of the algoithms fo solving this kind of poblems. The efficiency of one of the new algoithms which is based on a constant density is highe than the efficiency of Sipin s algoithm. The main eason fo this is the fact that the algoithm based on the constant density has the lowest computational complexity. One can conclude that algoithm A is pefeable, because of its high algoithmic efficiency and the fact that the speed of convegence of this algoithm is pactically the same as the speed of convegence of algoithm B o the Sipin s algoithm. Acknowledgements Suppoted by BIS Cente of Excellence (Contact No. ICA1-CT ) and by the Ministy of Education and Science of Bulgaia unde Gant #MM 902/99. Refeences [1] A.V. Bitzadze, Equations of the Mathematical Physics, Nauka, Moscow, [2] J.H. Cutiss, Monte Calo methods fo the iteation of the linea opeatos, J. Math. Phys. 32 (4) (1954) [3] I.T. Dimov, Efficiency estimatos the Monte Calo algoithms, in: Bl. Sendov, R. Lazaov, I.T. Dimov (Eds.), Poceedings of the Intenational Confeence on Numeical Methods and Applications, Publishing House of the Bulgaian Academy of Sciences, Sofia, 1989, pp [4] I.T. Dimov, T. Guov, Estimates of the computational complexity of iteative Monte Calo algoithm based on Geen s function appoach, Math. Comput. Simul. 47 (1998) [5] I.T. Dimov, A. Kaaivanova, H. Kuchen, H. Stoltze, Monte Calo algoithms fo elliptic diffeential equations a data paallel functional appoach, J. Paallel Algoithms Appl. 9 (1996) [6] I.T. Dimov, O. Tonev, Monte Calo algoithms: pefomance analysis fo some compute achitectues, J. Comput. Appl. Math. 48 (1993)
18 604 I.T. Dimov, R.Y. Papancheva / Mathematics and Computes in Simulation 63 (2003) [7] B.S. Elepov, A.A. Konbeg, G.A. Mikhailov, K.K. Sabelfeld, Solving the Bounday Value Poblems by the Monte Calo Method, Nauka, Novosibisk, [8] S.M. Emakov, G.A. Mikhailov, Statistical Simulation, Nauka, Moscow, [9] S.M. Emakov, V.V. Nekutkin, A.S. Sipin, Random Pocesses fo Solving Classical Equations of the Mathematical Physics, Nauka, Moscow, [10] G.A. Mikhailov, New Monte Calo Methods with Estimating Deivatives, VSP BV, Utecht, The Nethelands, [11] V.P. Mikhailov, Patial Diffeential Equations, Nauka, Moscow, [12] C. Mianda, Patial Diffeential Equations of Elliptic Type, Spinge-Velag, Belin, [13] M. Mulle, Some continuous Monte Calo methods fo the Diichlet poblem, Ann. Math. Stat. 27 (3) (1956) [14] K.K. Sabelfeld, Monte Calo Methods in Bounday Value Poblems, Spinge-Velag, Belin, [15] I.M. Sobol, Monte Calo Numeical Methods, Nauka, Moscow, [16] Z. Zlatev, I.T. Dimov, K. Geogiev, Thee-dimensional vesion of the Danish Euleian model, Z. Angew. Math. Mech. 76 (S4) (1996)
MONTE CARLO STUDY OF PARTICLE TRANSPORT PROBLEM IN AIR POLLUTION. R. J. Papancheva, T. V. Gurov, I. T. Dimov
Pliska Stud. Math. Bulga. 14 (23), 17 116 STUDIA MATHEMATICA BULGARICA MOTE CARLO STUDY OF PARTICLE TRASPORT PROBLEM I AIR POLLUTIO R. J. Papancheva, T. V. Guov, I. T. Dimov Abstact. The actual tanspot
More informationOn 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 informationPearson 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 informationAbsorption 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 informationTemporal-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 informationON 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 information15 Solving the Laplace equation by Fourier method
5 Solving the Laplace equation by Fouie method I aleady intoduced two o thee dimensional heat equation, when I deived it, ecall that it taes the fom u t = α 2 u + F, (5.) whee u: [0, ) D R, D R is the
More informationCompactly 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 informationNumerical 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 informationKOEBE DOMAINS FOR THE CLASSES OF FUNCTIONS WITH RANGES INCLUDED IN GIVEN SETS
Jounal of Applied Analysis Vol. 14, No. 1 2008), pp. 43 52 KOEBE DOMAINS FOR THE CLASSES OF FUNCTIONS WITH RANGES INCLUDED IN GIVEN SETS L. KOCZAN and P. ZAPRAWA Received Mach 12, 2007 and, in evised fom,
More informationMultiple 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 informationThe 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 informationSurveillance 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 informationNew 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 information3.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 informationPhysics 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 informationOn 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 informationCentral 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 informationMONTE CARLO SIMULATION OF FLUID FLOW
MONTE CARLO SIMULATION OF FLUID FLOW M. Ragheb 3/7/3 INTRODUCTION We conside the situation of Fee Molecula Collisionless and Reflective Flow. Collisionless flows occu in the field of aefied gas dynamics.
More informationMeasure Estimates of Nodal Sets of Polyharmonic Functions
Chin. Ann. Math. Se. B 39(5), 08, 97 93 DOI: 0.007/s40-08-004-6 Chinese Annals of Mathematics, Seies B c The Editoial Office of CAM and Spinge-Velag Belin Heidelbeg 08 Measue Estimates of Nodal Sets of
More informationGoodness-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 informationMath 124B February 02, 2012
Math 24B Febuay 02, 202 Vikto Gigoyan 8 Laplace s equation: popeties We have aleady encounteed Laplace s equation in the context of stationay heat conduction and wave phenomena. Recall that in two spatial
More informationResearch Article On Alzer and Qiu s Conjecture for Complete Elliptic Integral and Inverse Hyperbolic Tangent Function
Abstact and Applied Analysis Volume 011, Aticle ID 697547, 7 pages doi:10.1155/011/697547 Reseach Aticle On Alze and Qiu s Conjectue fo Complete Elliptic Integal and Invese Hypebolic Tangent Function Yu-Ming
More information1D2G - Numerical solution of the neutron diffusion equation
DG - Numeical solution of the neuton diffusion equation Y. Danon Daft: /6/09 Oveview A simple numeical solution of the neuton diffusion equation in one dimension and two enegy goups was implemented. Both
More informationAnalytical Solutions for Confined Aquifers with non constant Pumping using Computer Algebra
Poceedings of the 006 IASME/SEAS Int. Conf. on ate Resouces, Hydaulics & Hydology, Chalkida, Geece, May -3, 006 (pp7-) Analytical Solutions fo Confined Aquifes with non constant Pumping using Compute Algeba
More informationTHE LAPLACE EQUATION. The Laplace (or potential) equation is the equation. u = 0. = 2 x 2. x y 2 in R 2
THE LAPLACE EQUATION The Laplace (o potential) equation is the equation whee is the Laplace opeato = 2 x 2 u = 0. in R = 2 x 2 + 2 y 2 in R 2 = 2 x 2 + 2 y 2 + 2 z 2 in R 3 The solutions u of the Laplace
More informationConservative Averaging Method and its Application for One Heat Conduction Problem
Poceedings of the 4th WSEAS Int. Conf. on HEAT TRANSFER THERMAL ENGINEERING and ENVIRONMENT Elounda Geece August - 6 (pp6-) Consevative Aveaging Method and its Application fo One Heat Conduction Poblem
More informationDuality 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 informationA THREE CRITICAL POINTS THEOREM AND ITS APPLICATIONS TO THE ORDINARY DIRICHLET PROBLEM
A THREE CRITICAL POINTS THEOREM AND ITS APPLICATIONS TO THE ORDINARY DIRICHLET PROBLEM DIEGO AVERNA AND GABRIELE BONANNO Abstact. The aim of this pape is twofold. On one hand we establish a thee citical
More informationComputers and Mathematics with Applications
Computes and Mathematics with Applications 58 (009) 9 7 Contents lists available at ScienceDiect Computes and Mathematics with Applications jounal homepage: www.elsevie.com/locate/camwa Bi-citeia single
More informationAnalytical evaluation of 3D BEM integral representations using complex analysis
BIR Wokshop 5w5052 Moden Applications of Complex Vaiables: Modeling, Theoy and Computation Analytical evaluation of 3D BEM integal epesentations using complex analysis onia Mogilevskaya Depatment of Civil,
More informationMathematical Model of Magnetometric Resistivity. Sounding for a Conductive Host. with a Bulge Overburden
Applied Mathematical Sciences, Vol. 7, 13, no. 7, 335-348 Mathematical Model of Magnetometic Resistivity Sounding fo a Conductive Host with a Bulge Ovebuden Teeasak Chaladgan Depatment of Mathematics Faculty
More informationNuclear size corrections to the energy levels of single-electron atoms
Nuclea size coections to the enegy levels of single-electon atoms Babak Nadii Nii a eseach Institute fo Astonomy and Astophysics of Maagha (IAAM IAN P. O. Box: 554-44. Abstact A study is made of nuclea
More informationChapter 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 informationSOLVING THE VISCOUS COMPOSITE CYLINDER PROBLEM BY SOKOLOV S METHOD
1 SOLVING THE VISCOUS COMPOSITE CYLINDER PROBLEM BY SOKOLOV S METHOD NGO By CO. H. TRAN, PHONG. T. Faculty of Mathematics & Infomatics, Univesity of Natual Sciences VNU-HCM Abstact : The pape pesents some
More informationMath 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 informationSOME GENERAL NUMERICAL RADIUS INEQUALITIES FOR THE OFF-DIAGONAL PARTS OF 2 2 OPERATOR MATRICES
italian jounal of pue and applied mathematics n. 35 015 (433 44) 433 SOME GENERAL NUMERICAL RADIUS INEQUALITIES FOR THE OFF-DIAGONAL PARTS OF OPERATOR MATRICES Watheq Bani-Domi Depatment of Mathematics
More informationGreen s Identities and Green s Functions
LECTURE 7 Geen s Identities and Geen s Functions Let us ecall The ivegence Theoem in n-dimensions Theoem 7 Let F : R n R n be a vecto field ove R n that is of class C on some closed, connected, simply
More informationEM 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 informationRelating 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 informationGradient-based Neural Network for Online Solution of Lyapunov Matrix Equation with Li Activation Function
Intenational Confeence on Infomation echnology and Management Innovation (ICIMI 05) Gadient-based Neual Netwok fo Online Solution of Lyapunov Matix Equation with Li Activation unction Shiheng Wang, Shidong
More informationStanford 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 informationSOME SOLVABILITY THEOREMS FOR NONLINEAR EQUATIONS
Fixed Point Theoy, Volume 5, No. 1, 2004, 71-80 http://www.math.ubbcluj.o/ nodeacj/sfptcj.htm SOME SOLVABILITY THEOREMS FOR NONLINEAR EQUATIONS G. ISAC 1 AND C. AVRAMESCU 2 1 Depatment of Mathematics Royal
More informationContact impedance of grounded and capacitive electrodes
Abstact Contact impedance of gounded and capacitive electodes Andeas Hödt Institut fü Geophysik und extateestische Physik, TU Baunschweig The contact impedance of electodes detemines how much cuent can
More informationTESTING THE VALIDITY OF THE EXPONENTIAL MODEL BASED ON TYPE II CENSORED DATA USING TRANSFORMED SAMPLE DATA
STATISTICA, anno LXXVI, n. 3, 2016 TESTING THE VALIDITY OF THE EXPONENTIAL MODEL BASED ON TYPE II CENSORED DATA USING TRANSFORMED SAMPLE DATA Hadi Alizadeh Noughabi 1 Depatment of Statistics, Univesity
More informationOn 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 informationAsymptotically 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 informationOn Computing Optimal (Q, r) Replenishment Policies under Quantity Discounts
Annals of Opeations Reseach manuscipt No. will be inseted by the edito) On Computing Optimal, ) Replenishment Policies unde uantity Discounts The all - units and incemental discount cases Michael N. Katehakis
More informationLecture 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 informationHomework 7 Solutions
Homewok 7 olutions Phys 4 Octobe 3, 208. Let s talk about a space monkey. As the space monkey is oiginally obiting in a cicula obit and is massive, its tajectoy satisfies m mon 2 G m mon + L 2 2m mon 2
More information7.2. Coulomb s Law. The Electric Force
Coulomb s aw Recall that chaged objects attact some objects and epel othes at a distance, without making any contact with those objects Electic foce,, o the foce acting between two chaged objects, is somewhat
More informationMAGNETIC FIELD AROUND TWO SEPARATED MAGNETIZING COILS
The 8 th Intenational Confeence of the Slovenian Society fo Non-Destuctive Testing»pplication of Contempoay Non-Destuctive Testing in Engineeing«Septembe 1-3, 5, Potoož, Slovenia, pp. 17-1 MGNETIC FIELD
More informationDonnishJournals
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 informationMitscherlich s Law: Sum of two exponential Processes; Conclusions 2009, 1 st July
Mitschelich s Law: Sum of two exponential Pocesses; Conclusions 29, st July Hans Schneebege Institute of Statistics, Univesity of Elangen-Nünbeg, Gemany Summay It will be shown, that Mitschelich s fomula,
More informationApplication of homotopy perturbation method to the Navier-Stokes equations in cylindrical coordinates
Computational Ecology and Softwae 5 5(): 9-5 Aticle Application of homotopy petubation method to the Navie-Stokes equations in cylindical coodinates H. A. Wahab Anwa Jamal Saia Bhatti Muhammad Naeem Muhammad
More informationCALCULATING THE NUMBER OF TWIN PRIMES WITH SPECIFIED DISTANCE BETWEEN THEM BASED ON THE SIMPLEST PROBABILISTIC MODEL
U.P.B. Sci. Bull. Seies A, Vol. 80, Iss.3, 018 ISSN 13-707 CALCULATING THE NUMBER OF TWIN PRIMES WITH SPECIFIED DISTANCE BETWEEN THEM BASED ON THE SIMPLEST PROBABILISTIC MODEL Sasengali ABDYMANAPOV 1,
More informationChem 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 informationReview: 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 informationJournal 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 informationA Backward Identification Problem for an Axis-Symmetric Fractional Diffusion Equation
Mathematical Modelling and Analysis Publishe: Taylo&Fancis and VGTU Volume 22 Numbe 3, May 27, 3 32 http://www.tandfonline.com/tmma https://doi.og/.3846/3926292.27.39329 ISSN: 392-6292 c Vilnius Gediminas
More informationOn absence of solutions of a semi-linear elliptic equation with biharmonic operator in the exterior of a ball
Tansactions of NAS of Azebaijan, Issue Mathematics, 36, 63-69 016. Seies of Physical-Technical and Mathematical Sciences. On absence of solutions of a semi-linea elliptic euation with bihamonic opeato
More informationELASTIC ANALYSIS OF CIRCULAR SANDWICH PLATES WITH FGM FACE-SHEETS
THE 9 TH INTERNATIONAL CONFERENCE ON COMPOSITE MATERIALS ELASTIC ANALYSIS OF CIRCULAR SANDWICH PLATES WITH FGM FACE-SHEETS R. Sbulati *, S. R. Atashipou Depatment of Civil, Chemical and Envionmental Engineeing,
More informationScattering in Three Dimensions
Scatteing in Thee Dimensions Scatteing expeiments ae an impotant souce of infomation about quantum systems, anging in enegy fom vey low enegy chemical eactions to the highest possible enegies at the LHC.
More informationNumerical Inversion of the Abel Integral Equation using Homotopy Perturbation Method
Numeical Invesion of the Abel Integal Equation using Homotopy Petubation Method Sunil Kuma and Om P Singh Depatment of Applied Mathematics Institute of Technology Banaas Hindu Univesity Vaanasi -15 India
More informationGROWTH ESTIMATES THROUGH SCALING FOR QUASILINEAR PARTIAL DIFFERENTIAL EQUATIONS
Annales Academiæ Scientiaum Fennicæ Mathematica Volumen 32, 2007, 595 599 GROWTH ESTIMATES THROUGH SCALING FOR QUASILINEAR PARTIAL DIFFERENTIAL EQUATIONS Teo Kilpeläinen, Henik Shahgholian and Xiao Zhong
More informationExceptional regular singular points of second-order ODEs. 1. Solving second-order ODEs
(May 14, 2011 Exceptional egula singula points of second-ode ODEs Paul Gaett gaett@math.umn.edu http://www.math.umn.edu/ gaett/ 1. Solving second-ode ODEs 2. Examples 3. Convegence Fobenius method fo solving
More information12th WSEAS Int. Conf. on APPLIED MATHEMATICS, Cairo, Egypt, December 29-31,
th WSEAS Int. Conf. on APPLIED MATHEMATICS, Caio, Egypt, Decembe 9-3, 7 5 Magnetostatic Field calculations associated with thick Solenoids in the Pesence of Ion using a Powe Seies expansion and the Complete
More informationON THE TWO-BODY PROBLEM IN QUANTUM MECHANICS
ON THE TWO-BODY PROBLEM IN QUANTUM MECHANICS L. MICU Hoia Hulubei National Institute fo Physics and Nuclea Engineeing, P.O. Box MG-6, RO-0775 Buchaest-Maguele, Romania, E-mail: lmicu@theoy.nipne.o (Received
More informationImplicit Constraint Enforcement for Rigid Body Dynamic Simulation
Implicit Constaint Enfocement fo Rigid Body Dynamic Simulation Min Hong 1, Samuel Welch, John app, and Min-Hyung Choi 3 1 Division of Compute Science and Engineeing, Soonchunhyang Univesity, 646 Eupnae-i
More informationMASSACHUSETTS INSTITUTE OF TECHNOLOGY Physics Department Physics 8.07: Electromagnetism II September 15, 2012 Prof. Alan Guth PROBLEM SET 2
MASSACHUSETTS INSTITUTE OF TECHNOLOGY Physics Depatment Physics 8.07: Electomagnetism II Septembe 5, 202 Pof. Alan Guth PROBLEM SET 2 DUE DATE: Monday, Septembe 24, 202. Eithe hand it in at the lectue,
More informationPerturbation to Symmetries and Adiabatic Invariants of Nonholonomic Dynamical System of Relative Motion
Commun. Theo. Phys. Beijing, China) 43 25) pp. 577 581 c Intenational Academic Publishes Vol. 43, No. 4, Apil 15, 25 Petubation to Symmeties and Adiabatic Invaiants of Nonholonomic Dynamical System of
More informationPhysics 2212 GH Quiz #2 Solutions Spring 2016
Physics 2212 GH Quiz #2 Solutions Sping 216 I. 17 points) Thee point chages, each caying a chage Q = +6. nc, ae placed on an equilateal tiangle of side length = 3. mm. An additional point chage, caying
More information4/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 information6 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 informationDo not turn over until you are told to do so by the Invigilator.
UNIVERSITY OF EAST ANGLIA School of Mathematics Main Seies UG Examination 2015 16 FLUID DYNAMICS WITH ADVANCED TOPICS MTH-MD59 Time allowed: 3 Hous Attempt QUESTIONS 1 and 2, and THREE othe questions.
More informationFUSE Fusion Utility Sequence Estimator
FUSE Fusion Utility Sequence Estimato Belu V. Dasaathy Dynetics, Inc. P. O. Box 5500 Huntsville, AL 3584-5500 belu.d@dynetics.com Sean D. Townsend Dynetics, Inc. P. O. Box 5500 Huntsville, AL 3584-5500
More informationThis aticle was oiginally published in a jounal published by Elsevie, the attached copy is povided by Elsevie fo the autho s benefit fo the benefit of the autho s institution, fo non-commecial eseach educational
More informationChapter 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 informationChapter 13 Gravitation
Chapte 13 Gavitation In this chapte we will exploe the following topics: -Newton s law of gavitation, which descibes the attactive foce between two point masses and its application to extended objects
More information1 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 informationworking 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 informationA matrix method based on the Fibonacci polynomials to the generalized pantograph equations with functional arguments
A mati method based on the Fibonacci polynomials to the genealized pantogaph equations with functional aguments Ayşe Betül Koç*,a, Musa Çama b, Aydın Kunaz a * Coespondence: aysebetuloc @ selcu.edu.t a
More informationAlternative 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 informationA 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 informationRECTIFYING THE CIRCUMFERENCE WITH GEOGEBRA
ECTIFYING THE CICUMFEENCE WITH GEOGEBA A. Matín Dinnbie, G. Matín González and Anthony C.M. O 1 Intoducction The elation between the cicumfeence and the adius of a cicle is one of the most impotant concepts
More informationNumerical approximation to ζ(2n+1)
Illinois Wesleyan Univesity Fom the SelectedWoks of Tian-Xiao He 6 Numeical appoximation to ζ(n+1) Tian-Xiao He, Illinois Wesleyan Univesity Michael J. Dancs Available at: https://woks.bepess.com/tian_xiao_he/6/
More informationFresnel 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 informationON LACUNARY INVARIANT SEQUENCE SPACES DEFINED BY A SEQUENCE OF MODULUS FUNCTIONS
STUDIA UNIV BABEŞ BOLYAI, MATHEMATICA, Volume XLVIII, Numbe 4, Decembe 2003 ON LACUNARY INVARIANT SEQUENCE SPACES DEFINED BY A SEQUENCE OF MODULUS FUNCTIONS VATAN KARAKAYA AND NECIP SIMSEK Abstact The
More informationq 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 informationInformation 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 informationarxiv: v1 [math.na] 8 Feb 2013
A mixed method fo Diichlet poblems with adial basis functions axiv:1302.2079v1 [math.na] 8 Feb 2013 Nobet Heue Thanh Tan Abstact We pesent a simple discetization by adial basis functions fo the Poisson
More informationMASSACHUSETTS INSTITUTE OF TECHNOLOGY Physics Department. Problem Set 10 Solutions. r s
MASSACHUSETTS INSTITUTE OF TECHNOLOGY Physics Depatment Physics 8.033 Decembe 5, 003 Poblem Set 10 Solutions Poblem 1 M s y x test paticle The figue above depicts the geomety of the poblem. The position
More informationA method for solving dynamic problems for cylindrical domains
Tansactions of NAS of Azebaijan, Issue Mechanics, 35 (7), 68-75 (016). Seies of Physical-Technical and Mathematical Sciences. A method fo solving dynamic poblems fo cylindical domains N.B. Rassoulova G.R.
More informationA 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 informationTheWaveandHelmholtzEquations
TheWaveandHelmholtzEquations Ramani Duaiswami The Univesity of Mayland, College Pak Febuay 3, 2006 Abstact CMSC828D notes (adapted fom mateial witten with Nail Gumeov). Wok in pogess 1 Acoustic Waves 1.1
More informationLocalization 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 informationDIRECT INTERBAND LIGHT ABSORPTION IN A SPHERICAL QUANTUM DOT WITH THE MODIFIED PÖSCHEL-TELLER POTENTIAL
Lase Physics Intenational Jounal of Moden Physics: Confeence Seies Vol. 5 () 4 Wold Scientific Publishing Company DOI:.4/S945767 DIRECT INTERBAND LIGHT ABSORPTION IN A SPHERICAL QANTM DOT WITH THE MODIFIED
More informationME 210 Applied Mathematics for Mechanical Engineers
Tangent and Ac Length of a Cuve The tangent to a cuve C at a point A on it is defined as the limiting position of the staight line L though A and B, as B appoaches A along the cuve as illustated in the
More informationQuestion 1: The dipole
Septembe, 08 Conell Univesity, Depatment of Physics PHYS 337, Advance E&M, HW #, due: 9/5/08, :5 AM Question : The dipole Conside a system as discussed in class and shown in Fig.. in Heald & Maion.. Wite
More informationThe Schwartzchild Geometry
UNIVERSITY OF ROCHESTER The Schwatzchild Geomety Byon Osteweil Decembe 21, 2018 1 INTRODUCTION In ou study of geneal elativity, we ae inteested in the geomety of cuved spacetime in cetain special cases
More information