Skip to content

Commit 18d0512

Browse files
Fabio VerboFabio Verbo
authored andcommitted
Header changed
1 parent 69cacf5 commit 18d0512

File tree

5 files changed

+15
-2
lines changed

5 files changed

+15
-2
lines changed

diffreactstencil.m

Lines changed: 3 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -1,5 +1,8 @@
11
function S = diffreactstencil(U, x_old, nx, ny, bnd, alpha, dxs)
22
% DIFFREACTSTENCIL computes stencil as inner + bounds + corners
3+
%
4+
% Fabio VERBOSIO, Universita` della Svizzera italiana, November 2018
5+
%
36

47
bndN = bnd.bndN;
58
bndS = bnd.bndS;

diffusereact.m

Lines changed: 3 additions & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -1,6 +1,8 @@
11
function B = diffusereact(X, Xold, alpha, dxs)
22
% DIFFUSEREACT computes Diffusion/Reaction using Laplacian matrix.
3-
3+
%
4+
% Fabio VERBOSIO, Universita` della Svizzera italiana, November 2018
5+
%
46
global A;
57

68
aI = alpha*speye(size(A));

drstencil.m

Lines changed: 3 additions & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -1,6 +1,8 @@
11
function S = drstencil(X, Xold, nx, ny, alpha, dxs)
22
% DRSTENCIL computes the stencil for a grid including BCs
3-
3+
%
4+
% Fabio VERBOSIO, Universita` della Svizzera italiana, November 2018
5+
%
46
X = reshape(X,nx,ny);
57
Xold = reshape(Xold,nx,ny);
68
S = zeros(size(X));

reactdiff.m

Lines changed: 3 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -1,5 +1,8 @@
11
% REACTDIFF script for testing stencil-based CG for Diffusion/Reaction problems.
22
%
3+
%
4+
% Fabio VERBOSIO, Universita` della Svizzera italiana, November 2018
5+
%
36

47
%% SET PARAMS
58

stencil_pcg.m

Lines changed: 3 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -1,5 +1,8 @@
11
function [x, flg, resNew, its] = stencil_pcg(f, b, tol, maxit)
22
% STENCIL_PCG computes CG with 1st order approximation for f
3+
%
4+
% Fabio VERBOSIO, Universita` della Svizzera italiana, November 2018
5+
%
36

47
its = 0;
58
flg = 1;

0 commit comments

Comments
 (0)