Mercurial > fem-fenics-eugenio
view inst/ResidualForm.m @ 213:d4eef566511b
Extend feval function to support also passing the coordinates separately
as multiple arguments and evaluating the function at multiple
points at once. Author: Daniel Kraft
author | gedeone-octave <marcovass89@hotmail.it> |
---|---|
date | Tue, 04 Mar 2014 12:19:04 +0000 |
parents | 0e6d1d79244a |
children |
line wrap: on
line source
## Copyright (C) 2013 Marco Vassallo <gedeone-octave@users.sourceforge.net> ## ## 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 3 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/>. ## -*- texinfo -*- ## @deftypefn {Function File} {[@var{L}]} = @ ## LinearForm (@var{my_problem}, @var{U}, @var{coefficient_1}, @ ## @var{coefficient_2},...) ## ## Construct a ResidualForm previously imported from a ufl file with the ## function import_ufl_LinearForm. ## ## The compulsory arguments are: ## @itemize @bullet ## @item @var{my_problem} the name of the problem to solve. ## @item the FunctionSpace @var{U} where the problem is defined. ## @end itemize ## ## The optional arguments are the @var{coefficient_1}, @var{coefficient_2} ## which specify the parameters for the ResidualForm with the same name which ## was used in the ufl file. ## They can be either a Constant, a Function or an Expression. ## ## @seealso{import_ufl_LinearForm, import_ufl_Problem, BilinearForm, ## ResidualForm, BilinearForm} ## @end deftypefn function a = ResidualForm (name, V, varargin) if nargin < 2 error ("ResidualForm: wrong number of input parameters."); elseif ! ischar (name) error ("ResidualForm: first argument is not a valid string"); endif program = sprintf ("%s_LinearForm(V", name); for k = 1:length (varargin) eval(['f_' num2str(k) '=varargin{k};']); program = strjoin ({ program, strcat('f_',num2str(k))}, ','); end program = strjoin ({program, ');'}); a = eval (program); endfunction