forum.alglib.nethttp://forum.alglib.net/ CMatrixGEMM does not work proeprly???http://forum.alglib.net/viewtopic.php?f=2&t=2455 Page 1 of 1

 Author: jonmorga [ Fri Jun 19, 2015 8:08 am ] Post subject: CMatrixGEMM does not work proeprly??? Hello,I am using the CmatrixGemm to multiply two matrixes but it has debugging problems . I think everything is ok defined but it does not work.Sub main()M = 3K = 3N = 1IA = 0JA = 0OptypeA = 0IB = 0JB = 0OptypeB = 0IC = 0JC = 0Dim alpha As Complexalpha.x = 1alpha.y = 0Dim beta As Complexbeta.x = 0beta.y = 0Dim A(1 To 3, 1 To 3) As ComplexDim B(1 To 3, 1 To 1) As ComplexDim C(1 To 3, 1 To 1) As ComplexDim a11 As Complex, a12 As Complex, a13 As Complex, a21 As Complex, a22 As Complex, a23 As Complex, a31 As Complex, a32 As Complex, a33 As Complexa11.x = 1a11.y = 0a12.x = -0.5a12.y = 0a13.x = 0a13.y = 0a21.x = -0.5a21.y = 0a22.x = 1a22.y = 0a23.x = 0a23.y = 0a31.x = 0a31.y = 0a32.x = 0a32.y = 0a33.x = 0a33.y = 3Dim b11 As Complex, b21 As Complex, b31 As Complexb11.x = 1b11.y = 0b21.x = -0.5b21.y = 0b31.x = 0b31.y = 0A(1, 1).x = a11.x: A(1, 1).y = a11.y: A(1, 2).x = a12.x: A(1, 2).y = a12.y: A(1, 3).x = a13.x: A(1, 3).y = a13.yA(2, 1).x = a21.x: A(2, 1).y = a21.y: A(2, 2).x = a22.x: A(2, 2).y = a22.y: A(2, 3).x = a23.x: A(2, 3).y = a23.yA(3, 1).x = a31.x: A(3, 1).y = a31.y: A(3, 2).x = a32.x: A(3, 2).y = a32.y: A(3, 3).x = a33.x: A(3, 3).y = a33.yB(1, 1).x = b11.x: B(1, 1).y = b11.yB(2, 1).x = b21.x: B(2, 1).y = b21.yB(3, 1).x = b31.x: B(3, 1).y = b31.yCall CMatrixGEMM(M, N, K, alpha, A, IA, JA, 0, B, IB, JB, 0, beta, C, IC, JC)End Sub

 Page 1 of 1 All times are UTC Powered by phpBB © 2000, 2002, 2005, 2007 phpBB Grouphttp://www.phpbb.com/