Changeset 2773

Timestamp:

01/06/10 15:05:10 (15 years ago)

Author:

Mathieu Morlighem

Message:

forgot to normalize vector

File:

• issm/trunk/src/c/parallel/control_core.cpp (modified) (1 diff)

issm/trunk/src/c/parallel/control_core.cpp

@@ -116 +116 @@
                 /*Normalize if last gradient not satisfying (search_scalar==0)*/
                 if (n>0 && search_scalar==0){
-                        _printf_("%s","      normalizing directions...");
+                        _printf_("%s","      orthogonalization...");
                         Orthx(&new_grad_g,grad_g,grad_g_old);
                         _printf_("%s\n"," done.");
                 }
                 else{ 
-                        VecDuplicate(grad_g,&new_grad_g);
-                        VecCopy(grad_g,new_grad_g);
+                        _printf_("%s","      normalizing directions...");
+                        Orthx(&new_grad_g,grad_g,NULL);
+                        _printf_("%s\n"," done.");
                 }
                 VecFree(&grad_g); VecFree(&grad_g_old);