CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Main CFD Forum

the solution for illcondizioned matrix

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Display Modes
Old   December 9, 2008, 07:35
Default the solution for illcondizioned matrix
  #1
giuseppe casillo
Guest
 
Posts: n/a
%the A matrix end the b know terms must be i memory

%programma per la risoluzione di matrici malcondizionate e %p.i.m % matrice in A terminenoto in b %risolve pure le matrici a determinante nullo % Ainvertallok

AA=A*A;

sp=max(eig(AA))/1.9;

n=size(A);

m=n(1);

xold=zeros(m,1);

Ab=A*b/sp;

AA=AA/sp;

for i=1:2000

xn=Ab+xold-AA*xold ;

s=0.0;

for l=1:m

if (xn(l)==0.0)

else

s=s+abs((xn(l)-xold(l))/xn(l));

end;

end;

xold=xn;

if ( s <(0.0015*m))

break,

end;

end;

i

xn

%

%questo metodo risolve in media i problemi %con determinante nullo:in media perchè è % A*(A*x-b)=0
  Reply With Quote

Reply

Thread Tools
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
grid dependancy gueynard a. Main CFD Forum 19 June 27, 2014 21:22
Neumann Boundary Condition for Poisson Equation solution in Polar Coordinates prapanj Main CFD Forum 2 July 30, 2011 19:07
Steady solution from Transient simulations wawa FLUENT 2 November 9, 2010 18:44
Unsteady solution Christophe FLUENT 0 August 11, 2006 11:13
Wall functions Abhijit Tilak Main CFD Forum 6 February 5, 1999 02:16


All times are GMT -4. The time now is 16:02.