]> Generate Jacobian from FullBorn forward matrix

Generate Jacobian from FullBorn Forward Matrix

genJacobianFrom2pts() calculate the Jacobian matrix

J_{i j} = \frac{\partial φ_{i}}{\partial μ (r_{j})}

from the FullBorn forward matrix $A$ . The derivative is taken with respect to the optical properties $μ (r_{j})$ at the voxel $r_{j} .$

Function Summary

Syntax:	J = genJacobianFrom2pts(SD, Medium, MeasList, Phi2pt, PhiTot, Method, muvec);
Inputs:	SD	SD structure
	Medium	Medium structure
	MeasList	Measurement List that corresponds to experimental data
	Phi2pt	Returned by genBornData
	PhiTot	Returned by genBornData
	Method	How the Jacobian should be calculated, "Born" or "Rytov"
	muvec	Optical properties used to generate Phi2pt
Outputs:	J	The Jacobian

Detailed Descriptions

While the FullBorn calculation is the full solution to the Dyson equation (and thus represents perturbation theory to all orders), the Jacobian is calculated for small variations at a single voxel and thus still involves a first-Born (or Rytov-like) approximation.