ed2b80836fb92142143d8b7702d5519761671cf5,dipy/reconst/dki.py,,_ols_iter,#Any#Any#Any#Any#,700

Before Change

    log_s = np.log(sig)
    result = np.dot(inv_design, log_s)
    D=result[:6]
    tensor = from_lower_triangular(D)
    MD_square = ((tensor[0,0] + tensor[1,1] + tensor[2,2])/3.)**2  
    K_elements=result[6:21]/MD_square
    
    return decompose_tensors(tensor, K_elements,

After Change

    evals, evecs = decompose_tensor(from_lower_triangular(DT_elements))

    // Extracting kurtosis tensor
    MD_square = (evals.mean(0))**2  
    KT_elements = result[6:21] / MD_square
    
    dki_params = np.concatenate((evals, evecs[0], evecs[1], evecs[2], 
                                 KT_elements), axis=0)

    return dki_params

Instances