Fix division by zero problem

bernhardbachmann · OpenModelica-Hudson · commit cd21fa9e4162 · 2017-10-09T18:24:18.000+02:00
diff --git a/SimulationRuntime/c/simulation/solver/nonlinearSolverHomotopy.c b/SimulationRuntime/c/simulation/solver/nonlinearSolverHomotopy.c
@@ -613,23 +613,22 @@ void vecMultScaling(int n, double *a, double *b, double *c)
 {
   int i;
   for (i=0;i<n;i++)
-    c[i] = a[i]*fabs(b[i]);
+    c[i] = (fabs(b[i])>0 ? a[i]*fabs(b[i]):a[i]);
 }
 
 void vecDivScaling(int n, double *a, double *b, double *c)
 {
   int i;
   for (i=0;i<n;i++)
-    c[i] = a[i]/fabs(b[i]);
-    // c[i] = a[i]/fmax(1.0,fabs(b[i]));
+    c[i] = (fabs(b[i])>0 ? a[i]/fabs(b[i]):a[i]);
 }
 
 void vecNormalize(int n, double *a, double *b)
 {
   int i;
   double norm = vec2Norm(n,a);
   for (i=0;i<n;i++)
-    b[i] = a[i]/norm;
+    b[i] = (norm>0 ? a[i]/norm:a[i]);
 }
 
 void vecConst(int n, double value, double *a)

Original file line number	Diff line number	Diff line change
`@@ -613,23 +613,22 @@ void vecMultScaling(int n, double a, double b, double *c)`
`613`	`613`	`{`
`614`	`614`	`int i;`
`615`	`615`	`for (i=0;i<n;i++)`
`616`		`- c[i] = a[i]*fabs(b[i]);`
	`616`	`+ c[i] = (fabs(b[i])>0 ? a[i]*fabs(b[i]):a[i]);`
`617`	`617`	`}`
`618`	`618`
`619`	`619`	`void vecDivScaling(int n, double a, double b, double *c)`
`620`	`620`	`{`
`621`	`621`	`int i;`
`622`	`622`	`for (i=0;i<n;i++)`
`623`		`- c[i] = a[i]/fabs(b[i]);`
`624`		`- // c[i] = a[i]/fmax(1.0,fabs(b[i]));`
	`623`	`+ c[i] = (fabs(b[i])>0 ? a[i]/fabs(b[i]):a[i]);`
`625`	`624`	`}`
`626`	`625`
`627`	`626`	`void vecNormalize(int n, double a, double b)`
`628`	`627`	`{`
`629`	`628`	`int i;`
`630`	`629`	`double norm = vec2Norm(n,a);`
`631`	`630`	`for (i=0;i<n;i++)`
`632`		`- b[i] = a[i]/norm;`
	`631`	`+ b[i] = (norm>0 ? a[i]/norm:a[i]);`
`633`	`632`	`}`
`634`	`633`
`635`	`634`	`void vecConst(int n, double value, double *a)`