Compare commits

..

2 commits

Author SHA1 Message Date
Martin Reinecke
004b3aea9a fix typo 2022-09-29 14:59:32 +02:00
Martin Reinecke
3d376d9e8c update configure.ac 2022-05-22 09:38:27 +02:00
3 changed files with 4 additions and 30 deletions

View file

@ -1,4 +1,4 @@
AC_INIT([libsharp2], [1.0.0])
AC_INIT([libsharp2],[1.0.0])
AM_INIT_AUTOMAKE([foreign subdir-objects -Wall -Werror])
AM_MAINTAINER_MODE([enable])
@ -20,13 +20,13 @@ dnl
m4_ifdef([AM_SILENT_RULES], [AM_SILENT_RULES([yes])])
AC_PROG_CC_C99
m4_version_prereq(2.70, [AC_PROG_CC], [AC_PROG_CC_C99])
AC_OPENMP
# add math library
LIBS="-lm"
AC_PROG_LIBTOOL
LT_INIT
tmpval=`echo $CFLAGS | grep -c '\-DMULTIARCH'`
AM_CONDITIONAL([HAVE_MULTIARCH], [test $tmpval -gt 0])

View file

@ -94,7 +94,6 @@ typedef union
sxdata_s s;
} sxdata_u;
#if 0
static inline void Tvnormalize (Tv * restrict val, Tv * restrict scale,
double maxval)
{
@ -115,31 +114,6 @@ static inline void Tvnormalize (Tv * restrict val, Tv * restrict scale,
mask = vand_mask(vlt(vabs(*val),vfmin),vne(*val,vzero));
}
}
#else
static inline void Tvnormalize (Tv * restrict val_, Tv * restrict scale_,
double maxval)
{
Tv val=*val_, scale=*scale_;
const Tv vfmin=vload(sharp_fsmall*maxval), vfmax=vload(maxval);
const Tv vfsmall=vload(sharp_fsmall), vfbig=vload(sharp_fbig);
Tm mask = vgt(vabs(val),vfmax);
while (vanyTrue(mask))
{
vmuleq_mask(mask,val,vfsmall);
vaddeq_mask(mask,scale,vone);
mask = vgt(vabs(val),vfmax);
}
mask = vand_mask(vlt(vabs(val),vfmin),vne(val,vzero));
while (vanyTrue(mask))
{
vmuleq_mask(mask,val,vfbig);
vsubeq_mask(mask,scale,vone);
mask = vand_mask(vlt(vabs(val),vfmin),vne(val,vzero));
}
*val_ = val;
*scale_ = scale;
}
#endif
static void mypow(Tv val, int npow, const double * restrict powlimit,
Tv * restrict resd, Tv * restrict ress)

View file

@ -520,7 +520,7 @@ static void do_sht (sharp_geom_info *ginfo, sharp_alm_info *ainfo,
{
#ifdef USE_MPI
sharp_execute_mpi(MPI_COMM_WORLD,SHARP_MAP2ALM,spin,&alm[itrans*ncomp],&map[itrans*ncomp],ginfo,
ainfo,SHARP_DP|SHARP_ADD,&ttm2a,op_&tom2a);
ainfo,SHARP_DP|SHARP_ADD,&ttm2a,&tom2a);
#else
sharp_execute(SHARP_MAP2ALM,spin,&alm[itrans*ncomp],&map[itrans*ncomp],ginfo,ainfo,
SHARP_DP|SHARP_ADD,&ttm2a,&tom2a);