elementary_functions module.
[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  * This file must be used under the terms of the CeCILL.
6  * This source file is licensed as described in the file COPYING, which
7  * you should have received as part of this distribution.  The terms
8  * are also available at
9  * http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
10  *
11  */
12 /*--------------------------------------------------------------------------*/
13 #include "elem_func_gw.hxx"
14 #include "function.hxx"
15 #include "double.hxx"
16 #include "overload.hxx"
17 #include "execvisitor.hxx"
18
19
20 extern "C"
21 {
22 #include "Scierror.h"
23 #include "localization.h"
24 #include "cos.h"
25 }
26
27 /*--------------------------------------------------------------------------*/
28 types::Function::ReturnValue sci_acos(types::typed_list &in, int _iRetCount, types::typed_list &out)
29 {
30     types::Double* pDblIn   = NULL;
31     types::Double* pDblOut  = NULL;
32
33     if (in.size() != 1)
34     {
35         ScierrorW(77, _W("%ls: Wrong number of input argument(s): %d expected.\n"), L"acos", 1);
36         return types::Function::Error;
37     }
38
39     if (_iRetCount > 1)
40     {
41         ScierrorW(78, _W("%ls: Wrong number of output argument(s): %d expected.\n"), L"acos", 1);
42         return types::Function::Error;
43     }
44
45     if (in[0]->isDouble() == false)
46     {
47         std::wstring wstFuncName = L"%"  + in[0]->getShortTypeStr() + L"_acos";
48         return Overload::call(wstFuncName, in, _iRetCount, out, new ExecVisitor());
49     }
50
51     pDblIn = in[0]->getAs<types::Double>();
52
53     if (pDblIn->isComplex())
54     {
55         pDblOut = new types::Double(pDblIn->getDims(), pDblIn->getDimsArray(), true);
56         for (int i = 0 ; i < pDblIn->getSize() ; i++)
57         {
58             wacos(pDblIn->get(i), pDblIn->getImg(i), pDblOut->get() + i, pDblOut->getImg() + i);
59         }
60     }
61     else
62     {
63         bool bOutSide = 0;
64         //check if all variables are between [-1,1]
65         for (int i = 0 ; i < pDblIn->getSize() ; i++)
66         {
67             if (fabs(pDblIn->get(i)) > 1)
68             {
69                 bOutSide = 1;
70                 break;
71             }
72         }
73
74         if (bOutSide) // Values outside [-1,1]
75         {
76             pDblOut = new types::Double(pDblIn->getDims(), pDblIn->getDimsArray(), true);
77             for (int i = 0 ; i < pDblIn->getSize() ; i++)
78             {
79                 wacos(pDblIn->get(i), 0, pDblOut->get() + i, pDblOut->getImg() + i);
80             }
81         }
82         else //all values are in [-1,1]
83         {
84             pDblOut = new types::Double(pDblIn->getDims(), pDblIn->getDimsArray(), false);
85             for (int i = 0 ; i < pDblIn->getSize() ; i++)
86             {
87                 pDblOut->set(i, dacoss(pDblIn->get(i)));
88             }
89         }
90     }
91
92     out.push_back(pDblOut);
93     return types::Function::OK;
94 }
95 /*--------------------------------------------------------------------------*/