Lines Matching refs:n_y
3761 PetscReal n_y[3][3] = { in TPSNearestPointResJac() local
3767 feval(yreal, &f, grad, n_y); in TPSNearestPointResJac()
3771 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()
3792 …for (PetscInt i = 0; i < 3; i++) nd_y[i] = n[i] + n_y[0][i] * d[0] + n_y[1][i] * d[1] + n_y[2][i] … in TPSNearestPointResJac()
3800 J[1 + j * 3] = (j == 1) * 1. - 2 * (n_y[1][j] * nd + n[1] * nd_y[j]); in TPSNearestPointResJac()
3801 J[2 + j * 3] = (j == 2) * 1. - 2 * (n_y[2][j] * nd + n[2] * nd_y[j]); in TPSNearestPointResJac()