% /*-----------------------------------------------------------------*/
% /*! 
%   \file mcomputarray.m
%   \brief Check the output of the function calceph_compute with an INPOP file
%          for an array of date.
% 
%   \author  M. Gastineau 
%            Astronomie et Systemes Dynamiques, IMCCE, CNRS, Observatoire de Paris. 
% 
%    Copyright, 2018, CNRS
%    email of the author : Mickael.Gastineau@obspm.fr
% */
% /*-----------------------------------------------------------------*/
%  
% /*-----------------------------------------------------------------*/
% /* License  of this file :
%  This file is 'triple-licensed', you have to choose one  of the three licenses 
%  below to apply on this file.
%  
%     CeCILL-C
%     	The CeCILL-C license is close to the GNU LGPL.
%     	( http://www.cecill.info/licences/Licence_CeCILL-C_V1-en.html )
%    
%  or CeCILL-B
%         The CeCILL-B license is close to the BSD.
%         (http://www.cecill.info/licences/Licence_CeCILL-B_V1-en.txt)
%   
%  or CeCILL v2.1
%       The CeCILL license is compatible with the GNU GPL.
%       ( http://www.cecill.info/licences/Licence_CeCILL_V2.1-en.html )
%  
% 
%  This library is governed by the CeCILL-C, CeCILL-B or the CeCILL license under 
%  French law and abiding by the rules of distribution of free software.  
%  You can  use, modify and/ or redistribute the software under the terms 
%  of the CeCILL-C,CeCILL-B or CeCILL license as circulated by CEA, CNRS and INRIA  
%  at the following URL 'http://www.cecill.info'. 
%  
%  As a counterpart to the access to the source code and  rights to copy,
%  modify and redistribute granted by the license, users are provided only
%  with a limited warranty  and the software's author,  the holder of the
%  economic rights,  and the successive licensors  have only  limited
%  liability. 
%  
%  In this respect, the user's attention is drawn to the risks associated
%  with loading,  using,  modifying and/or developing or reproducing the
%  software by the user in light of its specific status of free software,
%  that may mean  that it is complicated to manipulate,  and  that  also
%  therefore means  that it is reserved for developers  and  experienced
%  professionals having in-depth computer knowledge. Users are therefore
%  encouraged to load and test the software's suitability as regards their
%  requirements in conditions enabling the security of their systems and/or 
%  data to be ensured and,  more generally, to use and operate it in the 
%  same conditions as regards security. 
%  
%  The fact that you are presently reading this means that you have had
%  knowledge of the CeCILL-C,CeCILL-B or CeCILL license and that you accept its terms.
%  */
%  /*-----------------------------------------------------------------*/

% /*-----------------------------------------------------------------*/
% /* main program */
% /*-----------------------------------------------------------------*/
function res = mcomputearray(prefetch)
        res = 1;
        peph = CalcephBin.open(openfiles('../../examples/example1.dat'));
        if (prefetch==0)
            peph.prefetch();
        end    
        
        target = 11;
        center = 4;
        jd0 = linspace(1,10, 10)+ 2442457;
        dt = linspace(1,10, 10)/10;
        PVarray = peph.compute(jd0, dt, target, center);
        if (size(PVarray,2)~=6)
            size(PVarray)
            size(jd0)
            error('Bad size array')
        end
        if (size(jd0,2)~=size(PVarray,1))
            size(PVarray)
            size(jd0)
            error('Bad size array')
        end
        
        for p=1:size(jd0,2)
            PV = peph.compute(jd0(p), dt(p), target, center);
            aline = PVarray(p,:);
            diff = PV - aline;
            if (max(abs(diff))~=0)
                    format longEng
                    printf('failure for %d %d at time %23.16E %23.16E\n',target, center, jd0(p), dt(p));
                    aline
                    PV
                    diff
                    printf('prefetch=%d\n', prefetch);
                    error('Bad array computation')
            end
        end    
           
        peph.close();
end
     
%!assert(mcomputearray(0)==1)
