One-dimensional I test and direction vector I test with array references by induction variable

Size: px
Start display at page:

Download "One-dimensional I test and direction vector I test with array references by induction variable"

Transcription

1 Int. J. High Performance Computing an Networking, Vol. 3, No. 4, One-imensional I test an irection vector I test with array references by inuction variable Minyi Guo School of Computer Science an Engineering, University of Aizu, Aizu-Wakamatsu City, Fukushima, Japan State Key Lab for Novel Software Technology, Nanjing University, PRC Fax: minyi@u-aizu.ac.jp Weng-Long Chang* Department of Computer Science an Information Engineering, National Kaohsiung University of Applie Sciences, 415 Chien Kung Roa, 807 Kaohsiung, Taiwan changwl@ mail.stut.eu.tw *Corresponing author Jian Lu State Key Lab for Novel Software Technology, Nanjing University, PRC lj@nju.eu.cn Minglu Li Department of Computer Science an Engineering, Shanghai JiaoTong University, Shanghai, China li-ml@cs.stut.eu.cn Abstract: In this paper, theoretical aspects to emonstrate the accuracy of the Interval Test (the I test an the irection vector I test) to be applie for resolving the problem state above is presente. Also, it is prove from the propose theoretical aspects that uner a specific irection vector θ = ( = 1,, = ) there are integer-value solutions for one-imensional arrays with subscripts forme by inuction variable an uner other specific irection vectors there are no integer-value solutions. Experiments with benchmarks, cite from Parallel loop, Vector loop an TRFD (Perfect benchmark), reveal that our framework can properly enhance the precision of ata epenence analysis for one-imensional arrays with subscripts mentione above. Keywors: parallelising/vectorising compilers; ata epenence analysis; loop parallelisation; loop vectorisation; automatic loop transformation. Reference to this paper shoul be mae as follows: Guo, M., Chang, W-L., Lu, J. an Li, M. (2005) One-imensional I test an irection vector I test with array references by inuction variable, Int. J. High Performance Computing an Networking, Vol. 3, No. 4, pp Copyright 2005 Inerscience Enterprises Lt.

2 220 M. GUO, W-L. CHANG, J. LU AND M. LI Biographical notes: Minyi Guo is currently an Associate Professor at the Department of Computer Software, at the University of Aizu, Japan. Guo has serve as programme committee or organising committee member for many international conferences. He is the Eitor-in-Chief of the Journal of Embee Systems. He is also on the eitorial boar of the International Journal of High Performance Computing an Networking, the Journal of Embee Computing, the Journal of Parallel an Distribute Scientific an Engineering Computing an the International Journal of Computer an Applications. Guo s research interests inclue parallel an istribute processing, parallelising compilers, ata parallel languages, ata mining, molecular computing an software engineering. He is a member of the ACM, IEEE, an IEICE. Weng-Long Chang receive his PhD egree in computer science an information engineering from the National Cheng Kung University, Taiwan, Republic of China, in He is currently an Assistant Professor at the Southern Taiwan University of Technology. His research interests inclue molecular computing an languages an compilers for parallel computing. Jian Lu receive his BS egree, MS egree, an PhD egree from Nanjing University of PR China in 1982, 1984, 1988, respectively. Currently, he is a Professor at Department of Computer Science an Technology, Vice Director of the Institute of Computer Software at Nanjing University an Director of State Key Laboratory of Novel Software Technology of P.R. China. His research interests inclue formal methos, mobile agent technology an software architecture. Dr. Minglu Li is a Full Professor an Vice Char of Department of Computer Science an Engineering of Shanghai Jiao Tong University (SJTU). Now, he is also Director of Gri Computing Center of SJTU, a member of Expert Committee of ChinaGri Program of Ministry of Eucation of China, PI of ShanghaiGri of Science an Technology Commission of Shanghai Municipality(STCSM), an Associate Eitor of International Journal of Gri an Utility Computing, an member of eitorial boar of International Journal of Web Services Research, a member of executive committee of Technical Committee on Services Computing of IEEE. 1 INTRODUCTION Achieving a goo ata epenence analysis is a critical, issue in orer to reuce the communication overhea an to exploit parallelism of applications as much as possible. This task becomes even more important in istribute memory systems where, in aition to accomplishing a high parallelism of computation an low communication overhea among the processors, it is essential to evelop a new analysis technique for ata epenence. There are several well-known ata epenence analysis algorithms applicable for one-imensional arrays uner constant bouns or variable bouns: the GCD test (Banerjee, 1988, 1993, 1997), Banerjee s metho (Banerjee, 1988, 1993, 1997), the I test an the irection vector I test (Kong et al., 1991; Niezielski an Psarris, 1999; Psarris et al., 1991, 1993; Psarris an Kyriakopoulos, 1999), the extene I test (Chang an Chu, 1998), the generalise irection vector I test (Chang an Chu, 2001) an the interval reuction test (Huang an Yang, 2000). There are also several well-known ata epenence analysis algorithms applicable for multi-imensional couple arrays uner constant bouns or variable bouns: the generalise GCD test (Banerjee, 1988, 1993, 1997), the Lamba test (Li et al., 1990), the generalise Lamba test (Chang et al., 1999), the multi-imensional I test (Chang et al., 2001), the multi-imensional irection vector I test (Chang et al., 2002), the Power test (Wolfe an Tseng, 1992) an the Omega test (Pugh, 1992). There are several well-known ata epenence analysis algorithms applicable for arrays with linear subscripts with symbolic coefficients, or with non-linear subscripts uner symbolic bouns: the infinity Banerjee test (Petersen, 1993), the Range test (Blume an Eigenmann, 1998), the infinity Lamba test (Chang an Chu, 2000), the access range test (Paek, 1997; Hoeflinger, 1998) an one analysis metho for pointers an inuction variables (Wu, 2001). In this paper, we propose a sophisticate technique of ata epenence analysis. This approach is to test epenence if there are integer-value solutions for one-imensional arrays with subscripts forme by inuction variable. Without irection vectors, there are integer-value solutions for one-imensional arrays with subscripts forme by inuction variable. Furthermore, it is also shown that uner a specific irection vector θ = ( = 1,, = ) there are integer-value solutions for one-imensional arrays with subscripts forme by inuction variable. Shen et al. (1992) ha inicate that in real programmes one-imensional arrays with subscripts forme by inuction variable occur quite frequently. A -neste loop accessing a one-imensional array with subscripts forme by inuction variable is shown in Figure 1.

3 ONE-DIMENSIONAL I TEST AND DIRECTION VECTOR I TEST WITH ARRAY REFERENCES BY INDUCTION VARIABLE 221 Figure 1 An example of a neste loop with inuction variable Inuction variable is a one scalar integer variable, which is use in a loop, to simulate loop s inex variables: it is incremente or ecremente by a constant amount in each iteration. Every inuction variable can be replace by a linear function in loop s inex variables. The transformation, which oes so, is calle inuction variable substitution. Since the variable K in Figure 1 is one inuction variable, it can be replace z (( I1 L1) ( p1= 2( U p L 1)) ( 1 p + + I 1 2 L2) ( p2= 3( U p L 2 p2 + 1) + + ( i L where is the number of common loops an z is one integer variable (Banerjee, 1988, 1993, 1997; Paek, 1997; Hoeflinger, 1998). Therefore, the coe in Figure 1 is transforme into the coe in Figure 2, after finishing the processing of inuction variable substitution for the variable K. (where n = 2). Furthermore, variables X 2k 1 an X 2k (1 k ) are ifferent instances of the same loop iteration variable, I k. Assume that L 2k 1, L 2k, U 2k 1 an U 2k are, respectively, lower bouns an upper bouns for X 2k 1 an X 2k. Because variables X 2k 1 an X 2k (1 k ) are ifferent instances of the same loop iteration variable, I k, lower bouns for X 2k 1 an X 2k are the same an upper bouns for X 2k 1 an X 2k also are the same. The problem of etermining whether there exists epenence for the array A between S 1 an S 2 in Figure 2 can be reuce to that of checking whether one system of a linear equation with n unknown variables has a simultaneous integer solution, which satisfies the constraints for each variable in the system. It is assume that one linear equation in a system is written as: z ( + 1) ( U L + 1) ( U L ( U L + + ( X2 3 X2 2) ( U L + 1) + ( X X )) = 0, where each L k an each U k are an integer variable an are, respectively, one lower boun an one upper boun for the k-the loop for (1 k ). Because z is the greatest common ivisor for all the coefficients in the left-han sie of equation (1), all the coefficients in equation (1) are ivie by z an equation (1) is rewritten as + 1) ( U L + 1) ( U L ( U L + + ( X2 3 X2 2) ( U L + 1) + ( X2 1 X2) = 0. (2) (1) It is postulate that the constraints to each variable in equation (2) are represente as L X an X U, (3) k 2k 1 2k k Figure 2 The transforme loop after inuction variable substitution for the inuction variable K Because epenence between S 1 an S 2 in Figure 2 may arise in ifferent iterations of the common loops, we eal with the loop iteration variables reference in S 1 as being ifferent variables from those reference in S 2 subject to common loop bouns. Therefore, when analysing epenence arising from a statement pair neste in common loops, the problem will involve n unique variable where (1 k ). Let us use an example to make clear the illustrations state above. Consier the neste o-loop in Figure 3. The variable K in Figure 3(a) is incremente by one in each iteration in the neste loop. So it is a one inuction variable. After finishing the processing of inuction variable substitution for the inuction variable K, the result is shown in Figure 3(b). In Figure 3(b), the lower an upper bouns of the first (outer) loop an the secon (inner) loop are, respectively, 1 an 10. Therefore, the bouns of the o-loop are constants. This o-loop executes 100 iterations by consecutively assigning the values 1, 2,, 10 to J an I an executing the boy (the main statement S) exactly once in each iteration. The net effect of the o-loop execution is then the orere execution of the statements:

4 222 M. GUO, W-L. CHANG, J. LU AND M. LI A(1) = B(1) A(2) = B(2) A(100) = B(100). (a) Figure 3 A neste o-loop in Fortran (a) a neste o-loop with inuction variable an (b) the transforme loop after inuction variable substitution for the inuction variable K To ascertain whether two references to the one-imensional array A may refer to the same element of A, it has to be checke if the following linear equation 10 X 1 10 X 2 + X 3 X 4 = 0 has a simultaneous integer solution uner the constant bouns 1 X 1, X 2 10, an 1 X 3, X It is well known that the problem of fining integer-value solutions to a system of linear equations is NP-har. Therefore, in practice, most well-known ata epenence analysis algorithms are use to solve as many particular cases of this problem as possible. The rest of this paper proffers the following: In Section 2, the summary accounts of ata epenence an interval equation are presente. In Section 3, the theoretical aspects of etermining whether there are integer-value solutions for one linear equation (2) with symbolic coefficients uner the bouns of equation (3) are escribe. Experimental results are given in Section 4. Finally, brief conclusions are rawn in Section 5. 2 BACKGROUND The summary accounts of ata epenence an interval equation are introuce briefly in this section. 2.1 Data epenence It is assume that S 1 an S 2 are two statements within the loop in Figure 2. The loop is presume to contain common loops. Statements S 1 an S 2 are postulate to be embee in common loops. An array A is suppose to appear simultaneously within statements S 1 an S 2. If S 2 uses the element of the array A efine first by S 1, then S 2 is true-epenent on S 1. If S 2 efines the element of the array A use first by S 1, then S 2 is anti-epenent on S 1. If S 2 (b) reefines the element of the array A efine first by S 1, then S 2 is output-epenent on S 1. Each iteration of a loop neste is ientifie by an iteration vector, whose elements are the values of the iteration variables for that iteration. For example, the instance of the statement S 1 uring iteration i = ( i1,, i ) is enote S1( i ); the instance of the statement S 2 uring iteration j = ( j1,, j ) is enote S 2 ( j). If (i 1,, i ) is ientical to (j 1,, j ) or (i 1,, i ) precees (j 1,, j ) lexicographically, then S 1 () i is sai to precee S 2 ( j), enote S1( i ) < S2 ( j). Otherwise, S2 ( j) is sai to precee S1( i ), enote S 1 () i > S2 ( j). In the following, Definition 1 efines irection vectors. Definition 1: A vector of the form θ = ( θ1,, θ ) is terme as a irection vector. The irection vector ( θ1,, θ ) is sai to be the irection vector from S1( i ) to S 2 ( j) if for (1 k ) i k θ k j k, i.e., the relation θ k is one element in a set {<,=,>}. 2.1 Interval equation A linear equation with symbolic coefficients uner the bouns of equation (3) will be sai to be integer solvable if the equation has an integer-value solution satisfying the bouns of each variable. Definitions 2 an 3, respectively, efine integer interval an interval equation (Kong et al., 1991; Psarris et al., 1993). Definition 2: Let [α 1,α 2 ] represent the integer intervals from α 1 to α 2, i.e., the set of all integers between α 1 an α 2. Definition 3: Let a 1,, a n 1, a n, L an U be integers. A linear equation ax 1 1+ ax an 1Xn 1 + ax n n = [ LU, ], (4) which is referre to as an interval equation will be use to enote the set of orinary equations consisting of: ax 1 1+ ax an 1Xn 1+ anxn = L ax 1 1+ ax an 1Xn 1+ ax n n = L+ 1 ax + ax + + a X + a X = U n 1 n 1 n n An interval equation (4) will be sai to be integer solvable if one of the equations in the set, which it efines, is integer solvable. The immeiate way to etermine this is to test if an integer in between L an U is ivisible by the GCD of the coefficients of the left-han-sie terms. If L > U in an interval equation (4), then there are no integer solutions for this interval equation. If the expression on the left-han sie of an interval equation (4) is reuce to zero items, in the processing of testing, then the interval equation (4) will be sai to be integer solvable if an only if L 0 U.

5 ONE-DIMENSIONAL I TEST AND DIRECTION VECTOR I TEST WITH ARRAY REFERENCES BY INDUCTION VARIABLE DETERMINING INTEGER-VALUED SOLUTIONS FOR ONE-DIMENSIONAL ARRAYS WITH SYMBOLIC COEFFICIENTS UNDER SYMBOLIC BOUNDS Given the ata epenence problem of one-imensional arrays with symbolic coefficients an symbolic bouns, we want to fin that uner what conitions there are integer-value solutions uner what conitions there are no integer-value solutions. We start to iscuss from the case without irection vectors for convenience of presentation. 3.1 Determine integer-value solutions for one-imensional arrays with symbolic coefficients uner symbolic bouns without irection vector One linear equation (2) with symbolic coefficients is euce from etermining whether there exists epenence for the array A between S 1 an S 2 in Figure 2. One boun (3) is erive from lower an upper bouns of loop inex variables in Figure 2. If there is no integer-value solution for one linear equation (2) with symbolic coefficients uner the limits of equation (3), then there is no epenence. Otherwise, there is epenence for one linear equation (2) with symbolic coefficients uner the limits of equation (3). Theorem 1 is presente to show that there are integer-value solutions for one linear equation (2) with symbolic coefficients uner the limits of equation (3). Theorem 1: There are integer-value solutions for one linear equation (2) with symbolic coefficients uner the limits of equation (3). Proof: 1 From the I test in Kong et al. (1991), one linear equation (2) with symbolic coefficients can be rewritten as the following interval equation: + 1) ( U L + 1) ( U L + + ( X2j 1 X2j) (( U ( U L + + ( X2 3 X2 2) ( U L + 1) + ( X X ) = [0,0] Accoring to the I test, the term X 2 in the interval equation (5) is move to the right-han sie to gain the new interval equation: (5) 4 Because (U L + 1) is the greatest common ivisor for all the coefficients in the left-han sie of the interval equation (7), all the coefficients in the interval equation (7) are ivie by (U L + 1) an the new interval equation is obtaine: + 1) ( U L + 1) ( U L ( U L + + ( X2 3 X2 2) ( U L + 1) = [ L U, U L ]. (7) + 1) ( U L + 1) ( U L + + ( X2j 1 X2j) (( U ( U L + + ( X2 3 X2 2) ( U L + 1) L U U L =, = [0,0]. (8) U L + 1 U L Repeat the processing of step (2) (4) until the term in left-han sie of the interval equation is reuce to zero item. Therefore, we will obtain the new interval equation 0 = [L 1 U 1,U 1 L 1 ]. Because U 1 L 1, L 1 U 1 0 U 1 L 1 is true. Thus, it is at once erive that there are integer-value solutions for one linear equation (2) with symbolic coefficients uner the limits of equation (3). The following example is use to explain the power of Theorem 1. Consier the o-loop in Figure 4(a). The front-en of one parallel compiler can recognise that the variable K is one inuction variable an the subscript of the array A is forme by the inuction variable K. After finishing the processing of inuction variable substitution for the inuction variable K, the result is shown in Figure 4(b). Because the coefficient of the subscript of the array A is an unknown variable (at compile time), the Power test an the Omega test cannot be applie to eal with the problem. However, in light of Theorem 1, it is at once conclue that there exists output epenence for the array A uner without consiering any irection vector. Therefore, it is inicate from the result that the precision of Theorem 1 is superior to that of the Omega test an the Power test. ( X X ) (( U L + 1) ( U L ( X X ) (( U L + 1) ( U L ( X X ) (( U L + 1) 2j 1 2j j+ 1 j+ 1 ( U L + + ( X X ) ( U L + 1) X = [ L, U ]. 2 1 (6) (a) (b) 3 In light of the I test, the term X 2 1 in the interval equation (6) is move to the right-han sie to gain the new interval equation. Figure 4 An example of o-loop in Fortran program (a) a neste o-loop with inuction variable an (b) the transforme loop after inuction variable substitution for the inuction variable K

6 224 M. GUO, W-L. CHANG, J. LU AND M. LI 3.2 Determine integer-value solutions for one-imensional arrays with symbolic coefficients uner symbolic bouns an irection vector From Theorem 1, it is very clear that there are integer-value solutions for one linear equation (2) with symbolic coefficients uner the limits of equation (3). Next, uner a specific given irection vector an the limits of equation (3), whether there are integer-value solutions for one linear equation (2) with symbolic coefficients will be iscusse. The following theorems are propose to show that uner what a specific irection vector there are no integer-value solutions uner what a specific irection vector there are integer-value solutions. Theorem 2: There are no integer-value solutions for one linear equation (2) with symbolic coefficients uner the limits of equation (3) an a specific irection vector (<, *,, *), where is the number of common loops an * is any one of {<, =, >}. Proof: 1 From the I test an the irection vector I test in Kong et al. (1991) an Psarris et al. (1993), one linear equation (2) with symbolic coefficients can be rewritten as the following interval equation: + 1) ( U L + 1) ( U L 1+ 1) ( U L + + ( X2 3 X2 2) ( U L + 1) + ( X X ) = [0,0] Repeat the processing of step (2) to step (4) in Theorem 1, until the term in left-han sie of the interval equation is reuce to two items. Therefore, we will obtain the new interval equation: 1 2 (9) ( X X ) = [0,0]. (10) 3 Because the irection vector is < for the terms X 1 an X 2, accoring to the irection vector I test in Psarris et al. (1993), the two terms X 1 an X 2 are move to the right-han sie. Therefore, we obtain the new interval equation 0 = [1, U 1 L 1 ]. Because 1 0 U 1 L 1 is false, it is, therefore, inferre that there are no integer-value solutions for one linear equation (2) with symbolic coefficients uner the limits of equation (3) an a specific irection vector (<, *,, *). Theorem 3: There are no integer-value solutions for one linear equation (2) with symbolic coefficients uner the limits of equation (3) an a specific irection vector (<, *,, *), where is the number of common loops an * is any one of {<, =, >}. Proof: Similar to Theorem 2. Theorem 4: There are no integer-value solutions for one linear equation (2) with symbolic coefficients uner the limits of equation (3) an a specific irection vector (=, θ 2,, θ 1, θ ), where is the number of common loops an θ k is any element of {<, =, >} for 2 k 1 an θ k is any element of {<, >}. Proof: Similar to Theorem 2. Theorem 5: There are integer-value solutions for one linear equation (2) with symbolic coefficients uner the limits of equation (3) an a specific irection vector (=, =,, =), where is the number of common loops. Proof: Similar to Theorem 2. We now use the o-loop in Figure 4 to explain the power of Theorem 2 to Theorem 5. Consier the o-loop in Figure 4(a) an Figure 4(b). Because the coefficient of the subscript of the array A is an unknown variable (at compile time), the Power test an the Omega test cannot be applie to eal with the problem uner any given irection vectors. However, in light of Theorem 5, it is right away inferre that there are integer-value solutions for the array A uner a specific irection vector (=, =). Accoring to Theorem 2 to Theorem 4, it is at once conclue that there are no integer-value solutions for the array A uner other specific irection vectors. Therefore, the front-en of one parallel compiler at once erive that there only exists loop-inepenence output epenence for the o-loop. Because the source statement an sink statement to the epenence relation is the same statement, the epenence relation can be ignore. The o-loop can be execute in parallel moe. Hence, it is inicate from the result that the precision of Theorem 5 is superior to that of the Omega test an the Power test. 3.3 Extening symbolic subscript forme by inuction variable The expression K + C for the array A in Figure 1 to a K + C, where a is an integer variable an is not equal to zero can be easily extene. Since the variable K in Figure 1 is one inuction variable, it can be replace by z (( I L ) ( ( U L + ( I L ) 1 1 p1 = 2 p1 p1 2 2 p2 = 3 U p L 2 p i 2 L ( ( + + ( where is the number of common loops an z is one integer variable (Banerjee, 1997; Paek, 1997; Hoeflinger, 1988).

7 ONE-DIMENSIONAL I TEST AND DIRECTION VECTOR I TEST WITH ARRAY REFERENCES BY INDUCTION VARIABLE 225 Therefore, the extene expression of the array A is transforme into the coe in Figure 5, after finishing the processing of inuction variable substitution. Eigenman et al., 1988). One hunre an sixty four pairs of one-imensional array references were observe to have subscripts forme by inuction variable. The propose theorems are only applie to test those one-imensional arrays with subscripts forme by inuction variable. It is very clear from the result shown in Table 1 that the propose theorems coul properly solve whether there are efinitive results for one-imensional arrays with subscripts forme by inuction variable. Table 1 The result of solving whether there are integer-value solutions for one-imensional arrays with subscripts forme by inuction variable The number The number of The number of no Benchmark of arrays teste integer-value solutions integer-value solutions Parallel loop Vector loop TRFD Figure 5 The result is obtaine after finishing the processing of inuction variable substitution for the extene expression of the array A The problem of etermining whether there exists epenence for the array A between S 1 an S 2 in Figure 5 is actually equal to that of checking whether there are integer-value solutions for one new linear equation (11) uner the constraints of equation (3). It is assume that the new linear equation (11) is written as a z (( X1 X2) (( U2 L2 + 1) ( U L + 1) ( U L ( U L + + ( X2 3 X2 2) ( U L + 1) + ( X X )) = 0, (11) where each L k an each U k are an integer variable an are, respectively, one lower boun an one upper boun for the k-th loop for 1 k. Because a z is not equal to zero, all the coefficients in equation (11) are ivie by a z an equation (11) is exactly equal to equation (2). Therefore, Theorem 5 also can be applie to eal with whether there is epenent relation for the array A between S 1 an S 2 in Figure 5. 4 EXPERIMENTAL RESULTS The propose theorems have been teste an experimente on the coes abstracte from three numerical packages: Parallel Loop, Vector Loop an TRFD (Perfect Benchmark) (Dongarra et al., 1991; Levine et al., 1991; The Omega test an the Power test have been implemente base on Wolfe an Tseng (1992) an Pugh (1992), to compare their effects with that of the propose theorems. The Omega test an the Power test were applie to resolve the same 164 pairs of one-imensional arrays with symbolic coefficients. It is very clear from the result shown in Table 2 that the Omega test an the Power test coul not be applie for etermining whether there are integer-value solutions for one-imensional arrays with symbolic coefficients. Therefore, it is inicate from the results shown in Tables 1 an 2 that the precision of the propose theorems is superior to that of the Omega test an the Power test. Table 2 The result of the Omega an Power tests for solving whether there are integer-value solutions for 164 pairs of one-imensional arrays with symbolic coefficients Depenence metho The Omega test The Power test The number of arrays teste The number of integer-value solution The number of no integer-value solutions : represents that the Omega an Power tests coul not eal with them. 5 CONCLUSIONS The front-en of a parallelising compiler can easily recognise inuction variable that forms subscripts of one-imensional arrays. Inuction variable can be replace by a linear function in form of loop s inex-variable after the front-en of a parallelising compiler finishes the processing of inuction variable substitution. This paper presents theoretical aspects to emonstrate the accuracy of

8 226 M. GUO, W-L. CHANG, J. LU AND M. LI the Interval Test (the I test an the irection vector I test) to be applie for resolving the problem state above. Also, it is prove from the propose theoretical aspects that uner a specific irection vector θ = ( = 1,, = ) there are integer-value solutions for one-imensional arrays with subscripts forme by inuction variable an uner other specific irection vectors there are no integer-value solutions. Experiments with benchmarks, cite from Parallel loop, Vector loop an TRFD (Perfect benchmark), reveal that this framework can properly enhance the precision of ata epenence analysis for one-imensional arrays with subscripts mentione above. The propose theorems can improve the precision of ata epenence analysis for one-imensional arrays with references forme by inuction variable. Depening on the application omains, it is suggeste that the propose theorems be applie together with the front-en of a parallelising compiler to provie ata epenence analysis for one-imensional arrays with references forme by inuction variable. ACKNOWLEDGEMENT This work was partially supporte by the 973 Program of China uner Grant 2002CB an NSFC of Jiangsu Province uner Grant BK REFERENCES Banerjee, U. (1988) Depenence Analysis for Supercomputing, Kluwer Acaemic Publishers, Norwell, Massachusetts. Banerjee, U. (1993) Loop Transformations for Restructuring Compilers: The Founations, Kluwer Acaemic Publishers, Norwell, Massachusetts. Banerjee, U. (1997) Depenence Analysis, Kluwer Acaemic Publishers, Norwell, Massachusetts. Blume, W. an Eigenmann, R. (1998) Nonlinear an symbolic ata epenence testing, IEEE Transaction on Parallel an Distribute Systems, Vol. 9, No. 12, pp Chang, Chu an Wu, J. (1999) The generalize lamba test: a multi-imensional version of banerjee s algorithm, International Journal of Parallel an Distribute Systems an Networks, Vol. 2, No. 2, pp Chang, W-L. an Chu, C-P. (1998) The extension of the I test, Parallel Computing, Vol. 24, pp Chang, W-L. an Chu, C-P. (2000) The infinity lamba test: a multi-imensional version of Banerjee s infinity test, Parallel Computing, Vol. 26, No. 10, August, pp Chang, W-L. an Chu, C-P. (2001) The generalize irection vector I test, Parallel Computing, Vol. 27, No. 8, July, pp Chang, W-L., Chu, C-P. an Wu, J-H. (2001a) A multi-imensional version of the I test, Parallel Computing, Vol. 27, No. 13, September, pp Chang, W-L., Chu, C-P. an Wu, J-H. (2002) A precise epenence analysis for multi-imensional arrays uner specific epenence irection, Journal of System an Software, Vol. 63, No. 2, pp Dongarra, J., Reinhart, M.F. an Russell, S.J. (1991) Parallel loops a test suite for parallelizing compilers: escription an example results, Parallel Computing, Vol. 17, pp Eigenman, R., Hoeflinger, J. an Paua, D. (1998) On the automatic parallelization of the perfect benchmarks, IEEE Transactions on Parallel an Distribute Systems, Vol. 9, No. 1, January, pp Hoeflinger, J.P. (1998) Inter-proceural Parallelization using Memory Classification Analysis, PhD Thesis, University of Illinois at Urbana-Champaign. Huang, T-C. an Yang, C-M. (2000) Data epenence analysis for array references, The Journal of System an Software, Vol. 52, pp Kong, X., Klappholz, D. an Psarris, K. (1991) The I test, IEEE Transaction on Parallel an Distribute Systems, Vol. 2, No. 3, July, pp Levine, D., Callahan, D. an Dongarra, J. (1991) A comparative stuy of automatic vectorizing compilers, Parallel Computing, Vol. 17, pp Li, Z., Yew, P-C. an Zhu, C-Q. (1990) An efficient ata epenence analysis for parallelizing compilers, IEEE Transaction on Parallel an Distribute Systems, Vol. 1, No. 1, January, pp Niezielski, D. an Psarris, K. (1999) An analytical comparison of the I-test an Omega test, LCPC 99: Twelfth International Workshop on Languages an Compilers for Parallel Computing, pp Paek, Y. (1997) Compiling for Distribute Memory Multiprocessors base on Access Region Analysis, PhD Thesis, University of Illinois at Urbana-Champaign. Petersen, P.M. (1993) Evaluation of Programs an Parallelizing Compilers Using Dynamic Analysis Techniques, PhD Thesis, University of Illinois at Urbana-Champaign, January. Psarris, K. an Kyriakopoulos, K. (1999) Data epenence testing in practice, Proceeings of the 1999 International Conference on Parallel Architectures an Compilation Techniques, pp Psarris, K., Klappholz, D. an Kong, X. (1991) On the accuracy of the banerjee test, Journal of Parallel an Distribute Computing, Vol. 12, No. 2, June, pp Psarris, K., Kong, X. an Klappholz, D. (1993) The irection vector I test, IEEE Transaction on Parallel an Distribute Systems, Vol. 4, No. 11, pp Pugh, W. (1992) A practical algorithm for exact array epenence analysis, Communication of the ACM, Vol. 35, No. 8, August, pp Shen, Z., Li, Z. an Yew, P-C. (1992) An empirical stuy of Fortran programs for parallelizing compilers, IEEE Transaction on Parallel an Distribute Systems, Vol. 1, No. 3, July, pp Wolfe, M. an Tseng, C-W. (1992) The power test for ata epenence, IEEE Transaction on Parallel an Distribute Systems, Vol. 3, No. 5, September, pp Wu, P. (2001) Analysis of Pointers an Inuction Variable an Container Objects for Depenence Tests, PhD Thesis, University of Illinois at Urbana-Champaign, UIUCDCS-R

Time-of-Arrival Estimation in Non-Line-Of-Sight Environments

Time-of-Arrival Estimation in Non-Line-Of-Sight Environments 2 Conference on Information Sciences an Systems, The Johns Hopkins University, March 2, 2 Time-of-Arrival Estimation in Non-Line-Of-Sight Environments Sinan Gezici, Hisashi Kobayashi an H. Vincent Poor

More information

'HVLJQ &RQVLGHUDWLRQ LQ 0DWHULDO 6HOHFWLRQ 'HVLJQ 6HQVLWLYLW\,1752'8&7,21

'HVLJQ &RQVLGHUDWLRQ LQ 0DWHULDO 6HOHFWLRQ 'HVLJQ 6HQVLWLYLW\,1752'8&7,21 Large amping in a structural material may be either esirable or unesirable, epening on the engineering application at han. For example, amping is a esirable property to the esigner concerne with limiting

More information

EVALUATING HIGHER DERIVATIVE TENSORS BY FORWARD PROPAGATION OF UNIVARIATE TAYLOR SERIES

EVALUATING HIGHER DERIVATIVE TENSORS BY FORWARD PROPAGATION OF UNIVARIATE TAYLOR SERIES MATHEMATICS OF COMPUTATION Volume 69, Number 231, Pages 1117 1130 S 0025-5718(00)01120-0 Article electronically publishe on February 17, 2000 EVALUATING HIGHER DERIVATIVE TENSORS BY FORWARD PROPAGATION

More information

Attribute Reduction and Information Granularity *

Attribute Reduction and Information Granularity * ttribute Reuction an nformation Granularity * Li-hong Wang School of Computer Engineering an Science, Shanghai University, Shanghai, 200072, P.R.C School of Computer Science an Technology, Yantai University,

More information

Lower Bounds for the Smoothed Number of Pareto optimal Solutions

Lower Bounds for the Smoothed Number of Pareto optimal Solutions Lower Bouns for the Smoothe Number of Pareto optimal Solutions Tobias Brunsch an Heiko Röglin Department of Computer Science, University of Bonn, Germany brunsch@cs.uni-bonn.e, heiko@roeglin.org Abstract.

More information

Least-Squares Regression on Sparse Spaces

Least-Squares Regression on Sparse Spaces Least-Squares Regression on Sparse Spaces Yuri Grinberg, Mahi Milani Far, Joelle Pineau School of Computer Science McGill University Montreal, Canaa {ygrinb,mmilan1,jpineau}@cs.mcgill.ca 1 Introuction

More information

Computing Exact Confidence Coefficients of Simultaneous Confidence Intervals for Multinomial Proportions and their Functions

Computing Exact Confidence Coefficients of Simultaneous Confidence Intervals for Multinomial Proportions and their Functions Working Paper 2013:5 Department of Statistics Computing Exact Confience Coefficients of Simultaneous Confience Intervals for Multinomial Proportions an their Functions Shaobo Jin Working Paper 2013:5

More information

Lecture 2 Lagrangian formulation of classical mechanics Mechanics

Lecture 2 Lagrangian formulation of classical mechanics Mechanics Lecture Lagrangian formulation of classical mechanics 70.00 Mechanics Principle of stationary action MATH-GA To specify a motion uniquely in classical mechanics, it suffices to give, at some time t 0,

More information

Optimization of Geometries by Energy Minimization

Optimization of Geometries by Energy Minimization Optimization of Geometries by Energy Minimization by Tracy P. Hamilton Department of Chemistry University of Alabama at Birmingham Birmingham, AL 3594-140 hamilton@uab.eu Copyright Tracy P. Hamilton, 1997.

More information

Chapter 6: Energy-Momentum Tensors

Chapter 6: Energy-Momentum Tensors 49 Chapter 6: Energy-Momentum Tensors This chapter outlines the general theory of energy an momentum conservation in terms of energy-momentum tensors, then applies these ieas to the case of Bohm's moel.

More information

05 The Continuum Limit and the Wave Equation

05 The Continuum Limit and the Wave Equation Utah State University DigitalCommons@USU Founations of Wave Phenomena Physics, Department of 1-1-2004 05 The Continuum Limit an the Wave Equation Charles G. Torre Department of Physics, Utah State University,

More information

Agmon Kolmogorov Inequalities on l 2 (Z d )

Agmon Kolmogorov Inequalities on l 2 (Z d ) Journal of Mathematics Research; Vol. 6, No. ; 04 ISSN 96-9795 E-ISSN 96-9809 Publishe by Canaian Center of Science an Eucation Agmon Kolmogorov Inequalities on l (Z ) Arman Sahovic Mathematics Department,

More information

Switching Time Optimization in Discretized Hybrid Dynamical Systems

Switching Time Optimization in Discretized Hybrid Dynamical Systems Switching Time Optimization in Discretize Hybri Dynamical Systems Kathrin Flaßkamp, To Murphey, an Sina Ober-Blöbaum Abstract Switching time optimization (STO) arises in systems that have a finite set

More information

Exponential Tracking Control of Nonlinear Systems with Actuator Nonlinearity

Exponential Tracking Control of Nonlinear Systems with Actuator Nonlinearity Preprints of the 9th Worl Congress The International Feeration of Automatic Control Cape Town, South Africa. August -9, Exponential Tracking Control of Nonlinear Systems with Actuator Nonlinearity Zhengqiang

More information

Lie symmetry and Mei conservation law of continuum system

Lie symmetry and Mei conservation law of continuum system Chin. Phys. B Vol. 20, No. 2 20 020 Lie symmetry an Mei conservation law of continuum system Shi Shen-Yang an Fu Jing-Li Department of Physics, Zhejiang Sci-Tech University, Hangzhou 3008, China Receive

More information

2Algebraic ONLINE PAGE PROOFS. foundations

2Algebraic ONLINE PAGE PROOFS. foundations Algebraic founations. Kick off with CAS. Algebraic skills.3 Pascal s triangle an binomial expansions.4 The binomial theorem.5 Sets of real numbers.6 Surs.7 Review . Kick off with CAS Playing lotto Using

More information

NOTES ON EULER-BOOLE SUMMATION (1) f (l 1) (n) f (l 1) (m) + ( 1)k 1 k! B k (y) f (k) (y) dy,

NOTES ON EULER-BOOLE SUMMATION (1) f (l 1) (n) f (l 1) (m) + ( 1)k 1 k! B k (y) f (k) (y) dy, NOTES ON EULER-BOOLE SUMMATION JONATHAN M BORWEIN, NEIL J CALKIN, AND DANTE MANNA Abstract We stuy a connection between Euler-MacLaurin Summation an Boole Summation suggeste in an AMM note from 196, which

More information

inflow outflow Part I. Regular tasks for MAE598/494 Task 1

inflow outflow Part I. Regular tasks for MAE598/494 Task 1 MAE 494/598, Fall 2016 Project #1 (Regular tasks = 20 points) Har copy of report is ue at the start of class on the ue ate. The rules on collaboration will be release separately. Please always follow the

More information

Math Notes on differentials, the Chain Rule, gradients, directional derivative, and normal vectors

Math Notes on differentials, the Chain Rule, gradients, directional derivative, and normal vectors Math 18.02 Notes on ifferentials, the Chain Rule, graients, irectional erivative, an normal vectors Tangent plane an linear approximation We efine the partial erivatives of f( xy, ) as follows: f f( x+

More information

Lectures - Week 10 Introduction to Ordinary Differential Equations (ODES) First Order Linear ODEs

Lectures - Week 10 Introduction to Ordinary Differential Equations (ODES) First Order Linear ODEs Lectures - Week 10 Introuction to Orinary Differential Equations (ODES) First Orer Linear ODEs When stuying ODEs we are consiering functions of one inepenent variable, e.g., f(x), where x is the inepenent

More information

ALGEBRAIC AND ANALYTIC PROPERTIES OF ARITHMETIC FUNCTIONS

ALGEBRAIC AND ANALYTIC PROPERTIES OF ARITHMETIC FUNCTIONS ALGEBRAIC AND ANALYTIC PROPERTIES OF ARITHMETIC FUNCTIONS MARK SCHACHNER Abstract. When consiere as an algebraic space, the set of arithmetic functions equippe with the operations of pointwise aition an

More information

Situation awareness of power system based on static voltage security region

Situation awareness of power system based on static voltage security region The 6th International Conference on Renewable Power Generation (RPG) 19 20 October 2017 Situation awareness of power system base on static voltage security region Fei Xiao, Zi-Qing Jiang, Qian Ai, Ran

More information

Thermal conductivity of graded composites: Numerical simulations and an effective medium approximation

Thermal conductivity of graded composites: Numerical simulations and an effective medium approximation JOURNAL OF MATERIALS SCIENCE 34 (999)5497 5503 Thermal conuctivity of grae composites: Numerical simulations an an effective meium approximation P. M. HUI Department of Physics, The Chinese University

More information

This module is part of the. Memobust Handbook. on Methodology of Modern Business Statistics

This module is part of the. Memobust Handbook. on Methodology of Modern Business Statistics This moule is part of the Memobust Hanbook on Methoology of Moern Business Statistics 26 March 2014 Metho: Balance Sampling for Multi-Way Stratification Contents General section... 3 1. Summary... 3 2.

More information

Schrödinger s equation.

Schrödinger s equation. Physics 342 Lecture 5 Schröinger s Equation Lecture 5 Physics 342 Quantum Mechanics I Wenesay, February 3r, 2010 Toay we iscuss Schröinger s equation an show that it supports the basic interpretation of

More information

Euler equations for multiple integrals

Euler equations for multiple integrals Euler equations for multiple integrals January 22, 2013 Contents 1 Reminer of multivariable calculus 2 1.1 Vector ifferentiation......................... 2 1.2 Matrix ifferentiation........................

More information

TEMPORAL AND TIME-FREQUENCY CORRELATION-BASED BLIND SOURCE SEPARATION METHODS. Yannick DEVILLE

TEMPORAL AND TIME-FREQUENCY CORRELATION-BASED BLIND SOURCE SEPARATION METHODS. Yannick DEVILLE TEMPORAL AND TIME-FREQUENCY CORRELATION-BASED BLIND SOURCE SEPARATION METHODS Yannick DEVILLE Université Paul Sabatier Laboratoire Acoustique, Métrologie, Instrumentation Bât. 3RB2, 8 Route e Narbonne,

More information

Make graph of g by adding c to the y-values. on the graph of f by c. multiplying the y-values. even-degree polynomial. graph goes up on both sides

Make graph of g by adding c to the y-values. on the graph of f by c. multiplying the y-values. even-degree polynomial. graph goes up on both sides Reference 1: Transformations of Graphs an En Behavior of Polynomial Graphs Transformations of graphs aitive constant constant on the outsie g(x) = + c Make graph of g by aing c to the y-values on the graph

More information

The Exact Form and General Integrating Factors

The Exact Form and General Integrating Factors 7 The Exact Form an General Integrating Factors In the previous chapters, we ve seen how separable an linear ifferential equations can be solve using methos for converting them to forms that can be easily

More information

The total derivative. Chapter Lagrangian and Eulerian approaches

The total derivative. Chapter Lagrangian and Eulerian approaches Chapter 5 The total erivative 51 Lagrangian an Eulerian approaches The representation of a flui through scalar or vector fiels means that each physical quantity uner consieration is escribe as a function

More information

A Novel Decoupled Iterative Method for Deep-Submicron MOSFET RF Circuit Simulation

A Novel Decoupled Iterative Method for Deep-Submicron MOSFET RF Circuit Simulation A Novel ecouple Iterative Metho for eep-submicron MOSFET RF Circuit Simulation CHUAN-SHENG WANG an YIMING LI epartment of Mathematics, National Tsing Hua University, National Nano evice Laboratories, an

More information

Separation of Variables

Separation of Variables Physics 342 Lecture 1 Separation of Variables Lecture 1 Physics 342 Quantum Mechanics I Monay, January 25th, 2010 There are three basic mathematical tools we nee, an then we can begin working on the physical

More information

TOEPLITZ AND POSITIVE SEMIDEFINITE COMPLETION PROBLEM FOR CYCLE GRAPH

TOEPLITZ AND POSITIVE SEMIDEFINITE COMPLETION PROBLEM FOR CYCLE GRAPH English NUMERICAL MATHEMATICS Vol14, No1 Series A Journal of Chinese Universities Feb 2005 TOEPLITZ AND POSITIVE SEMIDEFINITE COMPLETION PROBLEM FOR CYCLE GRAPH He Ming( Λ) Michael K Ng(Ξ ) Abstract We

More information

Linear First-Order Equations

Linear First-Order Equations 5 Linear First-Orer Equations Linear first-orer ifferential equations make up another important class of ifferential equations that commonly arise in applications an are relatively easy to solve (in theory)

More information

SYNCHRONOUS SEQUENTIAL CIRCUITS

SYNCHRONOUS SEQUENTIAL CIRCUITS CHAPTER SYNCHRONOUS SEUENTIAL CIRCUITS Registers an counters, two very common synchronous sequential circuits, are introuce in this chapter. Register is a igital circuit for storing information. Contents

More information

A Weak First Digit Law for a Class of Sequences

A Weak First Digit Law for a Class of Sequences International Mathematical Forum, Vol. 11, 2016, no. 15, 67-702 HIKARI Lt, www.m-hikari.com http://x.oi.org/10.1288/imf.2016.6562 A Weak First Digit Law for a Class of Sequences M. A. Nyblom School of

More information

LATTICE-BASED D-OPTIMUM DESIGN FOR FOURIER REGRESSION

LATTICE-BASED D-OPTIMUM DESIGN FOR FOURIER REGRESSION The Annals of Statistics 1997, Vol. 25, No. 6, 2313 2327 LATTICE-BASED D-OPTIMUM DESIGN FOR FOURIER REGRESSION By Eva Riccomagno, 1 Rainer Schwabe 2 an Henry P. Wynn 1 University of Warwick, Technische

More information

Classes of data dependence. Dependence analysis. Flow dependence (True dependence)

Classes of data dependence. Dependence analysis. Flow dependence (True dependence) Dependence analysis Classes of data dependence Pattern matching and replacement is all that is needed to apply many source-to-source transformations. For example, pattern matching can be used to determine

More information

Dependence analysis. However, it is often necessary to gather additional information to determine the correctness of a particular transformation.

Dependence analysis. However, it is often necessary to gather additional information to determine the correctness of a particular transformation. Dependence analysis Pattern matching and replacement is all that is needed to apply many source-to-source transformations. For example, pattern matching can be used to determine that the recursion removal

More information

Quantum Mechanics in Three Dimensions

Quantum Mechanics in Three Dimensions Physics 342 Lecture 20 Quantum Mechanics in Three Dimensions Lecture 20 Physics 342 Quantum Mechanics I Monay, March 24th, 2008 We begin our spherical solutions with the simplest possible case zero potential.

More information

Survey Sampling. 1 Design-based Inference. Kosuke Imai Department of Politics, Princeton University. February 19, 2013

Survey Sampling. 1 Design-based Inference. Kosuke Imai Department of Politics, Princeton University. February 19, 2013 Survey Sampling Kosuke Imai Department of Politics, Princeton University February 19, 2013 Survey sampling is one of the most commonly use ata collection methos for social scientists. We begin by escribing

More information

On colour-blind distinguishing colour pallets in regular graphs

On colour-blind distinguishing colour pallets in regular graphs J Comb Optim (2014 28:348 357 DOI 10.1007/s10878-012-9556-x On colour-blin istinguishing colour pallets in regular graphs Jakub Przybyło Publishe online: 25 October 2012 The Author(s 2012. This article

More information

Calculus of Variations

Calculus of Variations Calculus of Variations Lagrangian formalism is the main tool of theoretical classical mechanics. Calculus of Variations is a part of Mathematics which Lagrangian formalism is base on. In this section,

More information

Chapter 2 Lagrangian Modeling

Chapter 2 Lagrangian Modeling Chapter 2 Lagrangian Moeling The basic laws of physics are use to moel every system whether it is electrical, mechanical, hyraulic, or any other energy omain. In mechanics, Newton s laws of motion provie

More information

A Simple Model for the Calculation of Plasma Impedance in Atmospheric Radio Frequency Discharges

A Simple Model for the Calculation of Plasma Impedance in Atmospheric Radio Frequency Discharges Plasma Science an Technology, Vol.16, No.1, Oct. 214 A Simple Moel for the Calculation of Plasma Impeance in Atmospheric Raio Frequency Discharges GE Lei ( ) an ZHANG Yuantao ( ) Shanong Provincial Key

More information

Why Bernstein Polynomials Are Better: Fuzzy-Inspired Justification

Why Bernstein Polynomials Are Better: Fuzzy-Inspired Justification Why Bernstein Polynomials Are Better: Fuzzy-Inspire Justification Jaime Nava 1, Olga Kosheleva 2, an Vlaik Kreinovich 3 1,3 Department of Computer Science 2 Department of Teacher Eucation University of

More information

Analysis on a Localized Pruning Method for Connected Dominating Sets

Analysis on a Localized Pruning Method for Connected Dominating Sets JOURNAL OF INFORMATION SCIENCE AND ENGINEERING 23, 1073-1086 (2007) Analysis on a Localize Pruning Metho for Connecte Dominating Sets JOHN SUM 1, JIE WU 2 AND KEVIN HO 3 1 Department of Information Management

More information

Diagonalization of Matrices Dr. E. Jacobs

Diagonalization of Matrices Dr. E. Jacobs Diagonalization of Matrices Dr. E. Jacobs One of the very interesting lessons in this course is how certain algebraic techniques can be use to solve ifferential equations. The purpose of these notes is

More information

Topic Modeling: Beyond Bag-of-Words

Topic Modeling: Beyond Bag-of-Words Hanna M. Wallach Cavenish Laboratory, University of Cambrige, Cambrige CB3 0HE, UK hmw26@cam.ac.u Abstract Some moels of textual corpora employ text generation methos involving n-gram statistics, while

More information

Estimation of the Maximum Domination Value in Multi-Dimensional Data Sets

Estimation of the Maximum Domination Value in Multi-Dimensional Data Sets Proceeings of the 4th East-European Conference on Avances in Databases an Information Systems ADBIS) 200 Estimation of the Maximum Domination Value in Multi-Dimensional Data Sets Eleftherios Tiakas, Apostolos.

More information

Optimized Schwarz Methods with the Yin-Yang Grid for Shallow Water Equations

Optimized Schwarz Methods with the Yin-Yang Grid for Shallow Water Equations Optimize Schwarz Methos with the Yin-Yang Gri for Shallow Water Equations Abessama Qaouri Recherche en prévision numérique, Atmospheric Science an Technology Directorate, Environment Canaa, Dorval, Québec,

More information

Level Construction of Decision Trees in a Partition-based Framework for Classification

Level Construction of Decision Trees in a Partition-based Framework for Classification Level Construction of Decision Trees in a Partition-base Framework for Classification Y.Y. Yao, Y. Zhao an J.T. Yao Department of Computer Science, University of Regina Regina, Saskatchewan, Canaa S4S

More information

Integration Review. May 11, 2013

Integration Review. May 11, 2013 Integration Review May 11, 2013 Goals: Review the funamental theorem of calculus. Review u-substitution. Review integration by parts. Do lots of integration eamples. 1 Funamental Theorem of Calculus In

More information

Lower bounds on Locality Sensitive Hashing

Lower bounds on Locality Sensitive Hashing Lower bouns on Locality Sensitive Hashing Rajeev Motwani Assaf Naor Rina Panigrahy Abstract Given a metric space (X, X ), c 1, r > 0, an p, q [0, 1], a istribution over mappings H : X N is calle a (r,

More information

Final Exam Study Guide and Practice Problems Solutions

Final Exam Study Guide and Practice Problems Solutions Final Exam Stuy Guie an Practice Problems Solutions Note: These problems are just some of the types of problems that might appear on the exam. However, to fully prepare for the exam, in aition to making

More information

Hyperbolic Moment Equations Using Quadrature-Based Projection Methods

Hyperbolic Moment Equations Using Quadrature-Based Projection Methods Hyperbolic Moment Equations Using Quarature-Base Projection Methos J. Koellermeier an M. Torrilhon Department of Mathematics, RWTH Aachen University, Aachen, Germany Abstract. Kinetic equations like the

More information

The derivative of a function f(x) is another function, defined in terms of a limiting expression: f(x + δx) f(x)

The derivative of a function f(x) is another function, defined in terms of a limiting expression: f(x + δx) f(x) Y. D. Chong (2016) MH2801: Complex Methos for the Sciences 1. Derivatives The erivative of a function f(x) is another function, efine in terms of a limiting expression: f (x) f (x) lim x δx 0 f(x + δx)

More information

Construction of the Electronic Radial Wave Functions and Probability Distributions of Hydrogen-like Systems

Construction of the Electronic Radial Wave Functions and Probability Distributions of Hydrogen-like Systems Construction of the Electronic Raial Wave Functions an Probability Distributions of Hyrogen-like Systems Thomas S. Kuntzleman, Department of Chemistry Spring Arbor University, Spring Arbor MI 498 tkuntzle@arbor.eu

More information

Characterizing Real-Valued Multivariate Complex Polynomials and Their Symmetric Tensor Representations

Characterizing Real-Valued Multivariate Complex Polynomials and Their Symmetric Tensor Representations Characterizing Real-Value Multivariate Complex Polynomials an Their Symmetric Tensor Representations Bo JIANG Zhening LI Shuzhong ZHANG December 31, 2014 Abstract In this paper we stuy multivariate polynomial

More information

Similar Operators and a Functional Calculus for the First-Order Linear Differential Operator

Similar Operators and a Functional Calculus for the First-Order Linear Differential Operator Avances in Applie Mathematics, 9 47 999 Article ID aama.998.067, available online at http: www.iealibrary.com on Similar Operators an a Functional Calculus for the First-Orer Linear Differential Operator

More information

Perturbation Analysis and Optimization of Stochastic Flow Networks

Perturbation Analysis and Optimization of Stochastic Flow Networks IEEE TRANSACTIONS ON AUTOMATIC CONTROL, VOL. XX, NO. Y, MMM 2004 1 Perturbation Analysis an Optimization of Stochastic Flow Networks Gang Sun, Christos G. Cassanras, Yorai Wari, Christos G. Panayiotou,

More information

THE VAN KAMPEN EXPANSION FOR LINKED DUFFING LINEAR OSCILLATORS EXCITED BY COLORED NOISE

THE VAN KAMPEN EXPANSION FOR LINKED DUFFING LINEAR OSCILLATORS EXCITED BY COLORED NOISE Journal of Soun an Vibration (1996) 191(3), 397 414 THE VAN KAMPEN EXPANSION FOR LINKED DUFFING LINEAR OSCILLATORS EXCITED BY COLORED NOISE E. M. WEINSTEIN Galaxy Scientific Corporation, 2500 English Creek

More information

Closed and Open Loop Optimal Control of Buffer and Energy of a Wireless Device

Closed and Open Loop Optimal Control of Buffer and Energy of a Wireless Device Close an Open Loop Optimal Control of Buffer an Energy of a Wireless Device V. S. Borkar School of Technology an Computer Science TIFR, umbai, Inia. borkar@tifr.res.in A. A. Kherani B. J. Prabhu INRIA

More information

Fast image compression using matrix K-L transform

Fast image compression using matrix K-L transform Fast image compression using matrix K-L transform Daoqiang Zhang, Songcan Chen * Department of Computer Science an Engineering, Naning University of Aeronautics & Astronautics, Naning 2006, P.R. China.

More information

GLOBAL SOLUTIONS FOR 2D COUPLED BURGERS-COMPLEX-GINZBURG-LANDAU EQUATIONS

GLOBAL SOLUTIONS FOR 2D COUPLED BURGERS-COMPLEX-GINZBURG-LANDAU EQUATIONS Electronic Journal of Differential Equations, Vol. 015 015), No. 99, pp. 1 14. ISSN: 107-6691. URL: http://eje.math.txstate.eu or http://eje.math.unt.eu ftp eje.math.txstate.eu GLOBAL SOLUTIONS FOR D COUPLED

More information

arxiv: v4 [cs.ds] 7 Mar 2014

arxiv: v4 [cs.ds] 7 Mar 2014 Analysis of Agglomerative Clustering Marcel R. Ackermann Johannes Blömer Daniel Kuntze Christian Sohler arxiv:101.697v [cs.ds] 7 Mar 01 Abstract The iameter k-clustering problem is the problem of partitioning

More information

An inductance lookup table application for analysis of reluctance stepper motor model

An inductance lookup table application for analysis of reluctance stepper motor model ARCHIVES OF ELECTRICAL ENGINEERING VOL. 60(), pp. 5- (0) DOI 0.478/ v07-0-000-y An inuctance lookup table application for analysis of reluctance stepper motor moel JAKUB BERNAT, JAKUB KOŁOTA, SŁAWOMIR

More information

A Modification of the Jarque-Bera Test. for Normality

A Modification of the Jarque-Bera Test. for Normality Int. J. Contemp. Math. Sciences, Vol. 8, 01, no. 17, 84-85 HIKARI Lt, www.m-hikari.com http://x.oi.org/10.1988/ijcms.01.9106 A Moification of the Jarque-Bera Test for Normality Moawa El-Fallah Ab El-Salam

More information

Dependence Analysis. Dependence Examples. Last Time: Brief introduction to interprocedural analysis. do I = 2, 100 A(I) = A(I-1) + 1 enddo

Dependence Analysis. Dependence Examples. Last Time: Brief introduction to interprocedural analysis. do I = 2, 100 A(I) = A(I-1) + 1 enddo Dependence Analysis Dependence Examples Last Time: Brief introduction to interprocedural analysis Today: Optimization for parallel machines and memory hierarchies Dependence analysis Loop transformations

More information

Diophantine Approximations: Examining the Farey Process and its Method on Producing Best Approximations

Diophantine Approximations: Examining the Farey Process and its Method on Producing Best Approximations Diophantine Approximations: Examining the Farey Process an its Metho on Proucing Best Approximations Kelly Bowen Introuction When a person hears the phrase irrational number, one oes not think of anything

More information

6 General properties of an autonomous system of two first order ODE

6 General properties of an autonomous system of two first order ODE 6 General properties of an autonomous system of two first orer ODE Here we embark on stuying the autonomous system of two first orer ifferential equations of the form ẋ 1 = f 1 (, x 2 ), ẋ 2 = f 2 (, x

More information

Data Dependences and Parallelization. Stanford University CS243 Winter 2006 Wei Li 1

Data Dependences and Parallelization. Stanford University CS243 Winter 2006 Wei Li 1 Data Dependences and Parallelization Wei Li 1 Agenda Introduction Single Loop Nested Loops Data Dependence Analysis 2 Motivation DOALL loops: loops whose iterations can execute in parallel for i = 11,

More information

IERCU. Institute of Economic Research, Chuo University 50th Anniversary Special Issues. Discussion Paper No.210

IERCU. Institute of Economic Research, Chuo University 50th Anniversary Special Issues. Discussion Paper No.210 IERCU Institute of Economic Research, Chuo University 50th Anniversary Special Issues Discussion Paper No.210 Discrete an Continuous Dynamics in Nonlinear Monopolies Akio Matsumoto Chuo University Ferenc

More information

State-Space Model for a Multi-Machine System

State-Space Model for a Multi-Machine System State-Space Moel for a Multi-Machine System These notes parallel section.4 in the text. We are ealing with classically moele machines (IEEE Type.), constant impeance loas, an a network reuce to its internal

More information

THE ACCURATE ELEMENT METHOD: A NEW PARADIGM FOR NUMERICAL SOLUTION OF ORDINARY DIFFERENTIAL EQUATIONS

THE ACCURATE ELEMENT METHOD: A NEW PARADIGM FOR NUMERICAL SOLUTION OF ORDINARY DIFFERENTIAL EQUATIONS THE PUBISHING HOUSE PROCEEDINGS O THE ROMANIAN ACADEMY, Series A, O THE ROMANIAN ACADEMY Volume, Number /, pp. 6 THE ACCURATE EEMENT METHOD: A NEW PARADIGM OR NUMERICA SOUTION O ORDINARY DIERENTIA EQUATIONS

More information

A simplified macroscopic urban traffic network model for model-based predictive control

A simplified macroscopic urban traffic network model for model-based predictive control Delft University of Technology Delft Center for Systems an Control Technical report 9-28 A simplifie macroscopic urban traffic network moel for moel-base preictive control S. Lin, B. De Schutter, Y. Xi,

More information

Research Article When Inflation Causes No Increase in Claim Amounts

Research Article When Inflation Causes No Increase in Claim Amounts Probability an Statistics Volume 2009, Article ID 943926, 10 pages oi:10.1155/2009/943926 Research Article When Inflation Causes No Increase in Claim Amounts Vytaras Brazauskas, 1 Bruce L. Jones, 2 an

More information

II. First variation of functionals

II. First variation of functionals II. First variation of functionals The erivative of a function being zero is a necessary conition for the etremum of that function in orinary calculus. Let us now tackle the question of the equivalent

More information

Computing the Longest Common Subsequence of Multiple RLE Strings

Computing the Longest Common Subsequence of Multiple RLE Strings The 29th Workshop on Combinatorial Mathematics an Computation Theory Computing the Longest Common Subsequence of Multiple RLE Strings Ling-Chih Yao an Kuan-Yu Chen Grauate Institute of Networking an Multimeia

More information

Introduction to Markov Processes

Introduction to Markov Processes Introuction to Markov Processes Connexions moule m44014 Zzis law Gustav) Meglicki, Jr Office of the VP for Information Technology Iniana University RCS: Section-2.tex,v 1.24 2012/12/21 18:03:08 gustav

More information

Equilibrium in Queues Under Unknown Service Times and Service Value

Equilibrium in Queues Under Unknown Service Times and Service Value University of Pennsylvania ScholarlyCommons Finance Papers Wharton Faculty Research 1-2014 Equilibrium in Queues Uner Unknown Service Times an Service Value Laurens Debo Senthil K. Veeraraghavan University

More information

Robust Forward Algorithms via PAC-Bayes and Laplace Distributions. ω Q. Pr (y(ω x) < 0) = Pr A k

Robust Forward Algorithms via PAC-Bayes and Laplace Distributions. ω Q. Pr (y(ω x) < 0) = Pr A k A Proof of Lemma 2 B Proof of Lemma 3 Proof: Since the support of LL istributions is R, two such istributions are equivalent absolutely continuous with respect to each other an the ivergence is well-efine

More information

APPROXIMATE SOLUTION FOR TRANSIENT HEAT TRANSFER IN STATIC TURBULENT HE II. B. Baudouy. CEA/Saclay, DSM/DAPNIA/STCM Gif-sur-Yvette Cedex, France

APPROXIMATE SOLUTION FOR TRANSIENT HEAT TRANSFER IN STATIC TURBULENT HE II. B. Baudouy. CEA/Saclay, DSM/DAPNIA/STCM Gif-sur-Yvette Cedex, France APPROXIMAE SOLUION FOR RANSIEN HEA RANSFER IN SAIC URBULEN HE II B. Bauouy CEA/Saclay, DSM/DAPNIA/SCM 91191 Gif-sur-Yvette Ceex, France ABSRAC Analytical solution in one imension of the heat iffusion equation

More information

The Principle of Least Action

The Principle of Least Action Chapter 7. The Principle of Least Action 7.1 Force Methos vs. Energy Methos We have so far stuie two istinct ways of analyzing physics problems: force methos, basically consisting of the application of

More information

Exponential asymptotic property of a parallel repairable system with warm standby under common-cause failure

Exponential asymptotic property of a parallel repairable system with warm standby under common-cause failure J. Math. Anal. Appl. 341 (28) 457 466 www.elsevier.com/locate/jmaa Exponential asymptotic property of a parallel repairable system with warm stanby uner common-cause failure Zifei Shen, Xiaoxiao Hu, Weifeng

More information

Fractional Geometric Calculus: Toward A Unified Mathematical Language for Physics and Engineering

Fractional Geometric Calculus: Toward A Unified Mathematical Language for Physics and Engineering Fractional Geometric Calculus: Towar A Unifie Mathematical Language for Physics an Engineering Xiong Wang Center of Chaos an Complex Network, Department of Electronic Engineering, City University of Hong

More information

The Wiener Index of Trees with Prescribed Diameter

The Wiener Index of Trees with Prescribed Diameter 011 1 15 4 ± Dec., 011 Operations Research Transactions Vol.15 No.4 The Wiener Inex of Trees with Prescribe Diameter XING Baohua 1 CAI Gaixiang 1 Abstract The Wiener inex W(G) of a graph G is efine as

More information

Lecture Introduction. 2 Examples of Measure Concentration. 3 The Johnson-Lindenstrauss Lemma. CS-621 Theory Gems November 28, 2012

Lecture Introduction. 2 Examples of Measure Concentration. 3 The Johnson-Lindenstrauss Lemma. CS-621 Theory Gems November 28, 2012 CS-6 Theory Gems November 8, 0 Lecture Lecturer: Alesaner Mąry Scribes: Alhussein Fawzi, Dorina Thanou Introuction Toay, we will briefly iscuss an important technique in probability theory measure concentration

More information

On Characterizing the Delay-Performance of Wireless Scheduling Algorithms

On Characterizing the Delay-Performance of Wireless Scheduling Algorithms On Characterizing the Delay-Performance of Wireless Scheuling Algorithms Xiaojun Lin Center for Wireless Systems an Applications School of Electrical an Computer Engineering, Purue University West Lafayette,

More information

19 Eigenvalues, Eigenvectors, Ordinary Differential Equations, and Control

19 Eigenvalues, Eigenvectors, Ordinary Differential Equations, and Control 19 Eigenvalues, Eigenvectors, Orinary Differential Equations, an Control This section introuces eigenvalues an eigenvectors of a matrix, an iscusses the role of the eigenvalues in etermining the behavior

More information

3.7 Implicit Differentiation -- A Brief Introduction -- Student Notes

3.7 Implicit Differentiation -- A Brief Introduction -- Student Notes Fin these erivatives of these functions: y.7 Implicit Differentiation -- A Brief Introuction -- Stuent Notes tan y sin tan = sin y e = e = Write the inverses of these functions: y tan y sin How woul we

More information

JUST THE MATHS UNIT NUMBER DIFFERENTIATION 2 (Rates of change) A.J.Hobson

JUST THE MATHS UNIT NUMBER DIFFERENTIATION 2 (Rates of change) A.J.Hobson JUST THE MATHS UNIT NUMBER 10.2 DIFFERENTIATION 2 (Rates of change) by A.J.Hobson 10.2.1 Introuction 10.2.2 Average rates of change 10.2.3 Instantaneous rates of change 10.2.4 Derivatives 10.2.5 Exercises

More information

The maximum sustainable yield of Allee dynamic system

The maximum sustainable yield of Allee dynamic system Ecological Moelling 154 (2002) 1 7 www.elsevier.com/locate/ecolmoel The maximum sustainable yiel of Allee ynamic system Zhen-Shan Lin a, *, Bai-Lian Li b a Department of Geography, Nanjing Normal Uni ersity,

More information

Lecture XII. where Φ is called the potential function. Let us introduce spherical coordinates defined through the relations

Lecture XII. where Φ is called the potential function. Let us introduce spherical coordinates defined through the relations Lecture XII Abstract We introuce the Laplace equation in spherical coorinates an apply the metho of separation of variables to solve it. This will generate three linear orinary secon orer ifferential equations:

More information

Improving Estimation Accuracy in Nonrandomized Response Questioning Methods by Multiple Answers

Improving Estimation Accuracy in Nonrandomized Response Questioning Methods by Multiple Answers International Journal of Statistics an Probability; Vol 6, No 5; September 207 ISSN 927-7032 E-ISSN 927-7040 Publishe by Canaian Center of Science an Eucation Improving Estimation Accuracy in Nonranomize

More information

Nonlinear Dielectric Response of Periodic Composite Materials

Nonlinear Dielectric Response of Periodic Composite Materials onlinear Dielectric Response of Perioic Composite aterials A.G. KOLPAKOV 3, Bl.95, 9 th ovember str., ovosibirsk, 639 Russia the corresponing author e-mail: agk@neic.nsk.su, algk@ngs.ru A. K.TAGATSEV Ceramics

More information

IPA Derivatives for Make-to-Stock Production-Inventory Systems With Backorders Under the (R,r) Policy

IPA Derivatives for Make-to-Stock Production-Inventory Systems With Backorders Under the (R,r) Policy IPA Derivatives for Make-to-Stock Prouction-Inventory Systems With Backorers Uner the (Rr) Policy Yihong Fan a Benamin Melame b Yao Zhao c Yorai Wari Abstract This paper aresses Infinitesimal Perturbation

More information

Transmission Line Matrix (TLM) network analogues of reversible trapping processes Part B: scaling and consistency

Transmission Line Matrix (TLM) network analogues of reversible trapping processes Part B: scaling and consistency Transmission Line Matrix (TLM network analogues of reversible trapping processes Part B: scaling an consistency Donar e Cogan * ANC Eucation, 308-310.A. De Mel Mawatha, Colombo 3, Sri Lanka * onarecogan@gmail.com

More information

The effect of dissipation on solutions of the complex KdV equation

The effect of dissipation on solutions of the complex KdV equation Mathematics an Computers in Simulation 69 (25) 589 599 The effect of issipation on solutions of the complex KV equation Jiahong Wu a,, Juan-Ming Yuan a,b a Department of Mathematics, Oklahoma State University,

More information

d-dimensional Arrangement Revisited

d-dimensional Arrangement Revisited -Dimensional Arrangement Revisite Daniel Rotter Jens Vygen Research Institute for Discrete Mathematics University of Bonn Revise version: April 5, 013 Abstract We revisit the -imensional arrangement problem

More information