Click here to Skip to main content
15,947,785 members
Please Sign up or sign in to vote.
1.00/5 (2 votes)
See more:
Hello!
Can anybody help me to implement this piece of code in C++?

Python
import numpy as np

M = np.exp(-3j * k * n / 4)
return np.dot(M, x)


"k", "n" and "x" here are matrices.

What I have tried:

complex<double> z = exp(-3.0 * 3 * 2 / 4);

This code works, but I cant do this with complex number.

Thank u in advance!
Posted
Updated 23-Apr-23 3:59am
v5

 
Share this answer
 
Comments
bartello 21-Apr-23 3:48am    
I have already looked into this lib. I couldn't find anything to solve my problem.
The actual problem is to calculate the exponent of a complex number multiplied by matrix.
Found a bypass.
separately multiple matrices, then multiple matrix by a complex number, then use "exp" from "complex" library.
It seems to me that this works well.
 
Share this answer
 
The data type std::complex is available in C++, but for 2D matrices a template is probably needed. If the required arithmetic operations are all available with the datatype complex, only the matrix multiplication would have to be done. Matrix multiplication and, if necessary, other matrix arithmetic operations can be added quickly. However, this looks the same as for other data types.

First a general matrix template would have to be created:
C++
template<typename T>
class Matrix {
public:
   // default constructor (empty)
   Matrix() : _rows(0), _cols(0) {};

   // constructor with specified size
   Matrix(int rows, int cols): _rows(rows), _cols(cols), _data(rows* cols) { };

   // constructor with initializer_list
   Matrix(std::initializer_list<std::initializer_list<T>> elements) :
      _rows(elements.size()), _cols(elements.begin()->size()), _data(_rows* _cols)
   {
      int row = 0, col = 0;
      for (auto row_elem : elements) {
         col = 0;
         for (auto val : row_elem) {
            _data[row * _cols + col] = val;
            ++col;
         }
         ++row;
      }
   }
   // Access functions (with protection)
   T& operator()(int i, int j) {
      if (i < 0 || i >= _rows || j < 0 || j >= _cols) {
         throw std::out_of_range("Matrix index out of range");
      } return _data[i * _cols + j];
   }

   const T& operator()(int i, int j) const {
      if (i < 0 || i >= _rows || j < 0 || j >= _cols) {
         throw std::out_of_range("Matrix index out of range");
      } return _data[i * _cols + j];
   }

   // Getter functions
   int rows() const { return _rows; }
   int cols() const { return _cols; }

   // matrix multiplication operator
   Matrix<T> operator*(const Matrix<T>& other) const {
      if (_cols != other._rows) {
         throw std::invalid_argument("Number of columns of first matrix must be equal to number of rows of second matrix");
      }
 
      Matrix<T> result(_rows, other._cols);

      for (int i = 0; i < _rows; i++) {
        for (int j = 0; j < other._cols; j++) {
          // TODO: Matrix multiplication
          }
      }
      return result;
   }

   // Overload the output operator << as a friend function inside the matrix class. 
   friend std::ostream& operator<<(std::ostream& os, const Matrix<T>& m) {
      if (m._data.empty() || m._rows == 0 || m._cols == 0) {
         return os << "";
      }
      for (int i = 0; i < m._rows; ++i) {
         for (int j = 0; j < m._cols; ++j) {
            os << m(i, j) << " ";
         }
         os << std::endl;
      }
      return os;
   }

private:
   int _rows, _cols;
   std::vector<T> _data;
};


A test program could then look like this:
C++
using namespace std::complex_literals;

typedef std::complex<double> MatComplex;

// Defining the matrices
Matrix<MatComplex> A = { {1, 2i, 3}, {4i, 5, 6} };
Matrix<MatComplex> B = { {7, 8i}, {9i, 10}, {11, 12i} };
Matrix<MatComplex> C;

std::cout << "Matrix A:\n" << A << std::endl;
std::cout << "Matrix B:\n" << B << std::endl;

C = A * B; // multiply the matrices

std::cout << "Matrix C = A*B:\n" << C << std::endl;
 
Share this answer
 
v2

This content, along with any associated source code and files, is licensed under The Code Project Open License (CPOL)



CodeProject, 20 Bay Street, 11th Floor Toronto, Ontario, Canada M5J 2N8 +1 (416) 849-8900