Module jw_mf

Functions

calc_jw(data, params)
Spectral density function.

calc_S2_jw(data, params)
Spectral density function.

calc_S2_te_jw(data, params)
Spectral density function.

calc_S2f_S2_ts_jw(data, params)
Spectral density function.

calc_S2f_tf_S2_ts_jw(data, params)
Spectral density function.

calc_S2f_S2s_ts_jw(data, params)
Spectral density function.

calc_S2f_tf_S2s_ts_jw(data, params)
Spectral density function.

calc_diff_djw_dGj(data, params, j)
Spectral density gradient.

calc_ellipsoid_djw_dGj(data, params, j)
Spectral density gradient.

calc_diff_S2_djw_dGj(data, params, j)
Spectral density gradient.

calc_ellipsoid_S2_djw_dGj(data, params, j)
Spectral density gradient.

calc_diff_S2_te_djw_dGj(data, params, j)
Spectral density gradient.

calc_ellipsoid_S2_te_djw_dGj(data, params, j)
Spectral density gradient.

calc_diff_djw_dOj(data, params, j)
Spectral density gradient.

calc_diff_S2_djw_dOj(data, params, j)
Spectral density gradient.

calc_diff_S2_te_djw_dOj(data, params, j)
Spectral density gradient.

calc_S2_djw_dS2(data, params, j)
Spectral density gradient.

calc_S2_te_djw_dS2(data, params, j)
Spectral density gradient.

calc_S2_te_djw_dte(data, params, j)
Spectral density gradient.

calc_diff_S2f_S2_ts_djw_dGj(data, params, j)
Spectral density gradient.

calc_ellipsoid_S2f_S2_ts_djw_dGj(data, params, j)
Spectral density gradient.

calc_diff_S2f_tf_S2_ts_djw_dGj(data, params, j)
Spectral density gradient.

calc_ellipsoid_S2f_tf_S2_ts_djw_dGj(data, params, j)
Spectral density gradient.

calc_diff_S2f_S2_ts_djw_dOj(data, params, j)
Spectral density gradient.

calc_diff_S2f_tf_S2_ts_djw_dOj(data, params, j)
Spectral density gradient.

calc_S2f_S2_ts_djw_dS2(data, params, j)
Spectral density gradient.

calc_S2f_S2_ts_djw_dS2f(data, params, j)
Spectral density gradient.

calc_S2f_tf_S2_ts_djw_dS2f(data, params, j)
Spectral density gradient.

calc_S2f_tf_S2_ts_djw_dtf(data, params, j)
Spectral density gradient.

calc_S2f_S2_ts_djw_dts(data, params, j)
Spectral density gradient.

calc_diff_S2f_S2s_ts_djw_dGj(data, params, j)
Spectral density gradient.

calc_ellipsoid_S2f_S2s_ts_djw_dGj(data, params, j)
Spectral density gradient.

calc_diff_S2f_tf_S2s_ts_djw_dGj(data, params, j)
Spectral density gradient.

calc_ellipsoid_S2f_tf_S2s_ts_djw_dGj(data, params, j)
Spectral density gradient.

calc_diff_S2f_S2s_ts_djw_dOj(data, params, j)
Spectral density gradient.

calc_diff_S2f_tf_S2s_ts_djw_dOj(data, params, j)
Spectral density gradient.

calc_S2f_S2s_ts_djw_dS2f(data, params, j)
Spectral density gradient.

calc_S2f_tf_S2s_ts_djw_dS2f(data, params, j)
Spectral density gradient.

calc_S2f_tf_S2s_ts_djw_dS2s(data, params, j)
Spectral density gradient.

calc_S2f_tf_S2s_ts_djw_dtf(data, params, j)
Spectral density gradient.

calc_S2f_S2s_ts_djw_dts(data, params, j)
Spectral density gradient.

calc_diff_d2jw_dGjdGk(data, params, j, k)
Spectral density Hessian.

calc_ellipsoid_d2jw_dGjdGk(data, params, j, k)
Spectral density Hessian.

calc_diff_S2_d2jw_dGjdGk(data, params, j, k)
Spectral density Hessian.

calc_ellipsoid_S2_d2jw_dGjdGk(data, params, j, k)
Spectral density Hessian.

calc_diff_S2_te_d2jw_dGjdGk(data, params, j, k)
Spectral density Hessian.

calc_ellipsoid_S2_te_d2jw_dGjdGk(data, params, j, k)
Spectral density Hessian.

calc_diff_d2jw_dGjdOj(data, params, j, k)
Spectral density Hessian.

calc_ellipsoid_d2jw_dGjdOj(data, params, j, k)
Spectral density Hessian.

calc_diff_S2_d2jw_dGjdOj(data, params, j, k)
Spectral density Hessian.

calc_ellipsoid_S2_d2jw_dGjdOj(data, params, j, k)
Spectral density Hessian.

calc_diff_S2_te_d2jw_dGjdOj(data, params, j, k)
Spectral density Hessian.

calc_ellipsoid_S2_te_d2jw_dGjdOj(data, params, j, k)
Spectral density Hessian.

calc_diff_S2_d2jw_dGjdS2(data, params, j, k)
Spectral density Hessian.

calc_ellipsoid_S2_d2jw_dGjdS2(data, params, j, k)
Spectral density Hessian.

calc_diff_S2_te_d2jw_dGjdS2(data, params, j, k)
Spectral density Hessian.

calc_ellipsoid_S2_te_d2jw_dGjdS2(data, params, j, k)
Spectral density Hessian.

calc_diff_S2_te_d2jw_dGjdte(data, params, j, k)
Spectral density Hessian.

calc_ellipsoid_S2_te_d2jw_dGjdte(data, params, j, k)
Spectral density Hessian.

calc_diff_d2jw_dOjdOk(data, params, j, k)
Spectral density Hessian.

calc_diff_S2_d2jw_dOjdOk(data, params, j, k)
Spectral density Hessian.

calc_diff_S2_te_d2jw_dOjdOk(data, params, j, k)
Spectral density Hessian.

calc_diff_S2_d2jw_dOjdS2(data, params, j, k)
Spectral density Hessian.

calc_diff_S2_te_d2jw_dOjdS2(data, params, j, k)
Spectral density Hessian.

calc_diff_S2_te_d2jw_dOjdte(data, params, j, k)
Spectral density Hessian.

calc_S2_te_d2jw_dS2dte(data, params, j, k)
Spectral density Hessian.

calc_S2_te_d2jw_dte2(data, params, j, k)
Spectral density Hessian.

calc_diff_S2f_S2_ts_d2jw_dGjdGk(data, params, j, k)
Spectral density Hessian.

calc_ellipsoid_S2f_S2_ts_d2jw_dGjdGk(data, params, j, k)
Spectral density Hessian.

calc_diff_S2f_tf_S2_ts_d2jw_dGjdGk(data, params, j, k)
Spectral density Hessian.

calc_ellipsoid_S2f_tf_S2_ts_d2jw_dGjdGk(data, params, j, k)
Spectral density Hessian.

calc_diff_S2f_S2_ts_d2jw_dGjdOj(data, params, j, k)
Spectral density Hessian.

calc_ellipsoid_S2f_S2_ts_d2jw_dGjdOj(data, params, j, k)
Spectral density Hessian.

calc_diff_S2f_tf_S2_ts_d2jw_dGjdOj(data, params, j, k)
Spectral density Hessian.

calc_ellipsoid_S2f_tf_S2_ts_d2jw_dGjdOj(data, params, j, k)
Spectral density Hessian.

calc_diff_S2f_S2_ts_d2jw_dGjdS2(data, params, j, k)
Spectral density Hessian.

calc_ellipsoid_S2f_S2_ts_d2jw_dGjdS2(data, params, j, k)
Spectral density Hessian.

calc_diff_S2f_S2_ts_d2jw_dGjdS2f(data, params, j, k)
Spectral density Hessian.

calc_ellipsoid_S2f_S2_ts_d2jw_dGjdS2f(data, params, j, k)
Spectral density Hessian.

calc_diff_S2f_tf_S2_ts_d2jw_dGjdS2f(data, params, j, k)
Spectral density Hessian.

calc_ellipsoid_S2f_tf_S2_ts_d2jw_dGjdS2f(data, params, j, k)
Spectral density Hessian.

calc_diff_S2f_tf_S2_ts_d2jw_dGjdtf(data, params, j, k)
Spectral density Hessian.

calc_ellipsoid_S2f_tf_S2_ts_d2jw_dGjdtf(data, params, j, k)
Spectral density Hessian.

calc_diff_S2f_S2_ts_d2jw_dGjdts(data, params, j, k)
Spectral density Hessian.

calc_ellipsoid_S2f_S2_ts_d2jw_dGjdts(data, params, j, k)
Spectral density Hessian.

calc_diff_S2f_S2_ts_d2jw_dOjdOk(data, params, j, k)
Spectral density Hessian.

calc_diff_S2f_tf_S2_ts_d2jw_dOjdOk(data, params, j, k)
Spectral density Hessian.

calc_diff_S2f_S2_ts_d2jw_dOjdS2(data, params, j, k)
Spectral density Hessian.

calc_diff_S2f_S2_ts_d2jw_dOjdS2f(data, params, j, k)
Spectral density Hessian.

calc_diff_S2f_tf_S2_ts_d2jw_dOjdS2f(data, params, j, k)
Spectral density Hessian.

calc_diff_S2f_tf_S2_ts_d2jw_dOjdtf(data, params, j, k)
Spectral density Hessian.

calc_diff_S2f_S2_ts_d2jw_dOjdts(data, params, j, k)
Spectral density Hessian.

calc_S2f_S2_ts_d2jw_dS2dts(data, params, j, k)
Spectral density Hessian.

calc_S2f_tf_S2_ts_d2jw_dS2fdtf(data, params, j, k)
Spectral density Hessian.

calc_S2f_S2_ts_d2jw_dS2fdts(data, params, j, k)
Spectral density Hessian.

calc_S2f_tf_S2_ts_d2jw_dtf2(data, params, j, k)
Spectral density Hessian.

calc_S2f_S2_ts_d2jw_dts2(data, params, j, k)
Spectral density Hessian.

calc_diff_S2f_S2s_ts_d2jw_dGjdGk(data, params, j, k)
Spectral density Hessian.

calc_ellipsoid_S2f_S2s_ts_d2jw_dGjdGk(data, params, j, k)
Spectral density Hessian.

calc_diff_S2f_tf_S2s_ts_d2jw_dGjdGk(data, params, j, k)
Spectral density Hessian.

calc_ellipsoid_S2f_tf_S2s_ts_d2jw_dGjdGk(data, params, j, k)
Spectral density Hessian.

calc_diff_S2f_S2s_ts_d2jw_dGjdOj(data, params, j, k)
Spectral density Hessian.

calc_ellipsoid_S2f_S2s_ts_d2jw_dGjdOj(data, params, j, k)
Spectral density Hessian.

calc_diff_S2f_tf_S2s_ts_d2jw_dGjdOj(data, params, j, k)
Spectral density Hessian.

calc_ellipsoid_S2f_tf_S2s_ts_d2jw_dGjdOj(data, params, j, k)
Spectral density Hessian.

calc_diff_S2f_S2s_ts_d2jw_dGjdS2f(data, params, j, k)
Spectral density Hessian.

calc_ellipsoid_S2f_S2s_ts_d2jw_dGjdS2f(data, params, j, k)
Spectral density Hessian.

calc_diff_S2f_tf_S2s_ts_d2jw_dGjdS2f(data, params, j, k)
Spectral density Hessian.

calc_ellipsoid_S2f_tf_S2s_ts_d2jw_dGjdS2f(data, params, j, k)
Spectral density Hessian.

calc_diff_S2f_S2s_ts_d2jw_dGjdS2s(data, params, j, k)
Spectral density Hessian.

calc_ellipsoid_S2f_S2s_ts_d2jw_dGjdS2s(data, params, j, k)
Spectral density Hessian.

calc_diff_S2f_tf_S2s_ts_d2jw_dGjdtf(data, params, j, k)
Spectral density Hessian.

calc_ellipsoid_S2f_tf_S2s_ts_d2jw_dGjdtf(data, params, j, k)
Spectral density Hessian.

calc_diff_S2f_S2s_ts_d2jw_dGjdts(data, params, j, k)
Spectral density Hessian.

calc_ellipsoid_S2f_S2s_ts_d2jw_dGjdts(data, params, j, k)
Spectral density Hessian.

calc_diff_S2f_S2s_ts_d2jw_dOjdOk(data, params, j, k)
Spectral density Hessian.

calc_diff_S2f_tf_S2s_ts_d2jw_dOjdOk(data, params, j, k)
Spectral density Hessian.

calc_diff_S2f_S2s_ts_d2jw_dOjdS2f(data, params, j, k)
Spectral density Hessian.

calc_diff_S2f_tf_S2s_ts_d2jw_dOjdS2f(data, params, j, k)
Spectral density Hessian.

calc_diff_S2f_tf_S2s_ts_d2jw_dOjdS2s(data, params, j, k)
Spectral density Hessian.

calc_diff_S2f_tf_S2s_ts_d2jw_dOjdtf(data, params, j, k)
Spectral density Hessian.

calc_diff_S2f_S2s_ts_d2jw_dOjdts(data, params, j, k)
Spectral density Hessian.

calc_S2f_S2s_ts_d2jw_dS2fdS2s(data, params, j, k)
Spectral density Hessian.

calc_S2f_tf_S2s_ts_d2jw_dS2fdtf(data, params, j, k)
Spectral density Hessian.

calc_S2f_S2s_ts_d2jw_dS2fdts(data, params, j, k)
Spectral density Hessian.

calc_S2f_S2s_ts_d2jw_dS2sdts(data, params, j, k)
Spectral density Hessian.

calc_S2f_tf_S2s_ts_d2jw_dtf2(data, params, j, k)
Spectral density Hessian.

calc_S2f_S2s_ts_d2jw_dts2(data, params, j, k)
Spectral density Hessian.

Variables

__package__ = 'maths_fns'

Imports: sum

Function Details

calc_jw(data, params)

Spectral density function.

Calculate the spectral density values for the original model-free formula with no parameters {} with or without diffusion tensor parameters.

The formula is:

               _k_
            2  \                1
   J(w)  =  -   >  ci . ti ------------.
            5  /__         1 + (w.ti)^2
               i=-k

calc_S2_jw(data, params)

Spectral density function.

Calculate the spectral density values for the original model-free formula with the single parameter {S2} with or without diffusion tensor parameters.

The formula is:

                 _k_
            2    \                1
   J(w)  =  - S2  >  ci . ti ------------.
            5    /__         1 + (w.ti)^2
                 i=-k

calc_S2_te_jw(data, params)

Spectral density function.

Calculate the spectral density values for the original model-free formula with the parameters {S2, te} with or without diffusion tensor parameters.

The model-free formula is:

               _k_
            2  \           /      S2             (1 - S2)(te + ti)te    \ 
   J(w)  =  -   >  ci . ti | ------------  +  ------------------------- |.
            5  /__         \ 1 + (w.ti)^2     (te + ti)^2 + (w.te.ti)^2 /
               i=-k

calc_S2f_S2_ts_jw(data, params)

Spectral density function.

Calculate the spectral density values for the extended model-free formula with the parameters {S2f, S2, ts} with or without diffusion tensor parameters.

The model-free formula is:

               _k_
            2  \           /      S2            (S2f - S2)(ts + ti)ts   \ 
   J(w)  =  -   >  ci . ti | ------------  +  ------------------------- |.
            5  /__         \ 1 + (w.ti)^2     (ts + ti)^2 + (w.ts.ti)^2 /
               i=-k

calc_S2f_tf_S2_ts_jw(data, params)

Spectral density function.

Calculate the spectral density values for the extended model-free formula with the parameters {S2f, tf, S2, ts} with or without diffusion tensor parameters.

The model-free formula is:

               _k_
            2  \           /      S2            (1 - S2f)(tf + ti)tf          (S2f - S2)(ts + ti)ts   \ 
   J(w)  =  -   >  ci . ti | ------------  +  -------------------------  +  ------------------------- |.
            5  /__         \ 1 + (w.ti)^2     (tf + ti)^2 + (w.tf.ti)^2     (ts + ti)^2 + (w.ts.ti)^2 /
               i=-k

calc_S2f_S2s_ts_jw(data, params)

Spectral density function.

Calculate the spectral density values for the extended model-free formula with the parameters {S2f, S2s, ts} with or without diffusion tensor parameters.

The model-free formula is:

                  _k_
            2     \           /      S2s           (1 - S2s)(ts + ti)ts    \ 
   J(w)  =  - S2f  >  ci . ti | ------------  +  ------------------------- |.
            5     /__         \ 1 + (w.ti)^2     (ts + ti)^2 + (w.ts.ti)^2 /
                  i=-k

calc_S2f_tf_S2s_ts_jw(data, params)

Spectral density function.

Calculate the spectral density values for the extended model-free formula with the parameters {S2f, tf, S2s, ts} with or without diffusion tensor parameters.

The model-free formula is:

               _k_
            2  \           /   S2f . S2s        (1 - S2f)(tf + ti)tf         S2f(1 - S2s)(ts + ti)ts  \ 
   J(w)  =  -   >  ci . ti | ------------  +  -------------------------  +  ------------------------- |.
            5  /__         \ 1 + (w.ti)^2     (tf + ti)^2 + (w.tf.ti)^2     (ts + ti)^2 + (w.ts.ti)^2 /
               i=-k

calc_diff_djw_dGj(data, params, j)

Spectral density gradient.

Calculate the spectral desity values for the Gj partial derivative of the original model-free formula with no parameters {} together with diffusion tensor parameters.

The model-free gradient is:

                _k_
   dJ(w)     2  \        dti    1 - (w.ti)^2
   -----  =  -   >  ci . ---  ----------------.
    dGj      5  /__      dGj  (1 + (w.ti)^2)^2
                i=-k

calc_ellipsoid_djw_dGj(data, params, j)

Spectral density gradient.

Calculate the spectral desity values for the Gj partial derivative of the original model-free formula with no parameters {} together with diffusion tensor parameters.

The model-free gradient is:

                _k_
   dJ(w)     2  \   /      dti    1 - (w.ti)^2       dci           1       \ 
   -----  =  -   >  | ci . ---  ----------------  +  --- . ti ------------ |.
    dGj      5  /__ \      dGj  (1 + (w.ti)^2)^2     dGj      1 + (w.ti)^2 /
                i=-k

calc_diff_S2_djw_dGj(data, params, j)

Spectral density gradient.

Calculate the spectral desity values for the Gj partial derivative of the original model-free formula with the parameter {S2} together with diffusion tensor parameters.

The model-free gradient is:

                  _k_
   dJ(w)     2    \        dti    1 - (w.ti)^2
   -----  =  - S2  >  ci . ---  ----------------.
    dGj      5    /__      dGj  (1 + (w.ti)^2)^2
                  i=-k

calc_ellipsoid_S2_djw_dGj(data, params, j)

Spectral density gradient.

Calculate the spectral desity values for the Gj partial derivative of the original model-free formula with the parameter {S2} together with diffusion tensor parameters.

The model-free gradient is:

                  _k_
   dJ(w)     2    \   /      dti    1 - (w.ti)^2       dci           1       \ 
   -----  =  - S2  >  | ci . ---  ----------------  +  --- . ti ------------ |.
    dGj      5    /__ \      dGj  (1 + (w.ti)^2)^2     dGj      1 + (w.ti)^2 /
                  i=-k

calc_diff_S2_te_djw_dGj(data, params, j)

Spectral density gradient.

Calculate the spectral desity values for the Gj partial derivative of the original model-free formula with the parameters {S2, te} together with diffusion tensor parameters.

The model-free gradient is:

                _k_
   dJ(w)     2  \        dti  /      1 - (w.ti)^2                      (te + ti)^2 - (w.te.ti)^2   \ 
   -----  =  -   >  ci . ---  | S2 ----------------  +  (1 - S2)te^2 ----------------------------- |.
    dGj      5  /__      dGj  \    (1 + (w.ti)^2)^2                  ((te + ti)^2 + (w.te.ti)^2)^2 /
                i=-k

calc_ellipsoid_S2_te_djw_dGj(data, params, j)

Spectral density gradient.

Calculate the spectral desity values for the Gj partial derivative of the original model-free formula with the parameters {S2, te} together with diffusion tensor parameters.

The model-free gradient is:

                _k_
   dJ(w)     2  \   /      dti  /      1 - (w.ti)^2                      (te + ti)^2 - (w.te.ti)^2   \ 
   -----  =  -   >  | ci . ---  | S2 ----------------  +  (1 - S2)te^2 ----------------------------- |
    dGj      5  /__ \      dGj  \    (1 + (w.ti)^2)^2                  ((te + ti)^2 + (w.te.ti)^2)^2 /
                i=-k

                         dci      /      S2             (1 - S2)(te + ti)te    \ \ 
                      +  --- . ti | ------------  +  ------------------------- | |.
                         dGj      \ 1 + (w.ti)^2     (te + ti)^2 + (w.te.ti)^2 / /

calc_diff_djw_dOj(data, params, j)

Spectral density gradient.

Calculate the spectral desity values for the O partial derivative of the original model-free formula with no parameters {} together with diffusion tensor parameters.

The model-free gradient is:

               _k_
   dJ(w)     2 \   dci           1
   -----  =  -  >  --- . ti ------------.
    dOj      5 /__ dOj      1 + (w.ti)^2
               i=-k

calc_diff_S2_djw_dOj(data, params, j)

Spectral density gradient.

Calculate the spectral desity values for the O partial derivative of the original model-free formula with the parameter {S2} together with diffusion tensor parameters.

The model-free gradient is:

                  _k_
   dJ(w)     2    \   dci           1
   -----  =  - S2  >  --- . ti ------------.
    dOj      5    /__ dOj      1 + (w.ti)^2
                  i=-k

calc_diff_S2_te_djw_dOj(data, params, j)

Spectral density gradient.

Calculate the spectral desity values for the O partial derivative of the original model-free formula with the parameters {S2, te} together with diffusion tensor parameters.

The model-free gradient is:

                _k_
   dJ(w)     2  \   dci      /      S2             (1 - S2)(te + ti)te    \ 
   -----  =  -   >  --- . ti | ------------  +  ------------------------- |.
    dOj      5  /__ dOj      \ 1 + (w.ti)^2     (te + ti)^2 + (w.te.ti)^2 /
                i=-k

calc_S2_djw_dS2(data, params, j)

Spectral density gradient.

Calculate the spectral desity values for the S2 partial derivative of the original model-free formula with the single parameter {S2} with or without diffusion tensor parameters.

The model-free gradient is:

                _k_
   dJ(w)     2  \                1
   -----  =  -   >  ci . ti ------------.
    dS2      5  /__         1 + (w.ti)^2
                i=-k

calc_S2_te_djw_dS2(data, params, j)

Spectral density gradient.

Calculate the spectral desity values for the S2 partial derivative of the original model-free formula with the parameters {S2, te} with or without diffusion tensor parameters.

The model-free gradient is:

                _k_
   dJ(w)     2  \           /      1                 (te + ti)te         \ 
   -----  =  -   >  ci . ti | ------------  -  ------------------------- |.
    dS2      5  /__         \ 1 + (w.ti)^2     (te + ti)^2 + (w.te.ti)^2 /
                i=-k

calc_S2_te_djw_dte(data, params, j)

Spectral density gradient.

Calculate the spectral desity values for the te partial derivative of the original model-free formula with the parameters {S2, te} with or without diffusion tensor parameters.

The model-free gradient is:

                        _k_
   dJ(w)     2          \               (te + ti)^2 - (w.te.ti)^2
   -----  =  - (1 - S2)  >  ci . ti^2 -----------------------------.
    dte      5          /__           ((te + ti)^2 + (w.te.ti)^2)^2
                        i=-k

calc_diff_S2f_S2_ts_djw_dGj(data, params, j)

Spectral density gradient.

Calculate the spectral desity values for the Gj partial derivative of the extended model-free formula with the parameters {S2f, S2, ts} together with diffusion tensor parameters.

The formula is:

                _k_
   dJ(w)     2  \        dti  /      1 - (w.ti)^2                        (ts + ti)^2 - (w.ts.ti)^2   \ 
   -----  =  -   >  ci . ---  | S2 ----------------  +  (S2f - S2)ts^2 ----------------------------- |.
    dGj      5  /__      dGj  \    (1 + (w.ti)^2)^2                    ((ts + ti)^2 + (w.ts.ti)^2)^2 /
                i=-k

calc_ellipsoid_S2f_S2_ts_djw_dGj(data, params, j)

Spectral density gradient.

Calculate the spectral desity values for the Gj partial derivative of the extended model-free formula with the parameters {S2f, S2, ts} together with diffusion tensor parameters.

The formula is:

                _k_
   dJ(w)     2  \   /      dti  /      1 - (w.ti)^2                        (ts + ti)^2 - (w.ts.ti)^2   \ 
   -----  =  -   >  | ci . ---  | S2 ----------------  +  (S2f - S2)ts^2 ----------------------------- |
    dGj      5  /__ \      dGj  \    (1 + (w.ti)^2)^2                    ((ts + ti)^2 + (w.ts.ti)^2)^2 /
                i=-k

                         dci      /      S2            (S2f - S2)(ts + ti)ts   \ \ 
                      +  --- . ti | ------------  +  ------------------------- | |.
                         dGj      \ 1 + (w.ti)^2     (ts + ti)^2 + (w.ts.ti)^2 / /

calc_diff_S2f_tf_S2_ts_djw_dGj(data, params, j)

Spectral density gradient.

Calculate the spectral desity values for the Gj partial derivative of the extended model-free formula with the parameters {S2f, tf, S2, ts} together with diffusion tensor parameters.

The formula is:

                _k_
   dJ(w)     2  \        dti  /      1 - (w.ti)^2                       (tf + ti)^2 - (w.tf.ti)^2                        (ts + ti)^2 - (w.ts.ti)^2   \ 
   -----  =  -   >  ci . ---  | S2 ----------------  +  (1 - S2f)tf^2 -----------------------------  +  (S2f - S2)ts^2 ----------------------------- |.
    dGj      5  /__      dGj  \    (1 + (w.ti)^2)^2                   ((tf + ti)^2 + (w.tf.ti)^2)^2                    ((ts + ti)^2 + (w.ts.ti)^2)^2 /
                i=-k

calc_ellipsoid_S2f_tf_S2_ts_djw_dGj(data, params, j)

Spectral density gradient.

Calculate the spectral desity values for the Gj partial derivative of the extended model-free formula with the parameters {S2f, tf, S2, ts} together with diffusion tensor parameters.

The formula is:

                _k_
   dJ(w)     2  \   /      dti  /      1 - (w.ti)^2                       (tf + ti)^2 - (w.tf.ti)^2                        (ts + ti)^2 - (w.ts.ti)^2   \ 
   -----  =  -   >  | ci . ---  | S2 ----------------  +  (1 - S2f)tf^2 -----------------------------  +  (S2f - S2)ts^2 ----------------------------- |
    dGj      5  /__ \      dGj  \    (1 + (w.ti)^2)^2                   ((tf + ti)^2 + (w.tf.ti)^2)^2                    ((ts + ti)^2 + (w.ts.ti)^2)^2 /
                i=-k

                         dci      /      S2            (1 - S2f)(tf + ti)tf          (S2f - S2)(ts + ti)ts   \ \ 
                      +  --- . ti | ------------  +  -------------------------  +  ------------------------- | |.
                         dGj      \ 1 + (w.ti)^2     (tf + ti)^2 + (w.tf.ti)^2     (ts + ti)^2 + (w.ts.ti)^2 / /

calc_diff_S2f_S2_ts_djw_dOj(data, params, j)

Spectral density gradient.

Calculate the spectral desity values for the O partial derivative of the extended model-free formula with the parameters {S2f, S2, ts} together with diffusion tensor parameters.

The formula is:

                _k_
   dJ(w)     2  \   dci      /      S2            (S2f - S2)(ts + ti)ts   \ 
   -----  =  -   >  --- . ti | ------------  +  ------------------------- |.
    dOj      5  /__ dOj      \ 1 + (w.ti)^2     (ts + ti)^2 + (w.ts.ti)^2 /
                i=-k

calc_diff_S2f_tf_S2_ts_djw_dOj(data, params, j)

Spectral density gradient.

Calculate the spectral desity values for the O partial derivative of the extended model-free formula with the parameters {S2f, tf, S2, ts} together with diffusion tensor parameters.

The formula is:

                _k_
   dJ(w)     2  \   dci      /      S2            (1 - S2f)(tf + ti)tf          (S2f - S2)(ts + ti)ts   \ 
   -----  =  -   >  --- . ti | ------------  +  -------------------------  +  ------------------------- |.
    dOj      5  /__ dOj      \ 1 + (w.ti)^2     (tf + ti)^2 + (w.tf.ti)^2     (ts + ti)^2 + (w.ts.ti)^2 /
                i=-k

calc_S2f_S2_ts_djw_dS2(data, params, j)

Spectral density gradient.

Calculate the spectral desity values for the S2 partial derivative of the extended model-free formula with the parameters {S2f, S2, ts} or {S2f, tf, S2, ts} with or without diffusion tensor parameters.

The formula is:

                _k_
   dJ(w)     2  \           /      1                 (ts + ti).ts        \ 
   -----  =  -   >  ci . ti | ------------  -  ------------------------- |.
    dS2      5  /__         \ 1 + (w.ti)^2     (ts + ti)^2 + (w.ts.ti)^2 /
                i=-k

calc_S2f_S2_ts_djw_dS2f(data, params, j)

Spectral density gradient.

Calculate the spectral desity values for the S2f partial derivative of the extended model-free formula with the parameters {S2f, S2, ts} with or without diffusion tensor parameters.

The formula is:

                _k_
   dJ(w)     2  \                 (ts + ti).ts
   -----  =  -   >  ci . ti -------------------------.
   dS2f      5  /__         (ts + ti)^2 + (w.ts.ti)^2
                i=-k

calc_S2f_tf_S2_ts_djw_dS2f(data, params, j)

Spectral density gradient.

Calculate the spectral desity values for the S2f partial derivative of the extended model-free formula with the parameters {S2f, tf, S2, ts} with or without diffusion tensor parameters.

The formula is:

                  _k_
   dJ(w)       2  \           /       (tf + ti).tf                  (ts + ti).ts        \ 
   -----  =  - -   >  ci . ti | -------------------------  -  ------------------------- |.
   dS2f        5  /__         \ (tf + ti)^2 + (w.tf.ti)^2     (ts + ti)^2 + (w.ts.ti)^2 /
                  i=-k

calc_S2f_tf_S2_ts_djw_dtf(data, params, j)

Spectral density gradient.

Calculate the spectral desity values for the tf partial derivative of the extended model-free formula with the parameters {S2f, tf, S2, ts} with or without diffusion tensor parameters.

The formula is:

                         _k_
   dJ(w)     2           \               (tf + ti)^2 - (w.tf.ti)^2
   -----  =  - (1 - S2f)  >  ci . ti^2 -----------------------------.
    dtf      5           /__           ((tf + ti)^2 + (w.tf.ti)^2)^2
                         i=-k

calc_S2f_S2_ts_djw_dts(data, params, j)

Spectral density gradient.

Calculate the spectral desity values for the ts partial derivative of the extended model-free formula with the parameters {S2f, S2, ts} or {S2f, tf, S2, ts} with or without diffusion tensor parameters.

The formula is:

                          _k_
   dJ(w)     2            \               (ts + ti)^2 - (w.ts.ti)^2
   -----  =  - (S2f - S2)  >  ci . ti^2 -----------------------------.
    dts      5            /__           ((ts + ti)^2 + (w.ts.ti)^2)^2
                          i=-k

calc_diff_S2f_S2s_ts_djw_dGj(data, params, j)

Spectral density gradient.

Calculate the spectral desity values for the Gj partial derivative of the extended model-free formula with the parameters {S2f, S2s, ts} together with diffusion tensor parameters.

The formula is:

                   _k_
   dJ(w)     2     \        dti  /       1 - (w.ti)^2                       (ts + ti)^2 - (w.ts.ti)^2   \ 
   -----  =  - S2f  >  ci . ---  | S2s ----------------  +  (1 - S2s)ts^2 ----------------------------- |.
    dGj      5     /__      dGj  \     (1 + (w.ti)^2)^2                   ((ts + ti)^2 + (w.ts.ti)^2)^2 /
                   i=-k

calc_ellipsoid_S2f_S2s_ts_djw_dGj(data, params, j)

Spectral density gradient.

Calculate the spectral desity values for the Gj partial derivative of the extended model-free formula with the parameters {S2f, S2s, ts} together with diffusion tensor parameters.

The formula is:

                   _k_
   dJ(w)     2     \   /      dti  /       1 - (w.ti)^2                       (ts + ti)^2 - (w.ts.ti)^2   \ 
   -----  =  - S2f  >  | ci . ---  | S2s ----------------  +  (1 - S2s)ts^2 ----------------------------- |
    dGj      5     /__ \      dGj  \     (1 + (w.ti)^2)^2                   ((ts + ti)^2 + (w.ts.ti)^2)^2 /
                   i=-k

                            dci      /     S2s             (1 - S2s)(ts + ti)ts   \ \ 
                         +  --- . ti | ------------  +  ------------------------- | |.
                            dGj      \ 1 + (w.ti)^2     (ts + ti)^2 + (w.ts.ti)^2 / /

calc_diff_S2f_tf_S2s_ts_djw_dGj(data, params, j)

Spectral density gradient.

Calculate the spectral desity values for the Gj partial derivative of the extended model-free formula with the parameters {S2f, tf, S2s, ts} together with diffusion tensor parameters.

The formula is:

                _k_
   dJ(w)     2  \        dti  /           1 - (w.ti)^2                       (tf + ti)^2 - (w.tf.ti)^2                          (ts + ti)^2 - (w.ts.ti)^2   \ 
   -----  =  -   >  ci . ---  | S2f.S2s ----------------  +  (1 - S2f)tf^2 -----------------------------  +  S2f(1 - S2s)ts^2 ----------------------------- |
    dGj      5  /__      dGj  \         (1 + (w.ti)^2)^2                   ((tf + ti)^2 + (w.tf.ti)^2)^2                      ((ts + ti)^2 + (w.ts.ti)^2)^2 /
                i=-k

calc_ellipsoid_S2f_tf_S2s_ts_djw_dGj(data, params, j)

Spectral density gradient.

Calculate the spectral desity values for the Gj partial derivative of the extended model-free formula with the parameters {S2f, tf, S2s, ts} together with diffusion tensor parameters.

The formula is:

                _k_
   dJ(w)     2  \   /      dti  /           1 - (w.ti)^2                        (tf + ti)^2 - (w.tf.ti)^2                           (ts + ti)^2 - (w.ts.ti)^2   \ 
   -----  =  -   >  | ci . ---  | S2f.S2s ----------------  +  (1 - S2f).tf^2 -----------------------------  +  S2f(1 - S2s).ts^2 ----------------------------- |
    dGj      5  /__ \      dGj  \         (1 + (w.ti)^2)^2                    ((tf + ti)^2 + (w.tf.ti)^2)^2                       ((ts + ti)^2 + (w.ts.ti)^2)^2 /
                i=-k

                         dci      /  S2f . S2s         (1 - S2f)(tf + ti)tf         S2f(1 - S2s)(ts + ti)ts  \ \ 
                      +  --- . ti | ------------  +  -------------------------  +  ------------------------- | |.
                         dGj      \ 1 + (w.ti)^2     (tf + ti)^2 + (w.tf.ti)^2     (ts + ti)^2 + (w.ts.ti)^2 / /

calc_diff_S2f_S2s_ts_djw_dOj(data, params, j)

Spectral density gradient.

Calculate the spectral desity values for the O partial derivative of the extended model-free formula with the parameters {S2f, S2s, ts} together with diffusion tensor parameters.

The formula is:

                   _k_
   dJ(w)     2     \   dci      /     S2s             (1 - S2s)(ts + ti)ts   \ 
   -----  =  - S2f  >  --- . ti | ------------  +  ------------------------- |.
    dOj      5     /__ dOj      \ 1 + (w.ti)^2     (ts + ti)^2 + (w.ts.ti)^2 /
                   i=-k

calc_diff_S2f_tf_S2s_ts_djw_dOj(data, params, j)

Spectral density gradient.

Calculate the spectral desity values for the O partial derivative of the extended model-free formula with the parameters {S2f, tf, S2s, ts} together with diffusion tensor parameters.

The formula is:

                _k_
   dJ(w)     2  \   dci      /  S2f . S2s         (1 - S2f)(tf + ti)tf         S2f(1 - S2s)(ts + ti)ts  \ 
   -----  =  -   >  --- . ti | ------------  +  -------------------------  +  ------------------------- |.
    dOj      5  /__ dOj      \ 1 + (w.ti)^2     (tf + ti)^2 + (w.tf.ti)^2     (ts + ti)^2 + (w.ts.ti)^2 /
                i=-k

calc_S2f_S2s_ts_djw_dS2f(data, params, j)

Spectral density gradient.

Calculate the spectral desity values for the S2f partial derivative of the extended model-free formula with the parameters {S2f, S2s, ts} with or without diffusion tensor parameters.

The formula is:

                _k_
   dJ(w)     2  \           /     S2s            (1 - S2s)(ts + ti).ts   \ 
   -----  =  -   >  ci . ti | ------------  +  ------------------------- |.
   dS2f      5  /__         \ 1 + (w.ti)^2     (ts + ti)^2 + (w.ts.ti)^2 /
                i=-k

calc_S2f_tf_S2s_ts_djw_dS2f(data, params, j)

Spectral density gradient.

Calculate the spectral desity values for the S2f partial derivative of the extended model-free formula with the parameters {S2f, tf, S2s, ts} with or without diffusion tensor parameters.

The formula is:

                _k_
   dJ(w)     2  \           /     S2s                (tf + ti).tf              (1 - S2s)(ts + ti).ts   \ 
   -----  =  -   >  ci . ti | ------------  -  -------------------------  +  ------------------------- |.
   dS2f      5  /__         \ 1 + (w.ti)^2     (tf + ti)^2 + (w.tf.ti)^2     (ts + ti)^2 + (w.ts.ti)^2 /
                i=-k

calc_S2f_tf_S2s_ts_djw_dS2s(data, params, j)

Spectral density gradient.

Calculate the spectral desity values for the S2s partial derivative of the extended model-free formula with the parameters {S2f, S2s, ts} or {S2f, tf, S2s, ts} with or without diffusion tensor parameters.

The formula is:

                   _k_
   dJ(w)     2     \           /      1                 (ts + ti).ts        \ 
   -----  =  - S2f  >  ci . ti | ------------  -  ------------------------- |.
   dS2s      5     /__         \ 1 + (w.ti)^2     (ts + ti)^2 + (w.ts.ti)^2 /
                   i=-k

calc_S2f_tf_S2s_ts_djw_dtf(data, params, j)

Spectral density gradient.

Calculate the spectral desity values for the tf partial derivative of the extended model-free formula with the parameters {S2f, tf, S2s, ts} with or without diffusion tensor parameters.

The formula is:

                         _k_
   dJ(w)     2           \               (tf + ti)^2 - (w.tf.ti)^2
   -----  =  - (1 - S2f)  >  ci . ti^2 -----------------------------.
    dtf      5           /__           ((tf + ti)^2 + (w.tf.ti)^2)^2
                         i=-k

calc_S2f_S2s_ts_djw_dts(data, params, j)

Spectral density gradient.

Calculate the spectral desity values for the ts partial derivative of the extended model-free formula with the parameters {S2f, S2s, ts} or {S2f, tf, S2s, ts} with or without diffusion tensor parameters.

The formula is:

                            _k_
   dJ(w)     2              \               (ts + ti)^2 - (w.ts.ti)^2
   -----  =  - S2f(1 - S2s)  >  ci . ti^2 -----------------------------.
    dts      5              /__           ((ts + ti)^2 + (w.ts.ti)^2)^2
                            i=-k

calc_diff_d2jw_dGjdGk(data, params, j, k)

Spectral density Hessian.

Calculate the spectral desity values for the Gj - Gk double partial derivative of the original model-free formula with no parameters {} together with diffusion tensor parameters.

The model-free Hessian is:

                  _k_
    d2J(w)     2  \      /    dti   dti           3 - (w.ti)^2        d2ti     1 - (w.ti)^2   \ 
   -------  =  -   >  ci | -2 --- . ---  w^2.ti ----------------  +  ------- ---------------- |.
   dGj.dGk     5  /__    \    dGj   dGk         (1 + (w.ti)^2)^3     dGj.dGk (1 + (w.ti)^2)^2 /
                  i=-k

calc_ellipsoid_d2jw_dGjdGk(data, params, j, k)

Spectral density Hessian.

Calculate the spectral desity values for the Gj - Gk double partial derivative of the original model-free formula with no parameters {} together with diffusion tensor parameters.

The model-free Hessian is:

                  _k_
    d2J(w)     2  \   /      dti   dti           3 - (w.ti)^2       / dti   dci     dti   dci         d2ti   \   1 - (w.ti)^2        d2ci           1       \ 
   -------  =  -   >  | -2ci --- . ---  w^2.ti ----------------  +  | --- . ---  +  --- . ---  +  ci ------- | ----------------  +  ------- ti ------------ |.
   dGj.dGk     5  /__ \      dGj   dGk         (1 + (w.ti)^2)^3     \ dGj   dGk     dGk   dGj        dGj.dGk / (1 + (w.ti)^2)^2     dGj.dGk    1 + (w.ti)^2 /
                  i=-k

calc_diff_S2_d2jw_dGjdGk(data, params, j, k)

Spectral density Hessian.

Calculate the spectral desity values for the Gj - Gk double partial derivative of the original model-free formula with the parameter {S2} together with diffusion tensor parameters.

The model-free Hessian is:

                    _k_
    d2J(w)     2    \      /    dti   dti           3 - (w.ti)^2        d2ti     1 - (w.ti)^2   \ 
   -------  =  - S2  >  ci | -2 --- . ---  w^2.ti ----------------  +  ------- ---------------- |.
   dGj.dGk     5    /__    \    dGj   dGk         (1 + (w.ti)^2)^3     dGj.dGk (1 + (w.ti)^2)^2 /
                    i=-k

calc_ellipsoid_S2_d2jw_dGjdGk(data, params, j, k)

Spectral density Hessian.

Calculate the spectral desity values for the Gj - Gk double partial derivative of the original model-free formula with the parameter {S2} together with diffusion tensor parameters.

The model-free Hessian is:

                    _k_
    d2J(w)     2    \   /      dti   dti           3 - (w.ti)^2       / dti   dci     dti   dci         d2ti   \   1 - (w.ti)^2        d2ci           1       \ 
   -------  =  - S2  >  | -2ci --- . ---  w^2.ti ----------------  +  | --- . ---  +  --- . ---  +  ci ------- | ----------------  +  ------- ti ------------ |.
   dGj.dGk     5    /__ \      dGj   dGk         (1 + (w.ti)^2)^3     \ dGj   dGk     dGk   dGj        dGj.dGk / (1 + (w.ti)^2)^2     dGj.dGk    1 + (w.ti)^2 /
                    i=-k

calc_diff_S2_te_d2jw_dGjdGk(data, params, j, k)

Spectral density Hessian.

Calculate the spectral desity values for the Gj - Gk double partial derivative of the original model-free formula with the parameters {S2, te} together with diffusion tensor parameters.

The model-free Hessian is:

                  _k_
    d2J(w)     2  \      /    dti   dti  /             3 - (w.ti)^2                    (te + ti)^3 + 3.w^2.te^3.ti(te + ti) - (w.te)^4.ti^3 \ 
   -------  =  -   >  ci | -2 --- . ---  | S2.w^2.ti ----------------  +  (1 - S2)te^2 ---------------------------------------------------- |
   dGj.dGk     5  /__    \    dGj   dGk  \           (1 + (w.ti)^2)^3                            ((te + ti)^2 + (w.te.ti)^2)^3              /
                  i=-k

                            d2ti   /      1 - (w.ti)^2                      (te + ti)^2 - (w.te.ti)^2   \ \ 
                        +  ------- | S2 ----------------  +  (1 - S2)te^2 ----------------------------- | |.
                           dGj.dGk \    (1 + (w.ti)^2)^2                  ((te + ti)^2 + (w.te.ti)^2)^2 / /

calc_ellipsoid_S2_te_d2jw_dGjdGk(data, params, j, k)

Spectral density Hessian.

Calculate the spectral desity values for the Gj - Gk double partial derivative of the original model-free formula with the parameters {S2, te} together with diffusion tensor parameters.

The model-free Hessian is:

                  _k_
    d2J(w)     2  \   /      dti   dti  /             3 - (w.ti)^2                    (te + ti)^3 + 3.w^2.te^3.ti(te + ti) - (w.te)^4.ti^3 \ 
   -------  =  -   >  | -2ci --- . ---  | S2.w^2.ti ----------------  +  (1 - S2)te^2 ---------------------------------------------------- |
   dGj.dGk     5  /__ \      dGj   dGk  \           (1 + (w.ti)^2)^3                            ((te + ti)^2 + (w.te.ti)^2)^3              /
                  i=-k

                           / dti   dci     dti   dci         d2ti   \ /      1 - (w.ti)^2                      (te + ti)^2 - (w.te.ti)^2   \ 
                        +  | --- . ---  +  --- . ---  +  ci ------- | | S2 ----------------  +  (1 - S2)te^2 ----------------------------- |
                           \ dGj   dGk     dGk   dGj        dGj.dGk / \    (1 + (w.ti)^2)^2                  ((te + ti)^2 + (w.te.ti)^2)^2 /


                            d2ci      /      S2             (1 - S2)(te + ti)te    \ \ 
                        +  ------- ti | ------------  +  ------------------------- | |.
                           dGj.dGk    \ 1 + (w.ti)^2     (te + ti)^2 + (w.te.ti)^2 / /

calc_diff_d2jw_dGjdOj(data, params, j, k)

Spectral density Hessian.

Calculate the spectral desity values for the Gj - Oj double partial derivative of the original model-free formula with no parameters {} together with diffusion tensor parameters.

The model-free Hessian is:

                  _k_
    d2J(w)     2  \   dci   dti     1 - (w.ti)^2
   -------  =  -   >  --- . --- . ----------------.
   dGj.dOj     5  /__ dOj   dGj   (1 + (w.ti)^2)^2
                  i=-k

calc_ellipsoid_d2jw_dGjdOj(data, params, j, k)

Spectral density Hessian.

Calculate the spectral desity values for the Gj - Oj double partial derivative of the original model-free formula with no parameters {} together with diffusion tensor parameters.

The model-free Hessian is:

                  _k_
    d2J(w)     2  \   / dci   dti     1 - (w.ti)^2        d2ci           1       \ 
   -------  =  -   >  | --- . --- . ----------------  +  ------- ti ------------ |.
   dGj.dOj     5  /__ \ dOj   dGj   (1 + (w.ti)^2)^2     dGj.dOj    1 + (w.ti)^2 /
                  i=-k

calc_diff_S2_d2jw_dGjdOj(data, params, j, k)

Spectral density Hessian.

Calculate the spectral desity values for the Gj - Oj double partial derivative of the original model-free formula with the parameter {S2} together with diffusion tensor parameters.

The model-free Hessian is:

                    _k_
    d2J(w)     2    \   dci   dti     1 - (w.ti)^2
   -------  =  - S2  >  --- . --- . ----------------.
   dGj.dOj     5    /__ dOj   dGj   (1 + (w.ti)^2)^2
                    i=-k

calc_ellipsoid_S2_d2jw_dGjdOj(data, params, j, k)

Spectral density Hessian.

Calculate the spectral desity values for the Gj - Oj double partial derivative of the original model-free formula with the parameter {S2} together with diffusion tensor parameters.

The model-free Hessian is:

                    _k_
    d2J(w)     2    \   / dci   dti     1 - (w.ti)^2        d2ci           1       \ 
   -------  =  - S2  >  | --- . --- . ----------------  +  ------- ti ------------ |.
   dGj.dOj     5    /__ \ dOj   dGj   (1 + (w.ti)^2)^2     dGj.dOj    1 + (w.ti)^2 /
                    i=-k

calc_diff_S2_te_d2jw_dGjdOj(data, params, j, k)

Spectral density Hessian.

Calculate the spectral desity values for the Gj - Oj double partial derivative of the original model-free formula with the parameters {S2, te} together with diffusion tensor parameters.

The model-free Hessian is:

                  _k_
    d2J(w)     2  \   dci   dti  /      1 - (w.ti)^2                      (te + ti)^2 - (w.te.ti)^2   \ 
   -------  =  -   >  --- . ---  | S2 ----------------  +  (1 - S2)te^2 ----------------------------- |.
   dGj.dOj     5  /__ dOj   dGj  \    (1 + (w.ti)^2)^2                  ((te + ti)^2 + (w.te.ti)^2)^2 /
                  i=-k

calc_ellipsoid_S2_te_d2jw_dGjdOj(data, params, j, k)

Spectral density Hessian.

Calculate the spectral desity values for the Gj - Oj double partial derivative of the original model-free formula with the parameters {S2, te} together with diffusion tensor parameters.

The model-free Hessian is:

                  _k_
    d2J(w)     2  \   / dci   dti  /      1 - (w.ti)^2                      (te + ti)^2 - (w.te.ti)^2   \ 
   -------  =  -   >  | --- . ---  | S2 ----------------  +  (1 - S2)te^2 ----------------------------- |
   dGj.dOj     5  /__ \ dOj   dGj  \    (1 + (w.ti)^2)^2                  ((te + ti)^2 + (w.te.ti)^2)^2 /
                  i=-k

                              d2ci      /      S2             (1 - S2)(te + ti)te    \ \ 
                          +  ------- ti | ------------  +  ------------------------- | |.
                             dGj.dOj    \ 1 + (w.ti)^2     (te + ti)^2 + (w.te.ti)^2 / /

calc_diff_S2_d2jw_dGjdS2(data, params, j, k)

Spectral density Hessian.

Calculate the spectral desity values for the Gj - S2 double partial derivative of the original model-free formula with the parameter {S2} together with diffusion tensor parameters.

The model-free Hessian is:

                  _k_
    d2J(w)     2  \        dti    1 - (w.ti)^2
   -------  =  -   >  ci . ---  ----------------.
   dGj.dS2     5  /__      dGj  (1 + (w.ti)^2)^2
                  i=-k

calc_ellipsoid_S2_d2jw_dGjdS2(data, params, j, k)

Spectral density Hessian.

Calculate the spectral desity values for the Gj - S2 double partial derivative of the original model-free formula with the parameter {S2} together with diffusion tensor parameters.

The model-free Hessian is:

                  _k_
    d2J(w)     2  \   /      dti    1 - (w.ti)^2       dci           1       \ 
   -------  =  -   >  | ci . ---  ----------------  +  --- . ti ------------ |.
   dGj.dS2     5  /__ \      dGj  (1 + (w.ti)^2)^2     dGj      1 + (w.ti)^2 /
                  i=-k

calc_diff_S2_te_d2jw_dGjdS2(data, params, j, k)

Spectral density Hessian.

Calculate the spectral desity values for the Gj - S2 double partial derivative of the original model-free formula with the parameters {S2, te} together with diffusion tensor parameters.

The model-free Hessian is:

                  _k_
    d2J(w)     2  \        dti  /   1 - (w.ti)^2              (te + ti)^2 - (w.te.ti)^2   \ 
   -------  =  -   >  ci . ---  | ----------------  -  te^2 ----------------------------- |.
   dGj.dS2     5  /__      dGj  \ (1 + (w.ti)^2)^2          ((te + ti)^2 + (w.te.ti)^2)^2 /
                  i=-k

calc_ellipsoid_S2_te_d2jw_dGjdS2(data, params, j, k)

Spectral density Hessian.

Calculate the spectral desity values for the Gj - S2 double partial derivative of the original model-free formula with the parameters {S2, te} together with diffusion tensor parameters.

The model-free Hessian is:

                  _k_
    d2J(w)     2  \   /      dti  /   1 - (w.ti)^2              (te + ti)^2 - (w.te.ti)^2   \ 
   -------  =  -   >  | ci . ---  | ----------------  -  te^2 ----------------------------- |
   dGj.dS2     5  /__ \      dGj  \ (1 + (w.ti)^2)^2          ((te + ti)^2 + (w.te.ti)^2)^2 /
                  i=-k

                           dci      /      1                 (te + ti)te         \ \ 
                        +  --- . ti | ------------  -  ------------------------- | |.
                           dGj      \ 1 + (w.ti)^2     (te + ti)^2 + (w.te.ti)^2 / /

calc_diff_S2_te_d2jw_dGjdte(data, params, j, k)

Spectral density Hessian.

Calculate the spectral desity values for the Gj - te double partial derivative of the original model-free formula with the parameters {S2, te} together with diffusion tensor parameters.

The model-free Hessian is:

                               _k_
    d2J(w)     4               \        dti                   (te + ti)^2 - 3(w.te.ti)^2
   -------  =  - (1 - S2) . te  >  ci . --- . ti . (te + ti) -----------------------------.
   dGj.dte     5               /__      dGj                  ((te + ti)^2 + (w.te.ti)^2)^3
                               i=-k

calc_ellipsoid_S2_te_d2jw_dGjdte(data, params, j, k)

Spectral density Hessian.

Calculate the spectral desity values for the Gj - te double partial derivative of the original model-free formula with the parameters {S2, te} together with diffusion tensor parameters.

The model-free Hessian is:

                         _k_
    d2J(w)     2         \   /       dti                        (te + ti)^2 - 3(w.te.ti)^2       dci          (te + ti)^2 - (w.te.ti)^2   \ 
   -------  =  - (1 - S2) >  | 2ci . --- . te . ti . (te + ti) -----------------------------  +  --- . ti^2 ----------------------------- |.
   dGj.dte     5         /__ \       dGj                       ((te + ti)^2 + (w.te.ti)^2)^3     dGj        ((te + ti)^2 + (w.te.ti)^2)^2 /
                         i=-k

calc_diff_d2jw_dOjdOk(data, params, j, k)

Spectral density Hessian.

Calculate the spectral desity values for the Oj - Ok double partial derivative of the original model-free formula with no parameters {} together with diffusion tensor parameters.

The model-free Hessian is:

                  _k_
    d2J(w)     2  \    d2ci          ti
   -------  =  -   >  ------- . ------------.
   dOj.dOk     5  /__ dOj.dOk   1 + (w.ti)^2
                  i=-k

calc_diff_S2_d2jw_dOjdOk(data, params, j, k)

Spectral density Hessian.

Calculate the spectral desity values for the Oj - Ok double partial derivative of the original model-free formula with the parameter {S2} together with diffusion tensor parameters.

The model-free Hessian is:

                    _k_
    d2J(w)     2    \    d2ci          ti
   -------  =  - S2  >  ------- . ------------.
   dOj.dOk     5    /__ dOj.dOk   1 + (w.ti)^2
                    i=-k

calc_diff_S2_te_d2jw_dOjdOk(data, params, j, k)

Spectral density Hessian.

Calculate the spectral desity values for the Oj - Ok double partial derivative of the original model-free formula with the parameters {S2, te} together with diffusion tensor parameters.

The model-free Hessian is:

                  _k_
    d2J(w)     2  \    d2ci        /      S2             (1 - S2)(te + ti)te    \ 
   -------  =  -   >  ------- . ti | ------------  +  ------------------------- |.
   dOj.dOk     5  /__ dOj.dOk      \ 1 + (w.ti)^2     (te + ti)^2 + (w.te.ti)^2 /
                  i=-k

calc_diff_S2_d2jw_dOjdS2(data, params, j, k)

Spectral density Hessian.

Calculate the spectral desity values for the Oj - S2 double partial derivative of the original model-free formula with the parameter {S2} together with diffusion tensor parameters.

The model-free Hessian is:

                  _k_
    d2J(w)     2  \   dci            1
   -------  =  -   >  --- . ti ------------.
   dOj.dS2     5  /__ dOj      1 + (w.ti)^2
                  i=-k

calc_diff_S2_te_d2jw_dOjdS2(data, params, j, k)

Spectral density Hessian.

Calculate the spectral desity values for the Oj - S2 double partial derivative of the original model-free formula with the parameters {S2, te} together with diffusion tensor parameters.

The model-free Hessian is:

                  _k_
    d2J(w)     2  \   dci      /      1                 (te + ti)te         \ 
   -------  =  -   >  --- . ti | ------------  -  ------------------------- |.
   dOj.dS2     5  /__ dOj      \ 1 + (w.ti)^2     (te + ti)^2 + (w.te.ti)^2 /
                  i=-k

calc_diff_S2_te_d2jw_dOjdte(data, params, j, k)

Spectral density Hessian.

Calculate the spectral desity values for the Oj - te double partial derivative of the original model-free formula with the parameters {S2, te} together with diffusion tensor parameters.

The model-free Hessian is:

                          _k_
    d2J(w)     2          \   dci          (te + ti)^2 - (w.te.ti)^2
   -------  =  - (1 - S2)  >  --- . ti^2 -----------------------------.
   dOj.dte     5          /__ dOj        ((te + ti)^2 + (w.te.ti)^2)^2
                          i=-k

calc_S2_te_d2jw_dS2dte(data, params, j, k)

Spectral density Hessian.

Calculate the spectral desity values for the S2 - te double partial derivative of the original model-free formula with the parameters {S2, te} with or without diffusion tensor parameters.

The model-free Hessian is:

                   _k_
    d2J(w)       2 \               (te + ti)^2 - (w.te.ti)^2
   -------  =  - -  >  ci . ti^2 -----------------------------.
   dS2.dte       5 /__           ((te + ti)^2 + (w.te.ti)^2)^2
                   i=-k

calc_S2_te_d2jw_dte2(data, params, j, k)

Spectral density Hessian.

Calculate the spectral desity values for the te - te double partial derivative of the original model-free formula with the parameters {S2, te} with or without diffusion tensor parameters.

The model-free Hessian is:

                           _k_
   d2J(w)       4          \             (te + ti)^3 + 3.w^2.ti^3.te.(te + ti) - (w.ti)^4.te^3
   ------  =  - - (1 - S2)  >  ci . ti^2 -----------------------------------------------------.
   dte**2       5          /__                        ((te + ti)^2 + (w.te.ti)^2)^3
                           i=-k

calc_diff_S2f_S2_ts_d2jw_dGjdGk(data, params, j, k)

Spectral density Hessian.

Calculate the spectral desity values for the Gj - Gk double partial derivative of the extended model-free formula with the parameters {S2f, S2, ts} together with diffusion tensor parameters.

The model-free Hessian is:

                  _k_
    d2J(w)     2  \      /    dti   dti  /             3 - (w.ti)^2                      (ts + ti)^3 + 3.w^2.ts^3.ti(ts + ti) - (w.ts)^4.ti^3 \ 
   -------  =  -   >  ci | -2 --- . ---  | S2.w^2.ti ----------------  +  (S2f - S2)ts^2 ---------------------------------------------------- |
   dGj.dGk     5  /__    \    dGj   dGk  \           (1 + (w.ti)^2)^3                              ((ts + ti)^2 + (w.ts.ti)^2)^3              /
                  i=-k

                            d2ti   /      1 - (w.ti)^2                        (ts + ti)^2 - (w.ts.ti)^2   \ \ 
                        +  ------- | S2 ----------------  +  (S2f - S2)ts^2 ----------------------------- | |.
                           dGj.dGk \    (1 + (w.ti)^2)^2                    ((ts + ti)^2 + (w.ts.ti)^2)^2 / /

calc_ellipsoid_S2f_S2_ts_d2jw_dGjdGk(data, params, j, k)

Spectral density Hessian.

Calculate the spectral desity values for the Gj - Gk double partial derivative of the extended model-free formula with the parameters {S2f, S2, ts} together with diffusion tensor parameters.

The model-free Hessian is:

                  _k_
    d2J(w)     2  \   /      dti   dti  /             3 - (w.ti)^2                      (ts + ti)^3 + 3.w^2.ts^3.ti(ts + ti) - (w.ts)^4.ti^3 \ 
   -------  =  -   >  | -2ci --- . ---  | S2.w^2.ti ----------------  +  (S2f - S2)ts^2 ---------------------------------------------------- |
   dGj.dGk     5  /__ \      dGj   dGk  \           (1 + (w.ti)^2)^3                              ((ts + ti)^2 + (w.ts.ti)^2)^3              /
                  i=-k

                           / dti   dci     dti   dci         d2ti   \ /      1 - (w.ti)^2                        (ts + ti)^2 - (w.ts.ti)^2   \ 
                        +  | --- . ---  +  --- . ---  +  ci ------- | | S2 ----------------  +  (S2f - S2)ts^2 ----------------------------- |
                           \ dGj   dGk     dGk   dGj        dGj.dGk / \    (1 + (w.ti)^2)^2                    ((ts + ti)^2 + (w.ts.ti)^2)^2 /


                            d2ci        /      S2            (S2f - S2)(ts + ti)ts   \ \ 
                        +  ------- . ti | ------------  +  ------------------------- | |.
                           dGj.dGk      \ 1 + (w.ti)^2     (ts + ti)^2 + (w.ts.ti)^2 / /

calc_diff_S2f_tf_S2_ts_d2jw_dGjdGk(data, params, j, k)

Spectral density Hessian.

Calculate the spectral desity values for the Gj - Gk double partial derivative of the extended model-free formula with the parameters {S2f, tf, S2, ts} together with diffusion tensor parameters.

The model-free Hessian is:

                  _k_
    d2J(w)     2  \      /    dti   dti  /             3 - (w.ti)^2                     (tf + ti)^3 + 3.w^2.tf^3.ti(tf + ti) - (w.tf)^4.ti^3
   -------  =  -   >  ci | -2 --- . ---  | S2.w^2.ti ----------------  +  (1 - S2f)tf^2 ----------------------------------------------------
   dGj.dGk     5  /__    \    dGj   dGk  \           (1 + (w.ti)^2)^3                             ((tf + ti)^2 + (w.tf.ti)^2)^3
                  i=-k

                                                            (ts + ti)^3 + 3.w^2.ts^3.ti(ts + ti) - (w.ts)^4.ti^3 \ 
                                          +  (S2f - S2)ts^2 ---------------------------------------------------- |
                                                                      ((ts + ti)^2 + (w.ts.ti)^2)^3              /


                            d2ti   /      1 - (w.ti)^2                       (tf + ti)^2 - (w.tf.ti)^2                        (ts + ti)^2 - (w.ts.ti)^2   \ \ 
                        +  ------- | S2 ----------------  +  (1 - S2f)tf^2 -----------------------------  +  (S2f - S2)ts^2 ----------------------------- | |.
                           dGj.dGk \    (1 + (w.ti)^2)^2                   ((tf + ti)^2 + (w.tf.ti)^2)^2                    ((ts + ti)^2 + (w.ts.ti)^2)^2 / /

calc_ellipsoid_S2f_tf_S2_ts_d2jw_dGjdGk(data, params, j, k)

Spectral density Hessian.

Calculate the spectral desity values for the Gj - Gk double partial derivative of the extended model-free formula with the parameters {S2f, tf, S2, ts} together with diffusion tensor parameters.

The model-free Hessian is:

                  _k_
    d2J(w)     2  \   /      dti   dti  /             3 - (w.ti)^2                     (tf + ti)^3 + 3.w^2.tf^3.ti(tf + ti) - (w.tf)^4.ti^3
   -------  =  -   >  | -2ci --- . ---  | S2.w^2.ti ----------------  +  (1 - S2f)tf^2 ----------------------------------------------------
   dGj.dGk     5  /__ \      dGj   dGk  \           (1 + (w.ti)^2)^3                             ((tf + ti)^2 + (w.tf.ti)^2)^3
                  i=-k

                                                            (ts + ti)^3 + 3.w^2.ts^3.ti(ts + ti) - (w.ts)^4.ti^3 \ 
                                          +  (S2f - S2)ts^2 ---------------------------------------------------- |
                                                                      ((ts + ti)^2 + (w.ts.ti)^2)^3              /


                           / dti   dci     dti   dci         d2ti   \ /      1 - (w.ti)^2                       (tf + ti)^2 - (w.tf.ti)^2
                        +  | --- . ---  +  --- . ---  +  ci ------- | | S2 ----------------  +  (1 - S2f)tf^2 -----------------------------
                           \ dGj   dGk     dGk   dGj        dGj.dGk / \    (1 + (w.ti)^2)^2                   ((tf + ti)^2 + (w.tf.ti)^2)^2


                                                                                            (ts + ti)^2 - (w.ts.ti)^2   \ 
                                                                        +  (S2f - S2)ts^2 ----------------------------- |
                                                                                          ((ts + ti)^2 + (w.ts.ti)^2)^2 /


                            d2ci        /      S2            (1 - S2f)(tf + ti)tf          (S2f - S2)(ts + ti)ts   \ \ 
                        +  ------- . ti | ------------  +  -------------------------  +  ------------------------- | |.
                           dGj.dGk      \ 1 + (w.ti)^2     (tf + ti)^2 + (w.tf.ti)^2     (ts + ti)^2 + (w.ts.ti)^2 / /

calc_diff_S2f_S2_ts_d2jw_dGjdOj(data, params, j, k)

Spectral density Hessian.

Calculate the spectral desity values for the Gj - Oj double partial derivative of the extended model-free formula with the parameters {S2f, S2, ts} together with diffusion tensor parameters.

The model-free Hessian is:

                  _k_
    d2J(w)     2  \   dci   dti  /      1 - (w.ti)^2                        (ts + ti)^2 - (w.ts.ti)^2   \ 
   -------  =  -   >  --- . ---  | S2 ----------------  +  (S2f - S2)ts^2 ----------------------------- |.
   dGj.dOj     5  /__ dOj   dGj  \    (1 + (w.ti)^2)^2                    ((ts + ti)^2 + (w.ts.ti)^2)^2 /
                  i=-k

calc_ellipsoid_S2f_S2_ts_d2jw_dGjdOj(data, params, j, k)

Spectral density Hessian.

Calculate the spectral desity values for the Gj - Oj double partial derivative of the extended model-free formula with the parameters {S2f, S2, ts} together with diffusion tensor parameters.

The model-free Hessian is:

                  _k_
    d2J(w)     2  \   / dci   dti  /      1 - (w.ti)^2                        (ts + ti)^2 - (w.ts.ti)^2   \ 
   -------  =  -   >  | --- . ---  | S2 ----------------  +  (S2f - S2)ts^2 ----------------------------- |
   dGj.dOj     5  /__ \ dOj   dGj  \    (1 + (w.ti)^2)^2                    ((ts + ti)^2 + (w.ts.ti)^2)^2 /
                  i=-k

                              d2ci        /      S2            (S2f - S2)(ts + ti)ts   \ \ 
                          +  ------- . ti | ------------  +  ------------------------- | |.
                             dGj.dOj      \ 1 + (w.ti)^2     (ts + ti)^2 + (w.ts.ti)^2 / /

calc_diff_S2f_tf_S2_ts_d2jw_dGjdOj(data, params, j, k)

Spectral density Hessian.

Calculate the spectral desity values for the Gj - Oj double partial derivative of the extended model-free formula with the parameters {S2f, tf, S2, ts} together with diffusion tensor parameters.

The model-free Hessian is:

                  _k_
    d2J(w)     2  \   dci   dti  /      1 - (w.ti)^2                       (tf + ti)^2 - (w.tf.ti)^2                        (ts + ti)^2 - (w.ts.ti)^2   \ 
   -------  =  -   >  --- . ---  | S2 ----------------  +  (1 - S2f)tf^2 -----------------------------  +  (S2f - S2)ts^2 ----------------------------- |.
   dGj.dOj     5  /__ dOj   dGj  \    (1 + (w.ti)^2)^2                   ((tf + ti)^2 + (w.tf.ti)^2)^2                    ((ts + ti)^2 + (w.ts.ti)^2)^2 /
                  i=-k

calc_ellipsoid_S2f_tf_S2_ts_d2jw_dGjdOj(data, params, j, k)

Spectral density Hessian.

Calculate the spectral desity values for the Gj - Oj double partial derivative of the extended model-free formula with the parameters {S2f, tf, S2, ts} together with diffusion tensor parameters.

The model-free Hessian is:

                  _k_
    d2J(w)     2  \   / dci   dti  /      1 - (w.ti)^2                       (tf + ti)^2 - (w.tf.ti)^2                        (ts + ti)^2 - (w.ts.ti)^2   \ 
   -------  =  -   >  | --- . ---  | S2 ----------------  +  (1 - S2f)tf^2 -----------------------------  +  (S2f - S2)ts^2 ----------------------------- |
   dGj.dOj     5  /__ \ dOj   dGj  \    (1 + (w.ti)^2)^2                   ((tf + ti)^2 + (w.tf.ti)^2)^2                    ((ts + ti)^2 + (w.ts.ti)^2)^2 /
                  i=-k

                              d2ci        /      S2            (1 - S2f)(tf + ti)tf          (S2f - S2)(ts + ti)ts   \ \ 
                          +  ------- . ti | ------------  +  -------------------------  +  ------------------------- | |.
                             dGj.dOj      \ 1 + (w.ti)^2     (tf + ti)^2 + (w.tf.ti)^2     (ts + ti)^2 + (w.ts.ti)^2 / /

calc_diff_S2f_S2_ts_d2jw_dGjdS2(data, params, j, k)

Spectral density Hessian.

Calculate the spectral desity values for the Gj - S2 double partial derivative of the extended model-free formula with the parameters {S2f, S2, ts} or {S2f, tf, S2, ts} together with diffusion tensor parameters.

The model-free Hessian is:

                  _k_
    d2J(w)     2  \        dti  /   1 - (w.ti)^2              (ts + ti)^2 - (w.ts.ti)^2   \ 
   -------  =  -   >  ci . ---  | ----------------  -  ts^2 ----------------------------- |.
   dGj.dS2     5  /__      dGj  \ (1 + (w.ti)^2)^2          ((ts + ti)^2 + (w.ts.ti)^2)^2 /
                  i=-k

calc_ellipsoid_S2f_S2_ts_d2jw_dGjdS2(data, params, j, k)

Spectral density Hessian.

Calculate the spectral desity values for the Gj - S2 double partial derivative of the extended model-free formula with the parameters {S2f, S2, ts} or {S2f, tf, S2, ts} together with diffusion tensor parameters.

The model-free Hessian is:

                  _k_
    d2J(w)     2  \   /      dti  /   1 - (w.ti)^2              (ts + ti)^2 - (w.ts.ti)^2   \ 
   -------  =  -   >  | ci . ---  | ----------------  -  ts^2 ----------------------------- |
   dGj.dS2     5  /__ \      dGj  \ (1 + (w.ti)^2)^2          ((ts + ti)^2 + (w.ts.ti)^2)^2 /
                  i=-k

                           dci      /      1                 (ts + ti)ts         \ \ 
                        +  --- . ti | ------------  -  ------------------------- | |.
                           dGj      \ 1 + (w.ti)^2     (ts + ti)^2 + (w.ts.ti)^2 / /

calc_diff_S2f_S2_ts_d2jw_dGjdS2f(data, params, j, k)

Spectral density Hessian.

Calculate the spectral desity values for the Gj - S2f double partial derivative of the extended model-free formula with the parameters {S2f, S2, ts} together with diffusion tensor parameters.

The model-free Hessian is:

                   _k_
    d2J(w)      2  \        dti         (ts + ti)^2 - (w.ts.ti)^2
   --------  =  -   >  ci . ---  ts^2 -----------------------------.
   dGj.dS2f     5  /__      dGj       ((ts + ti)^2 + (w.ts.ti)^2)^2
                   i=-k

calc_ellipsoid_S2f_S2_ts_d2jw_dGjdS2f(data, params, j, k)

Spectral density Hessian.

Calculate the spectral desity values for the Gj - S2f double partial derivative of the extended model-free formula with the parameters {S2f, S2, ts} together with diffusion tensor parameters.

The model-free Hessian is:

                   _k_
    d2J(w)      2  \   /      dti         (ts + ti)^2 - (w.ts.ti)^2       dci            (ts + ti)ts         \ 
   --------  =  -   >  | ci . ---  ts^2 -----------------------------  +  --- . ti ------------------------- |.
   dGj.dS2f     5  /__ \      dGj       ((ts + ti)^2 + (w.ts.ti)^2)^2     dGj      (ts + ti)^2 + (w.ts.ti)^2 /
                   i=-k

calc_diff_S2f_tf_S2_ts_d2jw_dGjdS2f(data, params, j, k)

Spectral density Hessian.

Calculate the spectral desity values for the Gj - S2f double partial derivative of the extended model-free formula with the parameters {S2f, tf, S2, ts} together with diffusion tensor parameters.

The model-free Hessian is:

                     _k_
    d2J(w)        2  \        dti  /        (tf + ti)^2 - (w.tf.ti)^2              (ts + ti)^2 - (w.ts.ti)^2   \ 
   --------  =  - -   >  ci . ---  | tf^2 -----------------------------  -  ts^2 ----------------------------- |.
   dGj.dS2f       5  /__      dGj  \      ((tf + ti)^2 + (w.tf.ti)^2)^2          ((ts + ti)^2 + (w.ts.ti)^2)^2 /
                     i=-k

calc_ellipsoid_S2f_tf_S2_ts_d2jw_dGjdS2f(data, params, j, k)

Spectral density Hessian.

Calculate the spectral desity values for the Gj - S2f double partial derivative of the extended model-free formula with the parameters {S2f, tf, S2, ts} together with diffusion tensor parameters.

The model-free Hessian is:

                     _k_
    d2J(w)        2  \   /      dti  /        (tf + ti)^2 - (w.tf.ti)^2              (ts + ti)^2 - (w.ts.ti)^2   \ 
   --------  =  - -   >  | ci . ---  | tf^2 -----------------------------  -  ts^2 ----------------------------- |
   dGj.dS2f       5  /__ \      dGj  \      ((tf + ti)^2 + (w.tf.ti)^2)^2          ((ts + ti)^2 + (w.ts.ti)^2)^2 /
                     i=-k

                              dci      /       (tf + ti)tf                   (ts + ti)ts         \ \ 
                           +  --- . ti | -------------------------  -  ------------------------- | |.
                              dGj      \ (tf + ti)^2 + (w.tf.ti)^2     (ts + ti)^2 + (w.ts.ti)^2 / /

calc_diff_S2f_tf_S2_ts_d2jw_dGjdtf(data, params, j, k)

Spectral density Hessian.

Calculate the spectral desity values for the Gj - tf double partial derivative of the extended model-free formula with the parameters {S2f, tf, S2, ts} together with diffusion tensor parameters.

The model-free Hessian is:

                                _k_
    d2J(w)     4                \        dti                   (tf + ti)^2 - 3(w.tf.ti)^2
   -------  =  - (1 - S2f) . tf  >  ci . --- . ti . (tf + ti) -----------------------------.
   dGj.dtf     5                /__      dGj                  ((tf + ti)^2 + (w.tf.ti)^2)^3
                                i=-k

calc_ellipsoid_S2f_tf_S2_ts_d2jw_dGjdtf(data, params, j, k)

Spectral density Hessian.

Calculate the spectral desity values for the Gj - tf double partial derivative of the extended model-free formula with the parameters {S2f, tf, S2, ts} together with diffusion tensor parameters.

The model-free Hessian is:

                           _k_
    d2J(w)     2           \   /       dti                        (tf + ti)^2 - 3(w.tf.ti)^2       dci          (tf + ti)^2 - (w.tf.ti)^2   \ 
   -------  =  - (1 - S2f)  >  | 2ci . --- . tf . ti . (tf + ti) -----------------------------  +  --- . ti^2 ----------------------------- |.
   dGj.dtf     5           /__ \       dGj                       ((tf + ti)^2 + (w.tf.ti)^2)^3     dGj        ((tf + ti)^2 + (w.tf.ti)^2)^2 /
                           i=-k

calc_diff_S2f_S2_ts_d2jw_dGjdts(data, params, j, k)

Spectral density Hessian.

Calculate the spectral desity values for the Gj - ts double partial derivative of the extended model-free formula with the parameters {S2f, S2, ts} or {S2f, tf, S2, ts} together with diffusion tensor parameters.

The model-free Hessian is:

                                 _k_
    d2J(w)     4                 \        dti                   (ts + ti)^2 - 3(w.ts.ti)^2
   -------  =  - (S2f - S2) . ts  >  ci . --- . ti . (ts + ti) -----------------------------.
   dGj.dts     5                 /__      dGj                  ((ts + ti)^2 + (w.ts.ti)^2)^3
                                 i=-k

calc_ellipsoid_S2f_S2_ts_d2jw_dGjdts(data, params, j, k)

Spectral density Hessian.

Calculate the spectral desity values for the Gj - ts double partial derivative of the extended model-free formula with the parameters {S2f, S2, ts} or {S2f, tf, S2, ts} together with diffusion tensor parameters.

The model-free Hessian is:

                            _k_
    d2J(w)     2            \   /       dti                        (ts + ti)^2 - 3(w.ts.ti)^2       dci          (ts + ti)^2 - (w.ts.ti)^2   \ 
   -------  =  - (S2f - S2)  >  | 2ci . --- . ts . ti . (ts + ti) -----------------------------  +  --- . ti^2 ----------------------------- |.
   dGj.dts     5            /__ \       dGj                       ((ts + ti)^2 + (w.ts.ti)^2)^3     dGj        ((ts + ti)^2 + (w.ts.ti)^2)^2 /
                            i=-k

calc_diff_S2f_S2_ts_d2jw_dOjdOk(data, params, j, k)

Spectral density Hessian.

Calculate the spectral desity values for the Oj - Ok double partial derivative of the extended model-free formula with the parameters {S2f, S2, ts} together with diffusion tensor parameters.

The model-free Hessian is:

                  _k_
    d2J(w)     2  \    d2ci        /      S2            (S2f - S2)(ts + ti)ts   \ 
   -------  =  -   >  ------- . ti | ------------  +  ------------------------- |.
   dOj.dOk     5  /__ dOj.dOk      \ 1 + (w.ti)^2     (ts + ti)^2 + (w.ts.ti)^2 /
                  i=-k

calc_diff_S2f_tf_S2_ts_d2jw_dOjdOk(data, params, j, k)

Spectral density Hessian.

Calculate the spectral desity values for the Oj - Ok double partial derivative of the extended model-free formula with the parameters {S2f, tf, S2, ts} together with diffusion tensor parameters.

The model-free Hessian is:

                  _k_
    d2J(w)     2  \    d2ci        /      S2            (1 - S2f)(tf + ti)tf          (S2f - S2)(ts + ti)ts   \ 
   -------  =  -   >  ------- . ti | ------------  +  -------------------------  +  ------------------------- |.
   dOj.dOk     5  /__ dOj.dOk      \ 1 + (w.ti)^2     (tf + ti)^2 + (w.tf.ti)^2     (ts + ti)^2 + (w.ts.ti)^2 /
                  i=-k

calc_diff_S2f_S2_ts_d2jw_dOjdS2(data, params, j, k)

Spectral density Hessian.

Calculate the spectral desity values for the Oj - S2 double partial derivative of the extended model-free formula with the parameters {S2f, S2, ts} and {S2f, tf, S2, ts} together with diffusion tensor parameters.

The model-free Hessian is:

                  _k_
    d2J(w)     2  \   dci      /      1                 (ts + ti)ts         \ 
   -------  =  -   >  --- . ti | ------------  -  ------------------------- |.
   dOj.dS2     5  /__ dOj      \ 1 + (w.ti)^2     (ts + ti)^2 + (w.ts.ti)^2 /
                  i=-k

calc_diff_S2f_S2_ts_d2jw_dOjdS2f(data, params, j, k)

Spectral density Hessian.

Calculate the spectral desity values for the Oj - S2f double partial derivative of the extended model-free formula with the parameters {S2f, S2, ts} together with diffusion tensor parameters.

The model-free Hessian is:

                   _k_
    d2J(w)      2  \   dci             (ts + ti)ts
   --------  =  -   >  --- . ti -------------------------.
   dOj.dS2f     5  /__ dOj      (ts + ti)^2 + (w.ts.ti)^2
                   i=-k

calc_diff_S2f_tf_S2_ts_d2jw_dOjdS2f(data, params, j, k)

Spectral density Hessian.

Calculate the spectral desity values for the Oj - S2f double partial derivative of the extended model-free formula with the parameters {S2f, tf, S2, ts} together with diffusion tensor parameters.

The model-free Hessian is:

                     _k_
    d2J(w)        2  \   dci      /       (tf + ti)tf                   (ts + ti)ts         \ 
   --------  =  - -   >  --- . ti | -------------------------  -  ------------------------- |.
   dOj.dS2f       5  /__ dOj      \ (tf + ti)^2 + (w.tf.ti)^2     (ts + ti)^2 + (w.ts.ti)^2 /
                     i=-k

calc_diff_S2f_tf_S2_ts_d2jw_dOjdtf(data, params, j, k)

Spectral density Hessian.

Calculate the spectral desity values for the Oj - tf double partial derivative of the extended model-free formula with the parameters {S2f, tf, S2, ts} together with diffusion tensor parameters.

The model-free Hessian is:

                           _k_
    d2J(w)     2           \   dci          (tf + ti)^2 - (w.tf.ti)^2
   -------  =  - (1 - S2f)  >  --- . ti^2 -----------------------------.
   dOj.dtf     5           /__ dOj        ((tf + ti)^2 + (w.tf.ti)^2)^2
                           i=-k

calc_diff_S2f_S2_ts_d2jw_dOjdts(data, params, j, k)

Spectral density Hessian.

Calculate the spectral desity values for the Oj - ts double partial derivative of the extended model-free formula with the parameters {S2f, S2, ts} or {S2f, tf, S2, ts} together with diffusion tensor parameters.

The model-free Hessian is:

                            _k_
    d2J(w)     2            \   dci          (ts + ti)^2 - (w.ts.ti)^2
   -------  =  - (S2f - S2)  >  --- . ti^2 -----------------------------.
   dOj.dts     5            /__ dOj        ((ts + ti)^2 + (w.ts.ti)^2)^2
                            i=-k

calc_S2f_S2_ts_d2jw_dS2dts(data, params, j, k)

Spectral density Hessian.

Calculate the spectral desity values for the S2 - ts double partial derivative of the extended model-free formula with the parameters {S2f, S2, ts} or {S2f, tf, S2, ts} with or without diffusion tensor parameters.

The model-free Hessian is:

                    _k_
    d2J(w)       2  \               (ts + ti)^2 - (w.ts.ti)^2
   -------  =  - -   >  ci . ti^2 -----------------------------.
   dS2.dts       5  /__           ((ts + ti)^2 + (w.ts.ti)^2)^2
                    i=-k

calc_S2f_tf_S2_ts_d2jw_dS2fdtf(data, params, j, k)

Spectral density Hessian.

Calculate the spectral desity values for the S2f - tf double partial derivative of the extended model-free formula with the parameters {S2f, tf, S2, ts} with or without diffusion tensor parameters.

The model-free Hessian is:

                     _k_
    d2J(w)        2  \               (tf + ti)^2 - (w.tf.ti)^2
   --------  =  - -   >  ci . ti^2 -----------------------------.
   dS2f.dtf       5  /__           ((tf + ti)^2 + (w.tf.ti)^2)^2
                     i=-k

calc_S2f_S2_ts_d2jw_dS2fdts(data, params, j, k)

Spectral density Hessian.

Calculate the spectral desity values for the S2f - ts double partial derivative of the extended model-free formula with the parameters {S2f, S2, ts} or {S2f, tf, S2, ts} with or without diffusion tensor parameters.

The model-free Hessian is:

                   _k_
    d2J(w)      2  \               (ts + ti)^2 - (w.ts.ti)^2
   --------  =  -   >  ci . ti^2 -----------------------------.
   dS2f.dts     5  /__           ((ts + ti)^2 + (w.ts.ti)^2)^2
                   i=-k

calc_S2f_tf_S2_ts_d2jw_dtf2(data, params, j, k)

Spectral density Hessian.

Calculate the spectral desity values for the tf - tf double partial derivative of the extended model-free formula with the parameters {S2f, tf, S2, ts} with or without diffusion tensor parameters.

The model-free Hessian is:

                            _k_
   d2J(w)       4           \             (tf + ti)^3 + 3.w^2.ti^3.tf.(tf + ti) - (w.ti)^4.tf^3
   ------  =  - - (1 - S2f)  >  ci . ti^2 -----------------------------------------------------.
   dtf**2       5           /__                        ((tf + ti)^2 + (w.tf.ti)^2)^3
                            i=-k

calc_S2f_S2_ts_d2jw_dts2(data, params, j, k)

Spectral density Hessian.

Calculate the spectral desity values for the ts - ts double partial derivative of the extended model-free formula with the parameters {S2f, S2, ts} or {S2f, tf, S2, ts} with or without diffusion tensor parameters.

The model-free Hessian is:

                             _k_
   d2J(w)       4            \             (ts + ti)^3 + 3.w^2.ti^3.ts.(ts + ti) - (w.ti)^4.ts^3
   ------  =  - - (S2f - S2)  >  ci . ti^2 -----------------------------------------------------.
   dts**2       5            /__                        ((ts + ti)^2 + (w.ts.ti)^2)^3
                             i=-k

calc_diff_S2f_S2s_ts_d2jw_dGjdGk(data, params, j, k)

Spectral density Hessian.

Calculate the spectral desity values for the Gj - Gk double partial derivative of the extended model-free formula with the parameters {S2f, S2s, ts} together with diffusion tensor parameters.

The model-free Hessian is:

                  _k_
    d2J(w)     2  \      /    dti   dti  /                  3 - (w.ti)^2                        (ts + ti)^3 + 3.w^2.ts^3.ti(ts + ti) - (w.ts)^4.ti^3 \ 
   -------  =  -   >  ci | -2 --- . ---  | S2f.S2s.w^2.ti ----------------  +  S2f(1 - S2s)ts^2 ---------------------------------------------------- |
   dGj.dGk     5  /__    \    dGj   dGk  \                (1 + (w.ti)^2)^3                                ((ts + ti)^2 + (w.ts.ti)^2)^3              /
                  i=-k

                            d2ti   /           1 - (w.ti)^2                          (ts + ti)^2 - (w.ts.ti)^2   \ \ 
                        +  ------- | S2f.S2s ----------------  +  S2f(1 - S2s)ts^2 ----------------------------- | |.
                           dGj.dGk \         (1 + (w.ti)^2)^2                      ((ts + ti)^2 + (w.ts.ti)^2)^2 / /

calc_ellipsoid_S2f_S2s_ts_d2jw_dGjdGk(data, params, j, k)

Spectral density Hessian.

Calculate the spectral desity values for the Gj - Gk double partial derivative of the extended model-free formula with the parameters {S2f, S2s, ts} together with diffusion tensor parameters.

The model-free Hessian is:

                  _k_
    d2J(w)     2  \   /      dti   dti  /                  3 - (w.ti)^2                        (ts + ti)^3 + 3.w^2.ts^3.ti(ts + ti) - (w.ts)^4.ti^3 \ 
   -------  =  -   >  | -2ci --- . ---  | S2f.S2s.w^2.ti ----------------  +  S2f(1 - S2s)ts^2 ---------------------------------------------------- |
   dGj.dGk     5  /__ \      dGj   dGk  \                (1 + (w.ti)^2)^3                                ((ts + ti)^2 + (w.ts.ti)^2)^3              /
                  i=-k

                           / dti   dci     dti   dci         d2ti   \ /           1 - (w.ti)^2                          (ts + ti)^2 - (w.ts.ti)^2   \ 
                        +  | --- . ---  +  --- . ---  +  ci ------- | | S2f.S2s ----------------  +  S2f(1 - S2s)ts^2 ----------------------------- |
                           \ dGj   dGk     dGk   dGj        dGj.dGk / \         (1 + (w.ti)^2)^2                      ((ts + ti)^2 + (w.ts.ti)^2)^2 /


                            d2ci        /   S2f.S2s         S2f(1 - S2s)(ts + ti)ts  \ \ 
                        +  ------- . ti | ------------  +  ------------------------- | |.
                           dGj.dGk      \ 1 + (w.ti)^2     (ts + ti)^2 + (w.ts.ti)^2 / /

calc_diff_S2f_tf_S2s_ts_d2jw_dGjdGk(data, params, j, k)

Spectral density Hessian.

Calculate the spectral desity values for the Gj - Gk double partial derivative of the extended model-free formula with the parameters {S2f, tf, S2s, ts} together with diffusion tensor parameters.

The model-free Hessian is:

                  _k_
    d2J(w)     2  \      /    dti   dti  /                  3 - (w.ti)^2                     (tf + ti)^3 + 3.w^2.tf^3.ti(tf + ti) - (w.tf)^4.ti^3
   -------  =  -   >  ci | -2 --- . ---  | S2f.S2s.w^2.ti ----------------  +  (1 - S2f)tf^2 ----------------------------------------------------
   dGj.dGk     5  /__    \    dGj   dGk  \                (1 + (w.ti)^2)^3                             ((tf + ti)^2 + (w.tf.ti)^2)^3
                  i=-k

                                                              (ts + ti)^3 + 3.w^2.ts^3.ti(ts + ti) - (w.ts)^4.ti^3 \ 
                                          +  S2f(1 - S2s)ts^2 ---------------------------------------------------- |
                                                                        ((ts + ti)^2 + (w.ts.ti)^2)^3              /


                   d2ti   /           1 - (w.ti)^2                       (tf + ti)^2 - (w.tf.ti)^2                          (ts + ti)^2 - (w.ts.ti)^2   \ 
               +  ------- | S2f.S2s ----------------  +  (1 - S2f)tf^2 -----------------------------  +  S2f(1 - S2s)ts^2 ----------------------------- |.
                  dGj.dGk \         (1 + (w.ti)^2)^2                   ((tf + ti)^2 + (w.tf.ti)^2)^2                      ((ts + ti)^2 + (w.ts.ti)^2)^2 /

calc_ellipsoid_S2f_tf_S2s_ts_d2jw_dGjdGk(data, params, j, k)

Spectral density Hessian.

Calculate the spectral desity values for the Gj - Gk double partial derivative of the extended model-free formula with the parameters {S2f, tf, S2s, ts} together with diffusion tensor parameters.

The model-free Hessian is:

                  _k_
    d2J(w)     2  \   /      dti   dti  /                  3 - (w.ti)^2                     (tf + ti)^3 + 3.w^2.tf^3.ti(tf + ti) - (w.tf)^4.ti^3
   -------  =  -   >  | -2ci --- . ---  | S2f.S2s.w^2.ti ----------------  +  (1 - S2f)tf^2 ----------------------------------------------------
   dGj.dGk     5  /__ \      dGj   dGk  \                (1 + (w.ti)^2)^3                             ((tf + ti)^2 + (w.tf.ti)^2)^3
                  i=-k

                                                              (ts + ti)^3 + 3.w^2.ts^3.ti(ts + ti) - (w.ts)^4.ti^3 \ 
                                          +  S2f(1 - S2s)ts^2 ---------------------------------------------------- |
                                                                        ((ts + ti)^2 + (w.ts.ti)^2)^3              /


                           / dti   dci     dti   dci         d2ti   \ /           1 - (w.ti)^2                       (tf + ti)^2 - (w.tf.ti)^2
                        +  | --- . ---  +  --- . ---  +  ci ------- | | S2f.S2s ----------------  +  (1 - S2f)tf^2 -----------------------------
                           \ dGj   dGk     dGk   dGj        dGj.dGk / \         (1 + (w.ti)^2)^2                   ((tf + ti)^2 + (w.tf.ti)^2)^2


                                                                                              (ts + ti)^2 - (w.ts.ti)^2   \ 
                                                                        +  S2f(1 - S2s)ts^2 ----------------------------- |
                                                                                            ((ts + ti)^2 + (w.ts.ti)^2)^2 /


                            d2ci        /   S2f.S2s          (1 - S2f)(tf + ti)tf         S2f(1 - S2s)(ts + ti)ts  \ \ 
                        +  ------- . ti | ------------  +  -------------------------  +  ------------------------- | |.
                           dGj.dGk      \ 1 + (w.ti)^2     (tf + ti)^2 + (w.tf.ti)^2     (ts + ti)^2 + (w.ts.ti)^2 / /

calc_diff_S2f_S2s_ts_d2jw_dGjdOj(data, params, j, k)

Spectral density Hessian.

Calculate the spectral desity values for the Gj - Oj double partial derivative of the extended model-free formula with the parameters {S2f, S2s, ts} together with diffusion tensor parameters.

The model-free Hessian is:

                  _k_
    d2J(w)     2  \   dci   dti  /           1 - (w.ti)^2                          (ts + ti)^2 - (w.ts.ti)^2   \ 
   -------  =  -   >  --- . ---  | S2f.S2s ----------------  +  S2f(1 - S2s)ts^2 ----------------------------- |.
   dGj.dOj     5  /__ dOj   dGj  \         (1 + (w.ti)^2)^2                      ((ts + ti)^2 + (w.ts.ti)^2)^2 /
                  i=-k

calc_ellipsoid_S2f_S2s_ts_d2jw_dGjdOj(data, params, j, k)

Spectral density Hessian.

Calculate the spectral desity values for the Gj - Oj double partial derivative of the extended model-free formula with the parameters {S2f, S2s, ts} together with diffusion tensor parameters.

The model-free Hessian is:

                  _k_
    d2J(w)     2  \   / dci   dti  /           1 - (w.ti)^2                          (ts + ti)^2 - (w.ts.ti)^2   \ 
   -------  =  -   >  | --- . ---  | S2f.S2s ----------------  +  S2f(1 - S2s)ts^2 ----------------------------- |
   dGj.dOj     5  /__ \ dOj   dGj  \         (1 + (w.ti)^2)^2                      ((ts + ti)^2 + (w.ts.ti)^2)^2 /
                  i=-k

                              d2ci        /   S2f.S2s         S2f(1 - S2s)(ts + ti)ts  \ \ 
                          +  ------- . ti | ------------  +  ------------------------- | |.
                             dGj.dOj      \ 1 + (w.ti)^2     (ts + ti)^2 + (w.ts.ti)^2 / /

calc_diff_S2f_tf_S2s_ts_d2jw_dGjdOj(data, params, j, k)

Spectral density Hessian.

Calculate the spectral desity values for the Gj - Oj double partial derivative of the extended model-free formula with the parameters {S2f, tf, S2s, ts} together with diffusion tensor parameters.

The model-free Hessian is:

                  _k_
    d2J(w)     2  \   dci   dti  /           1 - (w.ti)^2                       (tf + ti)^2 - (w.tf.ti)^2                          (ts + ti)^2 - (w.ts.ti)^2   \ 
   -------  =  -   >  --- . ---  | S2f.S2s ----------------  +  (1 - S2f)tf^2 -----------------------------  +  S2f(1 - S2s)ts^2 ----------------------------- |.
   dGj.dOj     5  /__ dOj   dGj  \         (1 + (w.ti)^2)^2                   ((tf + ti)^2 + (w.tf.ti)^2)^2                      ((ts + ti)^2 + (w.ts.ti)^2)^2 /
                  i=-k

calc_ellipsoid_S2f_tf_S2s_ts_d2jw_dGjdOj(data, params, j, k)

Spectral density Hessian.

Calculate the spectral desity values for the Gj - Oj double partial derivative of the extended model-free formula with the parameters {S2f, tf, S2s, ts} together with diffusion tensor parameters.

The model-free Hessian is:

                  _k_
    d2J(w)     2  \   / dci   dti  /           1 - (w.ti)^2                       (tf + ti)^2 - (w.tf.ti)^2                          (ts + ti)^2 - (w.ts.ti)^2   \ 
   -------  =  -   >  | --- . ---  | S2f.S2s ----------------  +  (1 - S2f)tf^2 -----------------------------  +  S2f(1 - S2s)ts^2 ----------------------------- |
   dGj.dOj     5  /__ \ dOj   dGj  \         (1 + (w.ti)^2)^2                   ((tf + ti)^2 + (w.tf.ti)^2)^2                      ((ts + ti)^2 + (w.ts.ti)^2)^2 /
                  i=-k

                              d2ci        /   S2f.S2s          (1 - S2f)(tf + ti)tf         S2f(1 - S2s)(ts + ti)ts  \ \ 
                          +  ------- . ti | ------------  +  -------------------------  +  ------------------------- | |.
                             dGj.dOj      \ 1 + (w.ti)^2     (tf + ti)^2 + (w.tf.ti)^2     (ts + ti)^2 + (w.ts.ti)^2 / /

calc_diff_S2f_S2s_ts_d2jw_dGjdS2f(data, params, j, k)

Spectral density Hessian.

Calculate the spectral desity values for the Gj - S2f double partial derivative of the extended model-free formula with the parameters {S2f, S2s, ts} together with diffusion tensor parameters.

The model-free Hessian is:

                   _k_
    d2J(w)      2  \        dti  /       1 - (w.ti)^2                       (ts + ti)^2 - (w.ts.ti)^2   \ 
   --------  =  -   >  ci . ---  | S2s ----------------  +  (1 - S2s)ts^2 ----------------------------- |.
   dGj.dS2f     5  /__      dGj  \     (1 + (w.ti)^2)^2                   ((ts + ti)^2 + (w.ts.ti)^2)^2 /
                   i=-k

calc_ellipsoid_S2f_S2s_ts_d2jw_dGjdS2f(data, params, j, k)

Spectral density Hessian.

Calculate the spectral desity values for the Gj - S2f double partial derivative of the extended model-free formula with the parameters {S2f, S2s, ts} together with diffusion tensor parameters.

The model-free Hessian is:

                   _k_
    d2J(w)      2  \   /      dti  /       1 - (w.ti)^2                       (ts + ti)^2 - (w.ts.ti)^2   \ 
   --------  =  -   >  | ci . ---  | S2s ----------------  +  (1 - S2s)ts^2 ----------------------------- |
   dGj.dS2f     5  /__ \      dGj  \     (1 + (w.ti)^2)^2                   ((ts + ti)^2 + (w.ts.ti)^2)^2 /
                   i=-k

                           dci      /     S2s            (1 - S2s)(ts + ti)ts    \ \ 
                        +  --- . ti | ------------  +  ------------------------- | |.
                           dGj      \ 1 + (w.ti)^2     (ts + ti)^2 + (w.ts.ti)^2 / /

calc_diff_S2f_tf_S2s_ts_d2jw_dGjdS2f(data, params, j, k)

Spectral density Hessian.

Calculate the spectral desity values for the Gj - S2f double partial derivative of the extended model-free formula with the parameters {S2f, tf, S2s, ts} together with diffusion tensor parameters.

The model-free Hessian is:

                   _k_
    d2J(w)      2  \        dti  /       1 - (w.ti)^2              (tf + ti)^2 - (w.tf.ti)^2                       (ts + ti)^2 - (w.ts.ti)^2   \ 
   --------  =  -   >  ci . ---  | S2s ----------------  -  tf^2 -----------------------------  +  (1 - S2s)ts^2 ----------------------------- |.
   dGj.dS2f     5  /__      dGj  \     (1 + (w.ti)^2)^2          ((tf + ti)^2 + (w.tf.ti)^2)^2                   ((ts + ti)^2 + (w.ts.ti)^2)^2 /
                   i=-k

calc_ellipsoid_S2f_tf_S2s_ts_d2jw_dGjdS2f(data, params, j, k)

Spectral density Hessian.

Calculate the spectral desity values for the Gj - S2f double partial derivative of the extended model-free formula with the parameters {S2f, tf, S2s, ts} together with diffusion tensor parameters.

The model-free Hessian is:

                   _k_
    d2J(w)      2  \   /      dti  /       1 - (w.ti)^2              (tf + ti)^2 - (w.tf.ti)^2                       (ts + ti)^2 - (w.ts.ti)^2   \ 
   --------  =  -   >  | ci . ---  | S2s ----------------  -  tf^2 -----------------------------  +  (1 - S2s)ts^2 ----------------------------- |
   dGj.dS2f     5  /__ \      dGj  \     (1 + (w.ti)^2)^2          ((tf + ti)^2 + (w.tf.ti)^2)^2                   ((ts + ti)^2 + (w.ts.ti)^2)^2 /
                   i=-k

                            dci      /     S2s                (tf + ti)tf               (1 - S2s)(ts + ti)ts    \ \ 
                         +  --- . ti | ------------  -  -------------------------  +  ------------------------- | |.
                            dGj      \ 1 + (w.ti)^2     (tf + ti)^2 + (w.tf.ti)^2     (ts + ti)^2 + (w.ts.ti)^2 / /

calc_diff_S2f_S2s_ts_d2jw_dGjdS2s(data, params, j, k)

Spectral density Hessian.

Calculate the spectral desity values for the Gj - S2s double partial derivative of the extended model-free formula with the parameters {S2f, S2s, ts} or {S2f, tf, S2s, ts} together with diffusion tensor parameters.

The model-free Hessian is:

                      _k_
    d2J(w)      2     \        dti  /   1 - (w.ti)^2              (ts + ti)^2 - (w.ts.ti)^2   \ 
   --------  =  - S2f  >  ci . ---  | ----------------  -  ts^2 ----------------------------- |.
   dGj.dS2s     5     /__      dGj  \ (1 + (w.ti)^2)^2          ((ts + ti)^2 + (w.ts.ti)^2)^2 /
                      i=-k

calc_ellipsoid_S2f_S2s_ts_d2jw_dGjdS2s(data, params, j, k)

Spectral density Hessian.

Calculate the spectral desity values for the Gj - S2s double partial derivative of the extended model-free formula with the parameters {S2f, S2s, ts} or {S2f, tf, S2s, ts} together with diffusion tensor parameters.

The model-free Hessian is:

                      _k_
    d2J(w)      2     \   /      dti  /   1 - (w.ti)^2              (ts + ti)^2 - (w.ts.ti)^2   \ 
   --------  =  - S2f  >  | ci . ---  | ----------------  -  ts^2 ----------------------------- |
   dGj.dS2s     5     /__ \      dGj  \ (1 + (w.ti)^2)^2          ((ts + ti)^2 + (w.ts.ti)^2)^2 /
                      i=-k

                           dci      /      1                 (ts + ti)ts         \ \ 
                        +  --- . ti | ------------  -  ------------------------- | |.
                           dGj      \ 1 + (w.ti)^2     (ts + ti)^2 + (w.ts.ti)^2 / /

calc_diff_S2f_tf_S2s_ts_d2jw_dGjdtf(data, params, j, k)

Spectral density Hessian.

Calculate the spectral desity values for the Gj - tf double partial derivative of the extended model-free formula with the parameters {S2f, tf, S2s, ts} together with diffusion tensor parameters.

The model-free Hessian is:

                                _k_
    d2J(w)     4                \        dti                   (tf + ti)^2 - 3(w.tf.ti)^2
   -------  =  - (1 - S2f) . tf  >  ci . --- . ti . (tf + ti) -----------------------------.
   dGj.dtf     5                /__      dGj                  ((tf + ti)^2 + (w.tf.ti)^2)^3
                                i=-k

calc_ellipsoid_S2f_tf_S2s_ts_d2jw_dGjdtf(data, params, j, k)

Spectral density Hessian.

Calculate the spectral desity values for the Gj - tf double partial derivative of the extended model-free formula with the parameters {S2f, tf, S2s, ts} together with diffusion tensor parameters.

The model-free Hessian is:

                           _k_
    d2J(w)     2           \   /       dti                        (tf + ti)^2 - 3(w.tf.ti)^2       dci          (tf + ti)^2 - (w.tf.ti)^2   \ 
   -------  =  - (1 - S2f)  >  | 2ci . --- . tf . ti . (tf + ti) -----------------------------  +  --- . ti^2 ----------------------------- |.
   dGj.dtf     5           /__ \       dGj                       ((tf + ti)^2 + (w.tf.ti)^2)^3     dGj        ((tf + ti)^2 + (w.tf.ti)^2)^2 /
                           i=-k

calc_diff_S2f_S2s_ts_d2jw_dGjdts(data, params, j, k)

Spectral density Hessian.

Calculate the spectral desity values for the Gj - ts double partial derivative of the extended model-free formula with the parameters {S2f, S2s, ts} or {S2f, tf, S2s, ts} together with diffusion tensor parameters.

The model-free Hessian is:

                                   _k_
    d2J(w)     4                   \        dti                   (ts + ti)^2 - 3(w.ts.ti)^2
   -------  =  - S2f(1 - S2s) . ts  >  ci . --- . ti . (ts + ti) -----------------------------.
   dGj.dts     5                   /__      dGj                  ((ts + ti)^2 + (w.ts.ti)^2)^3
                                   i=-k

calc_ellipsoid_S2f_S2s_ts_d2jw_dGjdts(data, params, j, k)

Spectral density Hessian.

Calculate the spectral desity values for the Gj - ts double partial derivative of the extended model-free formula with the parameters {S2f, S2s, ts} or {S2f, tf, S2s, ts} together with diffusion tensor parameters.

The model-free Hessian is:

                              _k_
    d2J(w)     2              \   /       dti                        (ts + ti)^2 - 3(w.ts.ti)^2       dci          (ts + ti)^2 - (w.ts.ti)^2   \ 
   -------  =  - S2f(1 - S2s)  >  | 2ci . --- . ts . ti . (ts + ti) -----------------------------  +  --- . ti^2 ----------------------------- |.
   dGj.dts     5              /__ \       dGj                       ((ts + ti)^2 + (w.ts.ti)^2)^3     dGj        ((ts + ti)^2 + (w.ts.ti)^2)^2 /
                              i=-k

calc_diff_S2f_S2s_ts_d2jw_dOjdOk(data, params, j, k)

Spectral density Hessian.

Calculate the spectral desity values for the Oj - Ok double partial derivative of the extended model-free formula with the parameters {S2f, S2s, ts} together with diffusion tensor parameters.

The model-free Hessian is:

                  _k_
    d2J(w)     2  \    d2ci        /  S2f . S2s        S2f(1 - S2s)(ts + ti)ts  \ 
   -------  =  -   >  ------- . ti | ------------  +  ------------------------- |.
   dOj.dOk     5  /__ dOj.dOk      \ 1 + (w.ti)^2     (ts + ti)^2 + (w.ts.ti)^2 /
                  i=-k

calc_diff_S2f_tf_S2s_ts_d2jw_dOjdOk(data, params, j, k)

Spectral density Hessian.

Calculate the spectral desity values for the Oj - Ok double partial derivative of the extended model-free formula with the parameters {S2f, tf, S2s, ts} together with diffusion tensor parameters.

The model-free Hessian is:

                  _k_
    d2J(w)     2  \    d2ci        /  S2f . S2s         (1 - S2f)(tf + ti)tf         S2f(1 - S2s)(ts + ti)ts  \ 
   -------  =  -   >  ------- . ti | ------------  +  -------------------------  +  ------------------------- |.
   dOj.dOk     5  /__ dOj.dOk      \ 1 + (w.ti)^2     (tf + ti)^2 + (w.tf.ti)^2     (ts + ti)^2 + (w.ts.ti)^2 /
                  i=-k

calc_diff_S2f_S2s_ts_d2jw_dOjdS2f(data, params, j, k)

Spectral density Hessian.

Calculate the spectral desity values for the Oj - S2f double partial derivative of the extended model-free formula with the parameters {S2f, S2s, ts} together with diffusion tensor parameters.

The model-free Hessian is:

                   _k_
    d2J(w)      2  \   dci      /      S2s           (1 - S2s)(ts + ti)ts    \ 
   --------  =  -   >  --- . ti | ------------  +  ------------------------- |.
   dOj.dS2f     5  /__ dOj      \ 1 + (w.ti)^2     (ts + ti)^2 + (w.ts.ti)^2 /
                   i=-k

calc_diff_S2f_tf_S2s_ts_d2jw_dOjdS2f(data, params, j, k)

Spectral density Hessian.

Calculate the spectral desity values for the Oj - S2f double partial derivative of the extended model-free formula with the parameters {S2f, tf, S2s, ts} together with diffusion tensor parameters.

The model-free Hessian is:

                   _k_
    d2J(w)      2  \   dci      /      S2s               (tf + ti)tf               (1 - S2s)(ts + ti)ts    \ 
   --------  =  -   >  --- . ti | ------------  -  -------------------------  +  ------------------------- |.
   dOj.dS2f     5  /__ dOj      \ 1 + (w.ti)^2     (tf + ti)^2 + (w.tf.ti)^2     (ts + ti)^2 + (w.ts.ti)^2 /
                   i=-k

calc_diff_S2f_tf_S2s_ts_d2jw_dOjdS2s(data, params, j, k)

Spectral density Hessian.

Calculate the spectral desity values for the Oj - S2 double partial derivative of the extended model-free formula with the parameters {S2f, S2s, ts} and {S2f, tf, S2s, ts} together with diffusion tensor parameters.

The model-free Hessian is:

                      _k_
    d2J(w)      2     \   dci      /      1                 (ts + ti)ts         \ 
   --------  =  - S2f  >  --- . ti | ------------  -  ------------------------- |.
   dOj.dS2s     5     /__ dOj      \ 1 + (w.ti)^2     (ts + ti)^2 + (w.ts.ti)^2 /
                      i=-k

calc_diff_S2f_tf_S2s_ts_d2jw_dOjdtf(data, params, j, k)

Spectral density Hessian.

Calculate the spectral desity values for the Oj - tf double partial derivative of the extended model-free formula with the parameters {S2f, tf, S2s, ts} together with diffusion tensor parameters.

The model-free Hessian is:

                           _k_
    d2J(w)     2           \   dci          (tf + ti)^2 - (w.tf.ti)^2
   -------  =  - (1 - S2f)  >  --- . ti^2 -----------------------------.
   dOj.dtf     5           /__ dOj        ((tf + ti)^2 + (w.tf.ti)^2)^2
                           i=-k

calc_diff_S2f_S2s_ts_d2jw_dOjdts(data, params, j, k)

Spectral density Hessian.

Calculate the spectral desity values for the Oj - ts double partial derivative of the extended model-free formula with the parameters {S2f, S2s, ts} or {S2f, tf, S2s, ts} together with diffusion tensor parameters.

The model-free Hessian is:

                              _k_
    d2J(w)     2              \   dci          (ts + ti)^2 - (w.ts.ti)^2
   -------  =  - S2f(1 - S2s)  >  --- . ti^2 -----------------------------.
   dOj.dts     5              /__ dOj        ((ts + ti)^2 + (w.ts.ti)^2)^2
                              i=-k

calc_S2f_S2s_ts_d2jw_dS2fdS2s(data, params, j, k)

Spectral density Hessian.

Calculate the spectral desity values for the S2f - S2s double partial derivative of the extended model-free formula with the parameters {S2f, S2s, ts} or {S2f, tf, S2s, ts} with or without diffusion tensor parameters.

The model-free Hessian is:

                    _k_
     d2J(w)      2  \           /      1                 (ts + ti).ts        \ 
   ---------  =  -   >  ci . ti | ------------  -  ------------------------- |.
   dS2f.dS2s     5  /__         \ 1 + (w.ti)^2     (ts + ti)^2 + (w.ts.ti)^2 /
                    i=-k

calc_S2f_tf_S2s_ts_d2jw_dS2fdtf(data, params, j, k)

Spectral density Hessian.

Calculate the spectral desity values for the S2f - tf double partial derivative of the extended model-free formula with the parameters {S2f, tf, S2s, ts} with or without diffusion tensor parameters.

The model-free Hessian is:

                     _k_
    d2J(w)        2  \               (tf + ti)^2 - (w.tf.ti)^2
   --------  =  - -   >  ci . ti^2 -----------------------------.
   dS2f.dtf       5  /__           ((tf + ti)^2 + (w.tf.ti)^2)^2
                     i=-k

calc_S2f_S2s_ts_d2jw_dS2fdts(data, params, j, k)

Spectral density Hessian.

Calculate the spectral desity values for the S2f - ts double partial derivative of the extended model-free formula with the parameters {S2f, S2s, ts} or {S2f, tf, S2s, ts} with or without diffusion tensor parameters.

The model-free Hessian is:

                            _k_
    d2J(w)      2           \               (ts + ti)^2 - (w.ts.ti)^2
   --------  =  - (1 - S2s)  >  ci . ti^2 -----------------------------.
   dS2f.dts     5           /__           ((ts + ti)^2 + (w.ts.ti)^2)^2
                            i=-k

calc_S2f_S2s_ts_d2jw_dS2sdts(data, params, j, k)

Spectral density Hessian.

Calculate the spectral desity values for the S2s - ts double partial derivative of the extended model-free formula with the parameters {S2f, S2s, ts} or {S2f, tf, S2s, ts} with or without diffusion tensor parameters.

The model-free Hessian is:

                        _k_
    d2J(w)        2     \               (ts + ti)^2 - (w.ts.ti)^2
   --------  =  - - S2f  >  ci . ti^2 -----------------------------.
   dS2s.dts       5     /__           ((ts + ti)^2 + (w.ts.ti)^2)^2
                        i=-k

calc_S2f_tf_S2s_ts_d2jw_dtf2(data, params, j, k)

Spectral density Hessian.

Calculate the spectral desity values for the tf - tf double partial derivative of the extended model-free formula with the parameters {S2f, tf, S2s, ts} with or without diffusion tensor parameters.

The model-free Hessian is:

                            _k_
   d2J(w)       4           \             (tf + ti)^3 + 3.w^2.ti^3.tf.(tf + ti) - (w.ti)^4.tf^3
   ------  =  - - (1 - S2f)  >  ci . ti^2 -----------------------------------------------------.
   dtf**2       5           /__                        ((tf + ti)^2 + (w.tf.ti)^2)^3
                            i=-k

calc_S2f_S2s_ts_d2jw_dts2(data, params, j, k)

Spectral density Hessian.

Calculate the spectral desity values for the ts - ts double partial derivative of the extended model-free formula with the parameters {S2f, S2s, ts} or {S2f, tf, S2s, ts} with or without diffusion tensor parameters.

The model-free Hessian is:

                               _k_
   d2J(w)       4              \             (ts + ti)^3 + 3.w^2.ti^3.ts.(ts + ti) - (w.ti)^4.ts^3
   ------  =  - - S2f(1 - S2s)  >  ci . ti^2 -----------------------------------------------------.
   dts**2       5              /__                        ((ts + ti)^2 + (w.ts.ti)^2)^3
                               i=-k