Problema de iteración

[mikelromero]

Hola a todo el mundo. Estoy desarollando en C el metodo de aproximacion de funciones de Gauss-Newton. Mi problema es que al crear el while para iterar el vector optimo que estamos buscando, este no funciona y entra en un bucle infinito...
Este el codigo que he escrito, todas las matrices lo he comprobado y estan bien calculadas, el unico problema es al hacer la iteracion v=v+S. Agradeceria mucho si alguien me pudiera ayudar!


do{
L=cholesky(A,longueur1);

Lt=transposee(L,longueur1,longueur1);


S=algo_solution(L,Lt,C,longueur1);

int e;

for(e=0;e<longueur1;e++){

printf("%f " ,S [e]);
}
printf("\n");

for(int k=0; k<longueur1;k++){


v[k]=v[k]+S[k];
}


}while(sqrt(B[0]*B[0]+B[1]*B[1]+B[2]*B[2]+B[3]*B[3])>(10^-6));