AST 911: Numerical Techniques in Astronomy

Spring 2008
Assignment Codes
Arithmetic at finite precision test_int.f [output]
test_machine.f [output when compiled with option -fdefaut-real-8 ]
Finding roots roots.f
test_roots.f [ output ]
ODE basics rk2.f [ rms error for fixed no. of steps
rms error for integration over fixed interval ]
Stiff ODE's

fEuler_stiff.f [plots for h=0.05, h=0.2 ]
bEuler_stiff.f [plots for h=0.2, h=0.8 ]

Performance loop_order.f
matmul.f
cache.f
Advection

CDadvect.f [movie]
UDadvect.f [movie]

Project Codes
Equation of state  
A white dwarf model lane_emden.f, rk4.f, zbrent.f and a plot showing the error in the computed radius for order n = 1. As advertised the error scales as h4.
Flame crank_nicholson.f, driver.f, print_matrix.f, reaction.f, test_cn.f, test_flame.f
movie of flame propagation