Description

SF_DISC0 Constant discontinuous shape function.

[ VBASE, NLDOF, XLDOF, SFUN ] = SF_DISC0( I_EVAL, N_SDIM, N_VERT, I_DOF, XI, AINVJAC, VBASE ) Evaluates discontinuous constant shape functions.

Input       Value/[Size]           Description
-----------------------------------------------------------------------------------
i_eval      scalar:  1             Evaluate function values
                    >1             Evaluate values of derivatives
n_sdim      scalar: 1-3            Number of space dimensions
n_vert      scalar: 2-8            Number of vertices per cell
i_dof       scalar: 1-n_ldof       Local basis function to evaluate
xi          [n_sdim(+1)]           Local coordinates of evaluation point
aInvJac     [n,n_sdim(+1)*n_sdim]  Inverse of transformation Jacobian
vBase       [n]                    Preallocated output vector
                                                                                  .
Output      Value/[Size]           Description
-----------------------------------------------------------------------------------
vBase       [n]                    Evaluated function values
nLDof       [4]                    Number of local degrees of freedom on
                                   vertices, edges, faces, and cell interiors
xLDof       [n_sdim,n_ldof]        Local coordinates of local dofs
sfun        string                 Function name of called shape function

See also: sf_disc1

Code listing

 nLDof = [0 0 0 1];
 if( n_vert==(n_sdim+1) )
   xLDof = ones(n_sdim+1,1)/(n_sdim+1);
 else
   xLDof = zeros(n_sdim,1);
 end
 sfun  = 'sf_disc0';


 switch i_eval   % Evaluation type flag.

   case 1        % Evaluation of function values.
     vBase = 1;

   otherwise
     vBase = 0;
 end