sffitcrs (2.0-git)
index
user/roman/Mfitcrs.c

 
Synopsis
        sffitcrs < in.rsf > out.rsf > out_tcrs_params.rsf tcrs=

Input: T[m][h][t] and its sf-file m[], Nm h[], Nh
Output: a(3) three CRS parameters: a(0) m + a(1)m^2 + a(2) h^2

H_11 = 2 Nx Sum(m^2)
H_12 = 2 Nx Sum(m^4)
H_13 = 2 Sum(x^2) Sum(m)

H_22 = 2 Nx Sum(m^4)
H_23 = 2 Sum(x^2) Sum(m^2)

H_33 = 2 Nm Sum(x^2)

da = - grad H^{-1}

 
Parameters
       
 
string tcrs=