# Register now to access 7 million high quality study materials (What's Course Hero?) Course Hero is the premier provider of high quality online educational resources. With millions of study documents, online tutors, digital flashcards and free courseware, Course Hero is helping students learn more efficiently and effectively. Whether you're interested in exploring new subjects or mastering key topics for your next exam, Course Hero has the tools you need to achieve your goals.

5 Pages

### class06

Course: ASTR 415, Fall 2008
School: Maryland
Rating:

Word Count: 905

#### Document Preview

6. Class Numerical Linear Algebra, Part 1 Probably the simplest kind of problem. Occurs in many contexts, often as part of larger problem. Symbolic manipulation packages can do linear algebra analytically (e.g., Mathematica, Maple, etc.). Numerical methods needed when: Number of equations very large. One or more coecients numerical. Linear Systems Write linear system as: a11 x1 + a12 x2 + + a1n xn = b1...

Register Now

#### Unformatted Document Excerpt

Coursehero >> Maryland >> Maryland >> ASTR 415

Course Hero has millions of student submitted documents similar to the one
below including study guides, practice problems, reference materials, practice exams, textbook help and tutor support.

Course Hero has millions of student submitted documents similar to the one below including study guides, practice problems, reference materials, practice exams, textbook help and tutor support.

Find millions of documents on Course Hero - Study Guides, Lecture Notes, Reference Materials, Practice Exams and more. Course Hero has millions of course specific materials providing students with the best way to expand their education.

Below is a small sample set of documents:

Maryland - ASTR - 07
Class 7. Numerical Linear Algebra, Part 2LU Decomposition Suppose we can write A as a product of two matrices: A = LU, where L is lower triangular and U is upper triangular: 00 L= 0 U= 0 00 Then Ax = (LU) x = L (Ux) = b, i.e.,
Maryland - ASTR - 415
Class 7. Numerical Linear Algebra, Part 2LU Decomposition Suppose we can write A as a product of two matrices: A = LU, where L is lower triangular and U is upper triangular: 00 L= 0 U= 0 00 Then Ax = (LU) x = L (Ux) = b, i.e.,
Maryland - ASTR - 08
Class 8. Root Finding in 1-DNonlinear Equations Often (most of the time?) the relevant system of equations is nonlinear in the unknowns. Then, cannot decompose as Ax = b. Oh well. Instead write as: 1. f (x) = 0 (for functions of one variable, i.e
Maryland - ASTR - 415
Class 8. Root Finding in 1-DNonlinear Equations Often (most of the time?) the relevant system of equations is nonlinear in the unknowns. Then, cannot decompose as Ax = b. Oh well. Instead write as: 1. f (x) = 0 (for functions of one variable, i.e
Maryland - ASTR - 09
Class 9. Root Finding in Multi-D, and Numerical DifferentiationNonlinear Systems of Equations Consider the system f (x, y) = 0, g(x, y) = 0. Plot zero contours of f and g:
Maryland - ASTR - 415
Class 9. Root Finding in Multi-D, and Numerical DifferentiationNonlinear Systems of Equations Consider the system f (x, y) = 0, g(x, y) = 0. Plot zero contours of f and g:
Maryland - ASTR - 10
Class 10. Statistics and the K-S TestStatistical Description of Data Cf. NRiC 14. Statistics provides tools for understanding data. In the wrong hands these tools can be dangerous! Heres a typical data analysis cycle: 1. Apply some formula to da
Maryland - ASTR - 415
Class 10. Statistics and the K-S TestStatistical Description of Data Cf. NRiC 14. Statistics provides tools for understanding data. In the wrong hands these tools can be dangerous! Heres a typical data analysis cycle: 1. Apply some formula to da
Maryland - ASTR - 11
Class 11. Modeling of Data NRiC 15. Model depends on adjustable parameters. Can be used for constrained interpolation. Basic approach: 1. Choose gure-of-merit function (e.g., 2 ). 2. Adjust best-t parameters: minimize merit function. 3. Compute g
Maryland - ASTR - 415
Class 11. Modeling of Data NRiC 15. Model depends on adjustable parameters. Can be used for constrained interpolation. Basic approach: 1. Choose gure-of-merit function (e.g., 2 ). 2. Adjust best-t parameters: minimize merit function. 3. Compute g
Maryland - ASTR - 12
Class 12. Random Numbers NRiC 7. Frequently needed to generate initial conditions. Often used to solve problems statistically. How can a computer generate a random number? It cant! Generators are pseudo-random. Generators are deterministic: its
Maryland - ASTR - 415
Class 12. Random Numbers NRiC 7. Frequently needed to generate initial conditions. Often used to solve problems statistically. How can a computer generate a random number? It cant! Generators are pseudo-random. Generators are deterministic: its
Maryland - ASTR - 13
Class 13. Numerical IntegrationSimple Monte Carlo Integration (NRiC 7.6) Can use RNGs to estimate integrals. Suppose we pick n random points x1 , ., xN uniformly in a multi-D volume V . Basic theorem of Monte Carlo integration: f dV 1 NNVV &lt;
Maryland - ASTR - 415
Class 13. Numerical IntegrationSimple Monte Carlo Integration (NRiC 7.6) Can use RNGs to estimate integrals. Suppose we pick n random points x1 , ., xN uniformly in a multi-D volume V . Basic theorem of Monte Carlo integration: f dV 1 NNVV &lt;
Maryland - ASTR - 14
Class 14. Ordinary Dierential Equations NRiC 16. ODEs involve derivatives with respect to one independent variable, e.g., time t. ODEs can always be reduced to a set of rst-order equations (i.e., involving only rst derivatives). E.g., d2 y dy + b(
Maryland - ASTR - 415
Class 14. Ordinary Dierential Equations NRiC 16. ODEs involve derivatives with respect to one independent variable, e.g., time t. ODEs can always be reduced to a set of rst-order equations (i.e., involving only rst derivatives). E.g., d2 y dy + b(
Maryland - ASTR - 15
Class 15. ODEs, Part 2The Leapfrog Integrator Very useful for second-order DEs in which d2 x/dt2 = f(x), e.g., SHM, N-body, etc. NOTE: Now dropping the prime ( ) from f. Suppose x is position, so d2 x/dt2 is acceleration. Procedure: dene v = dx/
Maryland - ASTR - 415
Class 15. ODEs, Part 2The Leapfrog Integrator Very useful for second-order DEs in which d2 x/dt2 = f(x), e.g., SHM, N-body, etc. NOTE: Now dropping the prime ( ) from f. Suppose x is position, so d2 x/dt2 is acceleration. Procedure: dene v = dx/
Maryland - ASTR - 16
Class 16. ODEs, Part 3Sti ODEs A system of more than one ODE is sti if solutions vary on two or more widely disparate lengthscales. E.g., y = 100y. General solution: y = Ae10x + Be+10x . Suppose BCs are y(0) = 1, y (0) = 10. Then B = 0, i.e., pur
Maryland - ASTR - 415
Class 16. ODEs, Part 3Sti ODEs A system of more than one ODE is sti if solutions vary on two or more widely disparate lengthscales. E.g., y = 100y. General solution: y = Ae10x + Be+10x . Suppose BCs are y(0) = 1, y (0) = 10. Then B = 0, i.e., pur
Maryland - ASTR - 17
Class 17. ODEs, Part 42-pt BVP: Shooting MethodProcedure: 1. At x1 , must specify N starting values for yi , i = 1, ., N. n1 values given by BC at x1 . n2 = N n1 values can be freely chosen. 2. Represent the free values as a vector V of dimensi
Maryland - ASTR - 415
Class 17. ODEs, Part 42-pt BVP: Shooting MethodProcedure: 1. At x1 , must specify N starting values for yi , i = 1, ., N. n1 values given by BC at x1 . n2 = N n1 values can be freely chosen. 2. Represent the free values as a vector V of dimensi
Maryland - ASTR - 18
Class 18. N -body Techniques, Part 1The N -body Problem Study of the dynamics of interacting particles, usually involving mutual forces. E.g., Application Mutual Force gravity stellar dynamics, planetesimals QM molecular dynamics, solid-state physi
Maryland - ASTR - 415
Class 18. N -body Techniques, Part 1The N -body Problem Study of the dynamics of interacting particles, usually involving mutual forces. E.g., Application Mutual Force gravity stellar dynamics, planetesimals QM molecular dynamics, solid-state physi
Maryland - ASTR - 19
Class 19. N -body Techniques, Part 2Time-integration Schemes Clearly, Newtons laws are IVP. Could use any method (Euler, RK4, etc.). But, issue is to balance accuracy vs. eciency. Typically need many particles to capture dynamics correctly (e.g.,
Maryland - ASTR - 415
Class 19. N -body Techniques, Part 2Time-integration Schemes Clearly, Newtons laws are IVP. Could use any method (Euler, RK4, etc.). But, issue is to balance accuracy vs. eciency. Typically need many particles to capture dynamics correctly (e.g.,
Maryland - ASTR - 20
Class 20. N -body Techniques, Part 3The PM Method, ContinuedThere are several distinct steps in PM process: 1. Assign particles to mesh to compute i . 2. Get boundary conditions for (0 and N +1 ). 3. Solve discretized version of Poissons equation.
Maryland - ASTR - 415
Class 20. N -body Techniques, Part 3The PM Method, ContinuedThere are several distinct steps in PM process: 1. Assign particles to mesh to compute i . 2. Get boundary conditions for (0 and N +1 ). 3. Solve discretized version of Poissons equation.
Maryland - ASTR - 21
Class 21. N -body Techniques, Part 4Tree CodesEciency can be increased by grouping particles together: Nearest particles exert greatest forces direct summation. Distant particles exert smallest forces treat in groups.Treat distant particles as
Maryland - ASTR - 415
Class 21. N -body Techniques, Part 4Tree CodesEciency can be increased by grouping particles together: Nearest particles exert greatest forces direct summation. Distant particles exert smallest forces treat in groups.Treat distant particles as
Maryland - ASTR - 22
Class 22. PDEs, Part 1 Cf. NRiC 19.Classication of PDEs A PDE is simply a dierential equation of more than one variable (so an ODE is a special case of a PDE). PDEs are usually classied into three types: 1. Hyperbolic (second or rst order in time
Maryland - ASTR - 415
Class 22. PDEs, Part 1 Cf. NRiC 19.Classication of PDEs A PDE is simply a dierential equation of more than one variable (so an ODE is a special case of a PDE). PDEs are usually classied into three types: 1. Hyperbolic (second or rst order in time
Maryland - ASTR - 23
Class 23. PDEs, Part 2Solving Hyperbolic PDEs, ContinuedUpwind dierencing In addition to amplitude errors (instability or damping), scheme may also have phase errors (dispersion) or transport errors (spurious transport of information). Upwind die
Maryland - ASTR - 415
Class 23. PDEs, Part 2Solving Hyperbolic PDEs, ContinuedUpwind dierencing In addition to amplitude errors (instability or damping), scheme may also have phase errors (dispersion) or transport errors (spurious transport of information). Upwind die
Maryland - ASTR - 24
Class 24. Fluid Dynamics, Part 1 The equations of uid dynamics are coupled PDEs that form an IVP (hyperbolic). Use the techniques described so far, plus additions.Fluid Dynamics in Astrophysics Whenever mean free path problem scale L in a plasm
Maryland - ASTR - 415
Class 24. Fluid Dynamics, Part 1 The equations of uid dynamics are coupled PDEs that form an IVP (hyperbolic). Use the techniques described so far, plus additions.Fluid Dynamics in Astrophysics Whenever mean free path problem scale L in a plasm
Maryland - ASTR - 415
ASTR415 Survey ResultsSpring 2007 11 respondents 1. Computer familiarity [1=Master, 5=None]: Avg = 2.7 [Skilled], Min = 1 [Master], Max = 4 [Novice] 2. Unix familiarity [1=Master, 5=None]: Avg = 3.6 [Skilled], Min = 2 [Expert], Max = 5 [None] 3. Uni
Maryland - ASTR - 415
ASTR415 Spring 2007Due May 08, 2007Term Project1) For your term project you will install and learn how to use a freely available (opensource) 3D visualization tool based on Open-GL. You will write a short report and present in class the results
Maryland - ASTR - 415
| nsttuustBj1d4tssBnn nn Xdss &amp; Xs T a ` C F U IE W P U Rs y ` EY ' W F I 3 w i w w Te T UEY TYE S I P T C W I W p m UE H U T ` y H T y UE XVeiGu(uo&amp;d(Gvtb6dQ&amp;bc(b(xc iqgGgGeb(uu6(xn(g(BxoVXcgGY FE C
Maryland - ASTR - 415
x i e r | z i g j r v fdV i n rv G} ofe } i im f XGl Xl G im hl m f X u e i} g G i} g im g f fe ov i m v i u IfX u e ofe G i m n m i Vn} } i Vn} G oVe } } X i ofw iml g ine m l g l e iV } ne } dFVG r v m v f } | zx
Maryland - ASTR - 415
4Xeex xe e7jj \$v )e `j B I Va U b V V U R I F YH Y W V U W R H c W V V F Y WH Q V I Pgp1eCd)wss`Xgpwhe`Xe1wPX)1T Ia R Y W R S U b V Y I W R U U R S Y I R U V F Y I V H b U U T I Va b Y r WH Raa Rt V F D
Maryland - ASTR - 415
ICadafWWaHH\$CHIaRp\$a\$wraWCawHa 9 i i S x i U bIxpU B @ e F k i 6 d B wP v q D A B FP PT B 8 6 qT Q v D 4 B v 4 @ D B F f D F B A @ A B F 6P k Uif'fbS7)UdRT Y rtuCaCdRR99R7\$17t\$ICdP Y `CaC9ECGR'\$i S
Maryland - ASTR - 415
A Crash Course on UNIXUNIX is an &quot;operating system&quot;.Interface between user and data stored on computer. A Windows-style interface is not required. Many flavors of UNIX (and windows interfaces).Solaris, Mandrake, RedHat (fvwm, Gnome, KDE)
Maryland - ASTR - 415
Data RepresentationsComputers store data as different variable types, e.g. integer, floating point, complex, etc. Different machines have different wordlengths, e.g. 4-byte ints on a 32-bit machine (Pentium), 8-byte ints on a 64-bit machine (Alpha).
Maryland - ASTR - 415
VisualizationVisualization is useful for:1) Data 2) Codeentry (initial conditions) debugging and performance analysis and display of results3) Interpretation Our focus will be #3. The computational astrophysicist can either:1) Develop 2) U
Maryland - ASTR - 415
Numerical Linear Algebra Probably the simplest kind of problem. Occurs in many contexts, often as part of larger problem. Symbolic manipulation packages can do linear algebra &quot;analytically&quot; (e.g. Mathematica, Maple). Numerical methods needed when:
Maryland - ASTR - 415
Nonlinear EquationsOften (most of the time?) the relevant system of equations is not linear in the unknowns. Then, cannot decompose as Ax = b. Oh well. Instead write as:(1) (2) f(x) = 0 f(x) = 0function of one variable (1-D) x = (x1,x2,.,xn
Maryland - ASTR - 415
Statistical Description of Data Cf. NRiC, Chapter 14. Statistics provides tools for understanding data.In the wrong hands these tools can be dangerous! Apply some formula to data to compute a &quot;statistic&quot;. Find where value falls in a probability
Maryland - ASTR - 415
Modeling of Data NRiC Chapter 15. Model depends on adjustable parameters. Can be used for &quot;constrained interpolation&quot;. Basic approach:1. 2. 3. 4.Choose figure-of-merit function (e.g. 2). Adjust best-fit parameters: minimize merit function. Co
Maryland - ASTR - 415
Random Numbers NRiC Chapter 7. Frequently needed to generate initial conditions. Often used to solve problems statistically. How can a computer generate a random number? It can't! Generators are pseudo-random. Generators are deterministic: i
Maryland - ASTR - 415
Numerical Integration (Quadrature) NRiC Chapter 4. Already seen Monte Carlo integration. Can cast problem as a differential equation (DE): = is equivalent to solving for I y(b) the DE dy/dx = f(x) with the boundary condition (BC) y(a) =
Maryland - ASTR - 415
Ordinary Differential Equations (ODEs) NRiC Chapter 16. ODEs involve derivatives wrt one independent variable, e.g. time t. ODEs can always be reduced to a set of firstorder equations (involving only first derivatives).e.g. = is
Maryland - MATH - 241
MATH 241 CALCULUS III FIRST MIDTERM EXAM Instructions. Answer each question on a separate answer sheet. Show all your work. Be sure your name, section number, and problem number are on each answer sheet, and that you have copied and signed the honor
Maryland - MATH - 241
MATH 241 CALCULUS III FIRST MIDTERM EXAM SOLUTIONS (1) For this problem, u = + k, and v = + 2 + 3 k. ij i j (a) u v = 1 + 2 3 = 2. k ij (b) u v = det 1 1 1 = 5 2 + 3 k. i j 1 2 3 (c) The symmetric form of the equations are: 2 x =
Maryland - ASTR - 601
Astronomy 601 - Fall 2005 Radiative ProcessesInstructor Prof. Massimo Ricotti Oce: CSS 0213 E-mail: ricotti@astro.umd.edu Phone: (301) 405 5097 Oce hours: by appointment Class web page: http:/www.astro.umd.edu/ricotti/NEWWEB/teaching/ASTR601.html Sc
Maryland - ECE - 2003
Center for Satellite and Hybrid Communication NetworksIntegrated Security Services for Dynamic Coalition ManagementHimanshu Khurana and Vijay Bharadwaj Electrical and Computer Engineering Department, University of Maryland College Park, Maryland 2
Maryland - ECE - 2003
Proceedings of the DARPA Information Survivability Conference and Exposition (DISCEX03) 0-7695-1897-4/03 \$17.00 2003 IEEEProceedings of the DARPA Information Survivability Conference and Exposition (DISCEX03) 0-7695-1897-4/03 \$17.00 2003 IEEEPr
Maryland - ECE - 2003
Integrated Security Services for Dynamic CoalitionsHimanshu Khurana1, Serban Gavrila1, Rakeshbabu Bobba, Radostina Koleva, Anuja Sonalker, Emilian Dinu, Virgil Gligor, and John Baras Electrical and Computer Engineering Department, University of Mary
Maryland - ECE - 2003
Towards Automated Negotiation of Access Control PoliciesVijay G. Bharadwaj and John S. Baras Institute for Systems Research, University of Maryland, College Park MD 20742, USA. vgb,baras @umd.eduAbstractWe examine the problem of negotiating acces
Maryland - ECE - 2003
DYNAMIC ADAPTATION OF ACCESS CONTROL POLICIESVijay Bharadwaj and John Baras Institute for Systems Research University of Maryland College Park MD 20742ABSTRACTWe describe an architecture and algorithms for deriving an access control policy by com
Maryland - ECE - 2003
Center for Satellite and Hybrid Communication NetworksIntegrated Security Services for Dynamic Coalition ManagementHimanshu Khurana Electrical and Computer Engineering Department, University of Maryland College Park, Maryland 20742 DARPA DC PI Mee