Mercurial > fem-fenics-eugenio
view src/fem_bc.cc @ 44:fca8c3d75036
register_type is called before every type is used.
* A new static variable for every class has been defined.
It is named class_type_loaded, and it specify if the corresponding type have
yet been loaded or not. If not, it proceed with the registration.
This should be only a temporary solution.
author | gedeone-octave <marco.vassallo@outlook.com> |
---|---|
date | Mon, 22 Jul 2013 11:33:50 +0200 |
parents | 74cbe97f55af |
children | fcfecdd3a9b5 |
line wrap: on
line source
/* 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 "boundarycondition.h" #include "functionspace.h" #include "expression.h" DEFUN_DLD (fem_bc, args, , "-*- texinfo -*-\n\ @deftypefn {Function File} {[@var{bc}]} = \ fem_get_mesh (@var{Functional Space}, @var{Functio handle}, \ @var{Boundary}) \n\ The input parameters are\n\ @itemize @bullet \n\ @item @var{Functional Space} is a fem-fenics functional space where\ we want to apply the BC\n\ @item @var{Functio handle} is a function handle which contains the expression\ that we want to apply on the BC\n\ @item @var{Boundary} is an Array which specify the label of the \ sides where the BC is applied\n\ The output @var{BC} is an object which contains the boundary conditions\n\ @seealso{fem_init_mesh, fem_fs}\n\ @end deftypefn") { int nargin = args.length (); octave_value retval=0; if (nargin < 3 || nargin > 3) print_usage (); else { if (! functionspace_type_loaded) { functionspace::register_type (); functionspace_type_loaded = true; mlock (); } if (args(0).type_id () == functionspace::static_type_id ()) { const functionspace & fspo = static_cast<const functionspace&> (args(0).get_rep ()); octave_fcn_handle * fh = args(1).fcn_handle_value (); Array<octave_idx_type> side = args(2).array_value (); if (!error_state) { const boost::shared_ptr <const dolfin::FunctionSpace> & V (fspo.get_pfsp ()); const expression * pf = new expression (*fh); boost::shared_ptr<const expression> f (pf); if (! boundarycondition_type_loaded) { boundarycondition::register_type (); boundarycondition_type_loaded = true; mlock (); } boundarycondition * pbc = new boundarycondition (); for (octave_idx_type i = 0; i < side.length (); ++i) pbc->add_bc (V, f, side(i)); retval = octave_value (pbc); } } } return retval; }