Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Release v5.0.1, fixing complex retstyle infinite loop bug #62

Merged
merged 2 commits into from
Feb 2, 2022
Merged
Show file tree
Hide file tree
Changes from all commits
Commits
File filter

Filter by extension

Filter by extension

Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
2 changes: 2 additions & 0 deletions README.md
Original file line number Diff line number Diff line change
Expand Up @@ -60,6 +60,8 @@ You can always tell if your system is limited in this fashion by calling `lbt_ge

### Version History

v5.0.1 - Fix complex return wrapper infinite loop bug.

v5.0.0 - Add complex return value wrappers and CBLAS workaround. The complex return value wrapper ensures that all symbols maintain a standard ABI for returning complex numbers, and the CBLAS workaround maps CBLAS symbols to FORTRAN symbols when properly-suffixed CBLAS symbols do not exist, as is the case in MKL `v2022.0`.

v4.1.0 - Add `LBT_STRICT` environment variable that causes calling missing symbols to kill the process.
Expand Down
2 changes: 1 addition & 1 deletion src/Make.inc
Original file line number Diff line number Diff line change
Expand Up @@ -24,7 +24,7 @@ endif

LBT_SOVERSION_MAJOR := 5
LBT_SOVERSION_MINOR := 0
LBT_SOVERSION_PATCH := 0
LBT_SOVERSION_PATCH := 1

ifeq ($(OS), WINNT)
SHLIB_EXT := dll
Expand Down
60 changes: 30 additions & 30 deletions src/complex_return_style_adapters.c
Original file line number Diff line number Diff line change
Expand Up @@ -39,85 +39,85 @@ double complex cmplxret_zdotc_64_(const int64_t * N,


// zdotu
extern void (*zdotu__addr)(double complex * z,
const int32_t *,
const double complex *, const int32_t *,
const double complex *, const int32_t *);
extern void (*cmplxret_zdotu__addr)(double complex * z,
const int32_t *,
const double complex *, const int32_t *,
const double complex *, const int32_t *);
double complex cmplxret_zdotu_(const int32_t * N,
const double complex *X, const int32_t * incX,
const double complex *Y, const int32_t * incY)
{
double complex c;
zdotu__addr(&c, N, X, incX, Y, incY);
cmplxret_zdotu__addr(&c, N, X, incX, Y, incY);
return c;
}

extern void (*zdotu_64__addr)(double complex * z,
const int64_t *,
const double complex *, const int64_t *,
const double complex *, const int64_t *);
extern void (*cmplxret_zdotu_64__addr)(double complex * z,
const int64_t *,
const double complex *, const int64_t *,
const double complex *, const int64_t *);
double complex cmplxret_zdotu_64_(const int64_t * N,
const double complex *X, const int64_t * incX,
const double complex *Y, const int64_t * incY)
{
double complex c;
zdotu_64__addr(&c, N, X, incX, Y, incY);
cmplxret_zdotu_64__addr(&c, N, X, incX, Y, incY);
return c;
}


// cdotc
extern void (*cdotc__addr)(float complex * z,
const int32_t *,
const float complex *, const int32_t *,
const float complex *, const int32_t *);
extern void (*cmplxret_cdotc__addr)(float complex * z,
const int32_t *,
const float complex *, const int32_t *,
const float complex *, const int32_t *);
float complex cmplxret_cdotc_(const int32_t * N,
const float complex *X, const int32_t * incX,
const float complex *Y, const int32_t * incY)
{
float complex c;
cdotc__addr(&c, N, X, incX, Y, incY);
cmplxret_cdotc__addr(&c, N, X, incX, Y, incY);
return c;
}

extern void cdotc_64__addr(float complex * z,
const int64_t *,
const float complex *, const int64_t *,
const float complex *, const int64_t *);
extern void (*cmplxret_cdotc_64__addr)(float complex * z,
const int64_t *,
const float complex *, const int64_t *,
const float complex *, const int64_t *);
float complex cmplxret_cdotc_64_(const int64_t * N,
const float complex *X, const int64_t * incX,
const float complex *Y, const int64_t * incY)
{
float complex c;
cdotc_64__addr(&c, N, X, incX, Y, incY);
cmplxret_cdotc_64__addr(&c, N, X, incX, Y, incY);
return c;
}


// cdotu
extern void (*cdotu__addr)(float complex * z,
const int32_t *,
const float complex *, const int32_t *,
const float complex *, const int32_t *);
extern void (*cmplxret_cdotu__addr)(float complex * z,
const int32_t *,
const float complex *, const int32_t *,
const float complex *, const int32_t *);
float complex cmplxret_cdotu_(const int32_t * N,
const float complex *X, const int32_t * incX,
const float complex *Y, const int32_t * incY)
{
float complex c;
cdotu__addr(&c, N, X, incX, Y, incY);
cmplxret_cdotu__addr(&c, N, X, incX, Y, incY);
return c;
}

extern void (*cdotu_64__addr)(float complex * z,
const int64_t *,
const float complex *, const int64_t *,
const float complex *, const int64_t *);
extern void (*cmplxret_cdotu_64__addr)(float complex * z,
const int64_t *,
const float complex *, const int64_t *,
const float complex *, const int64_t *);
float complex cmplxret_cdotu_64_(const int64_t * N,
const float complex *X, const int64_t * incX,
const float complex *Y, const int64_t * incY)
{
float complex c;
cdotu_64__addr(&c, N, X, incX, Y, incY);
cmplxret_cdotu_64__addr(&c, N, X, incX, Y, incY);
return c;
}