octave-nkf: src/xdiv.cc comparison

[project @ 2003-10-18 03:53:52 by jwe]

comparison

equal deleted inserted replaced

-:2633831cbeb5
+:79df15d4470c
 #include <cassert>
 #include "CMatrix.h"
 #include "dMatrix.h"
+#include "CNDArray.h"
+#include "dNDArray.h"
 #include "oct-cmplx.h"
 #include "quit.h"
 #include "error.h"
 #include "xdiv.h"
 }
 return result;
 }
+// Funny element by element division operations.
+//
+//          op2 \ op1:   s   cs
+//               +--   +---+----+
+//   N-d array         | 1 |  3 |
+//                     +---+----+
+//   complex N-d array | 2 |  4 |
+//                     +---+----+
+NDArray
+x_el_div (double a, const NDArray& b)
+{
+NDArray result (b.dims ());
+for (int i = 0; i < b.length (); i++)
+{
+OCTAVE_QUIT;
+result (i) = a / b (i);
+}
+return result;
+}
+ComplexNDArray
+x_el_div (double a, const ComplexNDArray& b)
+{
+ComplexNDArray result (b.dims ());
+for (int i = 0; i < b.length (); i++)
+{
+OCTAVE_QUIT;
+result (i) = a / b (i);
+}
+return result;
+}
+ComplexNDArray
+x_el_div (const Complex a, const NDArray& b)
+{
+ComplexNDArray result (b.dims ());
+for (int i = 0; i < b.length (); i++)
+{
+OCTAVE_QUIT;
+result (i) = a / b (i);
+}
+return result;
+}
+ComplexNDArray
+x_el_div (const Complex a, const ComplexNDArray& b)
+{
+ComplexNDArray result (b.dims ());
+for (int i = 0; i < b.length (); i++)
+{
+OCTAVE_QUIT;
+result (i) = a / b (i);
+}
+return result;
+}
 // Left division functions.
 //
 //       op2 \ op1:   m   cm
 //            +--   +---+----+
 //   matrix         | 1 |  3 |

Mercurial > hg > octave-nkf