(t, m, s)-nets and Maximized Minimum Distance, Part II
|
|
- Samuel Greene
- 5 years ago
- Views:
Transcription
1 (t,, s)-nets and Maxiized Miniu Distance, Part II Leonhard Grünschloß and Alexander Keller Abstract The quality paraeter t of (t,, s)-nets controls extensive stratification properties of the generated saple points. However, the definition allows for points that are arbitrarily close across strata boundaries. We continue the investigation of (t,, s)-nets under the constraint of axiizing the utual distance of the points on the unit torus and present two new constructions along with algoriths. The first approach is based on the fact that reordering (t, s)-sequences can result in (t,, s + 1)-nets with varying toroidal distance, while the second algorith generates points by perutations instead of atrices. 1 Introduction Iportant probles in iage synthesis like e.g. anti-aliasing, heispherical integration, or the illuination by area light sources can be considered as low-diensional nuerical integration probles. Aong the ost successful approaches to coputing this kind of integrals are quasi-monte Carlo and randoized quasi-monte Carlo ethods [9, 13, ] based on the two-diensional Larcher-Pillichshaer points [10], which expose a coparatively large iniu toroidal distance. These points belong to the class of (t,, s)-nets in base q (for an extensive overview of the topic we refer to [1, Ch. 4, especially p. 48]), which is given by Leonhard Grünschloß ental iages GbH, Fasanenstraße 81, 1063 Berlin, Gerany. eail: leonhard@gruenschloss.org Alexander Keller ental iages GbH, Fasanenstraße 81, 1063 Berlin, Gerany. eail: alex@ental.co 1
2 L. Grünschloß, A. Keller Definition 1. For integers 0 t, a (t,, s)-net in base q is a point set of q points in [0, 1) s such that there are exactly q t points in each eleentary interval E with volue q t. Figure 1 shows an instance of a (0, 4, )-net in base and illustrates the concept of eleentary intervals E = s [ ai q bi, (a i + 1)q bi) [0, 1) s i=1 as used in the definition, where a i, b i Z, b i 0 and 0 a i < q bi. The concept of (t,, s)-nets can be generalized for sequences of points as given by Definition. For an integer t 0, a sequence x 0, x 1,... of points in [0, 1) s is a (t, s)-sequence in base q if, for all integers k 0 and > t, the point set x kq,..., x (k+1)q 1 is a (t,, s)-net in base q. Obviously, the stratification properties of (t,, s)-nets and (t, s)-sequences are best for the quality paraeter t = 0, because then every eleentary interval contains exactly one point, as illustrated in Figure 1. The conception does not consider the utual distance of the points, which allows points to lie arbitrarily close together across shared interval boundaries. In this paper we continue previous work fro [6] that used the iniu toroidal distance d in ({x 0,..., x N 1 ) := in 0 u<v<n x u x v T to classify point sets {x 0,..., x N 1, where the toroidal distance of two points x = (x 1,..., x s ) [0, 1) s and y = (y 1,..., y s ) [0, 1) s is defined as x y T := s (in{ x i y i, 1 x i y i ). i=1 Maxiizing the shift-invariant easure of toroidal distance further increases unifority, allows one to tile the resulting point sets, and to consider periodic integrands, which is especially useful in the aforeentioned graphics applications. In Section we therefore extend the construction of Larcher and Pillichshaer to s = 3 diensions resulting in a (0,, 3)-net in base, which exhibits a large iniu toroidal distance. In Section 3 we then present a new perutation-based construction for (0,, )-nets in base, which often have a larger iniu toroidal distance than can be obtained by any digital net in base.
3 (t,, s)-nets and Maxiized Miniu Distance, Part II 3 Fig. 1 The Larcher-Pillichshaer points as an exaple of a (0, 4, )-net in base superiposed on all possible eleentary intervals. There is exactly one point in each eleentary interval. A new (0,, 3)-Net in Base with Large Miniu Toroidal Distance Although both the Haersley and Larcher-Pillichshaer points are constructions of (0,, )-nets in base, the latter have a uch larger iniu toroidal distance [8]. In fact already in [10, 11] the Haersley points have been identified to be the worst construction with respect to certain other easures, too. The extension of the Larcher-Pillichshaer points to s = 3 diensions has been an open proble for long. In the following, we present a construction of a (0,, 3)-net in base, where one two-diensional projection equals the Larcher-Pillichshaer points [9] and therefore benefits fro their large iniu toroidal distance..1 Digital Nets and Sequences We briefly suarize necessary notation and algorithic facts on digital nets and sequences fro [1] in a siplified anner. While in general digital nets and sequences are defined using a finite field F q, with q being a prie power, we only consider the case q = in the following. Given suitable -atrices C 1,..., C s over F, the i-th coponent of the n-th point for 1 i s and 0 n < can be generated by x (i) n = 1. d 0 (n) C i. [0, 1), d 1 (n) T where the atrix-vector ultiplication has to be perfored in F and the digits d k (n) are defined by the binary expansion of n = 1 k=0 d k (n) k.
4 4 L. Grünschloß, A. Keller The eleents of the i-th generator atrix are denoted by c (i) j,r. The theoretical construction of (t, s)-sequences requires generator atrices of infinite size. However, in practice this does not pose a proble when enuerating the points, since d k (n) = 0 for all sufficiently large k and, in addition, we require the atrix entries c (i) j,r = 0 for all sufficiently large j [1, p. 7, (S6)]. Therefore, only finite upper left subatrices have to be considered when generating points. Finally, due to the finite precision of integer coputation, is finite in any case. Based on the results in [10, 11, 4], we will consider the (0, 1)-sequence in base generated by the non-singular upper triangular atrix ({ ) C 1 1 if j r, := (1) 0 otherwise j,r=0 in order to generate the Larcher-Pillichshaer points [9].. Reordering the Sobol -Sequence We now take a look at the first two coponents of the Sobol -sequence [15], which are a (0, )-sequence in base, and are defined by the infinite upper triangular atrices (( ) ) C 1 := (δ j,r ) r j,r=0 and C := od, j j,r=0 where δ j,r = 1 for j = r and zero otherwise. In the above definition and the reainder of this work, by a od we ean the coon residue. That is the nonnegative value b <, such that a b od. Suppose that the infinite atrices C 1,..., C s over F generate a (0, s)- sequence in base. Furtherore suppose that the infinite atrix C 1 over F generates a (0, 1)-sequence in base. This iplies that C 1 is a nonsingular upper triangular atrix. Then a (0, s)-sequence in base is generated by C 1, C D,..., C s D, where D := C1 1 C 1. This new sequence consists of the sae points as before, however, in a different order. This result is proven in ore general for in [5, Prop. 1]. Note that the first generator atrix C 1 is the one of Larcher-Pillichshaer. Since C 1 is the identity, we have D = C1 1 C 1 = C 1 and obtain
5 (t,, s)-nets and Maxiized Miniu Distance, Part II 5 (( ) ) r C D = od C 1 j j,r=0 ( r ( ) ) k = od j k=0 j,r=0 (( ) ) r + 1 = od =: C j + 1, () j,r=0 as the second generator atrix, where the last equality follows fro the Christas Stocking Theore. Note that the sae reordering can be applied to the Faure-sequence [3], which can be regarded as a generalization of the Sobol -sequence for any prie base q. An even ore general construction for any prie power base can be found in [1]..3 Construction n Cobining the coponent with the first points of a (0, )-sequence in base yields a (0,, 3)-net in base (see [1, Lea 4., p. 6]). The relationship of C 1 and C as expressed in Equation () in fact is a property of all (0, )-sequences in base generated by non-singular upper triangular atrices [7, Prop. 4]. Consequently fixing the first generator atrix uniquely deterines the second generator atrix and thus the (0,, 3)-net. In the previous section we reordered the Sobol -sequence in base such that one coponent atches the (0, 1)-sequence as defined in [10]. For the construction of the new (0,, 3)-net in base, it is sufficient to consider the atrices C 1, C fro Section., because for s = b = Sobol s construction is identical to Faure s construction. By construction two diensions of the resulting three-diensional net are the Larcher-Pillichshaer points [9]. In Figure we copared the iniu toroidal distance of the new construction to the one that uses the original Sobol generator atrices C 1 and C. Except for = 5 and = 6, the new construction is by far superior for all. We did not copute the iniu toroidal distance for >, though..3.1 Ipleentation The following code in C99 (the current ANSI standard of the C language) returns the n-th point of the (0,, 3)-net generated using C 1 and C (see Equations (1) and ()) for < 3 in O(). The vectorized ipleenta-
6 6 L. Grünschloß, A. Keller 1 (0,,3)-net fro Sobol -sequence (0,,3)-net fro LP-atched sequence 0.1 iniu toroidal distance e-05 1e Fig. A plot of the toroidal distance in [0, 1) 3 for both the Sobol (0,, 3)-net as well as the (0,, 3)-net constructed using C 1, C, where one of the two-diensional projections equals the Larcher-Pillichshaer point set. We would like to stress that the iniu distance is scaled logarithically, so the absolute difference is very significant for increasing. In contrast to the Sobol -net, the iniu distance decreases soothly for the new construction. tion is based on the fact that addition in F corresponds to the exclusive or operation. void x_n(unsigned int n, const unsigned int, float x[3]) { // first coponent: n / ^ x[0] = (float) n / (1U << ); // reaining coponents by atrix ultiplication unsigned int r1 = 0, r = 0; for (unsigned int v1 = 1U << 31, v = 3U << 30; n; n >>= 1) { if (n & 1) { // vector addition of atrix colun by XOR r1 ^= v1; r ^= v << 1; // update atrix coluns v1 = v1 >> 1; v ^= v >> 1; // ap to unit cube [0,1)^3 x[1] = r1 * (1.f / (1ULL << 3)); x[] = r * (1.f / (1ULL << 3));
7 (t,, s)-nets and Maxiized Miniu Distance, Part II 7 3 Perutation-Generated (0,, )-Nets in Base with Larger Miniu Toroidal Distance Taking the integer part of the coordinates of a (0,, s)-net in base q ultiplied by the nuber of points q results in each coponent being a perutation fro the syetric group S q. For s = diensions and base q = this is visualized in Figure 1, where each colun and each row (the leftost and rightost set of eleentary intervals in the figure) contain exactly one point. Using Heap s efficient perutation generation algorith [14] allows one to enuerate all perutations π : {0,..., 1 {0,..., 1 that 1 represent (t,, )-nets in base with the points (n, π(n)) [0, 1). We extended this basic backtracking algorith by pruning the search tree whenever the eleentary-interval property t = 0 was violated or already a net with larger iniu toroidal distance than the current one had been found. Verifying the t = 0 property is siple using the code published in [6, Section.]. Perforing such an exhaustive cobinatorial coputer search sees hopeless given the nuber of possible perutations S = ( )! and in fact we did not succeed in running an exhaustive search for > 5. In these cases we siply generated rando perutations, but kept the backtracking approach entioned above. However in [6], we were able to find very regular looking (0, 5, )-nets with a iniu toroidal distance of 3/ , which is larger than the iniu toroidal distance of all possible digital (0, 5, )-nets, where the axiu is 9/ In continuation of these findings we now present a first construction (see Figure 3) of such perutation-based nets. In Table 1 we copare the iniu toroidal distance of different (0,, )-nets in base. The new perutation construction clearly features the largest iniu toroidal distance. Due to the novelty of the approach, we present two different derivations and provide ore interpretations than usually necessary. 3.1 Iterative Construction by Quadrupling Point Sets Given the search results as displayed in Figure 3, an iterative construction procedure can be inferred. This procedure repeatedly quadruples an initial point set until the desired nuber of points is reached. We therefore need to distinguish the cases of even and odd as depicted by the exaples in Figure 4. We represent the (0, k, )-nets in base by
8 8 L. Grünschloß, A. Keller = 3 = 5 = 7 = 4 = 6 = 8 Fig. 3 Exaples of the perutation-generated (0,, )-nets in base for = 3,..., 8. Haersley Larcher-Pillichsh. Opt. atrices Perutation Table 1 Coparison of the Haersley points, the Larcher-Pillichshaer points, the points resulting fro the optiized atrices given in [6], and the new perutation construction with respect to iniu toroidal distance. Note that for easier coparison all distances have been ultiplied by the nuber of points of a (0,, )-net in base and squared afterwards. { 1 k (u od k, v) : (u, v) P k, (3) where the set P k contains k integer coordinates (u, v).
9 (t,, s)-nets and Maxiized Miniu Distance, Part II 9 For odd the initial point set P 1 consists of the points (0, 0) and (1, 1). The quadrupling rule to construct a point set P k+ := P k,0 P k,1 P k, P k,3 with k+ points fro a point set P k with P k = k points is as follows: Lower left: P k,0 := {(u, v) : (u, v) P k Lower right: P k,1 := { ( k+1 + u + 1, v + 1) : (u, v) P k Upper left: P k, := { ( k+1 + u, k+1 + v) : (u, v) P k Upper right: P k,3 := { ( k+ + u + 1, k+1 + v + 1) : (u, v) P k For even, the iterative construction is ore involved. We start out with four diagonals, each consisting of four points: P,0 := {(w + 1, 4w) : w = 0,..., 3, P,1 := {(w + 5, 4w + ) : w = 0,..., 3, P, := {(w + 9, 4w + 1) : w = 0,..., 3, P,3 := {(w + 13, 4w + 3) : w = 0,..., 3. The quadrupling rule for even k constructs P k+ := P k,0 P k,1 P k, P fro sets P k,l for l = 0,..., 3 as follows: 1. Multiply all points by two and add an offset for the points in P 1 and P 3 : P k,l := {(u + (l od ), v) : (u, v) P k,l for l = 0,..., 3.. Extend the diagonals: P k,l := P k,l { ( k 1 + u, k+1 + v) : (u, v) P k,l for l = 0,..., Cobine diagonals: 4. Append shifted copies: P k,0 := P k,0 P k,1, P k,1 := P k, P k,3. P k, := { ( k+1 + u, u + 1) : (u, v) P k,0, P k,3 := { ( k+1 + v, v + 1) : (u, v) P k,1. k,3 3. Direct Construction As these nets consist of points lying on parallel diagonals siilar to rank-1 lattices [], a direct construction is possible by coputing a single point on
10 10 L. Grünschloß, A. Keller a) Exaple for odd : Construction of a (0, 9, )-net in base. b) Exaple for even : Construction of a (0, 8, )-net in base. Fig. 4 Iterative construction of (0,, )-nets in base by quadrupling point sets, where the final net results fro wrapping the points such that they atch the unit square. This wrapping corresponds to the odulo operation in Equation (3). each diagonal. Since the points are evenly spaced along the diagonals, the reaining points follow iediately. These diagonals are to be understood with a odulo wrap-around for the first coordinate Construction for Odd The position of the point with the sallest second coordinate on the d-th diagonal is given by ( φ (d) + d, d), where 0 d < and φ (d) [0, 1) denotes the van der Corput radical inverse in base. Theore 1. For odd, the points { (u d,w, v d,w ) : 0 d <, 0 w <, where u d,w = ( φ (d) + d + w ) od, (4) v d,w = d + w, (5) constitute a (0,, )-net in base with a iniu toroidal distance of. This distance is easured using coponents ultiplied by, i.e. on integer scale.
11 (t,, s)-nets and Maxiized Miniu Distance, Part II 11 Case 1: width > height. Case : width < height. Fig. 5 The two cases in the proof corresponding to certain kinds of eleentary intervals, illustrated here for = 5. Proof. First, we will show that the eleentary interval property holds. For each kind 0 h of eleentary intervals, there ust be exactly one point in each integer-scaled eleentary interval [ x h, (x + 1) h) [ y h, (y + 1) h), where 0 x < h and 0 y < h. Partitioning the set of eleentary intervals as depicted in Figure 5, we need to consider two cases: 1. We first consider the case h. The width of these kinds of eleentary intervals is larger than their height. Looking at Equation (5), we can see that the least significant bits of vd,w are equal to the bits of d, while the reaining ost significant bits of vd,w are equal to the bits of w. Considering one horizontal strip of eleentary intervals, y deterines the h ost significant bits of v d,w. Since h, w is copletely deterined by y. What reains to show is that for this fixed w, the h ost significant bits of u d,w differ for suitable values for d. Analyzing Equation (4), we see that the ter w is constant, while the addition of d only odifies the least significant bits. However { φ (d) : 0 d < = { d : 0 d <, (6) thus due to the addition of φ (d) it is possible to guarantee that the h ost significant bits of ud,w are different by selecting suitable values for d. Thus the corresponding points fall into different eleentary intervals.. Now we consider the case 0 h. The width of these kinds of eleentary intervals is saller than their height. We consider one vertical strip of eleentary intervals of width h. Looking at Equation (4), we can see that the the only way to achieve consecutive values u d,w equal to x, x + 1,..., x + h 1 is by using consecutive values for d. That is because the ters φ (d) and w do not odify the h least significant bits of u d,w. However h consecutive values for d ean that the h ost significant bits of the values φ (d) are different for each d (cf. to Equation (6)). In order to stay in the vertical strip of eleentary intervals deterined by x, the h ost significant bits of w thus ust be
12 1 L. Grünschloß, A. Keller chosen accordingly for each point inside this strip in order to copensate using the ter w. As a consequence of the different values for w it follows fro Equation (5) that the h ost significant bits of vd,w are different for each point inside this strip, thus they fall into different eleentary intervals. We now consider the achieved iniu toroidal distance of( the point set. Points on the diagonals are placed with a ultiple of the offset, ), so their squared iniu toroidal distance to each other is ( ) =. The squared distance between the diagonals with slope 1 is ( ) =. In conclusion, the iniu toroidal distance is. As the diagonals can be tiled sealessly, the result is identical for the toroidal distance easure. On the unit square [0, 1), we need to divide the point coordinates by, so the iniu toroidal distance on the unit square is =, which concludes the proof. The theore allows for an additional interpretation of the structure: Using Equation (4), we can solve for d and w given u d,w : d = u d,w od, w = (u d,w φ (d) d) od. Inserting these equations into Equation (5) yields ( v d,w = d + (u d,w φ u d,w od ) ) d od ( = (u d,w φ u d,w od )) od, (7) which can be regarded as replicating a (0,, ) integer Haersley-net ties horizontally, scaling each one vertically by and finally adding a linear coponent to the cobination of nets. The resulting points are wrapped around the unit square.
13 (t,, s)-nets and Maxiized Miniu Distance, Part II 13 Fig. 6 On the left, the odified tiling for the perutation construction is shown: is even and each row of the tiling is shifted by (, 0) relative to its adjacent row below. Using an axis aligned tiling on the right reveals that all tiles in a row share the sae saple pattern, while the pattern differs by row by a odulo wrap-around in direction of the x-axis. 3.. Construction for Even For even 6 the perutation search did not reveal an iproveent of the iniu toroidal distance in coparison to the full atrix search. We would like to note that the perutation search for = 6 did not finish, though. However, when allowing a distance easure that accounts for a slightly irregular tiling, an approach very uch resebling the previous construction can be taken, yielding the nets shown in Figure 3. Using a tiling where each row is shifted by (, 0) relative to its adjacent row below (see Figure 6), we can get a sealess tiling using diagonals, each with points. Again, the points are spaced evenly on each diagonal, this tie by the offset vector (, ). The position of the point with the sallest second coordinate on the d-th diagonal, where d { 0,..., 1, is given by ( φ (d) + d 4 +, d). With the odified tiling described above, these diagonals are continued sealessly. These nets cannot be generated via the classical way of using generator atrices as the point (0, 0) is not included. While we believe that a siilar proof that t = 0 is possible as for the odd case, we only verified the t = 0 property using a coputer progra for all even. The iniu toroidal distance equals 41 for = 8 and 41 8 / for all even where 10. See Table 1 for iniu toroidal distance values for < 8. We would like to stress that the odified iniu toroidal distance easure that respects the shifts of the rows in the tiling has been used for these easureents. Considering a superiposed axis-aligned tiling as illustrated in Figure 6 reveals that in a row each tile contains the sae set of saples, while the set of saples varies by row. This in turn allows for constructing a larger pattern that consists of shifted copies of the original pattern as shown in Figure 7. Four rows and four coluns of such odulo-wrapped patterns were cobined
14 14 L. Grünschloß, A. Keller Fig. 7 By taking 4 4 odulo-wrapped copies of the sae pattern (see Figure 6), we obtain a larger pattern that is periodic and atches the pixel arrangeent. It can be tiled regularly without a decrease in iniu toroidal distance. to generate a larger pattern. This pattern can be tiled regularly without a decrease in iniu toroidal distance. 3.3 Ipleentation for General Exploiting the fact that we are generating nets in base, all operations needed to ipleent the perutation net constructions described above can be ipleented very efficiently without any ultiplications. A sall C++ class that generates the point sets directly for even and odd 4 is available at In ters of perforance, it is coparable to the generation of rank-1 lattice points and also includes code to generate points with the odified tiling explained above. As an exaple, we would like to show the ipleentation of Equation (7). Given the function phi to copute the van der Corput radical inverse in base [9], the ipleentation y has a one line body. inline unsigned int phi(unsigned int bits) { // for 3 bits bits = (bits << 16) (bits >> 16); bits = ((bits & 0x00ff00ff) << 8) ((bits & 0xff00ff00) >> 8); bits = ((bits & 0x0f0f0f0f) << 4) ((bits & 0xf0f0f0f0) >> 4); bits = ((bits & 0x ) << ) ((bits & 0xcccccccc) >> ); bits = ((bits & 0x ) << 1) ((bits & 0xaaaaaaaa) >> 1); return bits; inline unsigned int y(const unsigned int x, const unsigned int ) { return (x - (phi(x & ~-(1 << ( >> 1))) >> (3 - ))) & ~-(1U << ); For successive point requests this can be optiized by precoputing the bitasks and the different radical inverse values.
15 (t,, s)-nets and Maxiized Miniu Distance, Part II 15 4 Conclusion We constructed new (0,, s)-nets in base for s =, 3, which are constrained by axiizing the iniu toroidal distance. Especially the (0,, 3)-net has any applications in coputer graphics such as coputing anti-aliasing and otion blur in the Reyes architecture [1]. Acknowledgeents We would like to thank the reviewers and editors for their very helpful coents. We would also like to thank Sabrina Daertz, Johannes Hanika, Sehera Nawaz, Matthias Raab, and Daniel Seibert for helpful discussions. References 1. Cook, R., Carpenter, L., Catull, E.: The Reyes iage rendering architecture. SIG- GRAPH Coputer Graphics 1(4), (1987). Daertz, S., Keller, A.: Iage synthesis by rank-1 lattices. In: A. Keller, S. Heinrich, H. Niederreiter (eds.) Monte Carlo and Quasi-Monte Carlo Methods 006, pp Springer (008) 3. Faure, H.: Discrépance de suites associées à un systèe de nuération (en diension s). Acta Arithetica 41(4), (198) 4. Faure, H.: Discrepancy and diaphony of digital (0, 1)-sequences in prie base. Acta Arithetica 117, (005) 5. Faure, H., Tezuka, S.: Another rando scrabling of digital (t, s)-sequences. In: K. Fang, F. Hickernell, H. Niederreiter (eds.) Monte Carlo and Quasi-Monte Carlo Methods 000, pp Springer-Verlag, Berlin (00) 6. Grünschloß, L., Hanika, J., Schwede, R., Keller, A.: (t,, s)-nets and axiized iniu distance. In: A. Keller, S. Heinrich, H. Niederreiter (eds.) Monte Carlo and Quasi-Monte Carlo Methods 006, pp Springer (008) 7. Hofer, R., Larcher, G.: On existence and discrepancy of certain digital Niederreiter- Halton sequences (009). To appear in Acta Arithetica 8. Keller, A.: Myths of coputer graphics. In: H. Niederreiter (ed.) Monte Carlo and Quasi-Monte Carlo Methods 004, pp Springer (006) 9. Kollig, T., Keller, A.: Efficient ultidiensional sapling. Coputer Graphics Foru 1(3), (00) 10. Larcher, G., Pillichshaer, F.: Walsh series analysis of the L -discrepancy of syetrisized point sets. Monatsheft Matheatik 13, 1 18 (001) 11. Larcher, G., Pillichshaer, F.: Sus of distances to the nearest integer and the discrepancy of digital nets. Acta Arithetica 106, (003) 1. Niederreiter, H.: Rando Nuber Generation and Quasi-Monte Carlo Methods. SIAM, Philadelphia (199) 13. Pharr, M., Huphreys, G.: Physically Based Rendering: Fro Theory to Ipleentation. Morgan Kaufann Publishers Inc., San Francisco, CA, USA (004) 14. Sedgewick, R.: Perutation generation ethods. ACM Coputing Surveys 9(), (1977) 15. Sobol, I.: On the distribution of points in a cube and the approxiate evaluation of integrals. Zhurnal Vychislitelnoi Mateatiki i Mateaticheskoi Fiziki 7(4), (1967)
Necessity of low effective dimension
Necessity of low effective diension Art B. Owen Stanford University October 2002, Orig: July 2002 Abstract Practitioners have long noticed that quasi-monte Carlo ethods work very well on functions that
More informationFeature Extraction Techniques
Feature Extraction Techniques Unsupervised Learning II Feature Extraction Unsupervised ethods can also be used to find features which can be useful for categorization. There are unsupervised ethods that
More informationPolygonal Designs: Existence and Construction
Polygonal Designs: Existence and Construction John Hegean Departent of Matheatics, Stanford University, Stanford, CA 9405 Jeff Langford Departent of Matheatics, Drake University, Des Moines, IA 5011 G
More informationChapter 6 1-D Continuous Groups
Chapter 6 1-D Continuous Groups Continuous groups consist of group eleents labelled by one or ore continuous variables, say a 1, a 2,, a r, where each variable has a well- defined range. This chapter explores:
More informationBlock designs and statistics
Bloc designs and statistics Notes for Math 447 May 3, 2011 The ain paraeters of a bloc design are nuber of varieties v, bloc size, nuber of blocs b. A design is built on a set of v eleents. Each eleent
More informationBipartite subgraphs and the smallest eigenvalue
Bipartite subgraphs and the sallest eigenvalue Noga Alon Benny Sudaov Abstract Two results dealing with the relation between the sallest eigenvalue of a graph and its bipartite subgraphs are obtained.
More informationRandomized Recovery for Boolean Compressed Sensing
Randoized Recovery for Boolean Copressed Sensing Mitra Fatei and Martin Vetterli Laboratory of Audiovisual Counication École Polytechnique Fédéral de Lausanne (EPFL) Eail: {itra.fatei, artin.vetterli}@epfl.ch
More informationA1. Find all ordered pairs (a, b) of positive integers for which 1 a + 1 b = 3
A. Find all ordered pairs a, b) of positive integers for which a + b = 3 08. Answer. The six ordered pairs are 009, 08), 08, 009), 009 337, 674) = 35043, 674), 009 346, 673) = 3584, 673), 674, 009 337)
More information13.2 Fully Polynomial Randomized Approximation Scheme for Permanent of Random 0-1 Matrices
CS71 Randoness & Coputation Spring 018 Instructor: Alistair Sinclair Lecture 13: February 7 Disclaier: These notes have not been subjected to the usual scrutiny accorded to foral publications. They ay
More informationIntelligent Systems: Reasoning and Recognition. Perceptrons and Support Vector Machines
Intelligent Systes: Reasoning and Recognition Jaes L. Crowley osig 1 Winter Seester 2018 Lesson 6 27 February 2018 Outline Perceptrons and Support Vector achines Notation...2 Linear odels...3 Lines, Planes
More informationFast Montgomery-like Square Root Computation over GF(2 m ) for All Trinomials
Fast Montgoery-like Square Root Coputation over GF( ) for All Trinoials Yin Li a, Yu Zhang a, a Departent of Coputer Science and Technology, Xinyang Noral University, Henan, P.R.China Abstract This letter
More informationReed-Muller Codes. m r inductive definition. Later, we shall explain how to construct Reed-Muller codes using the Kronecker product.
Coding Theory Massoud Malek Reed-Muller Codes An iportant class of linear block codes rich in algebraic and geoetric structure is the class of Reed-Muller codes, which includes the Extended Haing code.
More informationarxiv: v1 [cs.ds] 17 Mar 2016
Tight Bounds for Single-Pass Streaing Coplexity of the Set Cover Proble Sepehr Assadi Sanjeev Khanna Yang Li Abstract arxiv:1603.05715v1 [cs.ds] 17 Mar 2016 We resolve the space coplexity of single-pass
More informationCh 12: Variations on Backpropagation
Ch 2: Variations on Backpropagation The basic backpropagation algorith is too slow for ost practical applications. It ay take days or weeks of coputer tie. We deonstrate why the backpropagation algorith
More informationarxiv: v1 [cs.ds] 3 Feb 2014
arxiv:40.043v [cs.ds] 3 Feb 04 A Bound on the Expected Optiality of Rando Feasible Solutions to Cobinatorial Optiization Probles Evan A. Sultani The Johns Hopins University APL evan@sultani.co http://www.sultani.co/
More informationA note on the multiplication of sparse matrices
Cent. Eur. J. Cop. Sci. 41) 2014 1-11 DOI: 10.2478/s13537-014-0201-x Central European Journal of Coputer Science A note on the ultiplication of sparse atrices Research Article Keivan Borna 12, Sohrab Aboozarkhani
More informationSharp Time Data Tradeoffs for Linear Inverse Problems
Sharp Tie Data Tradeoffs for Linear Inverse Probles Saet Oyak Benjain Recht Mahdi Soltanolkotabi January 016 Abstract In this paper we characterize sharp tie-data tradeoffs for optiization probles used
More informationarxiv: v1 [math.nt] 14 Sep 2014
ROTATION REMAINDERS P. JAMESON GRABER, WASHINGTON AND LEE UNIVERSITY 08 arxiv:1409.411v1 [ath.nt] 14 Sep 014 Abstract. We study properties of an array of nubers, called the triangle, in which each row
More informationASSUME a source over an alphabet size m, from which a sequence of n independent samples are drawn. The classical
IEEE TRANSACTIONS ON INFORMATION THEORY Large Alphabet Source Coding using Independent Coponent Analysis Aichai Painsky, Meber, IEEE, Saharon Rosset and Meir Feder, Fellow, IEEE arxiv:67.7v [cs.it] Jul
More informationCOS 424: Interacting with Data. Written Exercises
COS 424: Interacting with Data Hoework #4 Spring 2007 Regression Due: Wednesday, April 18 Written Exercises See the course website for iportant inforation about collaboration and late policies, as well
More informationList Scheduling and LPT Oliver Braun (09/05/2017)
List Scheduling and LPT Oliver Braun (09/05/207) We investigate the classical scheduling proble P ax where a set of n independent jobs has to be processed on 2 parallel and identical processors (achines)
More informationOcean 420 Physical Processes in the Ocean Project 1: Hydrostatic Balance, Advection and Diffusion Answers
Ocean 40 Physical Processes in the Ocean Project 1: Hydrostatic Balance, Advection and Diffusion Answers 1. Hydrostatic Balance a) Set all of the levels on one of the coluns to the lowest possible density.
More informationEfficient Filter Banks And Interpolators
Efficient Filter Banks And Interpolators A. G. DEMPSTER AND N. P. MURPHY Departent of Electronic Systes University of Westinster 115 New Cavendish St, London W1M 8JS United Kingdo Abstract: - Graphical
More informationCurious Bounds for Floor Function Sums
1 47 6 11 Journal of Integer Sequences, Vol. 1 (018), Article 18.1.8 Curious Bounds for Floor Function Sus Thotsaporn Thanatipanonda and Elaine Wong 1 Science Division Mahidol University International
More informationa a a a a a a m a b a b
Algebra / Trig Final Exa Study Guide (Fall Seester) Moncada/Dunphy Inforation About the Final Exa The final exa is cuulative, covering Appendix A (A.1-A.5) and Chapter 1. All probles will be ultiple choice
More informationtime time δ jobs jobs
Approxiating Total Flow Tie on Parallel Machines Stefano Leonardi Danny Raz y Abstract We consider the proble of optiizing the total ow tie of a strea of jobs that are released over tie in a ultiprocessor
More informationLATTICE POINT SOLUTION OF THE GENERALIZED PROBLEM OF TERQUEi. AND AN EXTENSION OF FIBONACCI NUMBERS.
i LATTICE POINT SOLUTION OF THE GENERALIZED PROBLEM OF TERQUEi. AND AN EXTENSION OF FIBONACCI NUMBERS. C. A. CHURCH, Jr. and H. W. GOULD, W. Virginia University, Morgantown, W. V a. In this paper we give
More informationTHE CONSTRUCTION OF GOOD EXTENSIBLE RANK-1 LATTICES. 1. Introduction We are interested in approximating a high dimensional integral [0,1]
MATHEMATICS OF COMPUTATION Volue 00, Nuber 0, Pages 000 000 S 0025-578(XX)0000-0 THE CONSTRUCTION OF GOOD EXTENSIBLE RANK- LATTICES JOSEF DICK, FRIEDRICH PILLICHSHAMMER, AND BENJAMIN J. WATERHOUSE Abstract.
More informationOn the Inapproximability of Vertex Cover on k-partite k-uniform Hypergraphs
On the Inapproxiability of Vertex Cover on k-partite k-unifor Hypergraphs Venkatesan Guruswai and Rishi Saket Coputer Science Departent Carnegie Mellon University Pittsburgh, PA 1513. Abstract. Coputing
More informationGenetic Algorithms using Low-Discrepancy Sequences
Genetic Algoriths using Low-Discrepancy Sequences ABSTRACT Shuhei Kiura Departent of Inforation and Knowledge Engineering, Faculty of Engineering, Tottori University 4-, Koyaa Minai, Tottori, JAPAN kiura@iketottori-uacjp
More informationConvex Programming for Scheduling Unrelated Parallel Machines
Convex Prograing for Scheduling Unrelated Parallel Machines Yossi Azar Air Epstein Abstract We consider the classical proble of scheduling parallel unrelated achines. Each job is to be processed by exactly
More informationKernel Methods and Support Vector Machines
Intelligent Systes: Reasoning and Recognition Jaes L. Crowley ENSIAG 2 / osig 1 Second Seester 2012/2013 Lesson 20 2 ay 2013 Kernel ethods and Support Vector achines Contents Kernel Functions...2 Quadratic
More informationLower Bounds for Quantized Matrix Completion
Lower Bounds for Quantized Matrix Copletion Mary Wootters and Yaniv Plan Departent of Matheatics University of Michigan Ann Arbor, MI Eail: wootters, yplan}@uich.edu Mark A. Davenport School of Elec. &
More informationA remark on a success rate model for DPA and CPA
A reark on a success rate odel for DPA and CPA A. Wieers, BSI Version 0.5 andreas.wieers@bsi.bund.de Septeber 5, 2018 Abstract The success rate is the ost coon evaluation etric for easuring the perforance
More informationNon-Parametric Non-Line-of-Sight Identification 1
Non-Paraetric Non-Line-of-Sight Identification Sinan Gezici, Hisashi Kobayashi and H. Vincent Poor Departent of Electrical Engineering School of Engineering and Applied Science Princeton University, Princeton,
More informationOBJECTIVES INTRODUCTION
M7 Chapter 3 Section 1 OBJECTIVES Suarize data using easures of central tendency, such as the ean, edian, ode, and idrange. Describe data using the easures of variation, such as the range, variance, and
More informationINNER CONSTRAINTS FOR A 3-D SURVEY NETWORK
eospatial Science INNER CONSRAINS FOR A 3-D SURVEY NEWORK hese notes follow closely the developent of inner constraint equations by Dr Willie an, Departent of Building, School of Design and Environent,
More informationA method to determine relative stroke detection efficiencies from multiplicity distributions
A ethod to deterine relative stroke detection eiciencies ro ultiplicity distributions Schulz W. and Cuins K. 2. Austrian Lightning Detection and Inoration Syste (ALDIS), Kahlenberger Str.2A, 90 Vienna,
More informationThe Weierstrass Approximation Theorem
36 The Weierstrass Approxiation Theore Recall that the fundaental idea underlying the construction of the real nubers is approxiation by the sipler rational nubers. Firstly, nubers are often deterined
More informationA Better Algorithm For an Ancient Scheduling Problem. David R. Karger Steven J. Phillips Eric Torng. Department of Computer Science
A Better Algorith For an Ancient Scheduling Proble David R. Karger Steven J. Phillips Eric Torng Departent of Coputer Science Stanford University Stanford, CA 9435-4 Abstract One of the oldest and siplest
More informationUpper bound on false alarm rate for landmine detection and classification using syntactic pattern recognition
Upper bound on false alar rate for landine detection and classification using syntactic pattern recognition Ahed O. Nasif, Brian L. Mark, Kenneth J. Hintz, and Nathalia Peixoto Dept. of Electrical and
More informationIn this chapter, we consider several graph-theoretic and probabilistic models
THREE ONE GRAPH-THEORETIC AND STATISTICAL MODELS 3.1 INTRODUCTION In this chapter, we consider several graph-theoretic and probabilistic odels for a social network, which we do under different assuptions
More informationMULTIPLAYER ROCK-PAPER-SCISSORS
MULTIPLAYER ROCK-PAPER-SCISSORS CHARLOTTE ATEN Contents 1. Introduction 1 2. RPS Magas 3 3. Ites as a Function of Players and Vice Versa 5 4. Algebraic Properties of RPS Magas 6 References 6 1. Introduction
More informationConstant-Space String-Matching. in Sublinear Average Time. (Extended Abstract) Wojciech Rytter z. Warsaw University. and. University of Liverpool
Constant-Space String-Matching in Sublinear Average Tie (Extended Abstract) Maxie Crocheore Universite de Marne-la-Vallee Leszek Gasieniec y Max-Planck Institut fur Inforatik Wojciech Rytter z Warsaw University
More informationThe Simplex Method is Strongly Polynomial for the Markov Decision Problem with a Fixed Discount Rate
The Siplex Method is Strongly Polynoial for the Markov Decision Proble with a Fixed Discount Rate Yinyu Ye April 20, 2010 Abstract In this note we prove that the classic siplex ethod with the ost-negativereduced-cost
More informationDistributed Subgradient Methods for Multi-agent Optimization
1 Distributed Subgradient Methods for Multi-agent Optiization Angelia Nedić and Asuan Ozdaglar October 29, 2007 Abstract We study a distributed coputation odel for optiizing a su of convex objective functions
More informationModel Fitting. CURM Background Material, Fall 2014 Dr. Doreen De Leon
Model Fitting CURM Background Material, Fall 014 Dr. Doreen De Leon 1 Introduction Given a set of data points, we often want to fit a selected odel or type to the data (e.g., we suspect an exponential
More informationarxiv: v1 [math.na] 10 Oct 2016
GREEDY GAUSS-NEWTON ALGORITHM FOR FINDING SPARSE SOLUTIONS TO NONLINEAR UNDERDETERMINED SYSTEMS OF EQUATIONS MÅRTEN GULLIKSSON AND ANNA OLEYNIK arxiv:6.395v [ath.na] Oct 26 Abstract. We consider the proble
More informationarxiv: v1 [cs.ds] 29 Jan 2012
A parallel approxiation algorith for ixed packing covering seidefinite progras arxiv:1201.6090v1 [cs.ds] 29 Jan 2012 Rahul Jain National U. Singapore January 28, 2012 Abstract Penghui Yao National U. Singapore
More informationProbability Distributions
Probability Distributions In Chapter, we ephasized the central role played by probability theory in the solution of pattern recognition probles. We turn now to an exploration of soe particular exaples
More informationPage 1 Lab 1 Elementary Matrix and Linear Algebra Spring 2011
Page Lab Eleentary Matri and Linear Algebra Spring 0 Nae Due /03/0 Score /5 Probles through 4 are each worth 4 points.. Go to the Linear Algebra oolkit site ransforing a atri to reduced row echelon for
More informationFixed-to-Variable Length Distribution Matching
Fixed-to-Variable Length Distribution Matching Rana Ali Ajad and Georg Böcherer Institute for Counications Engineering Technische Universität München, Gerany Eail: raa2463@gail.co,georg.boecherer@tu.de
More informationA Simple Regression Problem
A Siple Regression Proble R. M. Castro March 23, 2 In this brief note a siple regression proble will be introduced, illustrating clearly the bias-variance tradeoff. Let Y i f(x i ) + W i, i,..., n, where
More informationHomework 3 Solutions CSE 101 Summer 2017
Hoework 3 Solutions CSE 0 Suer 207. Scheduling algoriths The following n = 2 jobs with given processing ties have to be scheduled on = 3 parallel and identical processors with the objective of iniizing
More informationNBN Algorithm Introduction Computational Fundamentals. Bogdan M. Wilamoswki Auburn University. Hao Yu Auburn University
NBN Algorith Bogdan M. Wilaoswki Auburn University Hao Yu Auburn University Nicholas Cotton Auburn University. Introduction. -. Coputational Fundaentals - Definition of Basic Concepts in Neural Network
More informationLecture 21. Interior Point Methods Setup and Algorithm
Lecture 21 Interior Point Methods In 1984, Kararkar introduced a new weakly polynoial tie algorith for solving LPs [Kar84a], [Kar84b]. His algorith was theoretically faster than the ellipsoid ethod and
More informationDesign of Spatially Coupled LDPC Codes over GF(q) for Windowed Decoding
IEEE TRANSACTIONS ON INFORMATION THEORY (SUBMITTED PAPER) 1 Design of Spatially Coupled LDPC Codes over GF(q) for Windowed Decoding Lai Wei, Student Meber, IEEE, David G. M. Mitchell, Meber, IEEE, Thoas
More informationChaotic Coupled Map Lattices
Chaotic Coupled Map Lattices Author: Dustin Keys Advisors: Dr. Robert Indik, Dr. Kevin Lin 1 Introduction When a syste of chaotic aps is coupled in a way that allows the to share inforation about each
More informationON THE TWO-LEVEL PRECONDITIONING IN LEAST SQUARES METHOD
PROCEEDINGS OF THE YEREVAN STATE UNIVERSITY Physical and Matheatical Sciences 04,, p. 7 5 ON THE TWO-LEVEL PRECONDITIONING IN LEAST SQUARES METHOD M a t h e a t i c s Yu. A. HAKOPIAN, R. Z. HOVHANNISYAN
More informationAlgorithms for parallel processor scheduling with distinct due windows and unit-time jobs
BULLETIN OF THE POLISH ACADEMY OF SCIENCES TECHNICAL SCIENCES Vol. 57, No. 3, 2009 Algoriths for parallel processor scheduling with distinct due windows and unit-tie obs A. JANIAK 1, W.A. JANIAK 2, and
More informationSupport Vector Machine Classification of Uncertain and Imbalanced data using Robust Optimization
Recent Researches in Coputer Science Support Vector Machine Classification of Uncertain and Ibalanced data using Robust Optiization RAGHAV PAT, THEODORE B. TRAFALIS, KASH BARKER School of Industrial Engineering
More informationFairness via priority scheduling
Fairness via priority scheduling Veeraruna Kavitha, N Heachandra and Debayan Das IEOR, IIT Bobay, Mubai, 400076, India vavitha,nh,debayan}@iitbacin Abstract In the context of ulti-agent resource allocation
More informationVulnerability of MRD-Code-Based Universal Secure Error-Correcting Network Codes under Time-Varying Jamming Links
Vulnerability of MRD-Code-Based Universal Secure Error-Correcting Network Codes under Tie-Varying Jaing Links Jun Kurihara KDDI R&D Laboratories, Inc 2 5 Ohara, Fujiino, Saitaa, 356 8502 Japan Eail: kurihara@kddilabsjp
More informationQuantum algorithms (CO 781, Winter 2008) Prof. Andrew Childs, University of Waterloo LECTURE 15: Unstructured search and spatial search
Quantu algoriths (CO 781, Winter 2008) Prof Andrew Childs, University of Waterloo LECTURE 15: Unstructured search and spatial search ow we begin to discuss applications of quantu walks to search algoriths
More informationLecture 9 November 23, 2015
CSC244: Discrepancy Theory in Coputer Science Fall 25 Aleksandar Nikolov Lecture 9 Noveber 23, 25 Scribe: Nick Spooner Properties of γ 2 Recall that γ 2 (A) is defined for A R n as follows: γ 2 (A) = in{r(u)
More informationIterative Decoding of LDPC Codes over the q-ary Partial Erasure Channel
1 Iterative Decoding of LDPC Codes over the q-ary Partial Erasure Channel Rai Cohen, Graduate Student eber, IEEE, and Yuval Cassuto, Senior eber, IEEE arxiv:1510.05311v2 [cs.it] 24 ay 2016 Abstract In
More informationLeast Squares Fitting of Data
Least Squares Fitting of Data David Eberly, Geoetric Tools, Redond WA 98052 https://www.geoetrictools.co/ This work is licensed under the Creative Coons Attribution 4.0 International License. To view a
More informationMulti-Scale/Multi-Resolution: Wavelet Transform
Multi-Scale/Multi-Resolution: Wavelet Transfor Proble with Fourier Fourier analysis -- breaks down a signal into constituent sinusoids of different frequencies. A serious drawback in transforing to the
More informationCSE525: Randomized Algorithms and Probabilistic Analysis May 16, Lecture 13
CSE55: Randoied Algoriths and obabilistic Analysis May 6, Lecture Lecturer: Anna Karlin Scribe: Noah Siegel, Jonathan Shi Rando walks and Markov chains This lecture discusses Markov chains, which capture
More informationIdentical Maximum Likelihood State Estimation Based on Incremental Finite Mixture Model in PHD Filter
Identical Maxiu Lielihood State Estiation Based on Increental Finite Mixture Model in PHD Filter Gang Wu Eail: xjtuwugang@gail.co Jing Liu Eail: elelj20080730@ail.xjtu.edu.cn Chongzhao Han Eail: czhan@ail.xjtu.edu.cn
More informationCS Lecture 13. More Maximum Likelihood
CS 6347 Lecture 13 More Maxiu Likelihood Recap Last tie: Introduction to axiu likelihood estiation MLE for Bayesian networks Optial CPTs correspond to epirical counts Today: MLE for CRFs 2 Maxiu Likelihood
More informationInference in the Presence of Likelihood Monotonicity for Polytomous and Logistic Regression
Advances in Pure Matheatics, 206, 6, 33-34 Published Online April 206 in SciRes. http://www.scirp.org/journal/ap http://dx.doi.org/0.4236/ap.206.65024 Inference in the Presence of Likelihood Monotonicity
More informationA new type of lower bound for the largest eigenvalue of a symmetric matrix
Linear Algebra and its Applications 47 7 9 9 www.elsevier.co/locate/laa A new type of lower bound for the largest eigenvalue of a syetric atrix Piet Van Mieghe Delft University of Technology, P.O. Box
More informationRECOVERY OF A DENSITY FROM THE EIGENVALUES OF A NONHOMOGENEOUS MEMBRANE
Proceedings of ICIPE rd International Conference on Inverse Probles in Engineering: Theory and Practice June -8, 999, Port Ludlow, Washington, USA : RECOVERY OF A DENSITY FROM THE EIGENVALUES OF A NONHOMOGENEOUS
More informationThe Transactional Nature of Quantum Information
The Transactional Nature of Quantu Inforation Subhash Kak Departent of Coputer Science Oklahoa State University Stillwater, OK 7478 ABSTRACT Inforation, in its counications sense, is a transactional property.
More informationStatistical properties of contact maps
PHYSICAL REVIEW E VOLUME 59, NUMBER 1 JANUARY 1999 Statistical properties of contact aps Michele Vendruscolo, 1 Balakrishna Subraanian, 2 Ido Kanter, 3 Eytan Doany, 1 and Joel Lebowitz 2 1 Departent of
More informationA Note on the Applied Use of MDL Approximations
A Note on the Applied Use of MDL Approxiations Daniel J. Navarro Departent of Psychology Ohio State University Abstract An applied proble is discussed in which two nested psychological odels of retention
More informationGenetic Algorithm Search for Stent Design Improvements
Genetic Algorith Search for Stent Design Iproveents K. Tesch, M.A. Atherton & M.W. Collins, South Bank University, London, UK Abstract This paper presents an optiisation process for finding iproved stent
More informationInteractive Markov Models of Evolutionary Algorithms
Cleveland State University EngagedScholarship@CSU Electrical Engineering & Coputer Science Faculty Publications Electrical Engineering & Coputer Science Departent 2015 Interactive Markov Models of Evolutionary
More informationTesting Properties of Collections of Distributions
Testing Properties of Collections of Distributions Reut Levi Dana Ron Ronitt Rubinfeld April 9, 0 Abstract We propose a fraework for studying property testing of collections of distributions, where the
More informatione-companion ONLY AVAILABLE IN ELECTRONIC FORM
OPERATIONS RESEARCH doi 10.1287/opre.1070.0427ec pp. ec1 ec5 e-copanion ONLY AVAILABLE IN ELECTRONIC FORM infors 07 INFORMS Electronic Copanion A Learning Approach for Interactive Marketing to a Custoer
More informationLow-complexity, Low-memory EMS algorithm for non-binary LDPC codes
Low-coplexity, Low-eory EMS algorith for non-binary LDPC codes Adrian Voicila,David Declercq, François Verdier ETIS ENSEA/CP/CNRS MR-85 954 Cergy-Pontoise, (France) Marc Fossorier Dept. Electrical Engineering
More informationSoft Computing Techniques Help Assign Weights to Different Factors in Vulnerability Analysis
Soft Coputing Techniques Help Assign Weights to Different Factors in Vulnerability Analysis Beverly Rivera 1,2, Irbis Gallegos 1, and Vladik Kreinovich 2 1 Regional Cyber and Energy Security Center RCES
More informationAn Algorithm for Quantization of Discrete Probability Distributions
An Algorith for Quantization of Discrete Probability Distributions Yuriy A. Reznik Qualco Inc., San Diego, CA Eail: yreznik@ieee.org Abstract We study the proble of quantization of discrete probability
More informationLow complexity bit parallel multiplier for GF(2 m ) generated by equally-spaced trinomials
Inforation Processing Letters 107 008 11 15 www.elsevier.co/locate/ipl Low coplexity bit parallel ultiplier for GF generated by equally-spaced trinoials Haibin Shen a,, Yier Jin a,b a Institute of VLSI
More informationTopic 5a Introduction to Curve Fitting & Linear Regression
/7/08 Course Instructor Dr. Rayond C. Rup Oice: A 337 Phone: (95) 747 6958 E ail: rcrup@utep.edu opic 5a Introduction to Curve Fitting & Linear Regression EE 4386/530 Coputational ethods in EE Outline
More informationBoosting with log-loss
Boosting with log-loss Marco Cusuano-Towner Septeber 2, 202 The proble Suppose we have data exaples {x i, y i ) i =... } for a two-class proble with y i {, }. Let F x) be the predictor function with the
More informationUsing a De-Convolution Window for Operating Modal Analysis
Using a De-Convolution Window for Operating Modal Analysis Brian Schwarz Vibrant Technology, Inc. Scotts Valley, CA Mark Richardson Vibrant Technology, Inc. Scotts Valley, CA Abstract Operating Modal Analysis
More informationSymbolic Analysis as Universal Tool for Deriving Properties of Non-linear Algorithms Case study of EM Algorithm
Acta Polytechnica Hungarica Vol., No., 04 Sybolic Analysis as Universal Tool for Deriving Properties of Non-linear Algoriths Case study of EM Algorith Vladiir Mladenović, Miroslav Lutovac, Dana Porrat
More informationExponential sums and the distribution of inversive congruential pseudorandom numbers with prime-power modulus
ACTA ARITHMETICA XCII1 (2000) Exponential sus and the distribution of inversive congruential pseudorando nubers with prie-power odulus by Harald Niederreiter (Vienna) and Igor E Shparlinski (Sydney) 1
More informationFinite fields. and we ve used it in various examples and homework problems. In these notes I will introduce more finite fields
Finite fields I talked in class about the field with two eleents F 2 = {, } and we ve used it in various eaples and hoework probles. In these notes I will introduce ore finite fields F p = {,,...,p } for
More informationThe Fundamental Basis Theorem of Geometry from an algebraic point of view
Journal of Physics: Conference Series PAPER OPEN ACCESS The Fundaental Basis Theore of Geoetry fro an algebraic point of view To cite this article: U Bekbaev 2017 J Phys: Conf Ser 819 012013 View the article
More informationExact tensor completion with sum-of-squares
Proceedings of Machine Learning Research vol 65:1 54, 2017 30th Annual Conference on Learning Theory Exact tensor copletion with su-of-squares Aaron Potechin Institute for Advanced Study, Princeton David
More information. The univariate situation. It is well-known for a long tie that denoinators of Pade approxiants can be considered as orthogonal polynoials with respe
PROPERTIES OF MULTIVARIATE HOMOGENEOUS ORTHOGONAL POLYNOMIALS Brahi Benouahane y Annie Cuyt? Keywords Abstract It is well-known that the denoinators of Pade approxiants can be considered as orthogonal
More informationGenetic Quantum Algorithm and its Application to Combinatorial Optimization Problem
Genetic Quantu Algorith and its Application to Cobinatorial Optiization Proble Kuk-Hyun Han Dept. of Electrical Engineering, KAIST, 373-, Kusong-dong Yusong-gu Taejon, 305-70, Republic of Korea khhan@vivaldi.kaist.ac.kr
More informationA Self-Organizing Model for Logical Regression Jerry Farlow 1 University of Maine. (1900 words)
1 A Self-Organizing Model for Logical Regression Jerry Farlow 1 University of Maine (1900 words) Contact: Jerry Farlow Dept of Matheatics Univeristy of Maine Orono, ME 04469 Tel (07) 866-3540 Eail: farlow@ath.uaine.edu
More informationDeflation of the I-O Series Some Technical Aspects. Giorgio Rampa University of Genoa April 2007
Deflation of the I-O Series 1959-2. Soe Technical Aspects Giorgio Rapa University of Genoa g.rapa@unige.it April 27 1. Introduction The nuber of sectors is 42 for the period 1965-2 and 38 for the initial
More informationPrincipal Components Analysis
Principal Coponents Analysis Cheng Li, Bingyu Wang Noveber 3, 204 What s PCA Principal coponent analysis (PCA) is a statistical procedure that uses an orthogonal transforation to convert a set of observations
More informationGraphical Models in Local, Asymmetric Multi-Agent Markov Decision Processes
Graphical Models in Local, Asyetric Multi-Agent Markov Decision Processes Ditri Dolgov and Edund Durfee Departent of Electrical Engineering and Coputer Science University of Michigan Ann Arbor, MI 48109
More informationMulti-Dimensional Hegselmann-Krause Dynamics
Multi-Diensional Hegselann-Krause Dynaics A. Nedić Industrial and Enterprise Systes Engineering Dept. University of Illinois Urbana, IL 680 angelia@illinois.edu B. Touri Coordinated Science Laboratory
More information