{[ promptMessage ]}

Bookmark it

{[ promptMessage ]}

Biocompilerexercise

# Biocompilerexercise - Elia Ahadi 585.634.31 Due...

This preview shows pages 1–3. Sign up to view the full content.

Elia Ahadi 585.634.31 Due 03/05/09 Bio-Photonics Time Resolved diffusion Theory Compiler This exercise was to show how to compile the time resolved fluence rate referening http://omlc.ogi.edu/classroom/ece532/class5/trdt.c /******************************************** * trdt.c * Time-resolved diffusion theory * * by Steven L. Jacques, Dec. 28, 1998. * updated July 8, 1999. * * Computes the time-resolved fluence rate F(r,t) [W/cm^2] * in response to an impulse istropic point source Uo(t) [J] * deposited at time zero at the origin * within a uniform unbounded medium with optical properties of * absorption, reduced scattering, and refractive index. * * Creates output file "trdt.out" which records snapshots of * F(r,t) vs r at each of 3 timepoints. * **********/ #include <math.h> #include <stdio.h> #define PI 3.1415926 #define CC 2.997925e10 int main() { /* problem variables */ double mua, mus, g, nt; /* optical properties */ double D, delta; /* transport parameters */ double c; /* speed of light in medium */ double r; /* distance from source to observation point */ double t[3]; /* times of snapshots [s] */ double Uo; /* impulse energy of isotropic source point [J] */ double T; /* transport to observation point [1/(cm^2 s)] */ double F; /* fluence rate at observation point [W/cm^2] */ /* other variables */ double N; /* number of points in array */ double rmax; /* maximum position considered [cm] */ double dr; /* incremental step size [cm] */ short i, j; /* dummy indices */ double b, u; /* dummy variables */ FILE* target; /* pointer to output file */ /**** USER SPECIFIED VALUES *****/

This preview has intentionally blurred sections. Sign up to view the full version.

View Full Document
mua = 0.1; /* absorption coeff. [cm^-1] */ mus = 100.0; /* scattering coeff. [cm^-1] */ g = 0.90; /* anistoropy [-] */ nt = 1.33; /* refractive index of medium (eg., biological tissue) [-] */ Uo = 1; /* energy of impulse [J] */ rmax = 2; /* max distance between source and observation point [cm] */ N = 100; /* set number of points in array */ t[0] = 5e-12;
This is the end of the preview. Sign up to access the rest of the document.

{[ snackBarMessage ]}

### Page1 / 6

Biocompilerexercise - Elia Ahadi 585.634.31 Due...

This preview shows document pages 1 - 3. Sign up to view the full document.

View Full Document
Ask a homework question - tutors are online