common.h

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.
*/
#ifndef common_H
#define common_H

#define _USE_MATH_DEFINES
#include <cmath>
#define NOMINMAX
#include <algorithm>

#include <omp.h>
#include "QGDTypes.h"
#include "dot.h"
#include "matrix_sparse.h"
#include "matrix_real.h"

class Matrix_float;

#include <string>
#include <stdio.h>
#include <iostream>
#include <vector>
#include <cstring>
#include <sstream>

#ifdef __cplusplus
extern "C"
{
#endif

#if BLAS==0 // undefined blas
    void omp_set_num_threads(int num_threads);
    int omp_get_max_threads();
#elif BLAS==1 // MKL
    void MKL_Set_Num_Threads(int num_threads);
    int mkl_get_max_threads();
#elif BLAS==2 // OpenBLAS
    void openblas_set_num_threads(int num_threads);
    int openblas_get_num_threads();
#endif

#ifdef __cplusplus
}
#endif



double activation_function( double Phi, int limit );


void fread_wrapper( void * buffer, size_t size, size_t count, FILE * stream);


void* qgd_calloc( int element_num, int size, int alignment );

void* qgd_realloc(void* aligned_ptr, int element_num, int size, int alignment );

void qgd_free( void* ptr );

int Power_of_2(int n);


void add_unique_elelement( std::vector<int>& involved_qbits, int qbit );

Matrix create_identity( int matrix_size );

Matrix_float create_identity_float( int matrix_size );

Matrix reduce_zgemm( std::vector<Matrix>& mtxs );

void subtract_diag( Matrix& mtx,  QGD_Complex16 scalar );


QGD_Complex16 mult( QGD_Complex16& a, QGD_Complex16& b );

QGD_Complex8 mult( QGD_Complex8& a, QGD_Complex8& b );

QGD_Complex16 mult( double a, QGD_Complex16 b );

QGD_Complex8 mult( float a, QGD_Complex8 b );

void mult( QGD_Complex16 a, Matrix& b );

void mult( QGD_Complex8 a, Matrix_float& b );


Matrix mult(Matrix_sparse a, Matrix& b);

double arg( const QGD_Complex16& a );






void conjugate_gradient(Matrix_real A, Matrix_real b, Matrix_real& x0, double tol);


#endif