bdiv_mic.c

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.             *
 ***************************************************************************
*/
#include "bdiv_mic.h"


int dbdivg_mic(const int m, const int n, const int k, const double *A, double *W, double *H, const double expo, const int uType, const int nIter)
{
  int error=0;
  
  #ifdef With_MICOffload
    #pragma offload target(mic:0) \
    in(m,n,k,expo,nIter) \
    in(A:length(m*n)) \
    inout(W:length(m*k)) \
    inout(H:length(k*n)) \
    inout(error)
    {
  #else     
    mkl_mic_enable();
  #endif

    error =  dbdivg_cpu(m, n, k, A, W, H, expo, uType, nIter);
   
  #ifdef With_MICOffload
    }
  #endif

  return error;
}


int sbdivg_mic(const int m, const int n, const int k, const float *A, float *W, float *H, const float expo, const int uType, const int nIter)
{
  int error=0;
  
  #ifdef With_MICOffload
    #pragma offload target(mic:0) \
    in(m,n,k,expo,nIter) \
    in(A:length(m*n)) \
    inout(W:length(m*k)) \
    inout(H:length(k*n)) \
    inout(error)
    {
  #else     
    mkl_mic_enable();
  #endif

    error =  sbdivg_cpu(m, n, k, A, W, H, expo, uType, nIter);
    
  #ifdef With_MICOffload
    }
  #endif 

  return error;
}





int dbdivone_mic(const int m, const int n, const int k, const double *A, double *W, double *H, const int uType, const int nIter)
{
  int error=0;
  
  #ifdef With_MICOffload
    #pragma offload target(mic:0) \
    in(m,n,k,nIter) \
    in(A:length(m*n)) \
    inout(W:length(m*k)) \
    inout(H:length(k*n)) \
    inout(error)
    {
  #else     
    mkl_mic_enable();
  #endif

    error =  dbdivone_cpu(m, n, k, A, W, H, uType, nIter);
 
  #ifdef With_MICOffload
    }
  #endif

  return error;
}


int sbdivone_mic(const int m, const int n, const int k, const float *A, float *W, float *H, const int uType, const int nIter)
{
  int error=0;
  
  #ifdef With_MICOffload
    #pragma offload target(mic:0) \
    in(m,n,k,nIter) \
    in(A:length(m*n)) \
    inout(W:length(m*k)) \
    inout(H:length(k*n)) \
    inout(error)
    {
  #else     
    mkl_mic_enable();
  #endif

    error =  sbdivone_cpu(m, n, k, A, W, H, uType, nIter);
 
  #ifdef With_MICOffload
    }
  #endif

  return error;
}


int dbdiv_mic(const int m, const int n, const int k, const double *A, double *W, double *H, const double beta, const int uType, const int nIter)
{
  if ((beta < 0.0) || (nIter <= 0))
    return -1;

  if(beta>=2.0 && beta<=2.0)
    return dmlsa_mic(m, n, k, A, W, H, uType, nIter);
  else
  {
    if(beta>=1.0 && beta<=1.0)
      return dbdivone_mic(m, n, k, A, W, H, uType, nIter);
    else
      return dbdivg_mic(m, n, k, A, W, H, beta, uType, nIter);
  }
}


int sbdiv_mic(const int m, const int n, const int k, const float *A, float *W, float *H, const float beta, const int uType, const int nIter)
{
  if ((beta < 0.0) || (nIter <= 0))
    return -1;

  if(beta>=2.0 && beta<=2.0)
    return smlsa_mic(m, n, k, A, W, H, uType, nIter);
  else
  {
    if(beta>=1.0 && beta<=1.0)
      return sbdivone_mic(m, n, k, A, W, H, uType, nIter);
    else
      return sbdivg_mic(m, n, k, A, W, H, beta, uType, nIter);
  }
}