.\" -*- mode: troff; coding: utf-8 -*-
.\" Automatically generated by Pod::Man 5.01 (Pod::Simple 3.43)
.\"
.\" Standard preamble:
.\" ========================================================================
.de Sp \" Vertical space (when we can't use .PP)
.if t .sp .5v
.if n .sp
..
.de Vb \" Begin verbatim text
.ft CW
.nf
.ne \\$1
..
.de Ve \" End verbatim text
.ft R
.fi
..
.\" \*(C` and \*(C' are quotes in nroff, nothing in troff, for use with C<>.
.ie n \{\
. ds C` ""
. ds C' ""
'br\}
.el\{\
. ds C`
. ds C'
'br\}
.\"
.\" Escape single quotes in literal strings from groff's Unicode transform.
.ie \n(.g .ds Aq \(aq
.el .ds Aq '
.\"
.\" If the F register is >0, we'll generate index entries on stderr for
.\" titles (.TH), headers (.SH), subsections (.SS), items (.Ip), and index
.\" entries marked with X<> in POD. Of course, you'll have to process the
.\" output yourself in some meaningful fashion.
.\"
.\" Avoid warning from groff about undefined register 'F'.
.de IX
..
.nr rF 0
.if \n(.g .if rF .nr rF 1
.if (\n(rF:(\n(.g==0)) \{\
. if \nF \{\
. de IX
. tm Index:\\$1\t\\n%\t"\\$2"
..
. if !\nF==2 \{\
. nr % 0
. nr F 2
. \}
. \}
.\}
.rr rF
.\" ========================================================================
.\"
.IX Title "Math::GSL::BLAS 3pm"
.TH Math::GSL::BLAS 3pm 2024-03-07 "perl v5.38.2" "User Contributed Perl Documentation"
.\" For nroff, turn off justification. Always turn off hyphenation; it makes
.\" way too many mistakes in technical documents.
.if n .ad l
.nh
.SH NAME
Math::GSL::BLAS \- Basic Linear Algebra Subprograms
.SH SYNOPSIS
.IX Header "SYNOPSIS"
.Vb 2
\& use Math::GSL::BLAS qw/:all/;
\& use Math::GSL::Matrix qw/:all/;
\&
\& # matrix\-matrix product of double numbers
\& my $A = Math::GSL::Matrix\->new(2,2);
\& $A\->set_row(0, [1, 4]);
\& \->set_row(1, [3, 2]);
\& my $B = Math::GSL::Matrix\->new(2,2);
\& $B\->set_row(0, [2, 1]);
\& \->set_row(1, [5,3]);
\& my $C = Math::GSL::Matrix\->new(2,2);
\& gsl_matrix_set_zero($C\->raw);
\& gsl_blas_dgemm($CblasNoTrans, $CblasNoTrans, 1, $A\->raw, $B\->raw, 1, $C\->raw);
\& my @got = $C\->row(0)\->as_list;
\& print "The resulting matrix is: \en[";
\& print "$got[0] $got[1]\en";
\& @got = $C\->row(1)\->as_list;
\& print "$got[0] $got[1] ]\en";
\&
\& # compute the scalar product of two vectors :
\& use Math::GSL::Vector qw/:all/;
\& use Math::GSL::CBLAS qw/:all/;
\& my $vec1 = Math::GSL::Vector\->new([1,2,3,4,5]);
\& my $vec2 = Math::GSL::Vector\->new([5,4,3,2,1]);
\& my ($status, $result) = gsl_blas_ddot($vec1\->raw, $vec2\->raw);
\& if($status == 0) {
\& print "The function has succeeded. \en";
\& }
\& print "The result of the vector multiplication is $result.\en";
.Ve
.SH DESCRIPTION
.IX Header "DESCRIPTION"
The functions of this module are divised into 3 levels:
.SS "Level 1 \- Vector operations"
.IX Subsection "Level 1 - Vector operations"
.ie n .IP """gsl_blas_sdsdot""" 3
.el .IP \f(CWgsl_blas_sdsdot\fR 3
.IX Item "gsl_blas_sdsdot"
.PD 0
.ie n .IP """gsl_blas_dsdot""" 3
.el .IP \f(CWgsl_blas_dsdot\fR 3
.IX Item "gsl_blas_dsdot"
.ie n .IP """gsl_blas_sdot""" 3
.el .IP \f(CWgsl_blas_sdot\fR 3
.IX Item "gsl_blas_sdot"
.ie n .IP """gsl_blas_ddot($x, $y)""" 3
.el .IP "\f(CWgsl_blas_ddot($x, $y)\fR" 3
.IX Item "gsl_blas_ddot($x, $y)"
.PD
This function computes the scalar product x^T y for the vectors \f(CW$x\fR and \f(CW$y\fR. The
function returns two values, the first is 0 if the operation succeeded, 1
otherwise and the second value is the result of the computation.
.ie n .IP """gsl_blas_cdotu""" 3
.el .IP \f(CWgsl_blas_cdotu\fR 3
.IX Item "gsl_blas_cdotu"
.PD 0
.ie n .IP """gsl_blas_cdotc""" 3
.el .IP \f(CWgsl_blas_cdotc\fR 3
.IX Item "gsl_blas_cdotc"
.ie n .IP """gsl_blas_zdotu($x, $y, $dotu)""" 3
.el .IP "\f(CWgsl_blas_zdotu($x, $y, $dotu)\fR" 3
.IX Item "gsl_blas_zdotu($x, $y, $dotu)"
.PD
This function computes the complex scalar product x^T y for the complex vectors
\&\f(CW$x\fR and \f(CW$y\fR, returning the result in the complex number \f(CW$dotu\fR. The function
returns 0 if the operation succeeded, 1 otherwise.
.ie n .IP """gsl_blas_zdotc($x, $y, $dotc)""" 3
.el .IP "\f(CWgsl_blas_zdotc($x, $y, $dotc)\fR" 3
.IX Item "gsl_blas_zdotc($x, $y, $dotc)"
This function computes the complex conjugate scalar product x^H y for the
complex vectors \f(CW$x\fR and \f(CW$y\fR, returning the result in the complex number \f(CW$dotc\fR.
The function returns 0 if the operation succeeded, 1 otherwise.
.ie n .IP """gsl_blas_snrm2"" =item ""gsl_blas_sasum""" 3
.el .IP "\f(CWgsl_blas_snrm2\fR =item \f(CWgsl_blas_sasum\fR" 3
.IX Item "gsl_blas_snrm2 =item gsl_blas_sasum"
.PD 0
.ie n .IP gsl_blas_dnrm2($x) 3
.el .IP \f(CWgsl_blas_dnrm2($x)\fR 3
.IX Item "gsl_blas_dnrm2($x)"
.PD
This function computes the Euclidean norm
.Sp
.Vb 1
\& ||x||_2 = \esqrt {\esum x_i^2}
.Ve
.Sp
of the vector \f(CW$x\fR.
.ie n .IP gsl_blas_dasum($x) 3
.el .IP \f(CWgsl_blas_dasum($x)\fR 3
.IX Item "gsl_blas_dasum($x)"
This function computes the absolute sum \esum |x_i| of the elements of the vector \f(CW$x\fR.
.ie n .IP """gsl_blas_scnrm2""" 3
.el .IP \f(CWgsl_blas_scnrm2\fR 3
.IX Item "gsl_blas_scnrm2"
.PD 0
.ie n .IP """gsl_blas_scasum""" 3
.el .IP \f(CWgsl_blas_scasum\fR 3
.IX Item "gsl_blas_scasum"
.ie n .IP gsl_blas_dznrm2($x) 3
.el .IP \f(CWgsl_blas_dznrm2($x)\fR 3
.IX Item "gsl_blas_dznrm2($x)"
.PD
This function computes the Euclidean norm of the complex vector \f(CW$x\fR, ||x||_2 = \esqrt {\esum (\eRe(x_i)^2 + \eIm(x_i)^2)}.
.ie n .IP gsl_blas_dzasum($x) 3
.el .IP \f(CWgsl_blas_dzasum($x)\fR 3
.IX Item "gsl_blas_dzasum($x)"
This function computes the sum of the magnitudes of the real and imaginary parts of the complex vector \f(CW$x\fR, \esum |\eRe(x_i)| + |\eIm(x_i)|.
.ie n .IP """gsl_blas_isamax""" 3
.el .IP \f(CWgsl_blas_isamax\fR 3
.IX Item "gsl_blas_isamax"
.PD 0
.ie n .IP """gsl_blas_idamax""" 3
.el .IP \f(CWgsl_blas_idamax\fR 3
.IX Item "gsl_blas_idamax"
.ie n .IP """gsl_blas_icamax""" 3
.el .IP \f(CWgsl_blas_icamax\fR 3
.IX Item "gsl_blas_icamax"
.ie n .IP """gsl_blas_izamax """ 3
.el .IP "\f(CWgsl_blas_izamax \fR" 3
.IX Item "gsl_blas_izamax "
.ie n .IP """gsl_blas_sswap""" 3
.el .IP \f(CWgsl_blas_sswap\fR 3
.IX Item "gsl_blas_sswap"
.ie n .IP """gsl_blas_scopy""" 3
.el .IP \f(CWgsl_blas_scopy\fR 3
.IX Item "gsl_blas_scopy"
.ie n .IP """gsl_blas_saxpy""" 3
.el .IP \f(CWgsl_blas_saxpy\fR 3
.IX Item "gsl_blas_saxpy"
.ie n .IP """gsl_blas_dswap($x, $y)""" 3
.el .IP "\f(CWgsl_blas_dswap($x, $y)\fR" 3
.IX Item "gsl_blas_dswap($x, $y)"
.PD
This function exchanges the elements of the vectors \f(CW$x\fR and \f(CW$y\fR. The function returns 0 if the operation succeeded, 1 otherwise.
.ie n .IP """gsl_blas_dcopy($x, $y)""" 3
.el .IP "\f(CWgsl_blas_dcopy($x, $y)\fR" 3
.IX Item "gsl_blas_dcopy($x, $y)"
This function copies the elements of the vector \f(CW$x\fR into the vector \f(CW$y\fR. The function returns 0 if the operation succeeded, 1 otherwise.
.ie n .IP """gsl_blas_daxpy($alpha, $x, $y)""" 3
.el .IP "\f(CWgsl_blas_daxpy($alpha, $x, $y)\fR" 3
.IX Item "gsl_blas_daxpy($alpha, $x, $y)"
These functions compute the sum \f(CW$y\fR = \f(CW$alpha\fR * \f(CW$x\fR + \f(CW$y\fR for the vectors \f(CW$x\fR and \f(CW$y\fR.
.ie n .IP """gsl_blas_cswap""" 3
.el .IP \f(CWgsl_blas_cswap\fR 3
.IX Item "gsl_blas_cswap"
.PD 0
.ie n .IP """gsl_blas_ccopy """ 3
.el .IP "\f(CWgsl_blas_ccopy \fR" 3
.IX Item "gsl_blas_ccopy "
.ie n .IP """gsl_blas_caxpy""" 3
.el .IP \f(CWgsl_blas_caxpy\fR 3
.IX Item "gsl_blas_caxpy"
.ie n .IP """gsl_blas_zswap""" 3
.el .IP \f(CWgsl_blas_zswap\fR 3
.IX Item "gsl_blas_zswap"
.ie n .IP """gsl_blas_zcopy""" 3
.el .IP \f(CWgsl_blas_zcopy\fR 3
.IX Item "gsl_blas_zcopy"
.ie n .IP """gsl_blas_zaxpy """ 3
.el .IP "\f(CWgsl_blas_zaxpy \fR" 3
.IX Item "gsl_blas_zaxpy "
.ie n .IP """gsl_blas_srotg""" 3
.el .IP \f(CWgsl_blas_srotg\fR 3
.IX Item "gsl_blas_srotg"
.ie n .IP """gsl_blas_srotmg""" 3
.el .IP \f(CWgsl_blas_srotmg\fR 3
.IX Item "gsl_blas_srotmg"
.ie n .IP """gsl_blas_srot""" 3
.el .IP \f(CWgsl_blas_srot\fR 3
.IX Item "gsl_blas_srot"
.ie n .IP """gsl_blas_srotm """ 3
.el .IP "\f(CWgsl_blas_srotm \fR" 3
.IX Item "gsl_blas_srotm "
.ie n .IP """gsl_blas_drotg""" 3
.el .IP \f(CWgsl_blas_drotg\fR 3
.IX Item "gsl_blas_drotg"
.ie n .IP """gsl_blas_drotmg""" 3
.el .IP \f(CWgsl_blas_drotmg\fR 3
.IX Item "gsl_blas_drotmg"
.ie n .IP """gsl_blas_drot($x, $y, $c, $s)""" 3
.el .IP "\f(CWgsl_blas_drot($x, $y, $c, $s)\fR" 3
.IX Item "gsl_blas_drot($x, $y, $c, $s)"
.PD
This function applies a Givens rotation (x', y') = (c x + s y, \-s x + c y) to the vectors \f(CW$x\fR, \f(CW$y\fR.
.ie n .IP """gsl_blas_drotm """ 3
.el .IP "\f(CWgsl_blas_drotm \fR" 3
.IX Item "gsl_blas_drotm "
.PD 0
.ie n .IP """gsl_blas_sscal""" 3
.el .IP \f(CWgsl_blas_sscal\fR 3
.IX Item "gsl_blas_sscal"
.ie n .IP """gsl_blas_dscal($alpha, $x)""" 3
.el .IP "\f(CWgsl_blas_dscal($alpha, $x)\fR" 3
.IX Item "gsl_blas_dscal($alpha, $x)"
.PD
This function rescales the vector \f(CW$x\fR by the multiplicative factor \f(CW$alpha\fR.
.ie n .IP """gsl_blas_cscal""" 3
.el .IP \f(CWgsl_blas_cscal\fR 3
.IX Item "gsl_blas_cscal"
.PD 0
.ie n .IP """gsl_blas_zscal """ 3
.el .IP "\f(CWgsl_blas_zscal \fR" 3
.IX Item "gsl_blas_zscal "
.ie n .IP """gsl_blas_csscal""" 3
.el .IP \f(CWgsl_blas_csscal\fR 3
.IX Item "gsl_blas_csscal"
.ie n .IP """gsl_blas_zdscal""" 3
.el .IP \f(CWgsl_blas_zdscal\fR 3
.IX Item "gsl_blas_zdscal"
.PD
.SS "Level 2 \- Matrix-vector operations"
.IX Subsection "Level 2 - Matrix-vector operations"
.ie n .IP """gsl_blas_sgemv""" 3
.el .IP \f(CWgsl_blas_sgemv\fR 3
.IX Item "gsl_blas_sgemv"
.PD 0
.ie n .IP """gsl_blas_strmv """ 3
.el .IP "\f(CWgsl_blas_strmv \fR" 3
.IX Item "gsl_blas_strmv "
.ie n .IP """gsl_blas_strsv""" 3
.el .IP \f(CWgsl_blas_strsv\fR 3
.IX Item "gsl_blas_strsv"
.ie n .IP """gsl_blas_dgemv($TransA, $alpha, $A, $x, $beta, $y)"" \- This function computes the matrix-vector product and sum y = \ealpha op(A) x + \ebeta y, where op(A) = A, A^T, A^H for $TransA = $CblasNoTrans, $CblasTrans, $CblasConjTrans (constant values coming from the CBLAS module). $A is a matrix and $x and $y are vectors. The function returns 0 if the operation succeeded, 1 otherwise." 3
.el .IP "\f(CWgsl_blas_dgemv($TransA, $alpha, $A, $x, $beta, $y)\fR \- This function computes the matrix-vector product and sum y = \ealpha op(A) x + \ebeta y, where op(A) = A, A^T, A^H for \f(CW$TransA\fR = \f(CW$CblasNoTrans\fR, \f(CW$CblasTrans\fR, \f(CW$CblasConjTrans\fR (constant values coming from the CBLAS module). \f(CW$A\fR is a matrix and \f(CW$x\fR and \f(CW$y\fR are vectors. The function returns 0 if the operation succeeded, 1 otherwise." 3
.IX Item "gsl_blas_dgemv($TransA, $alpha, $A, $x, $beta, $y) - This function computes the matrix-vector product and sum y = alpha op(A) x + beta y, where op(A) = A, A^T, A^H for $TransA = $CblasNoTrans, $CblasTrans, $CblasConjTrans (constant values coming from the CBLAS module). $A is a matrix and $x and $y are vectors. The function returns 0 if the operation succeeded, 1 otherwise."
.ie n .IP """gsl_blas_dtrmv($Uplo, $TransA, $Diag, $A, $x)"" \- This function computes the matrix-vector product x = op(A) x for the triangular matrix $A, where op(A) = A, A^T, A^H for $TransA = $CblasNoTrans, $CblasTrans, $CblasConjTrans (constant values coming from the CBLAS module). When $Uplo is $CblasUpper then the upper triangle of $A is used, and when $Uplo is $CblasLower then the lower triangle of $A is used. If $Diag is $CblasNonUnit then the diagonal of the matrix is used, but if $Diag is $CblasUnit then the diagonal elements of the matrix $A are taken as unity and are not referenced. The function returns 0 if the operation succeeded, 1 otherwise." 3
.el .IP "\f(CWgsl_blas_dtrmv($Uplo, $TransA, $Diag, $A, $x)\fR \- This function computes the matrix-vector product x = op(A) x for the triangular matrix \f(CW$A\fR, where op(A) = A, A^T, A^H for \f(CW$TransA\fR = \f(CW$CblasNoTrans\fR, \f(CW$CblasTrans\fR, \f(CW$CblasConjTrans\fR (constant values coming from the CBLAS module). When \f(CW$Uplo\fR is \f(CW$CblasUpper\fR then the upper triangle of \f(CW$A\fR is used, and when \f(CW$Uplo\fR is \f(CW$CblasLower\fR then the lower triangle of \f(CW$A\fR is used. If \f(CW$Diag\fR is \f(CW$CblasNonUnit\fR then the diagonal of the matrix is used, but if \f(CW$Diag\fR is \f(CW$CblasUnit\fR then the diagonal elements of the matrix \f(CW$A\fR are taken as unity and are not referenced. The function returns 0 if the operation succeeded, 1 otherwise." 3
.IX Item "gsl_blas_dtrmv($Uplo, $TransA, $Diag, $A, $x) - This function computes the matrix-vector product x = op(A) x for the triangular matrix $A, where op(A) = A, A^T, A^H for $TransA = $CblasNoTrans, $CblasTrans, $CblasConjTrans (constant values coming from the CBLAS module). When $Uplo is $CblasUpper then the upper triangle of $A is used, and when $Uplo is $CblasLower then the lower triangle of $A is used. If $Diag is $CblasNonUnit then the diagonal of the matrix is used, but if $Diag is $CblasUnit then the diagonal elements of the matrix $A are taken as unity and are not referenced. The function returns 0 if the operation succeeded, 1 otherwise."
.ie n .IP """gsl_blas_dtrsv($Uplo, $TransA, $Diag, $A, $x)"" \- This function computes inv(op(A)) x for the vector $x, where op(A) = A, A^T, A^H for $TransA = $CblasNoTrans, $CblasTrans, $CblasConjTrans (constant values coming from the CBLAS module). When $Uplo is $CblasUpper then the upper triangle of $A is used, and when $Uplo is $CblasLower then the lower triangle of $A is used. If $Diag is $CblasNonUnit then the diagonal of the matrix is used, but if $Diag is $CblasUnit then the diagonal elements of the matrix $A are taken as unity and are not referenced. The function returns 0 if the operation succeeded, 1 otherwise." 3
.el .IP "\f(CWgsl_blas_dtrsv($Uplo, $TransA, $Diag, $A, $x)\fR \- This function computes inv(op(A)) x for the vector \f(CW$x\fR, where op(A) = A, A^T, A^H for \f(CW$TransA\fR = \f(CW$CblasNoTrans\fR, \f(CW$CblasTrans\fR, \f(CW$CblasConjTrans\fR (constant values coming from the CBLAS module). When \f(CW$Uplo\fR is \f(CW$CblasUpper\fR then the upper triangle of \f(CW$A\fR is used, and when \f(CW$Uplo\fR is \f(CW$CblasLower\fR then the lower triangle of \f(CW$A\fR is used. If \f(CW$Diag\fR is \f(CW$CblasNonUnit\fR then the diagonal of the matrix is used, but if \f(CW$Diag\fR is \f(CW$CblasUnit\fR then the diagonal elements of the matrix \f(CW$A\fR are taken as unity and are not referenced. The function returns 0 if the operation succeeded, 1 otherwise." 3
.IX Item "gsl_blas_dtrsv($Uplo, $TransA, $Diag, $A, $x) - This function computes inv(op(A)) x for the vector $x, where op(A) = A, A^T, A^H for $TransA = $CblasNoTrans, $CblasTrans, $CblasConjTrans (constant values coming from the CBLAS module). When $Uplo is $CblasUpper then the upper triangle of $A is used, and when $Uplo is $CblasLower then the lower triangle of $A is used. If $Diag is $CblasNonUnit then the diagonal of the matrix is used, but if $Diag is $CblasUnit then the diagonal elements of the matrix $A are taken as unity and are not referenced. The function returns 0 if the operation succeeded, 1 otherwise."
.ie n .IP """gsl_blas_cgemv """ 3
.el .IP "\f(CWgsl_blas_cgemv \fR" 3
.IX Item "gsl_blas_cgemv "
.ie n .IP """gsl_blas_ctrmv""" 3
.el .IP \f(CWgsl_blas_ctrmv\fR 3
.IX Item "gsl_blas_ctrmv"
.ie n .IP """gsl_blas_ctrsv""" 3
.el .IP \f(CWgsl_blas_ctrsv\fR 3
.IX Item "gsl_blas_ctrsv"
.ie n .IP """gsl_blas_zgemv """ 3
.el .IP "\f(CWgsl_blas_zgemv \fR" 3
.IX Item "gsl_blas_zgemv "
.ie n .IP """gsl_blas_ztrmv""" 3
.el .IP \f(CWgsl_blas_ztrmv\fR 3
.IX Item "gsl_blas_ztrmv"
.ie n .IP """gsl_blas_ztrsv""" 3
.el .IP \f(CWgsl_blas_ztrsv\fR 3
.IX Item "gsl_blas_ztrsv"
.ie n .IP """gsl_blas_ssymv""" 3
.el .IP \f(CWgsl_blas_ssymv\fR 3
.IX Item "gsl_blas_ssymv"
.ie n .IP """gsl_blas_sger """ 3
.el .IP "\f(CWgsl_blas_sger \fR" 3
.IX Item "gsl_blas_sger "
.ie n .IP """gsl_blas_ssyr""" 3
.el .IP \f(CWgsl_blas_ssyr\fR 3
.IX Item "gsl_blas_ssyr"
.ie n .IP """gsl_blas_ssyr2""" 3
.el .IP \f(CWgsl_blas_ssyr2\fR 3
.IX Item "gsl_blas_ssyr2"
.ie n .IP """gsl_blas_dsymv""" 3
.el .IP \f(CWgsl_blas_dsymv\fR 3
.IX Item "gsl_blas_dsymv"
.ie n .IP """gsl_blas_dger($alpha, $x, $y, $A)"" \- This function computes the rank\-1 update A = alpha x y^T + A of the matrix $A. $x and $y are vectors. The function returns 0 if the operation succeeded, 1 otherwise." 3
.el .IP "\f(CWgsl_blas_dger($alpha, $x, $y, $A)\fR \- This function computes the rank\-1 update A = alpha x y^T + A of the matrix \f(CW$A\fR. \f(CW$x\fR and \f(CW$y\fR are vectors. The function returns 0 if the operation succeeded, 1 otherwise." 3
.IX Item "gsl_blas_dger($alpha, $x, $y, $A) - This function computes the rank-1 update A = alpha x y^T + A of the matrix $A. $x and $y are vectors. The function returns 0 if the operation succeeded, 1 otherwise."
.ie n .IP """gsl_blas_dsyr($Uplo, $alpha, $x, $A)"" \- This function computes the symmetric rank\-1 update A = \ealpha x x^T + A of the symmetric matrix $A and the vector $x. Since the matrix $A is symmetric only its upper half or lower half need to be stored. When $Uplo is $CblasUpper then the upper triangle and diagonal of $A are used, and when $Uplo is $CblasLower then the lower triangle and diagonal of $A are used. The function returns 0 if the operation succeeded, 1 otherwise." 3
.el .IP "\f(CWgsl_blas_dsyr($Uplo, $alpha, $x, $A)\fR \- This function computes the symmetric rank\-1 update A = \ealpha x x^T + A of the symmetric matrix \f(CW$A\fR and the vector \f(CW$x\fR. Since the matrix \f(CW$A\fR is symmetric only its upper half or lower half need to be stored. When \f(CW$Uplo\fR is \f(CW$CblasUpper\fR then the upper triangle and diagonal of \f(CW$A\fR are used, and when \f(CW$Uplo\fR is \f(CW$CblasLower\fR then the lower triangle and diagonal of \f(CW$A\fR are used. The function returns 0 if the operation succeeded, 1 otherwise." 3
.IX Item "gsl_blas_dsyr($Uplo, $alpha, $x, $A) - This function computes the symmetric rank-1 update A = alpha x x^T + A of the symmetric matrix $A and the vector $x. Since the matrix $A is symmetric only its upper half or lower half need to be stored. When $Uplo is $CblasUpper then the upper triangle and diagonal of $A are used, and when $Uplo is $CblasLower then the lower triangle and diagonal of $A are used. The function returns 0 if the operation succeeded, 1 otherwise."
.ie n .IP """gsl_blas_dsyr2($Uplo, $alpha, $x, $y, $A)"" \- This function computes the symmetric rank\-2 update A = \ealpha x y^T + \ealpha y x^T + A of the symmetric matrix $A, the vector $x and vector $y. Since the matrix $A is symmetric only its upper half or lower half need to be stored. When $Uplo is $CblasUpper then the upper triangle and diagonal of $A are used, and when $Uplo is $CblasLower then the lower triangle and diagonal of $A are used." 3
.el .IP "\f(CWgsl_blas_dsyr2($Uplo, $alpha, $x, $y, $A)\fR \- This function computes the symmetric rank\-2 update A = \ealpha x y^T + \ealpha y x^T + A of the symmetric matrix \f(CW$A\fR, the vector \f(CW$x\fR and vector \f(CW$y\fR. Since the matrix \f(CW$A\fR is symmetric only its upper half or lower half need to be stored. When \f(CW$Uplo\fR is \f(CW$CblasUpper\fR then the upper triangle and diagonal of \f(CW$A\fR are used, and when \f(CW$Uplo\fR is \f(CW$CblasLower\fR then the lower triangle and diagonal of \f(CW$A\fR are used." 3
.IX Item "gsl_blas_dsyr2($Uplo, $alpha, $x, $y, $A) - This function computes the symmetric rank-2 update A = alpha x y^T + alpha y x^T + A of the symmetric matrix $A, the vector $x and vector $y. Since the matrix $A is symmetric only its upper half or lower half need to be stored. When $Uplo is $CblasUpper then the upper triangle and diagonal of $A are used, and when $Uplo is $CblasLower then the lower triangle and diagonal of $A are used."
.ie n .IP """gsl_blas_chemv""" 3
.el .IP \f(CWgsl_blas_chemv\fR 3
.IX Item "gsl_blas_chemv"
.ie n .IP """gsl_blas_cgeru """ 3
.el .IP "\f(CWgsl_blas_cgeru \fR" 3
.IX Item "gsl_blas_cgeru "
.ie n .IP """gsl_blas_cgerc""" 3
.el .IP \f(CWgsl_blas_cgerc\fR 3
.IX Item "gsl_blas_cgerc"
.ie n .IP """gsl_blas_cher""" 3
.el .IP \f(CWgsl_blas_cher\fR 3
.IX Item "gsl_blas_cher"
.ie n .IP """gsl_blas_cher2""" 3
.el .IP \f(CWgsl_blas_cher2\fR 3
.IX Item "gsl_blas_cher2"
.ie n .IP """gsl_blas_zhemv """ 3
.el .IP "\f(CWgsl_blas_zhemv \fR" 3
.IX Item "gsl_blas_zhemv "
.ie n .IP """gsl_blas_zgeru($alpha, $x, $y, $A)"" \- This function computes the rank\-1 update A = alpha x y^T + A of the complex matrix $A. $alpha is a complex number and $x and $y are complex vectors. The function returns 0 if the operation succeeded, 1 otherwise." 3
.el .IP "\f(CWgsl_blas_zgeru($alpha, $x, $y, $A)\fR \- This function computes the rank\-1 update A = alpha x y^T + A of the complex matrix \f(CW$A\fR. \f(CW$alpha\fR is a complex number and \f(CW$x\fR and \f(CW$y\fR are complex vectors. The function returns 0 if the operation succeeded, 1 otherwise." 3
.IX Item "gsl_blas_zgeru($alpha, $x, $y, $A) - This function computes the rank-1 update A = alpha x y^T + A of the complex matrix $A. $alpha is a complex number and $x and $y are complex vectors. The function returns 0 if the operation succeeded, 1 otherwise."
.ie n .IP """gsl_blas_zgerc""" 3
.el .IP \f(CWgsl_blas_zgerc\fR 3
.IX Item "gsl_blas_zgerc"
.ie n .IP """gsl_blas_zher($Uplo, $alpha, $x, $A)"" \- This function computes the hermitian rank\-1 update A = \ealpha x x^H + A of the hermitian matrix $A and of the complex vector $x. Since the matrix $A is hermitian only its upper half or lower half need to be stored. When $Uplo is $CblasUpper then the upper triangle and diagonal of $A are used, and when $Uplo is $CblasLower then the lower triangle and diagonal of $A are used. The imaginary elements of the diagonal are automatically set to zero. The function returns 0 if the operation succeeded, 1 otherwise." 3
.el .IP "\f(CWgsl_blas_zher($Uplo, $alpha, $x, $A)\fR \- This function computes the hermitian rank\-1 update A = \ealpha x x^H + A of the hermitian matrix \f(CW$A\fR and of the complex vector \f(CW$x\fR. Since the matrix \f(CW$A\fR is hermitian only its upper half or lower half need to be stored. When \f(CW$Uplo\fR is \f(CW$CblasUpper\fR then the upper triangle and diagonal of \f(CW$A\fR are used, and when \f(CW$Uplo\fR is \f(CW$CblasLower\fR then the lower triangle and diagonal of \f(CW$A\fR are used. The imaginary elements of the diagonal are automatically set to zero. The function returns 0 if the operation succeeded, 1 otherwise." 3
.IX Item "gsl_blas_zher($Uplo, $alpha, $x, $A) - This function computes the hermitian rank-1 update A = alpha x x^H + A of the hermitian matrix $A and of the complex vector $x. Since the matrix $A is hermitian only its upper half or lower half need to be stored. When $Uplo is $CblasUpper then the upper triangle and diagonal of $A are used, and when $Uplo is $CblasLower then the lower triangle and diagonal of $A are used. The imaginary elements of the diagonal are automatically set to zero. The function returns 0 if the operation succeeded, 1 otherwise."
.ie n .IP """gsl_blas_zher2 """ 3
.el .IP "\f(CWgsl_blas_zher2 \fR" 3
.IX Item "gsl_blas_zher2 "
.PD
.SS "Level 3 \- Matrix-matrix operations"
.IX Subsection "Level 3 - Matrix-matrix operations"
.ie n .IP """gsl_blas_sgemm""" 3
.el .IP \f(CWgsl_blas_sgemm\fR 3
.IX Item "gsl_blas_sgemm"
.PD 0
.ie n .IP """gsl_blas_ssymm""" 3
.el .IP \f(CWgsl_blas_ssymm\fR 3
.IX Item "gsl_blas_ssymm"
.ie n .IP """gsl_blas_ssyrk""" 3
.el .IP \f(CWgsl_blas_ssyrk\fR 3
.IX Item "gsl_blas_ssyrk"
.ie n .IP """gsl_blas_ssyr2k """ 3
.el .IP "\f(CWgsl_blas_ssyr2k \fR" 3
.IX Item "gsl_blas_ssyr2k "
.ie n .IP """gsl_blas_strmm""" 3
.el .IP \f(CWgsl_blas_strmm\fR 3
.IX Item "gsl_blas_strmm"
.ie n .IP """gsl_blas_strsm""" 3
.el .IP \f(CWgsl_blas_strsm\fR 3
.IX Item "gsl_blas_strsm"
.ie n .IP """gsl_blas_dgemm($TransA, $TransB, $alpha, $A, $B, $beta, $C)"" \- This function computes the matrix-matrix product and sum C = \ealpha op(A) op(B) + \ebeta C where op(A) = A, A^T, A^H for $TransA = $CblasNoTrans, $CblasTrans, $CblasConjTrans and similarly for the parameter $TransB. The function returns 0 if the operation succeeded, 1 otherwise." 3
.el .IP "\f(CWgsl_blas_dgemm($TransA, $TransB, $alpha, $A, $B, $beta, $C)\fR \- This function computes the matrix-matrix product and sum C = \ealpha op(A) op(B) + \ebeta C where op(A) = A, A^T, A^H for \f(CW$TransA\fR = \f(CW$CblasNoTrans\fR, \f(CW$CblasTrans\fR, \f(CW$CblasConjTrans\fR and similarly for the parameter \f(CW$TransB\fR. The function returns 0 if the operation succeeded, 1 otherwise." 3
.IX Item "gsl_blas_dgemm($TransA, $TransB, $alpha, $A, $B, $beta, $C) - This function computes the matrix-matrix product and sum C = alpha op(A) op(B) + beta C where op(A) = A, A^T, A^H for $TransA = $CblasNoTrans, $CblasTrans, $CblasConjTrans and similarly for the parameter $TransB. The function returns 0 if the operation succeeded, 1 otherwise."
.ie n .IP """gsl_blas_dsymm($Side, $Uplo, $alpha, $A, $B, $beta, $C)"" \- This function computes the matrix-matrix product and sum C = \ealpha A B + \ebeta C for $Side is $CblasLeft and C = \ealpha B A + \ebeta C for $Side is $CblasRight, where the matrix $A is symmetric. When $Uplo is $CblasUpper then the upper triangle and diagonal of $A are used, and when $Uplo is $CblasLower then the lower triangle and diagonal of $A are used. The function returns 0 if the operation succeeded, 1 otherwise." 3
.el .IP "\f(CWgsl_blas_dsymm($Side, $Uplo, $alpha, $A, $B, $beta, $C)\fR \- This function computes the matrix-matrix product and sum C = \ealpha A B + \ebeta C for \f(CW$Side\fR is \f(CW$CblasLeft\fR and C = \ealpha B A + \ebeta C for \f(CW$Side\fR is \f(CW$CblasRight\fR, where the matrix \f(CW$A\fR is symmetric. When \f(CW$Uplo\fR is \f(CW$CblasUpper\fR then the upper triangle and diagonal of \f(CW$A\fR are used, and when \f(CW$Uplo\fR is \f(CW$CblasLower\fR then the lower triangle and diagonal of \f(CW$A\fR are used. The function returns 0 if the operation succeeded, 1 otherwise." 3
.IX Item "gsl_blas_dsymm($Side, $Uplo, $alpha, $A, $B, $beta, $C) - This function computes the matrix-matrix product and sum C = alpha A B + beta C for $Side is $CblasLeft and C = alpha B A + beta C for $Side is $CblasRight, where the matrix $A is symmetric. When $Uplo is $CblasUpper then the upper triangle and diagonal of $A are used, and when $Uplo is $CblasLower then the lower triangle and diagonal of $A are used. The function returns 0 if the operation succeeded, 1 otherwise."
.ie n .IP """gsl_blas_dsyrk($Uplo, $Trans, $alpha, $A, $beta, $C)"" \- This function computes a rank-k update of the symmetric matrix $C, C = \ealpha A A^T + \ebeta C when $Trans is $CblasNoTrans and C = \ealpha A^T A + \ebeta C when $Trans is $CblasTrans. Since the matrix $C is symmetric only its upper half or lower half need to be stored. When $Uplo is $CblasUpper then the upper triangle and diagonal of $C are used, and when $Uplo is $CblasLower then the lower triangle and diagonal of $C are used. The function returns 0 if the operation succeeded, 1 otherwise." 3
.el .IP "\f(CWgsl_blas_dsyrk($Uplo, $Trans, $alpha, $A, $beta, $C)\fR \- This function computes a rank-k update of the symmetric matrix \f(CW$C\fR, C = \ealpha A A^T + \ebeta C when \f(CW$Trans\fR is \f(CW$CblasNoTrans\fR and C = \ealpha A^T A + \ebeta C when \f(CW$Trans\fR is \f(CW$CblasTrans\fR. Since the matrix \f(CW$C\fR is symmetric only its upper half or lower half need to be stored. When \f(CW$Uplo\fR is \f(CW$CblasUpper\fR then the upper triangle and diagonal of \f(CW$C\fR are used, and when \f(CW$Uplo\fR is \f(CW$CblasLower\fR then the lower triangle and diagonal of \f(CW$C\fR are used. The function returns 0 if the operation succeeded, 1 otherwise." 3
.IX Item "gsl_blas_dsyrk($Uplo, $Trans, $alpha, $A, $beta, $C) - This function computes a rank-k update of the symmetric matrix $C, C = alpha A A^T + beta C when $Trans is $CblasNoTrans and C = alpha A^T A + beta C when $Trans is $CblasTrans. Since the matrix $C is symmetric only its upper half or lower half need to be stored. When $Uplo is $CblasUpper then the upper triangle and diagonal of $C are used, and when $Uplo is $CblasLower then the lower triangle and diagonal of $C are used. The function returns 0 if the operation succeeded, 1 otherwise."
.ie n .IP """gsl_blas_dsyr2k($Uplo, $Trans, $alpha, $A, $B, $beta, $C)"" \- This function computes a rank\-2k update of the symmetric matrix $C, C = \ealpha A B^T + \ealpha B A^T + \ebeta C when $Trans is $CblasNoTrans and C = \ealpha A^T B + \ealpha B^T A + \ebeta C when $Trans is $CblasTrans. Since the matrix $C is symmetric only its upper half or lower half need to be stored. When $Uplo is $CblasUpper then the upper triangle and diagonal of $C are used, and when $Uplo is $CblasLower then the lower triangle and diagonal of $C are used. The function returns 0 if the operation succeeded, 1 otherwise." 3
.el .IP "\f(CWgsl_blas_dsyr2k($Uplo, $Trans, $alpha, $A, $B, $beta, $C)\fR \- This function computes a rank\-2k update of the symmetric matrix \f(CW$C\fR, C = \ealpha A B^T + \ealpha B A^T + \ebeta C when \f(CW$Trans\fR is \f(CW$CblasNoTrans\fR and C = \ealpha A^T B + \ealpha B^T A + \ebeta C when \f(CW$Trans\fR is \f(CW$CblasTrans\fR. Since the matrix \f(CW$C\fR is symmetric only its upper half or lower half need to be stored. When \f(CW$Uplo\fR is \f(CW$CblasUpper\fR then the upper triangle and diagonal of \f(CW$C\fR are used, and when \f(CW$Uplo\fR is \f(CW$CblasLower\fR then the lower triangle and diagonal of \f(CW$C\fR are used. The function returns 0 if the operation succeeded, 1 otherwise." 3
.IX Item "gsl_blas_dsyr2k($Uplo, $Trans, $alpha, $A, $B, $beta, $C) - This function computes a rank-2k update of the symmetric matrix $C, C = alpha A B^T + alpha B A^T + beta C when $Trans is $CblasNoTrans and C = alpha A^T B + alpha B^T A + beta C when $Trans is $CblasTrans. Since the matrix $C is symmetric only its upper half or lower half need to be stored. When $Uplo is $CblasUpper then the upper triangle and diagonal of $C are used, and when $Uplo is $CblasLower then the lower triangle and diagonal of $C are used. The function returns 0 if the operation succeeded, 1 otherwise."
.ie n .IP """gsl_blas_dtrmm($Side, $Uplo, $TransA, $Diag, $alpha, $A, $B)"" \- This function computes the matrix-matrix product B = \ealpha op(A) B for $Side is $CblasLeft and B = \ealpha B op(A) for $Side is $CblasRight. The matrix $A is triangular and op(A) = A, A^T, A^H for $TransA = $CblasNoTrans, $CblasTrans, $CblasConjTrans. When $Uplo is $CblasUpper then the upper triangle of $A is used, and when $Uplo is $CblasLower then the lower triangle of $A is used. If $Diag is $CblasNonUnit then the diagonal of $A is used, but if $Diag is $CblasUnit then the diagonal elements of the matrix $A are taken as unity and are not referenced. The function returns 0 if the operation succeeded, 1 otherwise." 3
.el .IP "\f(CWgsl_blas_dtrmm($Side, $Uplo, $TransA, $Diag, $alpha, $A, $B)\fR \- This function computes the matrix-matrix product B = \ealpha op(A) B for \f(CW$Side\fR is \f(CW$CblasLeft\fR and B = \ealpha B op(A) for \f(CW$Side\fR is \f(CW$CblasRight\fR. The matrix \f(CW$A\fR is triangular and op(A) = A, A^T, A^H for \f(CW$TransA\fR = \f(CW$CblasNoTrans\fR, \f(CW$CblasTrans\fR, \f(CW$CblasConjTrans\fR. When \f(CW$Uplo\fR is \f(CW$CblasUpper\fR then the upper triangle of \f(CW$A\fR is used, and when \f(CW$Uplo\fR is \f(CW$CblasLower\fR then the lower triangle of \f(CW$A\fR is used. If \f(CW$Diag\fR is \f(CW$CblasNonUnit\fR then the diagonal of \f(CW$A\fR is used, but if \f(CW$Diag\fR is \f(CW$CblasUnit\fR then the diagonal elements of the matrix \f(CW$A\fR are taken as unity and are not referenced. The function returns 0 if the operation succeeded, 1 otherwise." 3
.IX Item "gsl_blas_dtrmm($Side, $Uplo, $TransA, $Diag, $alpha, $A, $B) - This function computes the matrix-matrix product B = alpha op(A) B for $Side is $CblasLeft and B = alpha B op(A) for $Side is $CblasRight. The matrix $A is triangular and op(A) = A, A^T, A^H for $TransA = $CblasNoTrans, $CblasTrans, $CblasConjTrans. When $Uplo is $CblasUpper then the upper triangle of $A is used, and when $Uplo is $CblasLower then the lower triangle of $A is used. If $Diag is $CblasNonUnit then the diagonal of $A is used, but if $Diag is $CblasUnit then the diagonal elements of the matrix $A are taken as unity and are not referenced. The function returns 0 if the operation succeeded, 1 otherwise."
.ie n .IP """gsl_blas_dtrsm($Side, $Uplo, $TransA, $Diag, $alpha, $A, $B)"" \- This function computes the inverse-matrix matrix product B = \ealpha op(inv(A))B for $Side is $CblasLeft and B = \ealpha B op(inv(A)) for $Side is $CblasRight. The matrix $A is triangular and op(A) = A, A^T, A^H for $TransA = $CblasNoTrans, $CblasTrans, $CblasConjTrans. When $Uplo is $CblasUpper then the upper triangle of $A is used, and when $Uplo is $CblasLower then the lower triangle of $A is used. If $Diag is $CblasNonUnit then the diagonal of $A is used, but if $Diag is $CblasUnit then the diagonal elements of the matrix $A are taken as unity and are not referenced. The function returns 0 if the operation succeeded, 1 otherwise." 3
.el .IP "\f(CWgsl_blas_dtrsm($Side, $Uplo, $TransA, $Diag, $alpha, $A, $B)\fR \- This function computes the inverse-matrix matrix product B = \ealpha op(inv(A))B for \f(CW$Side\fR is \f(CW$CblasLeft\fR and B = \ealpha B op(inv(A)) for \f(CW$Side\fR is \f(CW$CblasRight\fR. The matrix \f(CW$A\fR is triangular and op(A) = A, A^T, A^H for \f(CW$TransA\fR = \f(CW$CblasNoTrans\fR, \f(CW$CblasTrans\fR, \f(CW$CblasConjTrans\fR. When \f(CW$Uplo\fR is \f(CW$CblasUpper\fR then the upper triangle of \f(CW$A\fR is used, and when \f(CW$Uplo\fR is \f(CW$CblasLower\fR then the lower triangle of \f(CW$A\fR is used. If \f(CW$Diag\fR is \f(CW$CblasNonUnit\fR then the diagonal of \f(CW$A\fR is used, but if \f(CW$Diag\fR is \f(CW$CblasUnit\fR then the diagonal elements of the matrix \f(CW$A\fR are taken as unity and are not referenced. The function returns 0 if the operation succeeded, 1 otherwise." 3
.IX Item "gsl_blas_dtrsm($Side, $Uplo, $TransA, $Diag, $alpha, $A, $B) - This function computes the inverse-matrix matrix product B = alpha op(inv(A))B for $Side is $CblasLeft and B = alpha B op(inv(A)) for $Side is $CblasRight. The matrix $A is triangular and op(A) = A, A^T, A^H for $TransA = $CblasNoTrans, $CblasTrans, $CblasConjTrans. When $Uplo is $CblasUpper then the upper triangle of $A is used, and when $Uplo is $CblasLower then the lower triangle of $A is used. If $Diag is $CblasNonUnit then the diagonal of $A is used, but if $Diag is $CblasUnit then the diagonal elements of the matrix $A are taken as unity and are not referenced. The function returns 0 if the operation succeeded, 1 otherwise."
.ie n .IP """gsl_blas_cgemm""" 3
.el .IP \f(CWgsl_blas_cgemm\fR 3
.IX Item "gsl_blas_cgemm"
.ie n .IP """gsl_blas_csymm""" 3
.el .IP \f(CWgsl_blas_csymm\fR 3
.IX Item "gsl_blas_csymm"
.ie n .IP """gsl_blas_csyrk""" 3
.el .IP \f(CWgsl_blas_csyrk\fR 3
.IX Item "gsl_blas_csyrk"
.ie n .IP """gsl_blas_csyr2k """ 3
.el .IP "\f(CWgsl_blas_csyr2k \fR" 3
.IX Item "gsl_blas_csyr2k "
.ie n .IP """gsl_blas_ctrmm""" 3
.el .IP \f(CWgsl_blas_ctrmm\fR 3
.IX Item "gsl_blas_ctrmm"
.ie n .IP """gsl_blas_ctrsm""" 3
.el .IP \f(CWgsl_blas_ctrsm\fR 3
.IX Item "gsl_blas_ctrsm"
.ie n .IP """gsl_blas_zgemm($TransA, $TransB, $alpha, $A, $B, $beta, $C)"" \- This function computes the matrix-matrix product and sum C = \ealpha op(A) op(B) + \ebeta C where op(A) = A, A^T, A^H for $TransA = $CblasNoTrans, $CblasTrans, $CblasConjTrans and similarly for the parameter $TransB. The function returns 0 if the operation succeeded, 1 otherwise. $A, $B and $C are complex matrices" 3
.el .IP "\f(CWgsl_blas_zgemm($TransA, $TransB, $alpha, $A, $B, $beta, $C)\fR \- This function computes the matrix-matrix product and sum C = \ealpha op(A) op(B) + \ebeta C where op(A) = A, A^T, A^H for \f(CW$TransA\fR = \f(CW$CblasNoTrans\fR, \f(CW$CblasTrans\fR, \f(CW$CblasConjTrans\fR and similarly for the parameter \f(CW$TransB\fR. The function returns 0 if the operation succeeded, 1 otherwise. \f(CW$A\fR, \f(CW$B\fR and \f(CW$C\fR are complex matrices" 3
.IX Item "gsl_blas_zgemm($TransA, $TransB, $alpha, $A, $B, $beta, $C) - This function computes the matrix-matrix product and sum C = alpha op(A) op(B) + beta C where op(A) = A, A^T, A^H for $TransA = $CblasNoTrans, $CblasTrans, $CblasConjTrans and similarly for the parameter $TransB. The function returns 0 if the operation succeeded, 1 otherwise. $A, $B and $C are complex matrices"
.ie n .IP """gsl_blas_zsymm($Side, $Uplo, $alpha, $A, $B, $beta, $C)"" \- This function computes the matrix-matrix product and sum C = \ealpha A B + \ebeta C for $Side is $CblasLeft and C = \ealpha B A + \ebeta C for $Side is $CblasRight, where the matrix $A is symmetric. When $Uplo is $CblasUpper then the upper triangle and diagonal of $A are used, and when $Uplo is $CblasLower then the lower triangle and diagonal of $A are used. $A, $B and $C are complex matrices. The function returns 0 if the operation succeeded, 1 otherwise." 3
.el .IP "\f(CWgsl_blas_zsymm($Side, $Uplo, $alpha, $A, $B, $beta, $C)\fR \- This function computes the matrix-matrix product and sum C = \ealpha A B + \ebeta C for \f(CW$Side\fR is \f(CW$CblasLeft\fR and C = \ealpha B A + \ebeta C for \f(CW$Side\fR is \f(CW$CblasRight\fR, where the matrix \f(CW$A\fR is symmetric. When \f(CW$Uplo\fR is \f(CW$CblasUpper\fR then the upper triangle and diagonal of \f(CW$A\fR are used, and when \f(CW$Uplo\fR is \f(CW$CblasLower\fR then the lower triangle and diagonal of \f(CW$A\fR are used. \f(CW$A\fR, \f(CW$B\fR and \f(CW$C\fR are complex matrices. The function returns 0 if the operation succeeded, 1 otherwise." 3
.IX Item "gsl_blas_zsymm($Side, $Uplo, $alpha, $A, $B, $beta, $C) - This function computes the matrix-matrix product and sum C = alpha A B + beta C for $Side is $CblasLeft and C = alpha B A + beta C for $Side is $CblasRight, where the matrix $A is symmetric. When $Uplo is $CblasUpper then the upper triangle and diagonal of $A are used, and when $Uplo is $CblasLower then the lower triangle and diagonal of $A are used. $A, $B and $C are complex matrices. The function returns 0 if the operation succeeded, 1 otherwise."
.ie n .IP """gsl_blas_zsyrk($Uplo, $Trans, $alpha, $A, $beta, $C)"" \- This function computes a rank-k update of the symmetric complex matrix $C, C = \ealpha A A^T + \ebeta C when $Trans is $CblasNoTrans and C = \ealpha A^T A + \ebeta C when $Trans is $CblasTrans. Since the matrix $C is symmetric only its upper half or lower half need to be stored. When $Uplo is $CblasUpper then the upper triangle and diagonal of $C are used, and when $Uplo is $CblasLower then the lower triangle and diagonal of $C are used. The function returns 0 if the operation succeeded, 1 otherwise." 3
.el .IP "\f(CWgsl_blas_zsyrk($Uplo, $Trans, $alpha, $A, $beta, $C)\fR \- This function computes a rank-k update of the symmetric complex matrix \f(CW$C\fR, C = \ealpha A A^T + \ebeta C when \f(CW$Trans\fR is \f(CW$CblasNoTrans\fR and C = \ealpha A^T A + \ebeta C when \f(CW$Trans\fR is \f(CW$CblasTrans\fR. Since the matrix \f(CW$C\fR is symmetric only its upper half or lower half need to be stored. When \f(CW$Uplo\fR is \f(CW$CblasUpper\fR then the upper triangle and diagonal of \f(CW$C\fR are used, and when \f(CW$Uplo\fR is \f(CW$CblasLower\fR then the lower triangle and diagonal of \f(CW$C\fR are used. The function returns 0 if the operation succeeded, 1 otherwise." 3
.IX Item "gsl_blas_zsyrk($Uplo, $Trans, $alpha, $A, $beta, $C) - This function computes a rank-k update of the symmetric complex matrix $C, C = alpha A A^T + beta C when $Trans is $CblasNoTrans and C = alpha A^T A + beta C when $Trans is $CblasTrans. Since the matrix $C is symmetric only its upper half or lower half need to be stored. When $Uplo is $CblasUpper then the upper triangle and diagonal of $C are used, and when $Uplo is $CblasLower then the lower triangle and diagonal of $C are used. The function returns 0 if the operation succeeded, 1 otherwise."
.ie n .IP """gsl_blas_zsyr2k($Uplo, $Trans, $alpha, $A, $B, $beta, $C)"" \- This function computes a rank\-2k update of the symmetric matrix $C, C = \ealpha A B^T + \ealpha B A^T + \ebeta C when $Trans is $CblasNoTrans and C = \ealpha A^T B + \ealpha B^T A + \ebeta C when $Trans is $CblasTrans. Since the matrix $C is symmetric only its upper half or lower half need to be stored. When $Uplo is $CblasUpper then the upper triangle and diagonal of $C are used, and when $Uplo is $CblasLower then the lower triangle and diagonal of $C are used. The function returns 0 if the operation succeeded, 1 otherwise. $A, $B and $C are complex matrices and $beta is a complex number." 3
.el .IP "\f(CWgsl_blas_zsyr2k($Uplo, $Trans, $alpha, $A, $B, $beta, $C)\fR \- This function computes a rank\-2k update of the symmetric matrix \f(CW$C\fR, C = \ealpha A B^T + \ealpha B A^T + \ebeta C when \f(CW$Trans\fR is \f(CW$CblasNoTrans\fR and C = \ealpha A^T B + \ealpha B^T A + \ebeta C when \f(CW$Trans\fR is \f(CW$CblasTrans\fR. Since the matrix \f(CW$C\fR is symmetric only its upper half or lower half need to be stored. When \f(CW$Uplo\fR is \f(CW$CblasUpper\fR then the upper triangle and diagonal of \f(CW$C\fR are used, and when \f(CW$Uplo\fR is \f(CW$CblasLower\fR then the lower triangle and diagonal of \f(CW$C\fR are used. The function returns 0 if the operation succeeded, 1 otherwise. \f(CW$A\fR, \f(CW$B\fR and \f(CW$C\fR are complex matrices and \f(CW$beta\fR is a complex number." 3
.IX Item "gsl_blas_zsyr2k($Uplo, $Trans, $alpha, $A, $B, $beta, $C) - This function computes a rank-2k update of the symmetric matrix $C, C = alpha A B^T + alpha B A^T + beta C when $Trans is $CblasNoTrans and C = alpha A^T B + alpha B^T A + beta C when $Trans is $CblasTrans. Since the matrix $C is symmetric only its upper half or lower half need to be stored. When $Uplo is $CblasUpper then the upper triangle and diagonal of $C are used, and when $Uplo is $CblasLower then the lower triangle and diagonal of $C are used. The function returns 0 if the operation succeeded, 1 otherwise. $A, $B and $C are complex matrices and $beta is a complex number."
.ie n .IP """gsl_blas_ztrmm($Side, $Uplo, $TransA, $Diag, $alpha, $A, $B)"" \- This function computes the matrix-matrix product B = \ealpha op(A) B for $Side is $CblasLeft and B = \ealpha B op(A) for $Side is $CblasRight. The matrix $A is triangular and op(A) = A, A^T, A^H for $TransA = $CblasNoTrans, $CblasTrans, $CblasConjTrans. When $Uplo is $CblasUpper then the upper triangle of $A is used, and when $Uplo is $CblasLower then the lower triangle of $A is used. If $Diag is $CblasNonUnit then the diagonal of $A is used, but if $Diag is $CblasUnit then the diagonal elements of the matrix $A are taken as unity and are not referenced. The function returns 0 if the operation succeeded, 1 otherwise. $A and $B are complex matrices and $alpha is a complex number." 3
.el .IP "\f(CWgsl_blas_ztrmm($Side, $Uplo, $TransA, $Diag, $alpha, $A, $B)\fR \- This function computes the matrix-matrix product B = \ealpha op(A) B for \f(CW$Side\fR is \f(CW$CblasLeft\fR and B = \ealpha B op(A) for \f(CW$Side\fR is \f(CW$CblasRight\fR. The matrix \f(CW$A\fR is triangular and op(A) = A, A^T, A^H for \f(CW$TransA\fR = \f(CW$CblasNoTrans\fR, \f(CW$CblasTrans\fR, \f(CW$CblasConjTrans\fR. When \f(CW$Uplo\fR is \f(CW$CblasUpper\fR then the upper triangle of \f(CW$A\fR is used, and when \f(CW$Uplo\fR is \f(CW$CblasLower\fR then the lower triangle of \f(CW$A\fR is used. If \f(CW$Diag\fR is \f(CW$CblasNonUnit\fR then the diagonal of \f(CW$A\fR is used, but if \f(CW$Diag\fR is \f(CW$CblasUnit\fR then the diagonal elements of the matrix \f(CW$A\fR are taken as unity and are not referenced. The function returns 0 if the operation succeeded, 1 otherwise. \f(CW$A\fR and \f(CW$B\fR are complex matrices and \f(CW$alpha\fR is a complex number." 3
.IX Item "gsl_blas_ztrmm($Side, $Uplo, $TransA, $Diag, $alpha, $A, $B) - This function computes the matrix-matrix product B = alpha op(A) B for $Side is $CblasLeft and B = alpha B op(A) for $Side is $CblasRight. The matrix $A is triangular and op(A) = A, A^T, A^H for $TransA = $CblasNoTrans, $CblasTrans, $CblasConjTrans. When $Uplo is $CblasUpper then the upper triangle of $A is used, and when $Uplo is $CblasLower then the lower triangle of $A is used. If $Diag is $CblasNonUnit then the diagonal of $A is used, but if $Diag is $CblasUnit then the diagonal elements of the matrix $A are taken as unity and are not referenced. The function returns 0 if the operation succeeded, 1 otherwise. $A and $B are complex matrices and $alpha is a complex number."
.ie n .IP """gsl_blas_ztrsm($Side, $Uplo, $TransA, $Diag, $alpha, $A, $B)"" \- This function computes the inverse-matrix matrix product B = \ealpha op(inv(A))B for $Side is $CblasLeft and B = \ealpha B op(inv(A)) for $Side is $CblasRight. The matrix $A is triangular and op(A) = A, A^T, A^H for $TransA = $CblasNoTrans, $CblasTrans, $CblasConjTrans. When $Uplo is $CblasUpper then the upper triangle of $A is used, and when $Uplo is $CblasLower then the lower triangle of $A is used. If $Diag is $CblasNonUnit then the diagonal of $A is used, but if $Diag is $CblasUnit then the diagonal elements of the matrix $A are taken as unity and are not referenced. The function returns 0 if the operation succeeded, 1 otherwise. $A and $B are complex matrices and $alpha is a complex number." 3
.el .IP "\f(CWgsl_blas_ztrsm($Side, $Uplo, $TransA, $Diag, $alpha, $A, $B)\fR \- This function computes the inverse-matrix matrix product B = \ealpha op(inv(A))B for \f(CW$Side\fR is \f(CW$CblasLeft\fR and B = \ealpha B op(inv(A)) for \f(CW$Side\fR is \f(CW$CblasRight\fR. The matrix \f(CW$A\fR is triangular and op(A) = A, A^T, A^H for \f(CW$TransA\fR = \f(CW$CblasNoTrans\fR, \f(CW$CblasTrans\fR, \f(CW$CblasConjTrans\fR. When \f(CW$Uplo\fR is \f(CW$CblasUpper\fR then the upper triangle of \f(CW$A\fR is used, and when \f(CW$Uplo\fR is \f(CW$CblasLower\fR then the lower triangle of \f(CW$A\fR is used. If \f(CW$Diag\fR is \f(CW$CblasNonUnit\fR then the diagonal of \f(CW$A\fR is used, but if \f(CW$Diag\fR is \f(CW$CblasUnit\fR then the diagonal elements of the matrix \f(CW$A\fR are taken as unity and are not referenced. The function returns 0 if the operation succeeded, 1 otherwise. \f(CW$A\fR and \f(CW$B\fR are complex matrices and \f(CW$alpha\fR is a complex number." 3
.IX Item "gsl_blas_ztrsm($Side, $Uplo, $TransA, $Diag, $alpha, $A, $B) - This function computes the inverse-matrix matrix product B = alpha op(inv(A))B for $Side is $CblasLeft and B = alpha B op(inv(A)) for $Side is $CblasRight. The matrix $A is triangular and op(A) = A, A^T, A^H for $TransA = $CblasNoTrans, $CblasTrans, $CblasConjTrans. When $Uplo is $CblasUpper then the upper triangle of $A is used, and when $Uplo is $CblasLower then the lower triangle of $A is used. If $Diag is $CblasNonUnit then the diagonal of $A is used, but if $Diag is $CblasUnit then the diagonal elements of the matrix $A are taken as unity and are not referenced. The function returns 0 if the operation succeeded, 1 otherwise. $A and $B are complex matrices and $alpha is a complex number."
.ie n .IP """gsl_blas_chemm""" 3
.el .IP \f(CWgsl_blas_chemm\fR 3
.IX Item "gsl_blas_chemm"
.ie n .IP """gsl_blas_cherk""" 3
.el .IP \f(CWgsl_blas_cherk\fR 3
.IX Item "gsl_blas_cherk"
.ie n .IP """gsl_blas_cher2k""" 3
.el .IP \f(CWgsl_blas_cher2k\fR 3
.IX Item "gsl_blas_cher2k"
.ie n .IP """gsl_blas_zhemm($Side, $Uplo, $alpha, $A, $B, $beta, $C)"" \- This function computes the matrix-matrix product and sum C = \ealpha A B + \ebeta C for $Side is $CblasLeft and C = \ealpha B A + \ebeta C for $Side is $CblasRight, where the matrix $A is hermitian. When Uplo is CblasUpper then the upper triangle and diagonal of A are used, and when Uplo is CblasLower then the lower triangle and diagonal of A are used. The imaginary elements of the diagonal are automatically set to zero." 3
.el .IP "\f(CWgsl_blas_zhemm($Side, $Uplo, $alpha, $A, $B, $beta, $C)\fR \- This function computes the matrix-matrix product and sum C = \ealpha A B + \ebeta C for \f(CW$Side\fR is \f(CW$CblasLeft\fR and C = \ealpha B A + \ebeta C for \f(CW$Side\fR is \f(CW$CblasRight\fR, where the matrix \f(CW$A\fR is hermitian. When Uplo is CblasUpper then the upper triangle and diagonal of A are used, and when Uplo is CblasLower then the lower triangle and diagonal of A are used. The imaginary elements of the diagonal are automatically set to zero." 3
.IX Item "gsl_blas_zhemm($Side, $Uplo, $alpha, $A, $B, $beta, $C) - This function computes the matrix-matrix product and sum C = alpha A B + beta C for $Side is $CblasLeft and C = alpha B A + beta C for $Side is $CblasRight, where the matrix $A is hermitian. When Uplo is CblasUpper then the upper triangle and diagonal of A are used, and when Uplo is CblasLower then the lower triangle and diagonal of A are used. The imaginary elements of the diagonal are automatically set to zero."
.ie n .IP """gsl_blas_zherk($Uplo, $Trans, $alpha, $A, $beta, $C)"" \- This function computes a rank-k update of the hermitian matrix $C, C = \ealpha A A^H + \ebeta C when $Trans is $CblasNoTrans and C = \ealpha A^H A + \ebeta C when $Trans is $CblasTrans. Since the matrix $C is hermitian only its upper half or lower half need to be stored. When $Uplo is $CblasUpper then the upper triangle and diagonal of $C are used, and when $Uplo is $CblasLower then the lower triangle and diagonal of $C are used. The imaginary elements of the diagonal are automatically set to zero. The function returns 0 if the operation succeeded, 1 otherwise. $A, $B and $C are complex matrices and $alpha and $beta are complex numbers." 3
.el .IP "\f(CWgsl_blas_zherk($Uplo, $Trans, $alpha, $A, $beta, $C)\fR \- This function computes a rank-k update of the hermitian matrix \f(CW$C\fR, C = \ealpha A A^H + \ebeta C when \f(CW$Trans\fR is \f(CW$CblasNoTrans\fR and C = \ealpha A^H A + \ebeta C when \f(CW$Trans\fR is \f(CW$CblasTrans\fR. Since the matrix \f(CW$C\fR is hermitian only its upper half or lower half need to be stored. When \f(CW$Uplo\fR is \f(CW$CblasUpper\fR then the upper triangle and diagonal of \f(CW$C\fR are used, and when \f(CW$Uplo\fR is \f(CW$CblasLower\fR then the lower triangle and diagonal of \f(CW$C\fR are used. The imaginary elements of the diagonal are automatically set to zero. The function returns 0 if the operation succeeded, 1 otherwise. \f(CW$A\fR, \f(CW$B\fR and \f(CW$C\fR are complex matrices and \f(CW$alpha\fR and \f(CW$beta\fR are complex numbers." 3
.IX Item "gsl_blas_zherk($Uplo, $Trans, $alpha, $A, $beta, $C) - This function computes a rank-k update of the hermitian matrix $C, C = alpha A A^H + beta C when $Trans is $CblasNoTrans and C = alpha A^H A + beta C when $Trans is $CblasTrans. Since the matrix $C is hermitian only its upper half or lower half need to be stored. When $Uplo is $CblasUpper then the upper triangle and diagonal of $C are used, and when $Uplo is $CblasLower then the lower triangle and diagonal of $C are used. The imaginary elements of the diagonal are automatically set to zero. The function returns 0 if the operation succeeded, 1 otherwise. $A, $B and $C are complex matrices and $alpha and $beta are complex numbers."
.ie n .IP """gsl_blas_zher2k($Uplo, $Trans, $alpha, $A, $B, $beta, $C)"" \- This function computes a rank\-2k update of the hermitian matrix $C, C = \ealpha A B^H + \ealpha^* B A^H + \ebeta C when $Trans is $CblasNoTrans and C = \ealpha A^H B + \ealpha^* B^H A + \ebeta C when $Trans is $CblasConjTrans. Since the matrix $C is hermitian only its upper half or lower half need to be stored. When $Uplo is $CblasUpper then the upper triangle and diagonal of $C are used, and when $Uplo is $CblasLower then the lower triangle and diagonal of $C are used. The imaginary elements of the diagonal are automatically set to zero. The function returns 0 if the operation succeeded, 1 otherwise." 3
.el .IP "\f(CWgsl_blas_zher2k($Uplo, $Trans, $alpha, $A, $B, $beta, $C)\fR \- This function computes a rank\-2k update of the hermitian matrix \f(CW$C\fR, C = \ealpha A B^H + \ealpha^* B A^H + \ebeta C when \f(CW$Trans\fR is \f(CW$CblasNoTrans\fR and C = \ealpha A^H B + \ealpha^* B^H A + \ebeta C when \f(CW$Trans\fR is \f(CW$CblasConjTrans\fR. Since the matrix \f(CW$C\fR is hermitian only its upper half or lower half need to be stored. When \f(CW$Uplo\fR is \f(CW$CblasUpper\fR then the upper triangle and diagonal of \f(CW$C\fR are used, and when \f(CW$Uplo\fR is \f(CW$CblasLower\fR then the lower triangle and diagonal of \f(CW$C\fR are used. The imaginary elements of the diagonal are automatically set to zero. The function returns 0 if the operation succeeded, 1 otherwise." 3
.IX Item "gsl_blas_zher2k($Uplo, $Trans, $alpha, $A, $B, $beta, $C) - This function computes a rank-2k update of the hermitian matrix $C, C = alpha A B^H + alpha^* B A^H + beta C when $Trans is $CblasNoTrans and C = alpha A^H B + alpha^* B^H A + beta C when $Trans is $CblasConjTrans. Since the matrix $C is hermitian only its upper half or lower half need to be stored. When $Uplo is $CblasUpper then the upper triangle and diagonal of $C are used, and when $Uplo is $CblasLower then the lower triangle and diagonal of $C are used. The imaginary elements of the diagonal are automatically set to zero. The function returns 0 if the operation succeeded, 1 otherwise."
.PD
.PP
You have to add the functions you want to use inside the qw /put_function_here /.
You can also write use Math::GSL::BLAS qw/:all/ to use all available functions of the module.
Other tags are also available, here is a complete list of all tags for this module :
.ie n .IP """level1""" 3
.el .IP \f(CWlevel1\fR 3
.IX Item "level1"
.PD 0
.ie n .IP """level2""" 3
.el .IP \f(CWlevel2\fR 3
.IX Item "level2"
.ie n .IP """level3""" 3
.el .IP \f(CWlevel3\fR 3
.IX Item "level3"
.PD
.PP
For more information on the functions, we refer you to the GSL official documentation:
.SH AUTHORS
.IX Header "AUTHORS"
Jonathan "Duke" Leto and Thierry Moisan
.SH "COPYRIGHT AND LICENSE"
.IX Header "COPYRIGHT AND LICENSE"
Copyright (C) 2008\-2023 Jonathan "Duke" Leto and Thierry Moisan
.PP
This program is free software; you can redistribute it and/or modify it
under the same terms as Perl itself.