Adaptive.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 "Adaptive.h"
#include "common.h"

namespace {

Matrix_real transform_adaptive_parameter(const Matrix_real& parameters, int limit) {
    Matrix_real phi_transformed(1, 1);
    phi_transformed[0] = activation_function(parameters[0], limit);
    return phi_transformed;
}

Matrix_real_float transform_adaptive_parameter(const Matrix_real_float& parameters, int limit) {
    Matrix_real_float phi_transformed(1, 1);
    phi_transformed[0] = static_cast<float>(activation_function(static_cast<double>(parameters[0]), limit));
    return phi_transformed;
}

}



//static tbb::spin_mutex my_mutex;
Adaptive::Adaptive() : CRY() {

    // A string labeling the gate operation
    name = "Adaptive";

    // A string describing the type of the gate
    type = ADAPTIVE_OPERATION;

    limit = 1;

}



Adaptive::Adaptive(int qbit_num_in, int target_qbit_in, int control_qbit_in) : CRY(qbit_num_in, target_qbit_in, control_qbit_in) {

    // A string labeling the gate operation
    name = "Adaptive";

    // A string describing the type of the gate
    type = ADAPTIVE_OPERATION;

    limit = 1;
}


Adaptive::Adaptive(int qbit_num_in, int target_qbit_in, int control_qbit_in, int limit_in) : CRY(qbit_num_in, target_qbit_in, control_qbit_in) {

    // A string labeling the gate operation
    name = "Adaptive";

    // A string describing the type of the gate
    type = ADAPTIVE_OPERATION;

    limit = limit_in;
}

Adaptive::~Adaptive() {

}




void 
Adaptive::apply_to( Matrix_real& parameters, Matrix& input, int parallel ) {


    if (input.rows != matrix_size ) {
        std::string err("Adaptive::apply_to: Wrong input size in Adaptive gate apply.");
        throw err;    
    }

    Matrix_real phi_transformed = transform_adaptive_parameter(parameters, limit);
    Matrix_real phi_precomputed_sincos = compute_precomputed_sincos(phi_transformed);
    CRY::apply_to_inner(phi_transformed, phi_precomputed_sincos, input, parallel);



}


void
Adaptive::apply_to( Matrix_real_float& parameters, Matrix_float& input, int parallel ) {

    if (input.rows != matrix_size ) {
        std::string err("Adaptive::apply_to: Wrong input size in Adaptive gate apply.");
        throw err;
    }

    Matrix_real_float phi_transformed = transform_adaptive_parameter(parameters, limit);
    Matrix_real_float phi_precomputed_sincos = compute_precomputed_sincos(phi_transformed);
    CRY::apply_to_inner(phi_transformed, phi_precomputed_sincos, input, parallel);
}


void
Adaptive::apply_to_inner( Matrix_real& parameters, const Matrix_real& /*precomputed_sincos*/, Matrix& input, int parallel ) {

    Matrix_real phi_transformed = transform_adaptive_parameter(parameters, limit);
    Matrix_real phi_precomputed_sincos = compute_precomputed_sincos(phi_transformed);
    CRY::apply_to_inner(phi_transformed, phi_precomputed_sincos, input, parallel);
}


void
Adaptive::apply_to_inner( Matrix_real_float& parameters, const Matrix_real_float& /*precomputed_sincos*/, Matrix_float& input, int parallel ) {

    Matrix_real_float phi_transformed = transform_adaptive_parameter(parameters, limit);
    Matrix_real_float phi_precomputed_sincos = compute_precomputed_sincos(phi_transformed);
    CRY::apply_to_inner(phi_transformed, phi_precomputed_sincos, input, parallel);
}



void 
Adaptive::apply_from_right( Matrix_real& parameters, Matrix& input ) {


    if (input.cols != matrix_size ) {
        std::stringstream sstream;
     sstream << "Wrong matrix size in Adaptive apply_from_right" << std::endl;
        print(sstream, 0);    
        exit(-1);
    }

    Matrix_real phi_transformed = transform_adaptive_parameter(parameters, limit);
    Matrix_real phi_precomputed_sincos = compute_precomputed_sincos(phi_transformed);
    CRY::apply_from_right_inner(phi_transformed, phi_precomputed_sincos, input);


}


void
Adaptive::apply_from_right( Matrix_real_float& parameters, Matrix_float& input ) {

    if (input.cols != matrix_size ) {
        std::string err("Wrong matrix size in Adaptive apply_from_right");
        throw err;
    }

    Matrix_real_float phi_transformed = transform_adaptive_parameter(parameters, limit);
    Matrix_real_float phi_precomputed_sincos = compute_precomputed_sincos(phi_transformed);
    CRY::apply_from_right_inner(phi_transformed, phi_precomputed_sincos, input);
}


void
Adaptive::apply_from_right_inner( Matrix_real& parameters, const Matrix_real& /*precomputed_sincos*/, Matrix& input ) {

    Matrix_real phi_transformed = transform_adaptive_parameter(parameters, limit);
    Matrix_real phi_precomputed_sincos = compute_precomputed_sincos(phi_transformed);
    CRY::apply_from_right_inner(phi_transformed, phi_precomputed_sincos, input);
}


void
Adaptive::apply_from_right_inner( Matrix_real_float& parameters, const Matrix_real_float& /*precomputed_sincos*/, Matrix_float& input ) {

    Matrix_real_float phi_transformed = transform_adaptive_parameter(parameters, limit);
    Matrix_real_float phi_precomputed_sincos = compute_precomputed_sincos(phi_transformed);
    CRY::apply_from_right_inner(phi_transformed, phi_precomputed_sincos, input);
}



std::vector<Matrix>
Adaptive::apply_derivate_to( Matrix_real& parameters, Matrix& input, int parallel ) {

    if (input.rows != matrix_size ) {
        std::stringstream sstream;
     sstream << "Wrong matrix size in Adaptive gate apply" << std::endl;
        print(sstream, 0);         
        exit(-1);
    }

    Matrix_real Phi_transformed = transform_adaptive_parameter(parameters, limit);
//    Phi_transformed[0] = Phi;
//    Phi_transformed[0] = 0.5*(1.0-std::cos(Phi))*M_PI;

/*
Phi = Phi + M_PI;
Phi = (1.0-std::cos(Phi/2))*M_PI;
Phi_transformed[0] = Phi - M_PI;
*/

    return CRY::apply_derivate_to( Phi_transformed, input, parallel );



}


std::vector<Matrix_float>
Adaptive::apply_derivate_to( Matrix_real_float& parameters, Matrix_float& input, int parallel ) {

    if (input.rows != matrix_size ) {
        std::stringstream sstream;
     sstream << "Wrong matrix size in Adaptive gate apply" << std::endl;
        print(sstream, 0);         
        exit(-1);
    }

    Matrix_real_float Phi_transformed = transform_adaptive_parameter(parameters, limit);
//    Phi_transformed[0] = Phi;
//    Phi_transformed[0] = 0.5*(1.0-std::cos(Phi))*M_PI;

/*
Phi = Phi + M_PI;
Phi = (1.0-std::cos(Phi/2))*M_PI;
Phi_transformed[0] = Phi - M_PI;
*/

    return CRY::apply_derivate_to( Phi_transformed, input, parallel );
}

Adaptive* Adaptive::clone() {

    Adaptive* ret = new Adaptive(qbit_num, target_qbit, control_qbit);
    
    ret->set_parameter_start_idx( get_parameter_start_idx() );
    ret->set_parents( parents );
    ret->set_children( children );

    ret->set_limit( limit );


    return ret;

}


void 
Adaptive::set_limit( int limit_in ) {

    limit = limit_in;

}


int 
Adaptive::get_limit() {

    return limit;
}