/*--------------------------------------------------------------------------------
QR分解を用いた最小二乗法
least square method using QR decomposition
main.h
--------------------------------------------------------------------------------*/
#include <stdio.h>
#include <stdlib.h>
#include <math.h>
#define ROW_NUM 3 // 行数, number of rows
#define COL_NUM 2 // 列数, number of number of columns
#define DELTA 1e-20 // 入力する行列の要素の値, value of the element of the matrix
/*------------------------------------------------------------------------*/
extern void householder_QR_decomp( double X[ROW_NUM][COL_NUM], double Q[ROW_NUM][ROW_NUM], double R[ROW_NUM][COL_NUM] );
extern void show_XQR( double X[ROW_NUM][COL_NUM], double Q[ROW_NUM][ROW_NUM], double R[ROW_NUM][COL_NUM] );
extern void solve_A_using_QR( double y[ROW_NUM], double Q[ROW_NUM][ROW_NUM], double R[ROW_NUM][COL_NUM], double A[COL_NUM] );
/*------------------------------------------------------------------------*/