4 id = 'L_pressure_mean';
7 arg_vars = {
'pressure'};
8 arg_vars_short = {
'P'};
12 function [INC, INCoff, J] = apply(
this, model, model_data, sim_data, NU_ind)
13 %
function INC = apply(
this, model, model_data, sim_data, NU_ind)
15 grid = model_data.grid;
19 if ~isfield(sim_data, 'm')
20 sim_data.m = length(P);
23 INC = sum(grid.A.* P);
31 Plength = grid.nelements;
32 if ~isfield(sim_data, 'Poff')
35 J = sparse(ones(1,Plength), (1:Plength) + sim_data.Poff, grid.A, n, sim_data.m);
41 function [LL, bb] = full_matrix(this, descr, model_data)
42 grid = model_data.grid;
45 LL = sparse(ones(1,m),1:m, grid.A, n, m);
49 function ipm = inner_product_matrix(this, model_data)
53 function [LL_comps, bb_comps] = components(this, descr, model_data)
54 [LL, bb] = full_matrix(this, descr, model_data);
59 function [LL_coeffs, bb_coeffs] = coefficients(
this, descr);