#include <gsl/gsl_test.h>
#include <gsl/gsl_ieee_utils.h>
#include <gsl/gsl_math.h>
#include <gsl/gsl_cblas.h>
#include "tests.h"
void
test_trsm (void) {
const double flteps = 1e-4, dbleps = 1e-6;
{
int order = 101;
int side = 141;
int uplo = 121;
int trans = 111;
int diag = 131;
int M = 2;
int N = 3;
float alpha = -0.3f;
float A[] = { -0.279f, 0.058f, 0.437f, 0.462f };
int lda = 2;
float B[] = { 0.578f, 0.473f, -0.34f, -0.128f, 0.503f, 0.2f };
int ldb = 3;
float B_expected[] = { 0.638784f, 0.440702f, -0.392589f, 0.0831169f, -0.326623f, -0.12987f };
cblas_strsm(order, side, uplo, trans, diag, M, N, alpha, A, lda, B, ldb);
{
int i;
for (i = 0; i < 6; i++) {
gsl_test_rel(B[i], B_expected[i], flteps, "strsm(case 1822)");
}
};
};
{
int order = 101;
int side = 141;
int uplo = 121;
int trans = 111;
int diag = 132;
int M = 2;
int N = 3;
float alpha = -0.3f;
float A[] = { -0.735f, -0.861f, 0.772f, -0.242f };
int lda = 2;
float B[] = { -0.793f, -0.162f, -0.844f, 0.143f, -0.379f, -0.46f };
int ldb = 3;
float B_expected[] = { 0.200963f, 0.146496f, 0.372018f, -0.0429f, 0.1137f, 0.138f };
cblas_strsm(order, side, uplo, trans, diag, M, N, alpha, A, lda, B, ldb);
{
int i;
for (i = 0; i < 6; i++) {
gsl_test_rel(B[i], B_expected[i], flteps, "strsm(case 1823)");
}
};
};
{
int order = 101;
int side = 141;
int uplo = 122;
int trans = 111;
int diag = 131;
int M = 2;
int N = 3;
float alpha = -0.3f;
float A[] = { -0.498f, 0.777f, -0.913f, 0.779f };
int lda = 2;
float B[] = { -0.831f, -0.663f, -0.098f, -0.894f, -0.059f, 0.468f };
int ldb = 3;
float B_expected[] = { -0.500602f, -0.399398f, -0.0590361f, -0.242426f, -0.445379f, -0.249422f };
cblas_strsm(order, side, uplo, trans, diag, M, N, alpha, A, lda, B, ldb);
{
int i;
for (i = 0; i < 6; i++) {
gsl_test_rel(B[i], B_expected[i], flteps, "strsm(case 1824)");
}
};
};
{
int order = 101;
int side = 141;
int uplo = 122;
int trans = 111;
int diag = 132;
int M = 2;
int N = 3;
float alpha = -0.3f;
float A[] = { -0.543f, 0.095f, -0.933f, -0.669f };
int lda = 2;
float B[] = { 0.068f, 0.715f, 0.012f, -0.785f, 0.378f, 0.251f };
int ldb = 3;
float B_expected[] = { -0.0204f, -0.2145f, -0.0036f, 0.216467f, -0.313528f, -0.0786588f };
cblas_strsm(order, side, uplo, trans, diag, M, N, alpha, A, lda, B, ldb);
{
int i;
for (i = 0; i < 6; i++) {
gsl_test_rel(B[i], B_expected[i], flteps, "strsm(case 1825)");
}
};
};
{
int order = 102;
int side = 141;
int uplo = 121;
int trans = 111;
int diag = 131;
int M = 2;
int N = 3;
float alpha = -0.3f;
float A[] = { 0.75f, 0.777f, -0.025f, 0.572f };
int lda = 2;
float B[] = { 0.03f, 0.392f, -0.056f, 0.399f, -0.489f, -0.167f };
int ldb = 2;
float B_expected[] = { -0.0188531f, -0.205594f, 0.0154245f, -0.209266f, 0.19852f, 0.0875874f };
cblas_strsm(order, side, uplo, trans, diag, M, N, alpha, A, lda, B, ldb);
{
int i;
for (i = 0; i < 6; i++) {
gsl_test_rel(B[i], B_expected[i], flteps, "strsm(case 1826)");
}
};
};
{
int order = 102;
int side = 141;
int uplo = 121;
int trans = 111;
int diag = 132;
int M = 2;
int N = 3;
float alpha = -0.3f;
float A[] = { 0.899f, -0.447f, 0.338f, -0.74f };
int lda = 2;
float B[] = { 0.964f, -0.104f, -0.199f, 0.503f, -0.386f, -0.764f };
int ldb = 2;
float B_expected[] = { -0.299746f, 0.0312f, 0.110704f, -0.1509f, 0.0383304f, 0.2292f };
cblas_strsm(order, side, uplo, trans, diag, M, N, alpha, A, lda, B, ldb);
{
int i;
for (i = 0; i < 6; i++) {
gsl_test_rel(B[i], B_expected[i], flteps, "strsm(case 1827)");
}
};
};
{
int order = 102;
int side = 141;
int uplo = 122;
int trans = 111;
int diag = 131;
int M = 2;
int N = 3;
float alpha = -0.3f;
float A[] = { 0.279f, 0.73f, -0.366f, 0.583f };
int lda = 2;
float B[] = { -0.572f, 0.75f, 0.603f, 0.697f, 0.908f, 0.119f };
int ldb = 2;
float B_expected[] = { 0.615054f, -1.15607f, -0.648387f, 0.453212f, -0.976344f, 1.16129f };
cblas_strsm(order, side, uplo, trans, diag, M, N, alpha, A, lda, B, ldb);
{
int i;
for (i = 0; i < 6; i++) {
gsl_test_rel(B[i], B_expected[i], flteps, "strsm(case 1828)");
}
};
};
{
int order = 102;
int side = 141;
int uplo = 122;
int trans = 111;
int diag = 132;
int M = 2;
int N = 3;
float alpha = -0.3f;
float A[] = { 0.581f, -0.911f, 0.438f, 0.731f };
int lda = 2;
float B[] = { 0.519f, 0.831f, 0.822f, 0.182f, 0.571f, -0.357f };
int ldb = 2;
float B_expected[] = { -0.1557f, -0.391143f, -0.2466f, -0.279253f, -0.1713f, -0.0489543f };
cblas_strsm(order, side, uplo, trans, diag, M, N, alpha, A, lda, B, ldb);
{
int i;
for (i = 0; i < 6; i++) {
gsl_test_rel(B[i], B_expected[i], flteps, "strsm(case 1829)");
}
};
};
{
int order = 101;
int side = 142;
int uplo = 121;
int trans = 111;
int diag = 131;
int M = 2;
int N = 3;
float alpha = -0.3f;
float A[] = { 0.065f, 0.099f, 0.48f, 0.746f, -0.739f, 0.695f, 0.197f, 0.621f, 0.063f };
int lda = 3;
float B[] = { 0.01f, -0.612f, 0.756f, -0.225f, 0.546f, 0.432f };
int ldb = 3;
float B_expected[] = { -0.0461538f, -0.254627f, -0.439373f, 1.03846f, 0.360768f, -13.9491f };
cblas_strsm(order, side, uplo, trans, diag, M, N, alpha, A, lda, B, ldb);
{
int i;
for (i = 0; i < 6; i++) {
gsl_test_rel(B[i], B_expected[i], flteps, "strsm(case 1830)");
}
};
};
{
int order = 101;
int side = 142;
int uplo = 121;
int trans = 111;
int diag = 132;
int M = 2;
int N = 3;
float alpha = -0.3f;
float A[] = { -0.86f, -0.653f, 0.87f, -0.037f, 0.788f, 0.015f, 0.028f, -0.804f, -0.357f };
int lda = 3;
float B[] = { -0.546f, 0.892f, -0.085f, -0.541f, -0.207f, 0.765f };
int ldb = 3;
float B_expected[] = { 0.1638f, -0.160639f, -0.114596f, 0.1623f, 0.168082f, -0.373222f };
cblas_strsm(order, side, uplo, trans, diag, M, N, alpha, A, lda, B, ldb);
{
int i;
for (i = 0; i < 6; i++) {
gsl_test_rel(B[i], B_expected[i], flteps, "strsm(case 1831)");
}
};
};
{
int order = 101;
int side = 142;
int uplo = 122;
int trans = 111;
int diag = 131;
int M = 2;
int N = 3;
float alpha = -0.3f;
float A[] = { 0.872f, -0.35f, 0.518f, -0.8f, -0.13f, -0.832f, 0.426f, 0.195f, -0.735f };
int lda = 3;
float B[] = { 0.773f, 0.069f, 0.45f, 0.189f, 0.504f, 0.996f };
int ldb = 3;
float B_expected[] = { 0.0431742f, 0.434741f, 0.183673f, 1.36286f, 1.77287f, 0.406531f };
cblas_strsm(order, side, uplo, trans, diag, M, N, alpha, A, lda, B, ldb);
{
int i;
for (i = 0; i < 6; i++) {
gsl_test_rel(B[i], B_expected[i], flteps, "strsm(case 1832)");
}
};
};
{
int order = 101;
int side = 142;
int uplo = 122;
int trans = 111;
int diag = 132;
int M = 2;
int N = 3;
float alpha = -0.3f;
float A[] = { -0.053f, -0.132f, -0.515f, -0.411f, 0.134f, 0.657f, 0.072f, -0.007f, -0.34f };
int lda = 3;
float B[] = { 0.494f, 0.072f, -0.882f, -0.112f, 0.904f, 0.755f };
int ldb = 3;
float B_expected[] = { -0.175368f, -0.0197478f, 0.2646f, -0.0622068f, -0.272786f, -0.2265f };
cblas_strsm(order, side, uplo, trans, diag, M, N, alpha, A, lda, B, ldb);
{
int i;
for (i = 0; i < 6; i++) {
gsl_test_rel(B[i], B_expected[i], flteps, "strsm(case 1833)");
}
};
};
{
int order = 102;
int side = 142;
int uplo = 121;
int trans = 111;
int diag = 131;
int M = 2;
int N = 3;
float alpha = -0.3f;
float A[] = { -0.154f, -0.54f, 0.146f, -0.106f, -0.478f, 0.938f, -0.731f, 0.25f, -0.4f }
没有合适的资源?快使用搜索试试~ 我知道了~
gsl-1.8.tar.gz
需积分: 0 6 下载量 185 浏览量
2008-07-06
09:37:40
上传
评论
收藏 2.31MB GZ 举报
温馨提示
共1637个文件
c:926个
h:328个
texi:77个
还不太会用,但知道很有用,分享一下,希望能跟大家共同交流学习。
资源推荐
资源详情
资源评论
收起资源包目录
gsl-1.8.tar.gz (1637个子文件)
gsl-histogram.1 1KB
gsl-config.1 1KB
gsl-randist.1 1KB
gsl.3 1KB
configure.ac 10KB
Makefile.am 4KB
Makefile.am 3KB
Makefile.am 2KB
Makefile.am 2KB
Makefile.am 1KB
Makefile.am 1KB
Makefile.am 1KB
Makefile.am 1KB
Makefile.am 1KB
Makefile.am 1KB
Makefile.am 1KB
Makefile.am 1KB
Makefile.am 981B
Makefile.am 953B
Makefile.am 901B
Makefile.am 862B
Makefile.am 805B
Makefile.am 777B
Makefile.am 774B
Makefile.am 759B
Makefile.am 732B
Makefile.am 729B
Makefile.am 718B
Makefile.am 710B
Makefile.am 679B
Makefile.am 670B
Makefile.am 636B
Makefile.am 634B
Makefile.am 555B
Makefile.am 549B
Makefile.am 539B
Makefile.am 518B
Makefile.am 509B
Makefile.am 486B
Makefile.am 479B
Makefile.am 433B
Makefile.am 428B
Makefile.am 387B
Makefile.am 382B
Makefile.am 378B
Makefile.am 373B
Makefile.am 366B
Makefile.am 352B
Makefile.am 346B
Makefile.am 306B
Makefile.am 241B
Makefile.am 193B
AUTHORS 1KB
fftalgorithms.bib 6KB
BUGS 11KB
test_trsm.c 123KB
test_trmm.c 123KB
test_auto.c 121KB
test_sf.c 108KB
test.c 99KB
test.c 82KB
test.c 68KB
gamma.c 59KB
hyperg_1F1.c 58KB
blas.c 56KB
test_tbsv.c 53KB
test_tbmv.c 53KB
test_rotg.c 47KB
hyperg_U.c 44KB
legendre_con.c 42KB
test_tpsv.c 42KB
test_tpmv.c 42KB
test_gemm.c 41KB
test_hyperg.c 40KB
test_bessel.c 40KB
test_trsv.c 39KB
test_trmv.c 39KB
coulomb.c 39KB
test.c 38KB
test_rotm.c 35KB
test_legendre.c 34KB
fermi_dirac.c 34KB
zeta.c 32KB
bessel_olver.c 31KB
bessel_zero.c 28KB
hyperg_2F1.c 27KB
bessel.c 26KB
test_nist.c 26KB
test.c 25KB
airy_der.c 25KB
test_gamma.c 23KB
airy.c 23KB
math.c 23KB
vegas.c 21KB
psi.c 21KB
landau.c 21KB
test_scal.c 20KB
test_symm.c 20KB
legendre_poly.c 20KB
test.c 20KB
共 1637 条
- 1
- 2
- 3
- 4
- 5
- 6
- 17
资源评论
yuyong57118194
- 粉丝: 0
- 资源: 7
上传资源 快速赚钱
- 我的内容管理 展开
- 我的资源 快来上传第一个资源
- 我的收益 登录查看自己的收益
- 我的积分 登录查看自己的积分
- 我的C币 登录后查看C币余额
- 我的收藏
- 我的下载
- 下载帮助
安全验证
文档复制为VIP权益,开通VIP直接复制
信息提交成功