matrix_real.cpp

Go to the documentation of this file.

/*
Created on Fri Jun 26 14:13:26 2020
Copyright 2020 Peter Rakyta, Ph.D.

Licensed under the Apache License, Version 2.0 (the "License");
you may not use this file except in compliance with the License.
You may obtain a copy of the License at

    http://www.apache.org/licenses/LICENSE-2.0

Unless required by applicable law or agreed to in writing, software
distributed under the License is distributed on an "AS IS" BASIS,
WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
See the License for the specific language governing permissions and
limitations under the License.


@author: Peter Rakyta, Ph.D.
*/
#include "matrix_real.h"
#include "matrix_real_float.h"
#include <cstring>
#include <iostream>
#include "tbb/tbb.h"
#include <math.h>


Matrix_real::Matrix_real() : matrix_base<double>() {

}

Matrix_real::Matrix_real( double* data_in, int rows_in, int cols_in) : matrix_base<double>(data_in, rows_in, cols_in) {

}


Matrix_real::Matrix_real( double* data_in, int rows_in, int cols_in, int stride_in) : matrix_base<double>(data_in, rows_in, cols_in, stride_in) {

}


Matrix_real::Matrix_real( int rows_in, int cols_in) : matrix_base<double>(rows_in, cols_in) {


}


Matrix_real::Matrix_real( int rows_in, int cols_in, int stride_in) : matrix_base<double>(rows_in, cols_in, stride_in) {

}



Matrix_real::Matrix_real(const Matrix_real &in) : matrix_base<double>(in)  {

}

Matrix_real&
Matrix_real::operator=(const Matrix_real &mtx) {
    matrix_base<double>::operator=(mtx);
    return *this;
}

Matrix_real
Matrix_real::copy() const {

  Matrix_real ret = Matrix_real(rows, cols, stride);

  // logical variable indicating whether the matrix needs to be conjugated in CBLAS operations
  ret.conjugated = conjugated;
  // logical variable indicating whether the matrix needs to be transposed in CBLAS operations
  ret.transposed = transposed;
  // logical value indicating whether the class instance is the owner of the stored data or not. (If true, the data array is released in the destructor)
  ret.owner = true;

  memcpy( ret.data, data, rows*stride*sizeof(double));

  return ret;

}

void
Matrix_real::copy_to(Matrix_real& target) const {

  matrix_base<double>::copy_to(target);

}

void
Matrix_real::copy_to(Matrix_real_float& target) const {

  if (target.rows != rows || target.cols != cols || target.stride != stride) {
    target = Matrix_real_float(rows, cols, stride);
  }

  for (int row=0; row<rows; row++) {
    for (int col=0; col<cols; col++) {
      const int idx = row*stride + col;
      target[idx] = static_cast<float>(data[idx]);
    }
  }

  if (target.is_conjugated() != conjugated) {
    target.conjugate();
  }
  if (target.is_transposed() != transposed) {
    target.transpose();
  }

}


bool
Matrix_real::isnan() {

    for (int idx=0; idx < rows*cols; idx++) {
        if ( std::isnan(data[idx]) ) {
            return true;
        }
    }

    return false;


}