[RESOLVED] NAG Libraries

[agmorgan]

Does anyone use the NAG (Numerical Algorithms Group) libraries?
I'm trying to do a least squares solving algorithm without much luck.
I'm trying to follow the example in the instructions but it crashes the IDE each time I call the routine.

Option Explicit
Option Base 1
 
Public Type Complex
    RE As Double
    IM As Double
End Type
 
 
    Declare Sub F04KLF Lib "DLL20DDS.DLL" (ByRef m As Integer, _
  ByRef n As Integer, ByRef p As Integer, ByRef a As Complex, _
  ByRef lda As Integer, ByRef b As Complex, ByRef ldb As Integer, _
  ByRef d As Complex, ByRef x As Complex, ByRef y As Complex, _
  ByRef work As Complex, ByRef lwork As Integer, ByRef ifail As Integer)
 
 
 
   Sub main()
 
 
 
        '        F04KLF Example Program Data
        '  4  3  4                                               :Values of M, N and P
        '( 0.96,-0.81) (-0.03, 0.96) (-0.91, 2.06)
        '(-0.98, 1.98) (-1.20, 0.19) (-0.66, 0.42)
        '( 0.62,-0.46) ( 1.01, 0.02) ( 0.63,-0.17)
        '( 1.08,-0.28) ( 0.20,-0.12) (-0.07, 1.23)               :End of matrix A
        '( 0.50,-1.00) ( 0.00, 0.00) ( 0.00, 0.00) ( 0.00, 0.00)
        '( 0.00, 0.00) ( 1.00,-2.00) ( 0.00, 0.00) ( 0.00, 0.00)
        '( 0.00, 0.00) ( 0.00, 0.00) ( 2.00,-3.00) ( 0.00, 0.00)
        '( 0.00, 0.00) ( 0.00, 0.00) ( 0.00, 0.00) ( 5.00,-4.00) :End of matrix B
        '( 6.01,-0.37)
        '(-5.27, 0.90)
        '( 2.72,-2.13)
        '(-1.34,-2.80)                                           :End of D
 
        '9.3  Program Results
        ' F04KLF Example Program Results
 
        ' Least-squares solution
        '    (-1.0000, 2.0000) ( 4.0000,-5.0000) (-3.0000, 1.0000)
        
        
        Dim m As Integer
        Dim n As Integer
        Dim p As Integer
        Dim a(4, 3) As Complex
        Dim lda As Integer
        Dim b(4, 4) As Complex
        Dim ldb As Integer
        Dim d(4) As Complex
        Dim x() As Complex
        Dim y() As Complex
        Dim work() As Complex
        Dim lwork As Integer
        Dim ifail As Integer
 
        ifail = 1
 
        m = 4
        n = 3
        p = 4
        ReDim x(p)
        ReDim y(p)
        lwork = n + m + 64 * (m + p)
        ReDim work(lwork)
        lda = m
        ldb = m
 
        '( 0.96,-0.81) (-0.03, 0.96) (-0.91, 2.06)
        '(-0.98, 1.98) (-1.20, 0.19) (-0.66, 0.42)
        '( 0.62,-0.46) ( 1.01, 0.02) ( 0.63,-0.17)
        '( 1.08,-0.28) ( 0.20,-0.12) (-0.07, 1.23)
 
        a(1, 1).RE = 0.96:   a(1, 1).IM = -0.81
        a(1, 2).RE = -0.03:  a(1, 2).IM = 0.96
        a(1, 3).RE = -0.91:  a(1, 3).IM = 2.06
 
        a(2, 1).RE = -0.98:  a(2, 1).IM = 1.98
        a(2, 2).RE = -1.2:   a(2, 2).IM = 0.19
        a(2, 3).RE = -0.66:  a(2, 3).IM = 0.42
        
        a(3, 1).RE = 0.62:   a(3, 1).IM = -0.46
        a(3, 2).RE = 1.01:   a(3, 2).IM = 0.02
        a(3, 3).RE = 0.63:   a(3, 3).IM = -0.17
        
        a(4, 1).RE = 1.08:   a(4, 1).IM = -0.28
        a(4, 2).RE = 0.2:    a(4, 2).IM = -0.12
        a(4, 3).RE = -0.07:  a(4, 3).IM = 1.23
 
 
        '( 0.50,-1.00) ( 0.00, 0.00) ( 0.00, 0.00) ( 0.00, 0.00)
        '( 0.00, 0.00) ( 1.00,-2.00) ( 0.00, 0.00) ( 0.00, 0.00)
        '( 0.00, 0.00) ( 0.00, 0.00) ( 2.00,-3.00) ( 0.00, 0.00)
        '( 0.00, 0.00) ( 0.00, 0.00) ( 0.00, 0.00) ( 5.00,-4.00)
        b(1, 1).RE = 0.5: b(1, 1).IM = -1#
        b(1, 2).RE = 0#: b(1, 2).IM = 0#
        b(1, 3).RE = 0#: b(1, 3).IM = 0#
        b(1, 4).RE = 0#: b(1, 4).IM = 0#
 
        b(2, 1).RE = 0#: b(2, 1).IM = 0#
        b(2, 2).RE = 1#: b(2, 2).IM = -2#
        b(2, 3).RE = 0#: b(2, 3).IM = 0#
        b(2, 4).RE = 0#: b(2, 4).IM = 0#
 
        b(3, 1).RE = 0#: b(3, 1).IM = 0#
        b(3, 2).RE = 0#: b(3, 2).IM = 0#
        b(3, 3).RE = 2#: b(3, 3).IM = -3#
        b(3, 4).RE = 0#: b(3, 4).IM = 0#
 
        b(4, 1).RE = 0#: b(4, 1).IM = 0#
        b(4, 2).RE = 0#: b(4, 2).IM = 0#
        b(4, 3).RE = 0#: b(4, 3).IM = 0#
        b(4, 4).RE = 5#: b(4, 4).IM = -4#
 
        '( 6.01,-0.37)
        '(-5.27, 0.90)
        '( 2.72,-2.13)
        '(-1.34,-2.80)
        d(1).RE = 6.01:  d(1).IM = -0.37
        d(2).RE = -5.27: d(2).IM = 0.9
        d(3).RE = 2.72:  d(3).IM = -2.13
        d(4).RE = -1.34: d(4).IM = -2.8
 
 
        Call F04KLF(m, n, p, a(1, 1), lda, b(1, 1), ldb, d(1), x(1), y(1), work(1), lwork, ifail)
        MsgBox "Done"
    End Sub

[krtxmrtz]

Isn't this library intended to be used with Fortran?

[agmorgan]

Its a DLL so you can call it from a variety of languages.
I think my problem was VBs definition of Integer.
I should have been using Longs.