Merge remote-tracking branch 'origin/master' into windows
[scilab.git] / scilab / modules / elementary_functions / sci_gateway / cpp / sci_acos.cpp
1 /*
2  * Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
3  * Copyright (C) 2012 - DIGITEO - Cedric DELAMARRE
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 "elem_func_gw.hxx"
17 #include "function.hxx"
18 #include "double.hxx"
19 #include "overload.hxx"
20
21 extern "C"
22 {
23 #include "Scierror.h"
24 #include "localization.h"
25 #include "cos.h"
26
27     int C2F(wacos)(double*, double*, double*, double*);
28 }
29
30 /*
31 clear a;nb = 2500;a = rand(nb, nb);tic();acos(a);toc
32 clear a;nb = 2500;a = rand(nb, nb) + 0.5;tic();acos(a);toc
33 clear a;nb = 2500;a = rand(nb, nb); a = a + a *%i;tic();acos(a);toc
34 */
35
36 /*--------------------------------------------------------------------------*/
37 types::Function::ReturnValue sci_acos(types::typed_list &in, int _iRetCount, types::typed_list &out)
38 {
39     types::Double* pDblIn   = NULL;
40     types::Double* pDblOut  = NULL;
41
42     if (in.size() != 1)
43     {
44         Scierror(77, _("%s: Wrong number of input argument(s): %d expected.\n"), "acos", 1);
45         return types::Function::Error;
46     }
47
48     if (_iRetCount > 1)
49     {
50         Scierror(78, _("%s: Wrong number of output argument(s): %d expected.\n"), "acos", 1);
51         return types::Function::Error;
52     }
53
54     if (in[0]->isDouble() == false)
55     {
56         std::string stFuncName = "%" + in[0]->getShortTypeStr() + "_acos";
57         return Overload::call(stFuncName, in, _iRetCount, out);
58     }
59
60     pDblIn = in[0]->getAs<types::Double>();
61
62     if (pDblIn->isComplex())
63     {
64         pDblOut = new types::Double(pDblIn->getDims(), pDblIn->getDimsArray(), true);
65         int size = pDblIn->getSize();
66
67         double* pInR = pDblIn->get();
68         double* pInI = pDblIn->getImg();
69         double* pOutR = pDblOut->get();
70         double* pOutI = pDblOut->getImg();
71
72         for (int i = 0 ; i < size ; i++)
73         {
74             C2F(wacos)(pInR + i, pInI + i, pOutR + i, pOutI + i);
75         }
76     }
77     else
78     {
79         bool bOutSide = 0;
80         //check if all variables are between [-1,1]
81         double* pInR = pDblIn->get();
82         int size = pDblIn->getSize();
83         for (int i = 0; i < size; i++)
84         {
85             if (std::abs(pInR[i]) > 1)
86             {
87                 bOutSide = 1;
88                 break;
89             }
90         }
91
92         if (bOutSide) // Values outside [-1,1]
93         {
94             pDblOut = new types::Double(pDblIn->getDims(), pDblIn->getDimsArray(), true);
95             double* pOutR = pDblOut->get();
96             double* pOutI = pDblOut->getImg();
97             double zero = 0;
98             for (int i = 0; i < size; i++)
99             {
100                 C2F(wacos)(pInR + i, &zero, pOutR + i, pOutI + i);
101             }
102         }
103         else //all values are in [-1,1]
104         {
105             pDblOut = new types::Double(pDblIn->getDims(), pDblIn->getDimsArray(), false);
106             double* pOutR = pDblOut->get();
107             for (int i = 0; i < size; i++)
108             {
109                 pOutR[i] = std::acos(pInR[i]);
110             }
111         }
112     }
113
114     out.push_back(pDblOut);
115     return types::Function::OK;
116 }
117 /*--------------------------------------------------------------------------*/