Merge branch 'scalar-multiplication' into next

c/src/operation.c

@@ -24,6 +24,24 @@
 #include "operation.h"
 #include "solve.h"
 
+void scalar_multiplication(env envelope, float scalar)
+{
+        int n = envelope->n;
+        int i;
+
+        if (scalar == 0) {
+                end_envelope(envelope);
+                envelope = init_eye_envelope(n, 0);
+                return;
+        }
+
+        for (i = 0; i < n; i++)
+                envelope->diagg[i] *= scalar;
+
+        for (i = 0; i < envelope->enveCol[n]; i++)
+                envelope->enve[i] *= scalar;
+}
+
 float *multiply_triangles(env triL, env triU)
 {
         int d, i, j, k;

c/src/operation.h

@@ -1 +1,2 @@
 float *multiply_triangles(env triL, env triU);
+void scalar_multiplication(env envelope, float scalar);

python/src/operation.py

@@ -20,6 +20,37 @@
 
 import solve
 
+def scalarMultiplication(enve, m):
+    envelope = enve['enve']
+    ENVcol = enve['enveCol']
+    ENVlin = enve['enveLin']
+    diagonal = enve['diagg']
+
+    p = 0
+
+    if m != 0:
+        n = len(envelope)
+        while(p < n):
+            envelope[p] = envelope[p] * m
+            p = p + 1
+    else:
+        n = len(diagonal)
+        enve['enveCol'] = [0] * (n + 1)
+        enve['enve'] = [0]
+
+
+    p = 0
+    n = len(diagonal)
+
+    if m == 0:
+        enve['enveLin'] = []
+
+    while(p < n):
+        if m == 0:
+            enve['enveLin'].append(p)
+        diagonal[p] = diagonal[p] * m
+        p = p + 1
+
 def multiplyTriangles(triL, triU):
     n = len(triL['diagg'])