eop_fd_operators.m File Reference

[L_I, L_E, b, L_I_adj, L_E_adj] = eop_fd_operators(model, model_data) More...

Go to the source code of this file.

Functions
function [ L_I , L_E , b , L_I_adj , L_E_adj ] =	eop_fd_operators (model, model_data)
	[L_I, L_E, b, L_I_adj, L_E_adj] = eop_fd_operators(model, model_data) More...

Detailed Description

[L_I, L_E, b, L_I_adj, L_E_adj] = eop_fd_operators(model, model_data)

Definition in file eop_fd_operators.m.

Function Documentation

function [ L_I , L_E , b , L_I_adj , L_E_adj ] = eop_fd_operators	(		model,
			model_data
	)

[L_I, L_E, b, L_I_adj, L_E_adj] = eop_fd_operators(model, model_data)

function which computes either L_I, L_E and b the components of the linear evolution scheme or L_I_adj and L_E_adj the components of the corresponding dual problem (controlled via model.want_dual). Supports affine decomposition via model.decomp_mode

decomp_mode: operation mode of the function

• 0=complete (default): no parameter dependence or decomposition is performed.
• 1=components: For each output argument a cell array of output arguments is returned representing the q-th component independent of the parameters given in mu_names
• 2=coefficients: For each output argument a cell array of output arguments is returned representing the q-th coefficient dependent of the parameters given in mu_names

the construction of L_I follows a 9-point-star described in my diploma thesis (chapter 4.2)

Parameters

model	model
model_data	model data

Return values

L_I	L I
L_E	L E
b	b
L_I_adj	L I adj
L_E_adj	L E adj

Required fields of model:

• decomp_mode —  flag indicating the operation mode of the function:
  • 0 (complete) : no affine parameter dependence or decomposition is performed.
  • 1 (components) : for each output argument a cell array of output matrices is returned representing the \(q\)-th component independent of the parameters given in mu_names.
  • 2 (coefficients) : returns a vector where each coordinate represents the \(q\)-the coefficient \(\sigma_{\cdot}^{q}(\mu)\) dependent on the parameters given in mu_names.
• sigma1 —  sigma1
• rho —  rho
• sigma2 —  sigma2
• N1 —  N1
• deltaT —  deltaT
• h1 —  h1
• r —  r
• N2 —  N2
• h2 —  h2
• want_dual —  want dual

Required fields of model_data:

• grid.nvertices —  grid.nvertices
• grid.X —  grid.X
• grid.Y —  grid.Y
• grid —  a structure containing geometry information of a mesh used for discretizations

Definition at line 17 of file eop_fd_operators.m.

Here is the caller graph for this function: