<html><!-- Created using the cpp_pretty_printer from the dlib C++ library. See http://dlib.net for updates. --><head><title>dlib C++ Library - cblas_zhbmv.c</title></head><body bgcolor='white'><pre>
<font color='#009900'>/*
* cblas_zhbmv.c
* The program is a C interface to zhbmv
*
* Keita Teranishi 5/18/98
*
*/</font>
<font color='#0000FF'>#include</font> "<a style='text-decoration:none' href='cblas.h.html'>cblas.h</a>"
<font color='#0000FF'>#include</font> "<a style='text-decoration:none' href='cblas_f77.h.html'>cblas_f77.h</a>"
<font color='#0000FF'>#include</font> <font color='#5555FF'><</font>stdio.h<font color='#5555FF'>></font>
<font color='#0000FF'>#include</font> <font color='#5555FF'><</font>stdlib.h<font color='#5555FF'>></font>
<font color='#0000FF'><u>void</u></font> <b><a name='cblas_zhbmv'></a>cblas_zhbmv</b><font face='Lucida Console'>(</font><font color='#0000FF'>const</font> <font color='#0000FF'>enum</font> CBLAS_ORDER order,
<font color='#0000FF'>const</font> <font color='#0000FF'>enum</font> CBLAS_UPLO Uplo,<font color='#0000FF'>const</font> <font color='#0000FF'><u>int</u></font> N,<font color='#0000FF'>const</font> <font color='#0000FF'><u>int</u></font> K,
<font color='#0000FF'>const</font> <font color='#0000FF'><u>void</u></font> <font color='#5555FF'>*</font>alpha, <font color='#0000FF'>const</font> <font color='#0000FF'><u>void</u></font> <font color='#5555FF'>*</font>A, <font color='#0000FF'>const</font> <font color='#0000FF'><u>int</u></font> lda,
<font color='#0000FF'>const</font> <font color='#0000FF'><u>void</u></font> <font color='#5555FF'>*</font>X, <font color='#0000FF'>const</font> <font color='#0000FF'><u>int</u></font> incX, <font color='#0000FF'>const</font> <font color='#0000FF'><u>void</u></font> <font color='#5555FF'>*</font>beta,
<font color='#0000FF'><u>void</u></font> <font color='#5555FF'>*</font>Y, <font color='#0000FF'>const</font> <font color='#0000FF'><u>int</u></font> incY<font face='Lucida Console'>)</font>
<b>{</b>
<font color='#0000FF'><u>char</u></font> UL;
<font color='#0000FF'>#ifdef</font> F77_CHAR
F77_CHAR F77_UL;
<font color='#0000FF'>#else</font>
<font color='#0000FF'>#define</font> F77_UL <font color='#5555FF'>&</font>UL
<font color='#0000FF'>#endif</font>
<font color='#0000FF'>#ifdef</font> F77_INT
F77_INT F77_N<font color='#5555FF'>=</font>N, F77_K<font color='#5555FF'>=</font>K, F77_lda<font color='#5555FF'>=</font>lda, F77_incX<font color='#5555FF'>=</font>incX, F77_incY<font color='#5555FF'>=</font>incY;
<font color='#0000FF'>#else</font>
<font color='#0000FF'>#define</font> F77_N N
<font color='#0000FF'>#define</font> F77_K K
<font color='#0000FF'>#define</font> F77_lda lda
<font color='#0000FF'>#define</font> F77_incX incx
<font color='#0000FF'>#define</font> F77_incY incY
<font color='#0000FF'>#endif</font>
<font color='#0000FF'><u>int</u></font> n, i<font color='#5555FF'>=</font><font color='#979000'>0</font>, incx<font color='#5555FF'>=</font>incX;
<font color='#0000FF'>const</font> <font color='#0000FF'><u>double</u></font> <font color='#5555FF'>*</font>xx<font color='#5555FF'>=</font> <font face='Lucida Console'>(</font><font color='#0000FF'><u>double</u></font> <font color='#5555FF'>*</font><font face='Lucida Console'>)</font>X, <font color='#5555FF'>*</font>alp<font color='#5555FF'>=</font> <font face='Lucida Console'>(</font><font color='#0000FF'><u>double</u></font> <font color='#5555FF'>*</font><font face='Lucida Console'>)</font>alpha, <font color='#5555FF'>*</font>bet <font color='#5555FF'>=</font> <font face='Lucida Console'>(</font><font color='#0000FF'><u>double</u></font> <font color='#5555FF'>*</font><font face='Lucida Console'>)</font>beta;
<font color='#0000FF'><u>double</u></font> ALPHA[<font color='#979000'>2</font>],BETA[<font color='#979000'>2</font>];
<font color='#0000FF'><u>int</u></font> tincY, tincx;
<font color='#0000FF'><u>double</u></font> <font color='#5555FF'>*</font>x<font color='#5555FF'>=</font><font face='Lucida Console'>(</font><font color='#0000FF'><u>double</u></font> <font color='#5555FF'>*</font><font face='Lucida Console'>)</font>X, <font color='#5555FF'>*</font>y<font color='#5555FF'>=</font><font face='Lucida Console'>(</font><font color='#0000FF'><u>double</u></font> <font color='#5555FF'>*</font><font face='Lucida Console'>)</font>Y, <font color='#5555FF'>*</font>st<font color='#5555FF'>=</font><font color='#979000'>0</font>, <font color='#5555FF'>*</font>tx;
<font color='#0000FF'>if</font> <font face='Lucida Console'>(</font>order <font color='#5555FF'>=</font><font color='#5555FF'>=</font> CblasColMajor<font face='Lucida Console'>)</font>
<b>{</b>
<font color='#0000FF'>if</font> <font face='Lucida Console'>(</font>Uplo <font color='#5555FF'>=</font><font color='#5555FF'>=</font> CblasLower<font face='Lucida Console'>)</font> UL <font color='#5555FF'>=</font> '<font color='#FF0000'>L</font>';
<font color='#0000FF'>else</font> <font color='#0000FF'>if</font> <font face='Lucida Console'>(</font>Uplo <font color='#5555FF'>=</font><font color='#5555FF'>=</font> CblasUpper<font face='Lucida Console'>)</font> UL <font color='#5555FF'>=</font> '<font color='#FF0000'>U</font>';
<font color='#0000FF'>else</font>
<b>{</b>
<font color='#BB00BB'>cblas_xerbla</font><font face='Lucida Console'>(</font><font color='#979000'>2</font>, "<font color='#CC0000'>cblas_zhbmv</font>","<font color='#CC0000'>Illegal Uplo setting, %d\n</font>",Uplo <font face='Lucida Console'>)</font>;
<font color='#0000FF'>return</font>;
<b>}</b>
<font color='#0000FF'>#ifdef</font> F77_CHAR
F77_UL <font color='#5555FF'>=</font> <font color='#BB00BB'>C2F_CHAR</font><font face='Lucida Console'>(</font><font color='#5555FF'>&</font>UL<font face='Lucida Console'>)</font>;
<font color='#0000FF'>#endif</font>
<font color='#BB00BB'>F77_zhbmv</font><font face='Lucida Console'>(</font>F77_UL, <font color='#5555FF'>&</font>F77_N, <font color='#5555FF'>&</font>F77_K, alpha, A, <font color='#5555FF'>&</font>F77_lda, X,
<font color='#5555FF'>&</font>F77_incX, beta, Y, <font color='#5555FF'>&</font>F77_incY<font face='Lucida Console'>)</font>;
<b>}</b>
<font color='#0000FF'>else</font> <font color='#0000FF'>if</font> <font face='Lucida Console'>(</font>order <font color='#5555FF'>=</font><font color='#5555FF'>=</font> CblasRowMajor<font face='Lucida Console'>)</font>
<b>{</b>
ALPHA[<font color='#979000'>0</font>]<font color='#5555FF'>=</font> <font color='#5555FF'>*</font>alp;
ALPHA[<font color='#979000'>1</font>]<font color='#5555FF'>=</font> <font color='#5555FF'>-</font>alp[<font color='#979000'>1</font>];
BETA[<font color='#979000'>0</font>]<font color='#5555FF'>=</font> <font color='#5555FF'>*</font>bet;
BETA[<font color='#979000'>1</font>]<font color='#5555FF'>=</font> <font color='#5555FF'>-</font>bet[<font color='#979000'>1</font>];
<font color='#0000FF'>if</font> <font face='Lucida Console'>(</font>N <font color='#5555FF'>></font> <font color='#979000'>0</font><font face='Lucida Console'>)</font>
<b>{</b>
n <font color='#5555FF'>=</font> N <font color='#5555FF'><</font><font color='#5555FF'><</font> <font color='#979000'>1</font>;
x <font color='#5555FF'>=</font> <font color='#BB00BB'>malloc</font><font face='Lucida Console'>(</font>n<font color='#5555FF'>*</font><font color='#0000FF'>sizeof</font><font face='Lucida Console'>(</font><font color='#0000FF'><u>double</u></font><font face='Lucida Console'>)</font><font face='Lucida Console'>)</font>;
tx <font color='#5555FF'>=</font> x;
<font color='#0000FF'>if</font><font face='Lucida Console'>(</font> incX <font color='#5555FF'>></font> <font color='#979000'>0</font> <font face='Lucida Console'>)</font> <b>{</b>
i <font color='#5555FF'>=</font> incX <font color='#5555FF'><</font><font color='#5555FF'><</font> <font color='#979000'>1</font> ;
tincx <font color='#5555FF'>=</font> <font color='#979000'>2</font>;
st<font color='#5555FF'>=</font> x<font color='#5555FF'>+</font>n;
<b>}</b> <font color='#0000FF'>else</font> <b>{</b>
i <font color='#5555FF'>=</font> incX <font color='#5555FF'>*</font><font face='Lucida Console'>(</font><font color='#5555FF'>-</font><font color='#979000'>2</font><font face='Lucida Console'>)</font>;
tincx <font color='#5555FF'>=</font> <font color='#5555FF'>-</font><font color='#979000'>2</font>;
st <font color='#5555FF'>=</font> x<font color='#5555FF'>-</font><font color='#979000'>2</font>;
x <font color='#5555FF'>+</font><font color='#5555FF'>=</font><font face='Lucida Console'>(</font>n<font color='#5555FF'>-</font><font color='#979000'>2</font><font face='Lucida Console'>)</font>;
<b>}</b>
<font color='#0000FF'>do</font>
<b>{</b>
<font color='#5555FF'>*</font>x <font color='#5555FF'>=</font> <font color='#5555FF'>*</font>xx;
x[<font color='#979000'>1</font>] <font color='#5555FF'>=</font> <font color='#5555FF'>-</font>xx[<font color='#979000'>1</font>];
x <font color='#5555FF'>+</font><font color='#5555FF'>=</font> tincx ;
xx <font color='#5555FF'>+</font><font color='#5555FF'>=</font> i;
<b>}</b>
<font color='#0000FF'>while</font> <font face='Lucida Console'>(</font>x <font color='#5555FF'>!</font><font color='#5555FF'>=</font> st<font face='Lucida Console'>)</font>;
x<font color='#5555FF'>=</font>tx;
<font color='#0000FF'>#ifdef</font> F77_INT
F77_incX <font color='#5555FF'>=</font> <font color='#979000'>1</font>;
<font color='#0000FF'>#else</font>
incx <font color='#5555FF'>=</font> <font color='#979000'>1</font>;
<font color='#0000FF'>#endif</font>
<font color='#0000FF'>if</font><font face='Lucida Console'>(</font>incY <font color='#5555FF'>></font> <font color='#979000'>0</font><font face='Lucida Console'>)</font>
tincY <font color='#5555FF'>=</font> incY;
<font color='#0000FF'>else</font>
tincY <font color='#5555FF'>=</font> <font color='#5555FF'>-</font>incY;
y<font color='#5555FF'>+</font><font color='#5555FF'>+</font>;
i <font color='#5555FF'>=</font> tincY <font color='#5555FF'><</font><font color='#5555FF'><</font> <font color='#979000'>1</font>;
n <font color='#5555FF'>=</font> i <font color='#5555FF'>*</font> N ;
st <font color='#5555FF'>=</font> y <font color='#5555FF'>+</font> n;
<font color='#0000FF'>do</font> <b>{</b>
<font color='#5555FF'>*</font>y <font color='#5555FF'>=</font> <font color='#5555FF'>-</font><font face='Lucida Console'>(</font><font color='#5555FF'>*</font>y<font face='Lucida Console'>)</font>;
y <font color='#5555FF'>+</font><font color='#5555FF'>=</font> i;
<b>}</b> <font color='#0000FF'>while</font><font face='Lucida Console'>(</font>y <font color='#5555FF'>!</font><font color='#5555FF'>=</font> st<font face='Lucida Console'>)</font>;
y <font color='#5555FF'>-</font><font color='#5555FF'>=</font> n;
<b>}</b> <font color='#0000FF'>else</font>
x <font color='#5555FF'>=</font> <font face='Lucida Console'>(</font><font color='#0000FF'><u>double</u></font> <font color='#5555FF'>*</font><font face='Lucida Console'>)</font> X;
<font color='#0000FF'>if</font> <font face='Lucida Console'>(</font>Uplo <font color='#5555FF'>=</font><font color='#5555FF'>=</font> CblasUpper<font face='Lucida Console'>)</font> UL <font color='#5555FF'>=</font> '<font color='#FF0000'>L</font>';
<font color='#0000FF'>else</font> <font color='#0000FF'>if</font> <font face='Lucida Console'>(</font>Uplo <font color='#5555FF'>=</font><font color='#5555FF'>=</font> CblasLower<font face='Lucida Console'>)</font> UL <font color='#5555FF'>=</font> '<font color='#FF0000'>U</font>';
<font color='#0000FF'>else</font>
<b>{</b>
<font color='#BB00BB'>cblas_xerbla</font><font face='Lucida Console'>(</font><font color='#979000'>2</font>, "<font color='#CC0000'>cblas_zhbmv</font>","<font color='#CC0000'>Illegal Uplo setting, %d\n</font>", Uplo<font face='Lucida Console'>)</font>;
<font color='#0000FF'>return</font>;
<b>}</b>
<font color='#0000FF'>#ifdef</font> F77_CHAR
F77_UL <font color='#5555FF'>=</font> <font color='#BB00BB'>C2F_CHAR</font><font face='Lucida Console'>(</font><font color='#5555FF'>&</font>UL<font face='Lucida Console'>)</font>;
<font color='#0000FF'>#endif</font>
<font color='#BB00BB'>F77_zhbmv</font><font face='Lucida Console'>(</font>F77_UL, <font color='#5555FF'>&</font>F77_N, <font color='#5555FF'>&</font>F77_K, ALPHA,
A ,<font color='#5555FF'>&</font>F77_lda, x,<font color='#5555FF'>&</font>F77_incX, BETA, Y, <font color='#5555FF'>&</font>F77_incY<font face='Lucida Console'>)</font>;
<b>}</b>
<font color='#0000FF'>else</font>
<b>{</b>
<font color='#BB00BB'>cblas_xerbla</font><font face='Lucida Console'>(</font><font color='#979000'>1</font>, "<font color='#CC0000'>cblas_zhbmv</font>","<font color='#CC0000'>Illegal Order setting, %d\n</font>", order<font face='Lucida Console'>)</font>;
<font color='#0000FF'>return</font>;
<b>}</b>
<font color='#0000FF'>if</font> <font face='Lucida Console'>(</font> order <font color='#5555FF'>=</font><font color='#5555FF'>=</font> CblasRowMajor <font face='Lucida Console'>)</font>
<b>{</b>
<font color='#0000FF'>if</font><font face='Lucida Console'>(</font>X<font color='#5555FF'>!</font><font color='#5555FF'>=</font>x<font face='Lucida Console'>)</font>
<font color='#BB00BB'>free</font><font face='Lucida Console'>(</font>x<font face='Lucida Console'>)</font>;
<font color='#0000FF'>if</font> <font face='Lucida Console'>(</font>N <font color='#5555FF'>></font> <font color='#979000'>0</font><font face='Lucida Console'>)</font>
<b>{</b>
<font color='#0000FF'>do</font>
<b>{</b>
<font color='#5555FF'>*</font>y <font color='#5555FF'>=</font> <font color='#5555FF'>-</font><font face='Lucida Console'>(</font><font color='#5555FF'>*</font>y<font face='Lucida Console'>)</font>;
y <font color='#5555FF'>+</font><font color='#5555FF'>=</font> i;
<b>}</b>
<font color='#0000FF'>while</font> <font face='Lucida Console'>(</font>y <font color='#5555FF'>!</font><font color='#5555FF'>=</font> st<font face='Lucida Console'>)</font>;
<b>}</b>
<b>}</b>
<font color='#0000FF'>return</font>;
<b>}</b>
</pre></body></html>