84         int i,j,k;  in arm_mat_cholesky_f64()  local
103         for(i=0 ;i < n ; i++)  in arm_mat_cholesky_f64()
105             for(j=i ; j+3 < n ; j+=4)  in arm_mat_cholesky_f64()
107                 pG[(j + 0) * n + i] = pA[(j + 0) * n + i];  in arm_mat_cholesky_f64()
108                 pG[(j + 1) * n + i] = pA[(j + 1) * n + i];  in arm_mat_cholesky_f64()
109                 pG[(j + 2) * n + i] = pA[(j + 2) * n + i];  in arm_mat_cholesky_f64()
110                 pG[(j + 3) * n + i] = pA[(j + 3) * n + i];  in arm_mat_cholesky_f64()
117                 kCnt = i >> 1U;  in arm_mat_cholesky_f64()
122                     vecGi=vld1q_f64(&pG[i * n + k]);  in arm_mat_cholesky_f64()
145                 kCnt = i & 1;  in arm_mat_cholesky_f64()
149                     sum0 = sum0 + pG[i * n + k] * pG[(j + 0) * n + k];  in arm_mat_cholesky_f64()
150                     sum1 = sum1 + pG[i * n + k] * pG[(j + 1) * n + k];  in arm_mat_cholesky_f64()
151                     sum2 = sum2 + pG[i * n + k] * pG[(j + 2) * n + k];  in arm_mat_cholesky_f64()
152                     sum3 = sum3 + pG[i * n + k] * pG[(j + 3) * n + k];  in arm_mat_cholesky_f64()
157                 pG[(j + 0) * n + i] -= sum0;  in arm_mat_cholesky_f64()
158                 pG[(j + 1) * n + i] -= sum1;  in arm_mat_cholesky_f64()
159                 pG[(j + 2) * n + i] -= sum2;  in arm_mat_cholesky_f64()
160                 pG[(j + 3) * n + i] -= sum3;  in arm_mat_cholesky_f64()
165                 pG[j * n + i] = pA[j * n + i];  in arm_mat_cholesky_f64()
169                 kCnt = i >> 1U;  in arm_mat_cholesky_f64()
174                     vecGi=vld1q_f64(&pG[i * n + k]);  in arm_mat_cholesky_f64()
186                 kCnt = i & 1;  in arm_mat_cholesky_f64()
189                     sum = sum + pG[i * n + k] * pG[(j + 0) * n + k];  in arm_mat_cholesky_f64()
196                 pG[j * n + i] -= sum;  in arm_mat_cholesky_f64()
199             if (pG[i * n + i] <= 0.0)  in arm_mat_cholesky_f64()
204             invSqrtVj = 1.0/sqrt(pG[i * n + i]);  in arm_mat_cholesky_f64()
205             SCALE_COL_F64(pDst,i,invSqrtVj,i);  in arm_mat_cholesky_f64()
240         int i,j,k;  in arm_mat_cholesky_f64()  local
249         for(i=0 ; i < n ; i++)  in arm_mat_cholesky_f64()
251             for(j=i ; j < n ; j++)  in arm_mat_cholesky_f64()
253                 pG[j * n + i] = pA[j * n + i];  in arm_mat_cholesky_f64()
255                 for(k=0; k < i ; k++)  in arm_mat_cholesky_f64()
257                     pG[j * n + i] = pG[j * n + i] - pG[i * n + k] * pG[j * n + k];  in arm_mat_cholesky_f64()
261             if (pG[i * n + i] <= 0.0)  in arm_mat_cholesky_f64()
266             invSqrtVj = 1.0/sqrt(pG[i * n + i]);  in arm_mat_cholesky_f64()
267             SCALE_COL_F64(pDst,i,invSqrtVj,i);  in arm_mat_cholesky_f64()