/*
 Copyright (C) 2013 Marco Vassallo

 This program is free software; you can redistribute it and/or modify it under
 the terms of the GNU General Public License as published by the Free Software
 Foundation; either version 2 of the License, or (at your option) any later
 version.

 This program is distributed in the hope that it will be useful, but WITHOUT
 ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
 FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License for more
 details.

 You should have received a copy of the GNU General Public License along with
 this program; if not, see <http://www.gnu.org/licenses/>.
*/

#include "mesh.h"

DEFUN_DLD (fem_init_mesh, args, ,"-*- texinfo -*-\n\
@deftypefn {Function File} {[@var{mesh_out}]} = \
fem_get_mesh (@var{mesh_in}) \n\
The @var{mesh_in} should be either\n\
@itemize @bullet \n\
@item a string containing the name of the file where the mesh is stored\
(compatible formats are the ones which are compatible with Fenics\n\
@item a a PDE-tool like structure with matrix fields (p,e,t)\n\
@end itemize\n\
fem_init_mesh().\n\
The output @var{mesh_out} is a representation of the\n\
@var{mesh_in} which is compatible with fem-fenics\n\
@seealso{fem_get_mesh}\n\
@end deftypefn")
{
  int nargin = args.length ();
  octave_value retval = 0;

  if (nargin < 1 || nargin > 1)
    print_usage ();
  else
    {
      if (!error_state)
        {
          if (args(0).is_string () == true)
            {
              std::string filename = args(0).string_value ();
              //if the filename is not valid, dolfin takes care of it

              if (! mesh_type_loaded)
                {
                  mesh::register_type ();
                  mesh_type_loaded = true;
                  mlock ();
                }
              retval = new mesh (filename);
            }

          else
            {
              octave_scalar_map a = args(0).scalar_map_value ();
              Array<double> p = a.contents ("p").matrix_value ();
              Array<octave_idx_type> t = a.contents ("t").matrix_value ();
              Array<octave_idx_type> e = a.contents ("e").matrix_value ();

              if (! error_state)
                {
                  if (! mesh_type_loaded)
                    {
                      mesh::register_type ();
                      mesh_type_loaded = true;
                      mlock ();
                    }

                  retval = new mesh (p, e, t);
                }
              else
               error ("fem_init_mesh: the argument you provide is invalid");
            }
        }
    }
  return retval;
}