Define - starting approximation for the parameters (p) - observational data (o) - solution criterion (e.g. number of iterations)
|
|
- Eleanore Rogers
- 5 years ago
- Views:
Transcription
1 Global Iterative Solution Distributed proessing of the attitude updating L. Lindegren (21 May 2001) SAG LL 37 Abstrat. The attitude updating algorithm given in GAIA LL 24 (v. 2) is modified to allow distributed proessing in whih eah proessor takes are of a subset of the soures aording to the sky partitioning. A detailed (but very simplified) numerial example is provided. The same priniple an be used to distribute the alibration and global proessing. 1 Introdution In GAIA LL 34 (v. 2), detailed algorithms in the form of fortran subroutines were provided for a number of proesses in the Global Iterative Solution (GIS) to be implemented in the on-going GAIA Data Aess and Analysis Study (GDAAS). It is realled that the GIS basially onsists of the following proesses to be exeuted ylially until onvergene: ffl soure proessing one soure (astronomial objet) at a time; ffl attitude proessing one attitude interval (some 12 hr) at a time; ffl alibration one CCD at a time; ffl global proessing needs to onsider the whole mission. It is envisaged that the GIS will be exeuted in a distributed enviroment with several proesses running in parallel. In order to minimize large-sale data shuffling it seems natural to assoiate a speifi area on the sky (following the HTM division) with eah loal proessor. This makes it very straightforward to distribute the soure proessing, sine that only requires aess to a very small area of the sky at a time. It is less obvious if and how the other proesses an be distributed without a dramati inrease in overhead. In this doument the attitude proessing is desribed in some detail, inluding a `worked example' using the subroutines supplied with GAIA LL 34 (v. 2). All the proesses mentioned above are instanes of non-linear least-squares estimation, implemented (as proposed for the GIS) by the lassial method using normal equations. The general idea for distributed proessing using normal equations is disussed in Set. 2 and then applied to the attitude proessing as an example involving fairly many unknowns (Set. 3). The same method an be used for the alibration and global proessing; here the normal equations only involve one or a few unknowns, so the situation is muh simpler than for the attitude proessing. 1
2 2 Least-squares by distributed proessing All the proesses desribed above solve a non-linear least-squares problem by the use of normal equations. The non-linearity arises primarily from the need to take are of outliers by means of downweighting. Let p be the parameter vetor to be estimated and m = dim(p) the number of unknowns (parameters). For given observations o (with dim(o) fl m) the quantity to minimize is Q(p) = X k ρ ok f(t k jp) ff k where f is the observation model, ρ the loss funtion, ff the standard error assigned to eah observation, and k the observation index. For ordinary least-squares the loss funtion is simply ρ(z) = z 2, but in order to aommodate outliers this should be modified by a fator whih redues the weight of observations with large jzj. Tentatively, the Huber metri has been proposed, (1) ρ(z) =z 2 ( 1 if jzj» (2jzj )z 2 otherwise (2) (note the typo in Eq. 2 of GAIA LL 34.2). Problem (1) an be solved iteratively from an initial approximation p by means of the equations: A ij = X k w k ff 2 j ; i; j =1:::m; (3) X w k b i = ff 2 (o k f k ) ; i =1:::m; (4) k i x = A 1 b ; (5) p ψ p + x : (6) This is iterated until jxj is negligible. Here w k is the weight redution fator alulated for eah observation based on the size of the normalized residual z k (o k f k )=ff k. 1 The linear system of equations Ax = b is alled the normal equations (A 1 is also an estimate of the ovariane of the resulting p). Sine A is symmetri of size m m and b is of size m 1 the total number of reals to be stored for the normal equations is (at most) m(m + 3)=2. (For the attitude updating A is band-diagonal whih means that the atual storage is muh smaller and inreases only linearly with m.) Figure 1 is a shemati illustration of the steps neessary to perform the non-linear leastsquares estimation aording to this method. In addition to the normal equations, some additional statistis are usually alulated in step 2, here exemplified by the total number of observations n and the mean weight redution fator w. 1 Atually, thedownweighting is not done quite as intended in the attahed subroutines, i.e. what is done does not orrespond to the Huber metri. This is of little onsequene at the present stage, but eventually the treatment ofoutliers should be onsidered more arefully, whih probably will lead to a modifiation at least of the funtion whuber. 2
3 Step 0: Initialize leastsquares problem Define - starting approximation for the parameters (p) - observational data (o) - solution riterion (e.g. number of iterations) Step 1: Initialize normal equations A := 0 b := 0 n := 0 wsum := 0 (et.) Step 2: Aumulate normal equations A := A + (new data) b := b + (new data) n := n + 1 wsum := wsum + (new data) (et.) loop over observations Iterate a few times, or until x is negligible Step 3: Solve n.e. and update parameters x := A -1 b p := p + x w := wsum/n (et.) After onvergene, output p, n, w and any other relevant statisti Figure 1. Flow diagram for a typial non-linear least-squares problem suh as the attitude updating. An important onsideration is that the aumulated arrays (A, b and various statistis) are relatively small and of onstant sizeforagiven problem. For instane, A and b together need at most m(m +3)=2 reals, if m is the dimension of p. For the attitude updating A is a sparse (band-diagonal) matrix whih an be stored muh more ompatly. The algorithm attupd takes full advantage of this property. The onveniene of the normal-equations approah derives partly from the fats (a) that their size is typially muh less than what is needed to store the observations thus it is more onvenient to send normal equations around than observations, and (b) that the sums in Eq. (2) an be aumulated independently for disjoint sets of observations (k), and then o-added before solving the normals and updating the parameter estimate. If the observational data are stored loally on N different nodes (eah with its own proessor or proessors), then an obvious modifiation of this sheme is as shown in Figure 2. 3 A worked example A numerial example demonstrates how the attitude updating subroutine attupd from GAIA LL 34.2 was modified to allow distributed proessing. 3
4 Step 0: Initialize leastsquares problem Define - starting approximation for the parameters (p) - observational data (o) - solution riterion (e.g. number of iterations) Distribute these to N proessors Pro 1 Pro 2 Pro N Step 1: Initialize loal normal equations A := 0 b := 0 n := 0 (et.) A := 0 b := 0 n := 0 (et.)... A := 0 b := 0 n := 0 (et.) Step 2: Aumulate loal normal equations A := A + (new data) b := b + (new data) n := n + 1 (et.) A := A + (new data) b := b + (new data) n := n + 1 (et.)... A := A + (new data) b := b + (new data) n := n + 1 (et.) loop over observations (N disjoint sets) Collet normal equations from loal proessors A := A 1 + A A N b := b 1 + b b N n := n 1 + n n N (et) Step 3: Solve n.e. and update parameters x := A -1 b p := p + x (et.) Iterate a few times, or until x is negligible After onvergene, output p, n, w and any other relevant statisti Figure 2. Flow diagram for the same non-linear least-squares problem as in Fig. 1, but with the initialization and aumulation of normal equations distributed on N loal proessors. Eah proessor takes are of a disjoint set of observations, e.g. orresponding to a ertain area on the sky. The data flow between proessors is minimal beause the dimensions of p, A, b, et are small ompared with the observations, whih are already stored loally. 3.1 The soure atalogue In order to avoid having to generate and store a huge soure atalogue, a seeded soure approah is used. This means that the soure identifiation number (an integer in the range from 1 to ' )isusedtoseedarandomnumber generator whih then provides the soure parameters aording to a given distribution model. The model used here is extremely simple: the position (ff; ffi) is drawn from a uniform distribution over the sky while the remaining four astrometri parameters (parallax, proper motion and radial veloity) are set to zero. The subroutine getstar returns the six astrometri parameters for a given soure number. Note that getstar uses a separate, very primitive random 4
5 number generator alled ran Observational data for input Simple geometri observations in the form of angular oordinates ( ; ) were generated for a time interval of one spin period (3 hr). The text file obs.prn is a list of 4253 simulated observations for the time interval 0:01 <t<0:135 days. The nominal sanning law was used as `true' attitude. Both Astro fields were used, eah of size 0:66 0:66 deg, and a basi angle of 106 deg. The path of eah soure aross the several CCDs was not simulated; instead seleted observations were taken at random time intervals. The total number of observations per time interval is muh smaller that for the real mission, but suffiient to fix the attitude for the hosen knot frequeny. The observational standard errors were taken to be onstant, viz. 0.1 mas along-san and 2 mas aross-san per observation. For eah observation, obs.prn gives the time (days), the soure number, the field used (1 or 2), the observed along-san oordinate (rad), the standard error ff (rad), the observed aross-san oordinate (rad), and the standard error ff (rad). The first five lines of obs.prn are shown by the following table: t ID Field ff ff Starting approximation The starting approximation for the attitude spline is given by the file q initial.prn, whih ontains one line per knot. The olumns are: time (days), S 1n, S 2n, S 3n, and S 4n, where S mn are the basi spline oeffiients defined in Set. 5.3 of GAIA LL Note that the first olumn ontains the adopted knot sequene, inluding three points on either side of the time interval over whih the attitude is defined. t S 1n S 2n S 3n S 4n Fitting the attitude spline (non-distributed proessing) The following main program implements the method in Figure 1 by suessive alls to the subroutine attupd (supplied with GAIA LL 34.2). 2 This routine should be kept separate from the more sophistiated routines used to generate observation noise and similar things where a good degree of randomness is important. never be reinitialized (given a new seed) within a simulation. These latter routines should 5
6 PROGRAM fit_attitude Fit an attitude spline to the observations in file_obs. An initial approximation of the attitude (as well as the knot sequene) is defined in file_ini. The final attitude estimate (after niter least-squares iterations) is given in file_out. N.B.: the parameter N must hosen suh that N+4 does not exeed the number of rows in file_ini; moreover, file_obs must ontain suffiient observations overing the time interval tau(2:n-1), where tau(-1:n+2) is the knot sequene defined by the first olumn of file_ini. L Lindegren (Lund Observatory) 19 May 2001 IMPLICIT REAL*8 (a-h,o-z) INCLUDE 'onst_math.h' PARAMETER (N = 273, NB = 4*N, LDA = 16) REAL*8 tau(-1:n+2), bs0(4,n), bs1(4,n), wa(lda,nb), wb(nb), + w(nb), wd(nb), glb(1), apar(6), + wr(2), fobs(2), ferr(2) CHARACTER file_obs*80, file_ini*80, file_out*80 lu0 = 10 lu1 = 11 file_obs = 'obs.prn' file_ini = 'q_initial.prn' file_out = 'qdata1.prn' glb(1) = 1d STEP Define knot sequene (tau), initial spline approximation (bs0) and number of iterations (niter): OPEN(lu0,file=file_ini) DO i = -1, N+2 IF (i.ge. 1.AND. i.le. N) THEN READ(lu0,*) tau(i), (bs0(j,i),j=1,4) ELSE READ(lu0,*) tau(i) ENDIF ENDDO CLOSE(lu0) tbeg = tau(2) tend = tau(n-1) niter = 5 This is the main iteration loop for the non-linear least-squares -- in eah loop, bs0 is improved to bs1: DO iter = 1, niter STEP Initialize normal equations et: mode = 0 CALL attupd(n, tau, bs0, mode, glb, apar, tobs, fobs, ferr, + wa, wb, w, wd, bs1, ierr, wr) Loop over the observations, aepting those in [tbeg,tend]: OPEN(lu0,file=file_obs) nobs = CONTINUE READ(lu0,*,END=200) tobs, id, ia, (fobs(i),ferr(i),i=1,2) 6
7 IF (tobs.lt. tbeg) GOTO 100 IF (tobs.gt. tend) GOTO 200 CALL getstar(id, apar) STEP Update normal equations: mode = 2 CALL attupd(n, tau, bs0, mode, glb, apar, tobs, fobs, ferr, + wa, wb, w, wd, bs1, ierr, wr) IF (ierr.ne. 0) print *, 'ierr,tobs=',ierr,tobs nobs = nobs + 1 GOTO 100 No more observations (or tobs exeeds tend): 200 CONTINUE CLOSE(lu0) STEP Solve normal equations and update spline oeffiients: mode = 3 CALL attupd(n, tau, bs0, mode, glb, apar, tobs, fobs, ferr, + wa, wb, w, wd, bs1, ierr, wr) WRITE(*,'(/x,a,i10)') 'iteration = ', iter WRITE(*,'(x,a,i10)') 'attupd ierr = ', ierr WRITE(*,'(x,a,2i10)') 'nobs = ', nobs, nobs WRITE(*,'(x,a,2f10.6)') 'mean weight redution = ', wr Calulate the RMS update and (in preparation for the next iteration) opy bs1 to bs0: rms = 0d0 DO i = 1, N DO j = 1, 4 rms = rms + (bs1(j,i) - bs0(j,i))**2 bs0(j,i) = bs1(j,i) ENDDO ENDDO rms = dsqrt(rms/dble(n*4)) WRITE(*,'(x,a,g12.6)') 'RMS update (uas) = ', rms/uas ENDDO OUTPUT Output updated spline data to file_out: OPEN(lu1,file=file_out) DO i = -1, N+2 IF (i.ge. 1.AND. i.le. N) THEN WRITE(lu1,'(f15.9,4f20.16)') tau(i), (bs1(j,i),j=1,4) ELSE WRITE(lu1,'(f15.9,4f20.16)') tau(i), (0d0,j=1,4) ENDIF ENDDO CLOSE(11) WRITE(*,'(/x,a,a)') 'spline data written to: ', file_out(1:20) END Exeuting fit attitude generates the following standard output (out fit attitude), where the onvergene of the five iterations an be seen. Between the last two iterations the RMS update of the spline oeffiients S mn is below 0:02 μas. 7
8 iteration = 1 attupd ierr = 0 nobs = mean weight redution = RMS update (uas) = iteration = 2 attupd ierr = 0 nobs = mean weight redution = RMS update (uas) = iteration = 3 attupd ierr = 0 nobs = mean weight redution = RMS update (uas) = iteration = 4 attupd ierr = 0 nobs = mean weight redution = RMS update (uas) = iteration = 5 attupd ierr = 0 nobs = mean weight redution = RMS update (uas) = E-01 spline data written to: qdata1.prn The first five lines of the output data file qdata1.prn are shown below. The meaning of the olumns is the same as for the starting approximation approximation. t S 1n S 2n S 3n S 4n Fitting the attitude spline (distributed proessing) In attupd the same subroutine is used for steps 1, 2 and 3, as ontrolled by the argument mode. In order to implement the sheme in Figure 2, the routine is broken up into three separate routines alled attupd ini (initialize), attupd a (aumulate) and attupd fin (finalize). These are ontained in the file attupd dis.f. It would be possible to run these routines from a single main program (as for attupd), but in order to bring out more learly the modularity and the interation between the different modules, they are here atually run by separate main programs alled att ini, att a and att fin. Standard input is used to ontrol e.g. the proessor number and data are exhanged between them by means of temporary binary files (*.uf). This struture in priniple allows the attitude fitting to be run on a luster without speial parallellization ode, although it would not be very effiient due to the use of disk files for data exhange. 8
9 Here the NPROC = 4 proesses are atually run in sequene rather than in parallel, using the following shell sript (att bat): #!/bin/sh # unix bath program to illustrate distributed proessing # of attitude determination. # NITER iterations are made, eah onsisting of three steps: # 1. initialize data (att_ini) - only needed for first iteration # 2. aumulate normals (att_a) # 3. finalize and solve normals (att_fin) # After the last iteration, results are written to a text file. # The seond step is divided into NPROC separate proesses, # here run sequentially, but they ould also run in parallel. # Step 3 has to wait until all NPROC proesses are ready. # define number of iterations and number of distributed proesses: set NITER = 5 set NPROC = 4 # Step 1 (initialize): eho $NPROC att_ini # this is the iteration loop: set ITER = 1 while ( $ITER <= $NITER ) # Step 2 (here written expliitly for IPROC = 1, NPROC): eho 1 att_a eho 2 att_a eho 3 att_a eho 4 att_a # Step 3: eho $ITER $NPROC ITER++ end # onvert final result (qdata.uf) to text file (qdata.prn) # for easier inspetion: print_qdata # leanup: rm work*.uf obs_*.uf qdata.uf Exeuting att bat generates standard output (out att bat) whih is idential to the output from fit attitude exept for the name of the output data file. The output data file (qdata.prn) is also virtually idential to the previous one (qdata1.prn). This merely shows that the distribution of steps 1 and 2 is matematially (if not quite numerially) equivalent to the original algorithm. If there is some slight numerial differene between qdata.prn and qdata1.prn it will be aused by rounding errors in the aumulation of the normal equations. 9
10 4 Files distributed with this doument The files available with this doument are as follows: --- program units for non-distributed proessing (Set. 3.4): fit_attitude.f out_fit_attitude qdata1.prn --- program units for distributed proessing (Set. 3.5): att_a.f main program to aumulate normals att_bat bath program to run att_ini, att_a, att_fin att_fin.f main program to finalize attitude updating att_ini.f main program to initialize attitude updating attupd_dis.f subroutines for attitude updating out_att_bat standard output from att_bat print_qdata.f onvert binary output data to printable qdata.prn output data file --- general routines: onst_astr.h inlude file for astronomial onstants onst_math.h inlude file for mathematial onstants getstar.f return soure parameters (Set. 3.1) LL34.f subroutine pakage (see below) makefile makefile for ompiling all programs obs.prn input observational data (Set. 3.2) q_initial.prn initial spline approximation (Set. 3.3) Remark: the subroutine pakage LL34.f (v. 2.1) used here is in priniple the same as distributed with GAIA LL 34.2, exept that some unneessary print statements used for debugging and then left there by mistake have been removed. 10
ICCGLU. A Fortran IV subroutine to solve large sparse general systems of linear equations. J.J. Dongarra, G.K. Leaf and M. Minkoff.
http://www.netlib.org/linalg/ig-do 1 of 8 12/7/2009 11:14 AM ICCGLU A Fortran IV subroutine to solve large sparse general systems of linear equations. J.J. Dongarra, G.K. Leaf and M. Minkoff July, 1982
More informationarxiv: v1 [cs.db] 13 Sep 2017
An effiient lustering algorithm from the measure of loal Gaussian distribution Yuan-Yen Tai (Dated: May 27, 2018) In this paper, I will introdue a fast and novel lustering algorithm based on Gaussian distribution
More informationOutline: Software Design
Outline: Software Design. Goals History of software design ideas Design priniples Design methods Life belt or leg iron? (Budgen) Copyright Nany Leveson, Sept. 1999 A Little History... At first, struggling
More informationExtracting Partition Statistics from Semistructured Data
Extrating Partition Statistis from Semistrutured Data John N. Wilson Rihard Gourlay Robert Japp Mathias Neumüller Department of Computer and Information Sienes University of Strathlyde, Glasgow, UK {jnw,rsg,rpj,mathias}@is.strath.a.uk
More informationA DYNAMIC ACCESS CONTROL WITH BINARY KEY-PAIR
Malaysian Journal of Computer Siene, Vol 10 No 1, June 1997, pp 36-41 A DYNAMIC ACCESS CONTROL WITH BINARY KEY-PAIR Md Rafiqul Islam, Harihodin Selamat and Mohd Noor Md Sap Faulty of Computer Siene and
More informationrepresent = as a finite deimal" either in base 0 or in base. We an imagine that the omputer first omputes the mathematial = then rounds the result to
Sientifi Computing Chapter I Computer Arithmeti Jonathan Goodman Courant Institute of Mathemaial Sienes Last revised January, 00 Introdution One of the many soures of error in sientifi omputing is inexat
More informationPipelined Multipliers for Reconfigurable Hardware
Pipelined Multipliers for Reonfigurable Hardware Mithell J. Myjak and José G. Delgado-Frias Shool of Eletrial Engineering and Computer Siene, Washington State University Pullman, WA 99164-2752 USA {mmyjak,
More informationLearning Convention Propagation in BeerAdvocate Reviews from a etwork Perspective. Abstract
CS 9 Projet Final Report: Learning Convention Propagation in BeerAdvoate Reviews from a etwork Perspetive Abstrat We look at the way onventions propagate between reviews on the BeerAdvoate dataset, and
More informationDrawing lines. Naïve line drawing algorithm. drawpixel(x, round(y)); double dy = y1 - y0; double dx = x1 - x0; double m = dy / dx; double y = y0;
Naïve line drawing algorithm // Connet to grid points(x0,y0) and // (x1,y1) by a line. void drawline(int x0, int y0, int x1, int y1) { int x; double dy = y1 - y0; double dx = x1 - x0; double m = dy / dx;
More informationCleanUp: Improving Quadrilateral Finite Element Meshes
CleanUp: Improving Quadrilateral Finite Element Meshes Paul Kinney MD-10 ECC P.O. Box 203 Ford Motor Company Dearborn, MI. 8121 (313) 28-1228 pkinney@ford.om Abstrat: Unless an all quadrilateral (quad)
More informationAutomatic Physical Design Tuning: Workload as a Sequence Sanjay Agrawal Microsoft Research One Microsoft Way Redmond, WA, USA +1-(425)
Automati Physial Design Tuning: Workload as a Sequene Sanjay Agrawal Mirosoft Researh One Mirosoft Way Redmond, WA, USA +1-(425) 75-357 sagrawal@mirosoft.om Eri Chu * Computer Sienes Department University
More informationOn - Line Path Delay Fault Testing of Omega MINs M. Bellos 1, E. Kalligeros 1, D. Nikolos 1,2 & H. T. Vergos 1,2
On - Line Path Delay Fault Testing of Omega MINs M. Bellos, E. Kalligeros, D. Nikolos,2 & H. T. Vergos,2 Dept. of Computer Engineering and Informatis 2 Computer Tehnology Institute University of Patras,
More information13.1 Numerical Evaluation of Integrals Over One Dimension
13.1 Numerial Evaluation of Integrals Over One Dimension A. Purpose This olletion of subprograms estimates the value of the integral b a f(x) dx where the integrand f(x) and the limits a and b are supplied
More informationAn Event Display for ATLAS H8 Pixel Test Beam Data
An Event Display for ATLAS H8 Pixel Test Beam Data George Gollin Centre de Physique des Partiules de Marseille and University of Illinois April 17, 1999 g-gollin@uiu.edu An event display program is now
More informationChapter 2: Introduction to Maple V
Chapter 2: Introdution to Maple V 2-1 Working with Maple Worksheets Try It! (p. 15) Start a Maple session with an empty worksheet. The name of the worksheet should be Untitled (1). Use one of the standard
More informationA Novel Validity Index for Determination of the Optimal Number of Clusters
IEICE TRANS. INF. & SYST., VOL.E84 D, NO.2 FEBRUARY 2001 281 LETTER A Novel Validity Index for Determination of the Optimal Number of Clusters Do-Jong KIM, Yong-Woon PARK, and Dong-Jo PARK, Nonmembers
More informationKERNEL SPARSE REPRESENTATION WITH LOCAL PATTERNS FOR FACE RECOGNITION
KERNEL SPARSE REPRESENTATION WITH LOCAL PATTERNS FOR FACE RECOGNITION Cuiui Kang 1, Shengai Liao, Shiming Xiang 1, Chunhong Pan 1 1 National Laboratory of Pattern Reognition, Institute of Automation, Chinese
More informationAdapting K-Medians to Generate Normalized Cluster Centers
Adapting -Medians to Generate Normalized Cluster Centers Benamin J. Anderson, Deborah S. Gross, David R. Musiant Anna M. Ritz, Thomas G. Smith, Leah E. Steinberg Carleton College andersbe@gmail.om, {dgross,
More informationQuery Evaluation Overview. Query Optimization: Chap. 15. Evaluation Example. Cost Estimation. Query Blocks. Query Blocks
Query Evaluation Overview Query Optimization: Chap. 15 CS634 Leture 12 SQL query first translated to relational algebra (RA) Atually, some additional operators needed for SQL Tree of RA operators, with
More informationExploring the Commonality in Feature Modeling Notations
Exploring the Commonality in Feature Modeling Notations Miloslav ŠÍPKA Slovak University of Tehnology Faulty of Informatis and Information Tehnologies Ilkovičova 3, 842 16 Bratislava, Slovakia miloslav.sipka@gmail.om
More informationReading Object Code. A Visible/Z Lesson
Reading Objet Code A Visible/Z Lesson The Idea: When programming in a high-level language, we rarely have to think about the speifi ode that is generated for eah instrution by a ompiler. But as an assembly
More informationthe data. Structured Principal Component Analysis (SPCA)
Strutured Prinipal Component Analysis Kristin M. Branson and Sameer Agarwal Department of Computer Siene and Engineering University of California, San Diego La Jolla, CA 9193-114 Abstrat Many tasks involving
More informationReading Object Code. A Visible/Z Lesson
Reading Objet Code A Visible/Z Lesson The Idea: When programming in a high-level language, we rarely have to think about the speifi ode that is generated for eah instrution by a ompiler. But as an assembly
More informationGraph-Based vs Depth-Based Data Representation for Multiview Images
Graph-Based vs Depth-Based Data Representation for Multiview Images Thomas Maugey, Antonio Ortega, Pasal Frossard Signal Proessing Laboratory (LTS), Eole Polytehnique Fédérale de Lausanne (EPFL) Email:
More informationSelf-Adaptive Parent to Mean-Centric Recombination for Real-Parameter Optimization
Self-Adaptive Parent to Mean-Centri Reombination for Real-Parameter Optimization Kalyanmoy Deb and Himanshu Jain Department of Mehanial Engineering Indian Institute of Tehnology Kanpur Kanpur, PIN 86 {deb,hjain}@iitk.a.in
More information1. The collection of the vowels in the word probability. 2. The collection of real numbers that satisfy the equation x 9 = 0.
C HPTER 1 SETS I. DEFINITION OF SET We begin our study of probability with the disussion of the basi onept of set. We assume that there is a ommon understanding of what is meant by the notion of a olletion
More informationHEXA: Compact Data Structures for Faster Packet Processing
Washington University in St. Louis Washington University Open Sholarship All Computer Siene and Engineering Researh Computer Siene and Engineering Report Number: 27-26 27 HEXA: Compat Data Strutures for
More informationAn Alternative Approach to the Fuzzifier in Fuzzy Clustering to Obtain Better Clustering Results
An Alternative Approah to the Fuzziier in Fuzzy Clustering to Obtain Better Clustering Results Frank Klawonn Department o Computer Siene University o Applied Sienes BS/WF Salzdahlumer Str. 46/48 D-38302
More informationAnalysis of input and output configurations for use in four-valued CCD programmable logic arrays
nalysis of input and output onfigurations for use in four-valued D programmable logi arrays J.T. utler H.G. Kerkhoff ndexing terms: Logi, iruit theory and design, harge-oupled devies bstrat: s in binary,
More informationIntroduction to Seismology Spring 2008
MIT OpenCourseWare http://ow.mit.edu 1.510 Introdution to Seismology Spring 008 For information about iting these materials or our Terms of Use, visit: http://ow.mit.edu/terms. 1.510 Leture Notes 3.3.007
More informationAbstract. Key Words: Image Filters, Fuzzy Filters, Order Statistics Filters, Rank Ordered Mean Filters, Channel Noise. 1.
Fuzzy Weighted Rank Ordered Mean (FWROM) Filters for Mixed Noise Suppression from Images S. Meher, G. Panda, B. Majhi 3, M.R. Meher 4,,4 Department of Eletronis and I.E., National Institute of Tehnology,
More informationCross-layer Resource Allocation on Broadband Power Line Based on Novel QoS-priority Scheduling Function in MAC Layer
Communiations and Networ, 2013, 5, 69-73 http://dx.doi.org/10.4236/n.2013.53b2014 Published Online September 2013 (http://www.sirp.org/journal/n) Cross-layer Resoure Alloation on Broadband Power Line Based
More informationSAND Unlimited Release Printed November 1995 Updated November 29, :26 PM EXODUS II: A Finite Element Data Model
SAND92-2137 Unlimited Release Printed November 1995 Updated November 29, 2006 12:26 PM EXODUS II: A Finite Element Data Model Gregory D. Sjaardema (updated version) Larry A. Shoof, Vitor R. Yarberry Computational
More informationUCSB Math TI-85 Tutorials: Basics
3 UCSB Math TI-85 Tutorials: Basis If your alulator sreen doesn t show anything, try adjusting the ontrast aording to the instrutions on page 3, or page I-3, of the alulator manual You should read the
More informationGray Codes for Reflectable Languages
Gray Codes for Refletable Languages Yue Li Joe Sawada Marh 8, 2008 Abstrat We lassify a type of language alled a refletable language. We then develop a generi algorithm that an be used to list all strings
More informationIntroductory Programming, IMM, DTU Systematic Software Test. Software test (afprøvning) Motivation. Structural test and functional test
Introdutory Programming, IMM, DTU Systemati Software Test Peter Sestoft a Programs often ontain unintended errors how do you find them? Strutural test Funtional test Notes: Systemati Software Test, http://www.dina.kvl.dk/
More informationSystem-Level Parallelism and Throughput Optimization in Designing Reconfigurable Computing Applications
System-Level Parallelism and hroughput Optimization in Designing Reonfigurable Computing Appliations Esam El-Araby 1, Mohamed aher 1, Kris Gaj 2, arek El-Ghazawi 1, David Caliga 3, and Nikitas Alexandridis
More informationUsing Augmented Measurements to Improve the Convergence of ICP
Using Augmented Measurements to Improve the onvergene of IP Jaopo Serafin, Giorgio Grisetti Dept. of omputer, ontrol and Management Engineering, Sapienza University of Rome, Via Ariosto 25, I-0085, Rome,
More informationBSPLND, A B-Spline N-Dimensional Package for Scattered Data Interpolation
BSPLND, A B-Spline N-Dimensional Pakage for Sattered Data Interpolation Mihael P. Weis Traker Business Systems 85 Terminal Drive, Suite Rihland, WA 995-59-946-544 mike@vidian.net Robert R. Lewis Washington
More informationThe recursive decoupling method for solving tridiagonal linear systems
Loughborough University Institutional Repository The reursive deoupling method for solving tridiagonal linear systems This item was submitted to Loughborough University's Institutional Repository by the/an
More informationA Load-Balanced Clustering Protocol for Hierarchical Wireless Sensor Networks
International Journal of Advanes in Computer Networks and Its Seurity IJCNS A Load-Balaned Clustering Protool for Hierarhial Wireless Sensor Networks Mehdi Tarhani, Yousef S. Kavian, Saman Siavoshi, Ali
More informationPath Sharing and Predicate Evaluation for High-Performance XML Filtering*
Path Sharing and Prediate Evaluation for High-Performane XML Filtering Yanlei Diao, Mihael J. Franklin, Hao Zhang, Peter Fisher EECS, University of California, Berkeley {diaoyl, franklin, nhz, fisherp}@s.erkeley.edu
More informationThe Minimum Redundancy Maximum Relevance Approach to Building Sparse Support Vector Machines
The Minimum Redundany Maximum Relevane Approah to Building Sparse Support Vetor Mahines Xiaoxing Yang, Ke Tang, and Xin Yao, Nature Inspired Computation and Appliations Laboratory (NICAL), Shool of Computer
More informationApproximate logic synthesis for error tolerant applications
Approximate logi synthesis for error tolerant appliations Doohul Shin and Sandeep K. Gupta Eletrial Engineering Department, University of Southern California, Los Angeles, CA 989 {doohuls, sandeep}@us.edu
More informationFast Rigid Motion Segmentation via Incrementally-Complex Local Models
Fast Rigid Motion Segmentation via Inrementally-Complex Loal Models Fernando Flores-Mangas Allan D. Jepson Department of Computer Siene, University of Toronto {mangas,jepson}@s.toronto.edu Abstrat The
More informationAnd, the (low-pass) Butterworth filter of order m is given in the frequency domain by
Problem Set no.3.a) The ideal low-pass filter is given in the frequeny domain by B ideal ( f ), f f; =, f > f. () And, the (low-pass) Butterworth filter of order m is given in the frequeny domain by B
More informationPerformance Benchmarks for an Interactive Video-on-Demand System
Performane Benhmarks for an Interative Video-on-Demand System. Guo,P.G.Taylor,E.W.M.Wong,S.Chan,M.Zukerman andk.s.tang ARC Speial Researh Centre for Ultra-Broadband Information Networks (CUBIN) Department
More informationTorpedo Trajectory Visual Simulation Based on Nonlinear Backstepping Control
orpedo rajetory Visual Simulation Based on Nonlinear Bakstepping Control Peng Hai-jun 1, Li Hui-zhou Chen Ye 1, 1. Depart. of Weaponry Eng, Naval Univ. of Engineering, Wuhan 400, China. Depart. of Aeronautial
More informationA {k, n}-secret Sharing Scheme for Color Images
A {k, n}-seret Sharing Sheme for Color Images Rastislav Luka, Konstantinos N. Plataniotis, and Anastasios N. Venetsanopoulos The Edward S. Rogers Sr. Dept. of Eletrial and Computer Engineering, University
More informationPlot-to-track correlation in A-SMGCS using the target images from a Surface Movement Radar
Plot-to-trak orrelation in A-SMGCS using the target images from a Surfae Movement Radar G. Golino Radar & ehnology Division AMS, Italy ggolino@amsjv.it Abstrat he main topi of this paper is the formulation
More informationGradient based progressive probabilistic Hough transform
Gradient based progressive probabilisti Hough transform C.Galambos, J.Kittler and J.Matas Abstrat: The authors look at the benefits of exploiting gradient information to enhane the progressive probabilisti
More informationPHYS 3437: Computational Methods in Physics, Assignment 2
PHYS 3437: Computational Methods in Physis, Assignment 2 Set January 27th due Feb 26th NOTE: This assignment is potentially quite lengthy if you are urrently developing your programming skills. If so,
More informationCluster-Based Cumulative Ensembles
Cluster-Based Cumulative Ensembles Hanan G. Ayad and Mohamed S. Kamel Pattern Analysis and Mahine Intelligene Lab, Eletrial and Computer Engineering, University of Waterloo, Waterloo, Ontario N2L 3G1,
More informationSupplementary Material: Geometric Calibration of Micro-Lens-Based Light-Field Cameras using Line Features
Supplementary Material: Geometri Calibration of Miro-Lens-Based Light-Field Cameras using Line Features Yunsu Bok, Hae-Gon Jeon and In So Kweon KAIST, Korea As the supplementary material, we provide detailed
More informationDetecting Outliers in High-Dimensional Datasets with Mixed Attributes
Deteting Outliers in High-Dimensional Datasets with Mixed Attributes A. Koufakou, M. Georgiopoulos, and G.C. Anagnostopoulos 2 Shool of EECS, University of Central Florida, Orlando, FL, USA 2 Dept. of
More informationA Partial Sorting Algorithm in Multi-Hop Wireless Sensor Networks
A Partial Sorting Algorithm in Multi-Hop Wireless Sensor Networks Abouberine Ould Cheikhna Department of Computer Siene University of Piardie Jules Verne 80039 Amiens Frane Ould.heikhna.abouberine @u-piardie.fr
More informationAlgorithms, Mechanisms and Procedures for the Computer-aided Project Generation System
Algorithms, Mehanisms and Proedures for the Computer-aided Projet Generation System Anton O. Butko 1*, Aleksandr P. Briukhovetskii 2, Dmitry E. Grigoriev 2# and Konstantin S. Kalashnikov 3 1 Department
More information特集 Road Border Recognition Using FIR Images and LIDAR Signal Processing
デンソーテクニカルレビュー Vol. 15 2010 特集 Road Border Reognition Using FIR Images and LIDAR Signal Proessing 高木聖和 バーゼル ファルディ Kiyokazu TAKAGI Basel Fardi ヘンドリック ヴァイゲル Hendrik Weigel ゲルド ヴァニーリック Gerd Wanielik This paper
More informationA Novel Bit Level Time Series Representation with Implication of Similarity Search and Clustering
A Novel Bit Level Time Series Representation with Impliation of Similarity Searh and lustering hotirat Ratanamahatana, Eamonn Keogh, Anthony J. Bagnall 2, and Stefano Lonardi Dept. of omputer Siene & Engineering,
More informationVideo Data and Sonar Data: Real World Data Fusion Example
14th International Conferene on Information Fusion Chiago, Illinois, USA, July 5-8, 2011 Video Data and Sonar Data: Real World Data Fusion Example David W. Krout Applied Physis Lab dkrout@apl.washington.edu
More informationMulti-Channel Wireless Networks: Capacity and Protocols
Multi-Channel Wireless Networks: Capaity and Protools Tehnial Report April 2005 Pradeep Kyasanur Dept. of Computer Siene, and Coordinated Siene Laboratory, University of Illinois at Urbana-Champaign Email:
More informationNONLINEAR BACK PROJECTION FOR TOMOGRAPHIC IMAGE RECONSTRUCTION. Ken Sauer and Charles A. Bouman
NONLINEAR BACK PROJECTION FOR TOMOGRAPHIC IMAGE RECONSTRUCTION Ken Sauer and Charles A. Bouman Department of Eletrial Engineering, University of Notre Dame Notre Dame, IN 46556, (219) 631-6999 Shool of
More informationAlleviating DFT cost using testability driven HLS
Alleviating DFT ost using testability driven HLS M.L.Flottes, R.Pires, B.Rouzeyre Laboratoire d Informatique, de Robotique et de Miroéletronique de Montpellier, U.M. CNRS 5506 6 rue Ada, 34392 Montpellier
More informationEffecting Parallel Graph Eigensolvers Through Library Composition
Effeting Parallel Graph Eigensolvers Through Library Composition Alex Breuer, Peter Gottshling, Douglas Gregor, Andrew Lumsdaine Open Systems Laboratory Indiana University Bloomington, IN 47405 {abreuer,pgottsh,dgregor,lums@osl.iu.edu
More informationPreliminary investigation of multi-wavelet denoising in partial discharge detection
Preliminary investigation of multi-wavelet denoising in partial disharge detetion QIAN YONG Department of Eletrial Engineering Shanghai Jiaotong University 8#, Donghuan Road, Minhang distrit, Shanghai
More informationA Numerical optimization technique for the design of airfoils in viscous flows
Rohester Institute of Tehnology RIT Sholar Works Theses Thesis/Dissertation Colletions 7-1-1995 A Numerial optimization tehnique for the design of airfoils in visous flows Robert MaNeill Follow this and
More informationDr.Hazeem Al-Khafaji Dept. of Computer Science, Thi-Qar University, College of Science, Iraq
Volume 4 Issue 6 June 014 ISSN: 77 18X International Journal of Advaned Researh in Computer Siene and Software Engineering Researh Paper Available online at: www.ijarsse.om Medial Image Compression using
More informationMultiple-Criteria Decision Analysis: A Novel Rank Aggregation Method
3537 Multiple-Criteria Deision Analysis: A Novel Rank Aggregation Method Derya Yiltas-Kaplan Department of Computer Engineering, Istanbul University, 34320, Avilar, Istanbul, Turkey Email: dyiltas@ istanbul.edu.tr
More informationOptimization of Two-Stage Cylindrical Gear Reducer with Adaptive Boundary Constraints
5 JOURNAL OF SOFTWARE VOL. 8 NO. 8 AUGUST Optimization of Two-Stage Cylindrial Gear Reduer with Adaptive Boundary Constraints Xueyi Li College of Mehanial and Eletroni Engineering Shandong University of
More informationNew Fuzzy Object Segmentation Algorithm for Video Sequences *
JOURNAL OF INFORMATION SCIENCE AND ENGINEERING 24, 521-537 (2008) New Fuzzy Obet Segmentation Algorithm for Video Sequenes * KUO-LIANG CHUNG, SHIH-WEI YU, HSUEH-JU YEH, YONG-HUAI HUANG AND TA-JEN YAO Department
More informationParallelization and Performance of 3D Ultrasound Imaging Beamforming Algorithms on Modern Clusters
Parallelization and Performane of 3D Ultrasound Imaging Beamforming Algorithms on Modern Clusters F. Zhang, A. Bilas, A. Dhanantwari, K.N. Plataniotis, R. Abiprojo, and S. Stergiopoulos Dept. of Eletrial
More informationRotation Invariant Spherical Harmonic Representation of 3D Shape Descriptors
Eurographis Symposium on Geometry Proessing (003) L. Kobbelt, P. Shröder, H. Hoppe (Editors) Rotation Invariant Spherial Harmoni Representation of 3D Shape Desriptors Mihael Kazhdan, Thomas Funkhouser,
More informationAlgorithms for External Memory Lecture 6 Graph Algorithms - Weighted List Ranking
Algorithms for External Memory Leture 6 Graph Algorithms - Weighted List Ranking Leturer: Nodari Sithinava Sribe: Andi Hellmund, Simon Ohsenreither 1 Introdution & Motivation After talking about I/O-effiient
More informationLocality-Sensitive Hashing Scheme Based on p-stable Distributions
Loality-Sensitive Hashing Sheme Based on p-stable Distributions Mayur Datar Department of Computer Siene, Stanford University datar@s.stanford.edu Piotr Indyk Laboratory for Computer Siene, MIT indyk@theory.ls.mit.edu
More informationSequential Incremental-Value Auctions
Sequential Inremental-Value Autions Xiaoming Zheng and Sven Koenig Department of Computer Siene University of Southern California Los Angeles, CA 90089-0781 {xiaominz,skoenig}@us.edu Abstrat We study the
More informationCapturing Large Intra-class Variations of Biometric Data by Template Co-updating
Capturing Large Intra-lass Variations of Biometri Data by Template Co-updating Ajita Rattani University of Cagliari Piazza d'armi, Cagliari, Italy ajita.rattani@diee.unia.it Gian Lua Marialis University
More informationThis fact makes it difficult to evaluate the cost function to be minimized
RSOURC LLOCTION N SSINMNT In the resoure alloation step the amount of resoures required to exeute the different types of proesses is determined. We will refer to the time interval during whih a proess
More informationA scheme for racquet sports video analysis with the combination of audio-visual information
A sheme for raquet sports video analysis with the ombination of audio-visual information Liyuan Xing a*, Qixiang Ye b, Weigang Zhang, Qingming Huang a and Hua Yu a a Graduate Shool of the Chinese Aadamy
More informationAcoustic Links. Maximizing Channel Utilization for Underwater
Maximizing Channel Utilization for Underwater Aousti Links Albert F Hairris III Davide G. B. Meneghetti Adihele Zorzi Department of Information Engineering University of Padova, Italy Email: {harris,davide.meneghetti,zorzi}@dei.unipd.it
More information3-D IMAGE MODELS AND COMPRESSION - SYNTHETIC HYBRID OR NATURAL FIT?
3-D IMAGE MODELS AND COMPRESSION - SYNTHETIC HYBRID OR NATURAL FIT? Bernd Girod, Peter Eisert, Marus Magnor, Ekehard Steinbah, Thomas Wiegand Te {girod eommuniations Laboratory, University of Erlangen-Nuremberg
More informationASSESSMENT OF TWO CHEAP CLOSE-RANGE FEATURE EXTRACTION SYSTEMS
ASSESSMENT OF TWO CHEAP CLOSE-RANGE FEATURE EXTRACTION SYSTEMS Ahmed Elaksher a, Mohammed Elghazali b, Ashraf Sayed b, and Yasser Elmanadilli b a Shool of Civil Engineering, Purdue University, West Lafayette,
More informationwith respect to the normal in each medium, respectively. The question is: How are θ
Prof. Raghuveer Parthasarathy University of Oregon Physis 35 Winter 8 3 R EFRACTION When light travels from one medium to another, it may hange diretion. This phenomenon familiar whenever we see the bent
More informationFinding the Equation of a Straight Line
Finding the Equation of a Straight Line You should have, before now, ome aross the equation of a straight line, perhaps at shool. Engineers use this equation to help determine how one quantity is related
More informationWeak Dependence on Initialization in Mixture of Linear Regressions
Proeedings of the International MultiConferene of Engineers and Computer Sientists 8 Vol I IMECS 8, Marh -6, 8, Hong Kong Weak Dependene on Initialization in Mixture of Linear Regressions Ryohei Nakano
More informationAllocating Rotating Registers by Scheduling
Alloating Rotating Registers by Sheduling Hongbo Rong Hyunhul Park Cheng Wang Youfeng Wu Programming Systems Lab Intel Labs {hongbo.rong,hyunhul.park,heng..wang,youfeng.wu}@intel.om ABSTRACT A rotating
More informationtimestamp, if silhouette(x, y) 0 0 if silhouette(x, y) = 0, mhi(x, y) = and mhi(x, y) < timestamp - duration mhi(x, y), else
3rd International Conferene on Multimedia Tehnolog(ICMT 013) An Effiient Moving Target Traking Strateg Based on OpenCV and CAMShift Theor Dongu Li 1 Abstrat Image movement involved bakground movement and
More informationWe P9 16 Eigenray Tracing in 3D Heterogeneous Media
We P9 Eigenray Traing in 3D Heterogeneous Media Z. Koren* (Emerson), I. Ravve (Emerson) Summary Conventional two-point ray traing in a general 3D heterogeneous medium is normally performed by a shooting
More informationFOREGROUND OBJECT EXTRACTION USING FUZZY C MEANS WITH BIT-PLANE SLICING AND OPTICAL FLOW
FOREGROUND OBJECT EXTRACTION USING FUZZY C EANS WITH BIT-PLANE SLICING AND OPTICAL FLOW SIVAGAI., REVATHI.T, JEGANATHAN.L 3 APSG, SCSE, VIT University, Chennai, India JRF, DST, Dehi, India. 3 Professor,
More informationContents Contents...I List of Tables...VIII List of Figures...IX 1. Introduction Information Retrieval... 8
Contents Contents...I List of Tables...VIII List of Figures...IX 1. Introdution... 1 1.1. Internet Information...2 1.2. Internet Information Retrieval...3 1.2.1. Doument Indexing...4 1.2.2. Doument Retrieval...4
More informationOne Against One or One Against All : Which One is Better for Handwriting Recognition with SVMs?
One Against One or One Against All : Whih One is Better for Handwriting Reognition with SVMs? Jonathan Milgram, Mohamed Cheriet, Robert Sabourin To ite this version: Jonathan Milgram, Mohamed Cheriet,
More informationChemical, Biological and Radiological Hazard Assessment: A New Model of a Plume in a Complex Urban Environment
hemial, Biologial and Radiologial Haard Assessment: A New Model of a Plume in a omplex Urban Environment Skvortsov, A.T., P.D. Dawson, M.D. Roberts and R.M. Gailis HPP Division, Defene Siene and Tehnology
More informationCOST PERFORMANCE ASPECTS OF CCD FAST AUXILIARY MEMORY
COST PERFORMANCE ASPECTS OF CCD FAST AUXILIARY MEMORY Dileep P, Bhondarkor Texas Instruments Inorporated Dallas, Texas ABSTRACT Charge oupled devies (CCD's) hove been mentioned as potential fast auxiliary
More informationA Fast Kernel-based Multilevel Algorithm for Graph Clustering
A Fast Kernel-based Multilevel Algorithm for Graph Clustering Inderjit Dhillon Dept. of Computer Sienes University of Texas at Austin Austin, TX 78712 inderjit@s.utexas.edu Yuqiang Guan Dept. of Computer
More informationFacility Location: Distributed Approximation
Faility Loation: Distributed Approximation Thomas Mosibroda Roger Wattenhofer Distributed Computing Group PODC 2005 Where to plae ahes in the Internet? A distributed appliation that has to dynamially plae
More informationFitting conics to paracatadioptric projections of lines
Computer Vision and Image Understanding 11 (6) 11 16 www.elsevier.om/loate/viu Fitting onis to paraatadioptri projetions of lines João P. Barreto *, Helder Araujo Institute for Systems and Robotis, Department
More informationCS:APP2e Web Aside ASM:X87: X87-Based Support for Floating Point
CS:APP2e Web Aside ASM:X87: X87-Based Support for Floating Point Randal E. Bryant David R. O Hallaron June 5, 2012 Notie The material in this doument is supplementary material to the book Computer Systems,
More informationConnection Guide. Installing the printer locally (Windows) What is local printing? Installing the printer using the Software and Documentation CD
Page 1 of 7 Connetion Guide Installing the printer loally (Windows) Note: When installing a loally attahed printer, if the operating system is not supported y the Software and Doumentation CD, then the
More informationMeasurement of the stereoscopic rangefinder beam angular velocity using the digital image processing method
Measurement of the stereosopi rangefinder beam angular veloity using the digital image proessing method ROMAN VÍTEK Department of weapons and ammunition University of defense Kouniova 65, 62 Brno CZECH
More informationA radiometric analysis of projected sinusoidal illumination for opaque surfaces
University of Virginia tehnial report CS-21-7 aompanying A Coaxial Optial Sanner for Synhronous Aquisition of 3D Geometry and Surfae Refletane A radiometri analysis of projeted sinusoidal illumination
More informationDivide-and-conquer algorithms 1
* 1 Multipliation Divide-and-onquer algorithms 1 The mathematiian Gauss one notied that although the produt of two omplex numbers seems to! involve four real-number multipliations it an in fat be done
More information