scnc
/
phono3py
mirror of https://github.com/phonopy/phono3py.git
1
0
Fork 0
Code Issues Packages Projects Releases Wiki Activity

Define complex for VC

This commit is contained in:
Atsushi Togo 2024-12-23 15:45:51 +09:00
parent 7b8f7a8d84
commit 71dcf84108
2 changed files with 14 additions and 21 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

5
c/lapack_wrapper.c
View File

@ -38,7 +38,7 @@
#define min(a, b) ((a) > (b) ? (b) : (a))
#if (defined(MKL_BLAS) || defined(SCIPY_MKL_H)) && !defined(NO_INCLUDE_LAPACKE)
#if defined(_MSC_VER) || defined(MKL_BLAS) || defined(SCIPY_MKL_H)
lapack_complex_double lapack_make_complex_double(double re, double im) {
lapack_complex_double z;
z.real = re;
@ -47,7 +47,8 @@ lapack_complex_double lapack_make_complex_double(double re, double im) {
}
#endif
#if defined(MKL_BLAS) || defined(SCIPY_MKL_H) || defined(NO_INCLUDE_LAPACKE)
#if defined(_MSC_VER) || defined(MKL_BLAS) || defined(SCIPY_MKL_H) || \
defined(NO_INCLUDE_LAPACKE)
#ifndef LAPACKE_malloc
#define LAPACKE_malloc(size) malloc(size)
#endif

30
c/lapack_wrapper.h
View File

@ -35,21 +35,21 @@
#ifndef __lapack_wrapper_H__
#define __lapack_wrapper_H__
#ifdef NO_INCLUDE_LAPACKE
#if defined(_MSC_VER) || defined(MKL_BLAS) || defined(SCIPY_MKL_H)
#if defined(_MSC_VER)
typedef struct {
double real;
double imag;
} lapack_complex_double;
lapack_complex_double lapack_make_complex_double(double re, double im) {
lapack_complex_double c;
c.real = re;
c.imag = im;
return c;
}
#else
#include <mkl.h>
#define lapack_complex_double MKL_Complex16
#endif
lapack_complex_double lapack_make_complex_double(double re, double im);
#define lapack_complex_double_real(z) ((z).real)
#define lapack_complex_double_imag(z) ((z).imag)
#else
#if defined(NO_INCLUDE_LAPACKE)
#include <complex.h>
#define lapack_complex_double double _Complex
#ifdef CMPLX
@ -59,20 +59,12 @@ lapack_complex_double lapack_make_complex_double(double re, double im) {
#endif
#define lapack_complex_double_real(z) (creal(z))
#define lapack_complex_double_imag(z) (cimag(z))
#endif
#endif
#if defined(MKL_BLAS) || defined(SCIPY_MKL_H)
#include <mkl.h>
#define lapack_complex_double MKL_Complex16
MKL_Complex16 lapack_make_complex_double(double re, double im);
#define lapack_complex_double_real(z) ((z).real)
#define lapack_complex_double_imag(z) ((z).imag)
#endif
#if !defined(MKL_BLAS) && !defined(SCIPY_MKL_H) && !defined(NO_INCLUDE_LAPACKE)
#else
#if !defined(MKL_BLAS) && !defined(SCIPY_MKL_H)
#include <lapacke.h>
#endif
#endif
#endif
lapack_complex_double phonoc_complex_prod(const lapack_complex_double a,
const lapack_complex_double b);