kronecker product

[Saman Kaviani]

Hi there

Can anybody tell me how I can do Kronecker product in fortran90?

NB: I know google is my friend and I found two subroutines and they did not work. If you know the answer please respond this email.

best, 

Saman 

[Youngung Jeong]

Hi,
Just for curiosity, when you said they don't work, you mean they are not compilable or what?

Saman 

--
You received this message because you are subscribed to the Google Groups "gg95" group.
To post to this group, send email to gg...@googlegroups.com.
To unsubscribe from this group, send email to gg95+uns...@googlegroups.com.
For more options, visit this group at http://groups.google.com/group/gg95?hl=en.

[Saman Kaviani]

On Thu, Jun 28, 2012 at 4:33 PM, Youngung Jeong <youngun...@gmail.com> wrote:

Hi,
Just for curiosity, when you said they don't work, you mean they are not compilable or what?

I mean there is always error in compiling!

here is the code maybe you can try and let me the solution :

!-----------------------------------------------------------------------------------------------

SUBROUTINE KRON(K, A, B)

    IMPLICIT NONE

    REAL, INTENT(IN)  :: A(:,:), B(:,:)

    REAL, INTENT(INOUT) :: K(:,:)

    INTEGER :: I, J, MA, NA, MB, NB

    MA = UBOUND(A, 1)

    NA = UBOUND(A, 2)

    MB = UBOUND(B, 1)

    NB = UBOUND(B, 2)

    IF (SIZE(K,1) /= MA*MB .OR. SIZE(K,2) /= NA*NB) THEN

        WRITE(*,*) 'K has invalid size'

        CALL ABORT

    END IF

    FORALL(I=1:MA, J=1:NA)

        K(MB*(I-1)+1:MB*I,NB*(J-1)+1:NB*J) = A(I,J)*B

    END FORALL

END SUBROUTINE KRON

!----------------------------------------------------------------------------------------------------------------

[Youngung Jeong]

What error messages you've got? What compiler did you use?

  Youngung Jeong, 정영웅

[evan]

Delete the following line

CALL ABORT