Howdy, Stranger!

It looks like you're new here. If you want to get involved, click one of these buttons!

Sign In with Facebook Sign In with Google Sign In with OpenID

Categories

We have migrated to a new platform! Please note that you will need to reset your password to log in (your credentials are still in-tact though). Please contact lee@programmersheaven.com if you have questions.
Welcome to the new platform of Programmer's Heaven! We apologize for the inconvenience caused, if you visited us from a broken link of the previous version. The main reason to move to a new platform is to provide more effective and collaborative experience to you all. Please feel free to experience the new platform and use its exciting features. Contact us for any issue that you need to get clarified. We are more than happy to help you.

Fortran code Segmentation error :(

Hi guys I've been stuck on this for the best part of a day now. I'm new to programming (3 days so far) and just can' seem to get this to work!!! My code is below:

Module powermethord
use tri
implicit none
contains

subroutine powerit(a,b,c,E,n)
implicit none

real, intent(in) :: E
real, intent(inout) :: a(:),b(:),c(:)
!real, intent(out) :: x(:)
real,allocatable :: x(:),d(:)
real :: phi_0,lambda_0, lambda, newsum, oldsum, phi
integer :: k, I, j
integer, intent(in) :: n

j = 10
lambda_0 = 1
phi_0 = 1
phi = phi_0
lambda = lambda_0
oldsum = phi_0
allocate(d(I))
d(1)=1
Do k = 1, j-1
print*, d
call tridiag(a,b,c,d,n,x)
d(k) = lambda*(E*x(i))
newsum = sum(x)
lambda = (newsum/oldsum)
print*, k, lambda, x(i)
oldsum = newsum
End do
end subroutine powerit
end module powermethord

the sub routine tridiag is given by:

Module tri
implicit none
contains

subroutine tridiag(a,b,c,d,n,x)
implicit none



integer, intent(in) :: n
real, intent(in) :: a(:),b(:),c(:),d(:)
real, intent(out) :: x(:)
real :: bp(n),dp(n)
real :: dx
integer j,l

l=10
bp(1) = b(1)
dp(1) = d(1)
firstpass: do j = 2,n
dx = a(j)/bp(j-1)
bp(j) = b(j) - dx*c(j-1)
dp(j) = d(j) - dx*dp(j-1)
end do firstpass
x(n) = dp(n)/bp(n)
backsub:do j = l-1, 1, -1
x(j) = (dp(j) - c(j)*x(j+1))/bp(j)
print*, j,x(j)
end do backsub
!print*, i,x(i)
end subroutine tridiag
end module tri

any ideas what's going wrong?! Thanks guys!!!


Sign In or Register to comment.