353    float grad[3] = {0.,0.,0.};  in Basic_Grads_to_B7()  local
359          grad[j] = *(fptr + 3*i + j);  in Basic_Grads_to_B7()
360       k = GradConv_BmatA_from_Bsign( bseven[i+1], grad );  in Basic_Grads_to_B7()
401 int GradConv_Gsign_from_GmatA( float *grad, float *matr )  in GradConv_Gsign_from_GmatA()  argument
416          grad[i] = (float) sqrt(matr[i]);  in GradConv_Gsign_from_GmatA()
417          grad[i]*= signum[i];  in GradConv_Gsign_from_GmatA()
421          grad[i] = 0;  in GradConv_Gsign_from_GmatA()
429 int GradConv_GmatA_from_Gsign( float *matr, float *grad )  in GradConv_GmatA_from_Gsign()  argument
434       matr[i] =  grad[i]*grad[i];  in GradConv_GmatA_from_Gsign()
435    matr[3] = grad[0]*grad[1];  in GradConv_GmatA_from_Gsign()
436    matr[4] = grad[0]*grad[2];  in GradConv_GmatA_from_Gsign()
437    matr[5] = grad[1]*grad[2];  in GradConv_GmatA_from_Gsign()
443 int GradConv_BmatA_from_Bsign( float *matr, float *grad )  in GradConv_BmatA_from_Bsign()  argument
449       matr[i] =  grad[i]*grad[i];  in GradConv_BmatA_from_Bsign()
452    matr[3] = grad[0]*grad[1];  in GradConv_BmatA_from_Bsign()
453    matr[4] = grad[0]*grad[2];  in GradConv_BmatA_from_Bsign()
454    matr[5] = grad[1]*grad[2];  in GradConv_BmatA_from_Bsign()
525                      THD_3dim_dataset *D,  in RicianNoiseDWIs()  argument
535    float *grad;  in RicianNoiseDWIs()  local
541    grad = MRI_FLOAT_PTR(g);  in RicianNoiseDWIs()
550             gscale = sqrt( grad[3*i]*grad[3*i] +  in RicianNoiseDWIs()
551                            grad[3*i+1]*grad[3*i+1] +  in RicianNoiseDWIs()
552                            grad[3*i+2]*grad[3*i+2] ); //bval;  apr,2016  in RicianNoiseDWIs()
556             sig+= THD_get_voxel(D,k,0)*grad[3*i]*grad[3*i];  in RicianNoiseDWIs()
557             sig+= THD_get_voxel(D,k,2)*grad[3*i+1]*grad[3*i+1];  in RicianNoiseDWIs()
558             sig+= THD_get_voxel(D,k,5)*grad[3*i+2]*grad[3*i+2];  in RicianNoiseDWIs()
559             sig+= 2*THD_get_voxel(D,k,1)*grad[3*i]*grad[3*i+1];  in RicianNoiseDWIs()
560             sig+= 2*THD_get_voxel(D,k,3)*grad[3*i]*grad[3*i+2];  in RicianNoiseDWIs()
561             sig+= 2*THD_get_voxel(D,k,4)*grad[3*i+1]*grad[3*i+2];  in RicianNoiseDWIs()
578 int DT_TORTOISEtoAFNI(float **D,   in DT_TORTOISEtoAFNI()  argument
587       D[0][k] = THD_get_voxel(DTS,k,0)/Lscale;  in DT_TORTOISEtoAFNI()
588       D[1][k] = INV[0]*INV[1]*THD_get_voxel(DTS,k,3)/Lscale;  in DT_TORTOISEtoAFNI()
589       D[2][k] = THD_get_voxel(DTS,k,1)/Lscale;  in DT_TORTOISEtoAFNI()
590       D[3][k] = INV[0]*INV[2]*THD_get_voxel(DTS,k,4)/Lscale;  in DT_TORTOISEtoAFNI()
591       D[4][k] = INV[1]*INV[2]*THD_get_voxel(DTS,k,5)/Lscale;  in DT_TORTOISEtoAFNI()
592       D[5][k] = THD_get_voxel(DTS,k,2)/Lscale;  in DT_TORTOISEtoAFNI()