sistema tridiagonale

[maxithron]

dovrei risolvere un sistema tridiagonale.

Suggerimenti per il linguaggio + consono da utilizzare?

[cionci]

Se ti va bene: http://www.prato.linux.it/~fcasadei/math/calcolo/html/node38.html

In generale qualsiasi linguaggio che permette di gestire i vettori...ad esempio quello sopra dovrebbe essere Matlab...

[maxithron]

che può andar bene il link

anche se con google sto trovando esempi con quasi tutti i linguaggi.....ora:

Se cerco di risolvere con un solo linguaggio, mi pentirò per non aver curiosato anche con gli altri ma......se comincio a curiosare anche con gli altri....mi pentirò perchè non riuscirò a risolvere nemmeno con uno

[a2000]

Originally posted by "maxithron"

dovrei risolvere un sistema tridiagonale.

Suggerimenti per il linguaggio + consono da utilizzare?

da volumi finiti 2D o ... colonne a stadi ?

è un classico: algoritmo di Thomas.

il Fortran è la morte sua (che vuol dire 7 righe).

[maxithron]

vramente dovrei utilizzare l'algoritmo di eliminazione di gauss per risolvere:

classico caso Ax = B e la matrice A è tridiagonale.

Dici che il fortran è il miglior metodo?

[a2000]

Nel caso di matrici tridiagonali il metodo più efficiente è quello di Thomas:
due passate una in avanti e una all'indietro e hai risolto.

il linguaggio più adatto è il Fortran poi trovi i listati (si fa per dire: < 10 righe) già bell'e pronti, oltre che in Fortran, anche in Pascal, Basic, C, forse anche Java

[maxithron]

Originally posted by "a2000"


il linguaggio più adatto è il Fortran poi trovi i listati (si fa per dire: < 10 righe) già bell'e pronti, oltre che in Fortran, anche in Pascal, Basic, C, forse anche Java

Nooooo....non voglio trovare i listati già fatti! Sono un masochista convinto

Ehm..cmq... sette righe dici?

[a2000]

9 con

Subroutine

e

End Subroutine

se vuoi te lo passo

[maxithron]

Originally posted by "a2000"


se vuoi te lo passo

Non c'è sfizio :o

hai detto proprio 9?

no perchè io

[a2000]

con Thomas sì ....

con Gauss ... 15

beh però con la ricerca pivot ... 20 !




te lo passo ????

[maxithron]

Originally posted by "a2000"



te lo passo ????

No......No......No........


Com'è?

[a2000]

ma da dove nasce il sistema ?

no, perchè se è a diagonale dominante è garantita la convergenza anche con metodi iterativi: 5 righe

[a2000]

Originally posted by "maxithron"



Com'è?

pezzettino (cioè la metà di andata)

Do i = 2, II
   DEN = P(i) - W(i) * EE(i - 1)
   EE(i) = E(i) / DEN
   BB(i) = (B(i) + W(i) * BB(i - 1)) / DEN
End Do

[maxithron]

Grrrrrrrr!!!!!

5 righe????

[maxithron]

ma non è così?:

program main
real, dimension (10):: a,d,c,b,x
integer ::n = 10

interface
subroutine tri(n,a,d,c,b,x)
integer, intent(in)::n
real, dimension(, intent(in):: a,d,c,b
real, dimension(, intent(out)::x
end subroutine tri
end interface
do i=1,n
d(i) = 2.0
a(i) = 0.5
c(i) = 0.5
b(i) = 3.0
end do
b(1) = 2.5
b(n) = 2.5
call tri(n,a,d,c,b,x)
print*, "subroutine tri"
print"(f22.14)",(x(i),i=1,n)
do i=1,n
d(i) = 2.0
c(i) = 0.5
b(i) = 3.0
end do
b(1) = 2.5
b(n) = 2.5
call tri(n,c,d,c,b,x)
print*, "subroutine tri again"
print "(f22.14)",(x(i),i=1,n)
end program main

subroutine tri(n,a,d,c,b,x)
integer, intent(in)::n
real, dimension(, intent(in)::a,d,c,b
real, dimension(, intent(out):: x
integer ::i
real :: xmult
do i = 2,n
xmult = a(i-1)/d(i-1)
d(i) = d(i) - xmult*c(i-1)
b(i) = b(i) - xmult*b(i-1)
end do
x(n) = b(n)/d(n)
do i = n-1,1,-1
x(i) = (b(i) - c(i)*x(i+1))/d(i)
end do
end subroutine tri

[a2000]

5 + 4 = 9

ah, ma ca++o non mi dai elementi ...

da dove viene sto' sistema o è un esercizio di calcolo numerico

cosa studi ?

[a2000]

Originally posted by "maxithron"

ma non è coì?:
....
subroutine tri(n,a,d,c,b,x)
integer, intent(in)::n
real, dimension(, intent(in)::a,d,c,b
real, dimension(, intent(out):: x
integer ::i
real :: xmult
do i = 2,n
xmult = a(i-1)/d(i-1)
d(i) = d(i) - xmult*c(i-1)
b(i) = b(i) - xmult*b(i-1)
end do
x(n) = b(n)/d(n)
do i = n-1,1,-1
x(i) = (b(i) - c(i)*x(i+1))/d(i)
end do
end subroutine tri[/size]

sì bravo, con qualche ridondanza, ma direi che ci siamo !

[maxithron]

E' un esercizio di calcolo numerico.

Cosa studio?

Quello che mi capita...

Scherzi a parte non sono uno studente universitario ma sono divorato da una strana passione schizo-frenetico-psicotica per la matematica e la programmazione

[a2000]

e figa ?

[maxithron]

ca@@o ci fanno quelle faccine nel list?