Searched refs:norm_y (Results 1 – 1 of 1) sorted by relevance
3760 PetscReal f, grad[3], n[3], norm, norm_y[3], nd, nd_y[3], sign; in TPSNearestPointResJac() local3771 for (PetscInt i = 0; i < 3; i++) norm_y[i] = 1. / norm * n[i] * n_y[i][i]; in TPSNearestPointResJac()3776 for (PetscInt i = 0; i < 3; i++) n_y[0][i] += norm_y[i] * sign; in TPSNearestPointResJac()3779 norm_y[0] = 1. / norm * (n[0] * n_y[0][0]); in TPSNearestPointResJac()3780 norm_y[1] = 1. / norm * (n[0] * n_y[0][1] + n[1] * n_y[1][1]); in TPSNearestPointResJac()3781 norm_y[2] = 1. / norm * (n[0] * n_y[0][2] + n[2] * n_y[2][2]); in TPSNearestPointResJac()3787 n_y[i][j] = n_y[i][j] / norm - n[i] / norm * norm_y[j]; in TPSNearestPointResJac()