Jacobian

Here we will calculate the log of the Jacobian [Graphics:../Images/index_gr_747.gif].

First, we obtain the expression of [Graphics:../Images/index_gr_748.gif]=[Graphics:../Images/index_gr_749.gif].

foo23 =[Graphics:../Images/index_gr_750.gif]

[Graphics:../Images/index_gr_751.gif]
[Graphics:../Images/index_gr_752.gif]

foo24 =[Graphics:../Images/index_gr_753.gif]

[Graphics:../Images/index_gr_754.gif]
[Graphics:../Images/index_gr_755.gif]

foo25 =[Graphics:../Images/index_gr_756.gif]

[Graphics:../Images/index_gr_757.gif]
[Graphics:../Images/index_gr_758.gif]

foo26 =[Graphics:../Images/index_gr_759.gif]

[Graphics:../Images/index_gr_760.gif]
[Graphics:../Images/index_gr_761.gif]

[Graphics:../Images/index_gr_762.gif]

[Graphics:../Images/index_gr_763.gif]
[Graphics:../Images/index_gr_764.gif]

foo28=foo27 foo24 =[Graphics:../Images/index_gr_765.gif]

[Graphics:../Images/index_gr_766.gif]
[Graphics:../Images/index_gr_767.gif]

foo29=foo23-foo28[Graphics:../Images/index_gr_768.gif]=foo23-[Graphics:../Images/index_gr_769.gif][Graphics:../Images/index_gr_770.gif].

[Graphics:../Images/index_gr_771.gif]
[Graphics:../Images/index_gr_772.gif]

Now, [Graphics:../Images/index_gr_773.gif] is transformed to [Graphics:../Images/index_gr_774.gif] by the simple conventions preserving the determinant.

Then, [Graphics:../Images/index_gr_775.gif]1+(foo26-1)) = [Graphics:../Images/index_gr_776.gif].

foo30=trace(foo29)

[Graphics:../Images/index_gr_777.gif]
[Graphics:../Images/index_gr_778.gif]

[Graphics:../Images/index_gr_779.gif]

[Graphics:../Images/index_gr_780.gif]
[Graphics:../Images/index_gr_781.gif]

[Graphics:../Images/index_gr_782.gif]

[Graphics:../Images/index_gr_783.gif]
[Graphics:../Images/index_gr_784.gif]

Finally,  logdetJ=[Graphics:../Images/index_gr_785.gif]foo31 [Graphics:../Images/index_gr_786.gif]

[Graphics:../Images/index_gr_787.gif]
[Graphics:../Images/index_gr_788.gif]
[Graphics:../Images/index_gr_789.gif]
[Graphics:../Images/index_gr_790.gif]


Converted by Mathematica      July 21, 2003