fbdc2c7f72d351b95ce041f93ace0f1fb367fc88
[scilab.git] / scilab / modules / elementary_functions / sci_gateway / cpp / sci_sign.cpp
1 /*
2  * Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
3  * Copyright (C) 2014 - Scilab Enterprises - Anais AUBERT
4  *
5  * Copyright (C) 2012 - 2016 - Scilab Enterprises
6  *
7  * This file is hereby licensed under the terms of the GNU GPL v2.0,
8  * pursuant to article 5.3.4 of the CeCILL v.2.1.
9  * This file was originally licensed under the terms of the CeCILL v2.1,
10  * and continues to be available under such terms.
11  * For more information, see the COPYING file which you should have received
12  * along with this program.
13  *
14  */
15 /*--------------------------------------------------------------------------*/
16 #include <string.h>
17 #include "elem_func_gw.hxx"
18 #include "function.hxx"
19 #include "double.hxx"
20 #include "overload.hxx"
21
22 extern "C"
23 {
24 #include "basic_functions.h"
25 #include "Scierror.h"
26 }
27
28 /*--------------------------------------------------------------------------*/
29 types::Function::ReturnValue sci_sign(types::typed_list &in, int _iRetCount, types::typed_list &out)
30 {
31     int iRows   = 0;
32     int iCols   = 0;
33     int iType   = 0;
34     int* piAddr = NULL;
35
36     double *pdblReal    = NULL;
37     double *pdblImg     = NULL;
38     double *pdblRealRet = NULL;
39     double *pdblImgRet  = NULL;
40
41     if (in.size() != 1)
42     {
43         Scierror(77, _("%s: Wrong number of input argument(s): %d expected.\n"), "sign", 1);
44         return types::Function::Error;
45     }
46
47     if (_iRetCount > 1)
48     {
49         Scierror(78, _("%s: Wrong number of output argument(s): %d expected.\n"), "sign", 1);
50         return types::Function::Error;
51     }
52
53     if (in[0]->isDouble())
54     {
55         types::Double* pDblIn = in[0]->getAs<types::Double>();
56         int dims = pDblIn->getDims();
57         int *dimsArray = pDblIn->getDimsArray();
58         int size = pDblIn->getSize();
59
60         if (pDblIn->isComplex())
61         {
62             types::Double* pOut = new types::Double(dims, dimsArray, true);
63             double *dOutR   = pOut->getReal();
64             double *dOutImg = pOut->getImg() ;
65             double *dInR    = pDblIn->getReal();
66             double *dInImg  = pDblIn->getImg() ;
67
68             for (int i = 0; i < size; i++)
69             {
70                 double dblTemp = dpythags(dInR[i], dInImg[i]);
71                 if (dblTemp == 0)
72                 {
73                     dOutR[i] = 0;
74                     dOutImg[i] = 0;
75                 }
76                 else
77                 {
78                     dOutR[i] =  dInR[i] / dblTemp;
79                     dOutImg[i] = dInImg[i] / dblTemp;
80                 }
81             }
82
83             out.push_back(pOut);
84         }
85         else
86         {
87             types::Double* pOut = new types::Double(dims, dimsArray);
88             double * dOutR  = pOut->getReal();
89             double * dInR   = pDblIn->getReal();
90
91             for (int i = 0; i < size; i++)
92             {
93                 dOutR[i] = dsignsEx(dInR[i]);
94             }
95
96             out.push_back(pOut);
97         }
98     }
99     else
100     {
101         std::wstring wstFuncName = L"%" + in[0]->getShortTypeStr() + L"_sign";
102         return Overload::call(wstFuncName, in, _iRetCount, out);
103     }
104
105     return types::Function::OK;
106 }
107 /*--------------------------------------------------------------------------*/