Sub_Matrix_Decomposition.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 SUB_MATRIX_DECOMPOSITION_H
#define SUB_MATRIX_DECOMPOSITION_H

#include "Decomposition_Base.h"



class Sub_Matrix_Decomposition : public Decomposition_Base {


public:

    bool subdisentaglement_done;

    Matrix subdecomposed_mtx;

    bool optimize_layer_num;

    std::map<int,int> identical_blocks;

    Gates_block* unit_gate_structure;

    long long max_inner_iterations;


public:

Sub_Matrix_Decomposition( );

Sub_Matrix_Decomposition( Matrix Umtx_in, int qbit_num_in, bool optimize_layer_num_in, std::map<std::string, Config_Element>& config_in, guess_type initial_guess_in );

~Sub_Matrix_Decomposition();


void disentangle_submatrices();

virtual void add_gate_layers();


void add_default_gate_layers();


void set_custom_gate_layers( Gates_block* block_in);



void solve_layer_optimization_problem( int num_of_parameters, Matrix_real solution_guess_gsl);




double optimization_problem( double* parameters);



static double optimization_problem( Matrix_real parameters, void* void_instance );


static void optimization_problem_grad( Matrix_real parameters, void* void_instance, Matrix_real& grad  );


static void optimization_problem_combined( Matrix_real parameters, void* void_instance, double* f0, Matrix_real& grad );

int set_identical_blocks( int qbit, int identical_blocks_in );


int set_identical_blocks( std::map<int, int> identical_blocks_in );


Sub_Matrix_Decomposition* clone();


};


#endif //SUB_MATRIX_DECOMPOSITION