/****************************************************************************
 *
 * bidiagonal.h
 *
 * Abstract
 *   Define prototypes of functions for bidiagonalization.
 *
 * Student:  Guohong Liu
 *
 * ID:       0385117
 *
 ****************************************************************************/

#ifndef __BIDIAGONAL_H
#define __BIDIAGONAL_H

#include "computing.h"

int bidiagonal (DoubleComplexMat *A);
int brd (int n, DoubleComplex *A_m, DoubleReal *d_m, DoubleReal *b_m,
         DoubleComplex *Q_m, DoubleComplex *P_m);
int brdErrChk (int n, DoubleComplex *A_m, DoubleReal *d_m, DoubleReal *b_m, 
               DoubleComplex *Q_m, DoubleComplex *P_m); 

#endif