utils_x86.h

Go to the documentation of this file.

/***************************************************************************
 *   Copyright (C) 2014 by PIR (University of Oviedo) and                  *
 *   INCO2 (Polytechnic University of Valencia) groups.                    *
 *   nnmfpack@gmail.com                                                    *
 *                                                                         *
 *   This program is free software; you can redistribute it and/or modify  *
 *   it under the terms of the GNU General Public License as published by  *
 *   the Free Software Foundation; either version 2 of the License, or     *
 *   (at your option) any later version.                                   *
 *                                                                         *
 *   This program is distributed in the hope that it will be useful,       *
 *   but WITHOUT ANY WARRANTY; without even the implied warranty of        *
 *   MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the         *
 *   GNU General Public License for more details.                          *
 *                                                                         *
 *   You should have received a copy of the GNU General Public License     *
 *   along with this program; if not, write to the                         *
 *   Free Software Foundation, Inc.,                                       *
 *   59 Temple Place - Suite 330, Boston, MA  02111-1307, USA.             *
 ***************************************************************************
*/
#ifndef UTILSX86_H
#define UTILSX86_H

#include <stdio.h>
#include <stdlib.h>
#include <math.h>
#include <assert.h>

#ifdef With_OMP
  #include <omp.h>
#endif

#ifdef With_MKL
  #include <mkl.h>
#else
  #include <cblas.h>
  #ifndef With_ARM
    #include <mm_malloc.h>
  #endif
  extern void dlarnv_(int *, int *, int *, double *);
  extern void slarnv_(int *, int *, int *, float  *);
#endif

#ifdef With_MIC
  #define WRDLEN 64
#else
  #define WRDLEN 32
#endif

#define max(a,b) (((a)>(b ))?( a):(b))
#define min(a,b) (((a)<(b ))?( a):(b))

void dmemset_x86(const int n, double *__restrict__ x, const double val);
void smemset_x86(const int n, float  *__restrict__ x, const float  val);

void ddiv_x86(const int n, const double *x, double *__restrict__ y);
void sdiv_x86(const int n, const float  *x, float  *__restrict__ y);

void dsub_x86(const int n, const double *x, double *__restrict__ y);
void ssub_x86(const int n, const float  *x, float  *__restrict__ y);

void dlarngenn_x86(const int m, const int n, const int seed, double *X);
void slarngenn_x86(const int m, const int n, const int seed, float  *X);

double derror_x86(const int m, const int n, const int k, const double *A, const double *W, const double *H);
float  serror_x86(const int m, const int n, const int k, const float  *A, const float  *W, const float  *H);


#endif