D:/demona/2.0/manual/doc/2-31set/matlab/cont.m File Reference

Go to the source code of this file.

Functions/Subroutines
*	norm (f0)%delta
	if (norm(f0)<=tol) disp('Newton converged') break end end s = zeros(N/unk,1)
	v (i)
	subplot (1, 2, 1) contourf(S')
	subplot (1, 2, 2) contourf(V')
	disp ('maximum temperature:')%max(x)%pause disp('re = myFp(x,n)
	delta (133:134) reOld
	x (133:134) f0
	surf (X')
	delta (67)%%lambdaOld
	x (67)%%%%f0
Variables
clear all clc	n = 65
global lambda global unk global re	lambda = 0.1
	re = 1000
	unk = 2
	N = n*n*unk
	maxNewton = 10
	tol = 1e-5
	maxit = 1000
initial guess for the solution vector	x = zeros(N,1)
	f0 = myF(x,n)
	initial_norm
for	i
	Jac = myJac(x,n)
	delta = gmres(Jac,-f0,min(100,N),tol,min(1000,N))
	v = zeros(N/unk,1)
end	S = reshape(s,n,n)
	V = reshape(v,n,n)

---

Function Documentation

delta

(

67

)

delta

(

133:134

)

disp

(

'maximum temperature:'

)

= myFp(x,n)

if

(

norm(f0)<=

tol

)

= zeros(N/unk,1)

* norm

(

f0

)

subplot	(	1	,
		2	,
		2
	)

subplot	(	1	,
		2	,
		1
	)

surf

(

X'

)

v

(

i

)

x

(

67

)

x

(

133:134

)

---

Variable Documentation

delta = gmres(Jac,-f0,min(100,N),tol,min(1000,N))

Definition at line 45 of file cont.m.

f0 = myF(x,n)

Definition at line 27 of file cont.m.

end i

Initial value:

1:maxNewton
    
    disp('Newton iter')
    i
    
    %initialize the update vector
    delta=zeros(N,1)

Definition at line 33 of file cont.m.

Referenced by generateCoordinates().

initial_norm

Initial value:

norm(f0)

jac=sparse(N,N)

Definition at line 29 of file cont.m.

Jac = myJac(x,n)

Definition at line 41 of file cont.m.

global lambda lambda = 0.1

Definition at line 9 of file cont.m.

maxit = 1000

Definition at line 21 of file cont.m.

maxNewton = 10

Definition at line 18 of file cont.m.

N = n*n*unk

Definition at line 13 of file cont.m.

clear all clc n = 65

Definition at line 3 of file cont.m.

re = 1000

Definition at line 10 of file cont.m.

end S = reshape(s,n,n)

Definition at line 69 of file cont.m.

tol = 1e-5

Definition at line 19 of file cont.m.

unk = 2

Definition at line 11 of file cont.m.

V = reshape(v,n,n)

Definition at line 72 of file cont.m.

v = zeros(N/unk,1)

Definition at line 61 of file cont.m.

initial guess for the solution vector x = zeros(N,1)

Definition at line 25 of file cont.m.