elementary_functions module.
[scilab.git] / scilab / modules / elementary_functions / sci_gateway / cpp / sci_expm.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 "expm.h"
25 }
26
27 /*--------------------------------------------------------------------------*/
28 types::Function::ReturnValue sci_expm(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 to %d expmected.\n"), L"expm", 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 expmected.\n"), L"expm", 1);
42         return types::Function::Error;
43     }
44
45     if (in[0]->isDouble() == false)
46     {
47         std::wstring wstFuncName = L"%"  + in[0]->getShortTypeStr() + L"_expm";
48         return Overload::call(wstFuncName, in, _iRetCount, out, new ExecVisitor());
49     }
50
51     pDblIn = in[0]->getAs<types::Double>();
52
53     if (pDblIn->getDims() > 2)
54     {
55         std::wstring wstFuncName = L"%hm_expm";
56         return Overload::call(wstFuncName, in, _iRetCount, out, new ExecVisitor());
57     }
58
59     if (pDblIn->getSize() == 0)
60     {
61         out.push_back(types::Double::Empty());
62         return types::Function::OK;
63     }
64
65     if (pDblIn->getCols() != pDblIn->getRows())
66     {
67         ScierrorW(999, _W("%ls: Wrong size for input argument #%d : A square matrix expected.\n"), L"expm", 1);
68         return types::Function::Error;
69     }
70
71     pDblOut = new types::Double(pDblIn->getDims(), pDblIn->getDimsArray(), pDblIn->isComplex());
72
73     if (pDblIn->isComplex())
74     {
75         zexpms2(pDblIn->get(), pDblIn->getImg(), pDblOut->get(), pDblOut->getImg(), pDblIn->getCols());
76     }
77     else
78     {
79         dexpms2(pDblIn->get(), pDblOut->get(), pDblIn->getCols());
80     }
81
82     out.push_back(pDblOut);
83     return types::Function::OK;
84 }
85 /*--------------------------------------------------------------------------*/