IdentifiantMot de passe
Loading...
Mot de passe oublié ?Je m'inscris ! (gratuit)
Navigation

Inscrivez-vous gratuitement
pour pouvoir participer, suivre les réponses en temps réel, voter pour les messages, poser vos propres questions et recevoir la newsletter

 .NET Discussion :

Equivalence syntaxe vba -> vb.net


Sujet :

.NET

  1. #1
    Candidat au Club
    Homme Profil pro
    Développeur .NET
    Inscrit en
    octobre 2019
    Messages
    3
    Détails du profil
    Informations personnelles :
    Sexe : Homme
    Localisation : Espagne

    Informations professionnelles :
    Activité : Développeur .NET
    Secteur : Bâtiment Travaux Publics

    Informations forums :
    Inscription : octobre 2019
    Messages : 3
    Points : 4
    Points
    4
    Par défaut Equivalence syntaxe vba -> vb.net
    Bonjour,

    Avant tout, desolé si c'est pas le forum correct pour cette question.

    J'essai de "Traudire" un logiciel excel vers vb.net.

    Le problemme est que la syntaxe n'est pas exactement la meme, et en plus l'option Base 1 de vba n'est pas disponible sur vb.net.

    J'ai essaie de la "repliquer" avec changemments comme:

    Dim nf as long (1 to ndof, 1 to nn)

    Remplace par:

    Dim nf As array = Array.CreateInstance(GetType(Long), {ndof,nn}, {1,1})

    Mais le code n'arrive a marcher bien.

    Est ce que quelqu'un pourrait m'orienter?

    Vous pouvez trouver un exemple du logiciel d'origin:

    https://newtonexcelbach.com/2009/09/...oads-on-piles/






    Les 2 codes sont ci desus:
    Code vba : Sélectionner tout - Visualiser dans une fenêtre à part
    1
    2
    3
    4
    5
    6
    7
    8
    9
    10
    11
    12
    13
    14
    15
    16
    17
    18
    19
    20
    21
    22
    23
    24
    25
    26
    27
    28
    29
    30
    31
    32
    33
    34
    35
    36
    37
    38
    39
    40
    41
    42
    43
    44
    45
    46
    47
    48
    49
    50
    51
    52
    53
    54
    55
    56
    57
    58
    59
    60
    61
    62
    63
    64
    65
    66
    67
    68
    69
    70
    71
    72
    73
    74
    75
    76
    77
    78
    79
    80
    81
    82
    83
    84
    85
    86
    87
    88
    89
    90
    91
    92
    93
    94
    95
    96
    97
    98
    99
    100
    101
    102
    103
    104
    105
    106
    107
    108
    109
    110
    111
    112
    113
    114
    115
    116
    117
    118
    119
    120
    121
    122
    123
    124
    125
    126
    127
    128
    129
    130
    131
    132
    133
    134
    135
    136
    137
    138
    139
    140
    141
    142
    143
    144
    145
    146
    147
    148
    149
    150
    151
    152
    153
    154
    155
    156
    157
    158
    159
    160
    161
    162
    163
    164
    165
    166
    167
    168
    169
    170
    171
    172
    173
    174
    175
    176
    177
    178
    179
    180
    181
    182
    183
    184
    185
    186
    187
    188
    189
    190
    191
    192
    193
    194
    195
    196
    197
    198
    199
    200
    201
    202
    203
    204
    205
    206
    207
    208
    209
    210
    211
    212
    213
    214
    215
    216
    217
    218
    219
    220
    221
    222
    223
    224
    225
    226
    227
    228
    229
    230
    231
    232
    233
    234
    235
    236
    237
    238
    239
    240
    241
    242
    243
    244
    245
    246
    247
    248
    249
    250
    251
    252
    253
    254
    255
    256
    257
    258
    259
    260
    261
    262
    263
    264
    265
    266
    267
    268
    269
    270
    271
    272
    273
    274
    275
    276
    277
    278
    279
    280
    281
    282
    283
    284
    285
    286
    287
    288
    289
    290
    291
    292
    293
    294
    295
    296
    297
    298
    299
    300
    301
    302
    303
    304
    305
    306
    307
    308
    309
    310
    311
    312
    313
    314
    315
    316
    317
    318
    319
    320
    321
    322
    323
    324
    325
    326
    327
    328
    329
    330
    331
    332
    333
    334
    335
    336
    337
    338
    339
    340
    341
    342
    343
    344
    345
    346
    347
    348
    349
    350
    351
    352
    353
    354
    355
    356
    357
    358
    359
    360
    361
    362
    363
    364
    365
    366
    367
    368
    369
    370
    371
    372
    373
    374
    375
    376
    377
    378
    379
    380
    381
    382
    383
    384
    385
    386
    387
    388
    389
    390
    391
    392
    393
    394
    395
    396
    397
    398
    399
    400
    401
    402
    403
    404
    405
    406
    407
    408
    409
    410
    411
    412
    413
    414
    415
    416
    417
    418
    419
    420
    421
    422
    423
    424
    425
    426
    427
    428
    429
    430
    431
    432
    433
    434
    435
    436
    437
    438
    439
    440
    441
    442
    443
    444
    445
    446
    447
    448
    449
    450
    451
    452
    453
    454
    455
    456
    457
    458
    459
    460
    461
    462
    463
    464
    465
    466
    467
    468
    469
    470
    471
    472
    473
    474
    475
    476
    477
    478
    479
    480
    481
    482
    483
    484
    485
    486
    487
    488
    489
    490
    491
    492
    493
    494
    495
    496
    497
    498
    499
    500
    501
    502
    503
    504
    505
    506
    507
    508
    509
    510
    511
    512
    513
    514
    515
    516
    517
    518
    519
    520
    521
    522
    523
    524
    525
    526
    527
    528
    529
    530
    531
    532
    533
    534
    535
    536
    537
    538
    539
    540
    541
    542
    543
    544
    545
    546
    547
    548
    549
    550
    551
    552
    553
    554
    555
    556
    557
    558
    559
    560
    561
    562
    563
    564
    565
    566
    567
    Option Explicit
     
    Option Base 1
    Function BoEFA(beamA As Variant, LoadA As Variant, Optional FixedNodes As Variant, Optional NIP As Long = 2) As Variant
        Dim Numnodes As Long, resa() As Double
     
        If TypeName(beamA) = "Range" Then beamA = beamA.Value2
        If TypeName(LoadA) = "Range" Then LoadA = LoadA.Value2
        If TypeName(FixedNodes) = "Range" Then FixedNodes = FixedNodes.Value2
     
        Numnodes = UBound(beamA)
        If IsMissing(FixedNodes) = True Then FixedNodes = 0
     
        ReDim resa(1 To Numnodes, 1 To 4)
        Call BoEF(beamA, LoadA, FixedNodes, NIP, resa)
     
        BoEFA = resa
     
    End Function
     
    ' Programs in this module are based on Fortran programs taken from "Prgramming The Finite Element Method"
    ' by IM Smith ans D.V. Griffiths, published by John Wiley & Sons
     
     
    Sub BoEF(beamA As Variant, LoadA As Variant, FixedNodeA As Variant, NIP As Long, resa() As Double)
     
        Dim nels As Long, neq As Long, nn As Long, nband As Long, nr As Long
        Dim iel As Long, i As Long, j As Long, k As Long, l As Long
        Dim loaded_nodes As Long, fixed_nodes As Long, FixedCols As Long, Sense As Long, Value As Double, Node As Long
        Dim fs As Double, fs0 As Double, fs1 As Double, x As Double, samp_pt As Double, EI0 As Double, EI1 As Double, EI As Double
     
        Dim km() As Double, mm() As Double, eld() As Double, kv() As Double, Loads() As Double, coord() As Double
        Dim action As Variant, g_coord() As Double, ftf() As Double, dtd() As Double
        Dim der2() As Double, fun() As Double, mom() As Double, store_km() As Double, points() As Double
        Dim weights() As Double, ell() As Double
        Dim nf() As Long, g() As Long, num() As Long, g_num() As Long, no() As Long, g_g() As Long
        Dim etype() As Long, n As Long
     
        Const nod As Long = 2, nodof As Long = 2, ndof As Long = 4, ndim As Long = 1
        Const element As String = "line"
     
        If TypeName(beamA) = "Range" Then beamA = beamA.Value2
        If TypeName(FixedNodeA) = "Range" Then FixedNodeA = FixedNodeA.Value2
        If TypeName(LoadA) = "Range" Then LoadA = LoadA.Value2
     
        loaded_nodes = UBound(LoadA)
        nels = UBound(beamA) - 1
        If NIP < 2 Then NIP = 2
        If NIP > 6 Then NIP = 6
     
        nn = nels + 1
     
     
        ReDim nf(nodof, nn)
        ReDim km(ndof, ndof)
        ReDim coord(nod, ndim)
        ReDim g_coord(ndim, nn)
        ReDim eld(ndof)
        ReDim g_num(nod, nels)
        ReDim num(nod)
        ReDim g(ndof)
        ReDim g_g(ndof, nels)
        ReDim mm(ndof, ndof)
        ReDim ftf(ndof, ndof)
        ReDim ell(nels)
        ReDim dtd(ndof, ndof)
        ReDim store_km(ndof, ndof, nels)
        ReDim der2(ndof)
        ReDim fun(ndof)
        ReDim mom(nn)
        ReDim points(NIP, ndim)
        ReDim weights(NIP)
        ReDim etype(nels)
     
        fs0 = beamA(1, 3)
        fs1 = beamA(2, 3)
        EI = beamA(1, 2)
     
        For i = 1 To nels
            ell(i) = beamA(i + 1, 1) - beamA(i, 1)
            etype(i) = 1
        Next i
        nr = 0
     
        For i = 1 To nodof
            For j = 1 To nn
                nf(i, j) = 1
            Next j
        Next i
     
        If nr > 0 Then
            'read(10,*)(k,nf(:,k),i=1,nr)
        End If
        Call formnf(nf)
        neq = WorksheetFunction.max(nf)
     
        '!--------------loop the elements to find global array sizes--------------------
        nband = 0
     
        For iel = 1 To nels
            Call geometry_2l(iel, ell(iel), coord, num)
            Call num_to_g(num, nf, g, ndof)
     
            For i = 1 To nod
                g_num(i, iel) = num(i)
            Next i
     
            For j = 1 To nod
                g_coord(1, num(j)) = coord(j, 1)
            Next j
     
            For i = 1 To ndof
                g_g(i, iel) = g(i)
            Next i
     
     
            If nband < bandwidth(g) Then nband = bandwidth(g)
        Next iel
     
        ReDim kv(neq * (nband + 1))
        ReDim Loads(0 To neq)
     
        Call sample(element, points, weights)
     
        x = 0#
     
        For iel = 1 To nels
            ReDim km(ndof, ndof)
            ReDim mm(ndof, ndof)
            For i = 1 To ndof
                g(i) = g_g(i, iel)
            Next i
     
            fs0 = beamA(iel, 3)
            fs1 = beamA(iel + 1, 3)
            EI0 = beamA(iel, 2)
            EI1 = beamA(iel + 1, 2)
            For i = 1 To NIP
                samp_pt = ell(iel) * 0.5 * (points(i, 1) + 1#)
                fs = samp_pt / (ell(iel)) * (fs1 - fs0) + fs0
                EI = samp_pt / (ell(iel)) * (EI1 - EI0) + EI0
                Call fmbeam(der2, fun, points, i, ell(iel))
                For k = 1 To ndof
                    For l = 1 To ndof
                        ftf(k, l) = fun(k) * fun(l) * weights(i) * 0.5 * ell(iel) * fs
                        dtd(k, l) = der2(k) * der2(l) * weights(i) * 8# * EI / (ell(iel) ^ 3)
                        mm(k, l) = mm(k, l) + ftf(k, l)
                        km(k, l) = km(k, l) + dtd(k, l)
     
                    Next l
                Next k
            Next i
     
            For i = 1 To ndof
                For j = 1 To ndof
                    km(i, j) = km(i, j) + mm(i, j)
                    store_km(i, j, iel) = km(i, j)
                Next j
            Next i
     
            x = beamA(iel, 1)
            Call formkv(kv, km, g, neq, iel)
        Next iel
     
        '!-----------------------------read loads---------------------------------------
     
        If loaded_nodes <> 0 Then
     
            For i = 1 To loaded_nodes
                k = LoadA(i, 1)
                For j = 1 To nodof
                    Loads(nf(j, k)) = LoadA(i, j + 1)
                Next j
            Next i
        End If
     
        'read (10,*)fixed_nodes
        If IsArray(FixedNodeA) = True Then
        fixed_nodes = UBound(FixedNodeA)
        FixedCols = UBound(FixedNodeA, 2)
        Else
        fixed_nodes = 0
        End If
     
        If fixed_nodes <> 0 Then
              For i = 1 To fixed_nodes
              Node = FixedNodeA(i, 1)
              Sense = 0
              If FixedNodeA(i, 2) = 1 Then Sense = 1
              If FixedNodeA(i, 3) = 1 Then Sense = Sense + 2
              If Sense = 1 Or Sense = 3 Then
              If FixedCols > 3 Then Value = FixedNodeA(i, 4) Else Value = 0
              n = nf(1, Node)
              kv(n) = kv(n) + 1E+20
              Loads(n) = kv(n) * Value
              Sense = Sense - 1
              End If
              If Sense = 2 Then
               If FixedCols > 4 Then Value = FixedNodeA(i, 5) Else Value = 0
              n = nf(2, Node)
              kv(n) = kv(n) + 1E+20
              Loads(n) = kv(n) * Value
              End If
              Next i
     
        End If
     
        '!-----------------------------equation solution -------------------------------
        Call banred(kv, neq)
        Call bacsub(kv, Loads)
     
        '!-----------------------retrieve element end actions---------------------------
        For iel = 1 To nels
            For i = 1 To ndof
                For j = 1 To ndof
                    km(i, j) = store_km(i, j, iel)
                Next j
                g(i) = g_g(i, iel)
                eld(i) = Loads(g(i))
            Next i
     
     
            action = MATMUL(km, eld)
            resa(iel, 3) = -action(2)
            resa(iel, 4) = action(1)
     
        Next iel
        resa(iel, 3) = -action(4)
        resa(iel, 4) = -action(3)
        For i = 1 To nn
            resa(i, 1) = beamA(i, 1)
            resa(i, 2) = Loads(2 * i - 1)
        Next i
     
    End Sub
     
     
    Sub formnf(nf)
     
    '  ! reform nf
        Dim i As Long, j As Long, M As Long
        M = 0
        For j = 1 To UBound(nf, 2)
            For i = 1 To UBound(nf, 1)
                If (nf(i, j) <> 0) Then
                    M = M + 1
                    nf(i, j) = M
                End If
            Next i
        Next j
     
    End Sub
     
     
    '!----------------Node to freedom number conversion ----------------------------
    Sub num_to_g(num, nf, g, ndof)
     
    ' !finds the g vector from num and nf
        Dim i As Long, j As Long, k As Long, nod As Long, nodof As Long
     
        nod = UBound(num, 1)
        nodof = UBound(nf, 1)
        For i = 1 To nod
            k = i * nodof
            For j = 1 To nodof
                g(k - nodof + j) = nf(j, num(i))
            Next j
        Next i
     
    End Sub
     
     
    ' !-------------------------------- Lines  --------------------------------------
     
    Sub geometry_2l(iel, ell, coord, num)
     
    ' ! node numbers, nodal coordinates and steering vectors for
    ' ! a line of (nonuniform) beam elements
        num(1) = iel
        num(2) = iel + 1
        If (iel = 1) Then
            coord(1, 1) = 0#
            coord(2, 1) = ell
        Else
            coord(1, 1) = coord(2, 1)
            coord(2, 1) = coord(2, 1) + ell
        End If
    End Sub
     
    Sub rod_km(km, ea, Length)
     
        Dim Stiff As Double
        Stiff = ea / Length
        km(1, 1) = Stiff
        km(2, 2) = Stiff
        km(1, 2) = -Stiff
        km(2, 1) = -Stiff
     
    End Sub
     
    Sub formkv(bk, km, g, n, iel)
     
    ' !global stiffness matrix stored as a vector (upper triangle)
        Dim idof As Long, i As Long, j As Long, icd As Long, ival As Long
        idof = UBound(km, 1) - LBound(km, 1) + 1
     
        For i = 1 To idof
            If (g(i) <> 0) Then
                For j = 1 To idof
                    If (g(j) <> 0) Then
                        icd = g(j) - g(i) + 1
                        If (icd - 1 >= 0) Then
                            ival = n * (icd - 1) + g(i)
                            bk(ival) = bk(ival) + km(i, j)
                        End If
                    End If
                Next j
     
            End If
     
        Next i
     
     
        ' Return
    End Sub
     
     
    Sub banred(bk, n)
     
    ' ! gaussian reduction on a vector stored as an upper triangle
        Dim i As Long, il1 As Long, kbl As Long, j As Long, ij As Long, nkb As Long, M As Long, ni As Long, nj As Long, iw As Long
        Dim sum As Double
     
        iw = UBound(bk, 1) / n - 1
        For i = 2 To n
            il1 = i - 1
            kbl = il1 + iw + 1
            If (kbl - n > 0) Then kbl = n
            For j = i To kbl
                ij = (j - i) * n + i
                sum = bk(ij)
                nkb = j - iw
                If (nkb <= 0) Then nkb = 1
                If (nkb - il1 <= 0) Then
                    For M = nkb To il1
                        ni = (i - M) * n + M
                        nj = (j - M) * n + M
                        sum = sum - bk(ni) * bk(nj) / bk(M)
                    Next M
                End If
                bk(ij) = sum
            Next j
        Next i
     
    End Sub
     
    Sub bacsub(bk, Loads)
     
     
    ' ! performs the complete gaussian backsubstitution
        Dim nkb As Long, k As Long, i As Long, jn As Long, jj As Long, i1 As Long, n As Long, iw As Long
        Dim sum As Double
        n = UBound(Loads, 1)
        iw = UBound(bk, 1) / n - 1
        Loads(1) = Loads(1) / bk(1)
        For i = 2 To n
            sum = Loads(i)
            i1 = i - 1
            nkb = i - iw
            If (nkb <= 0) Then nkb = 1
            For k = nkb To i1
                jn = (i - k) * n + k
                sum = sum - bk(jn) * Loads(k)
            Next k
            Loads(i) = sum / bk(i)
        Next i
        For jj = 2 To n
            i = n - jj + 1
            sum = 0#
            i1 = i + 1
            nkb = i + iw
            If (nkb - n > 0) Then nkb = n
            For k = i1 To nkb
                jn = (k - i) * n + i
                sum = sum + bk(jn) * Loads(k)
            Next k
            Loads(i) = Loads(i) - sum / bk(i)
        Next jj
        ' Return
    End Sub
     
     
    Function MATMUL(A, b)
    '
    '      PRODUCT OF TWO MATRICES
    '
        Dim i As Long, j As Long, k As Long, Mat, x As Double, l As Long, n As Long
        l = UBound(A, 2)
        n = UBound(A)
        ReDim Mat(l) As Double
        For i = 1 To l
            For j = 1 To n
                x = 0#
                For k = 1 To l
                    x = x + A(i, k) * b(k)
                Next k
                Mat(i) = x
            Next j
        Next i
        MATMUL = Mat
    End Function
     
     
    '
    Sub MATRAN(A, IA, b, IB, M, n)
    '
    '      FORMS THE TRANSPOSE OF A MATRIX
    '
    '     REAL A(IA,*),B(IB,*)
    '     DO 1 I = 1,M
    '         DO 1 J = 1,N
    '     A(j, i) = B(i, j)
     
    End Sub
     
    Sub MSMULT(A, IA, c, M, n)
    '
    '      MULTIPLIES A MATRIX BY A SCALAR
    '
    '      REAL A(IA,*)
    '      DO 1 I = 1,M
    '          DO 1 J = 1,N
    '     A(i, j) = A(i, j) * c
    End Sub
     
    Sub MVMULT(M, v, k, l, y)
     
    '      MULTIPLIES A MATRIX BY A VECTOR
     
     
        Dim i As Long, j As Long, x As Double
        ' ReDim Preserve M(k, l), V(l), Y(k)
        ReDim y(k)
        For i = 1 To k
            x = 0#
            For j = 1 To l
                x = x + M(i, j) * v(j)
            Next j
            y(i) = x
        Next i
    End Sub
     
    Function bandwidth(g)
     
     
    '! finds the element bandwidth from g
        If WorksheetFunction.min(g) > 0 Then
            bandwidth = WorksheetFunction.max(g) - WorksheetFunction.min(g)
        Else
            bandwidth = 0
        End If
    End Function
     
     
    Sub sample(element, s, wt)
     
    ' ! returns the local coordinates of the integrating points
     
    ' ****  Edited for line element only  *****
     
        Dim NIP As Long, root3 As Double, R15 As Double, w(1 To 3) As Double, v(1 To 9) As Double, b As Double, c As Double, i As Long
     
        root3 = 1# / (3#) ^ 0.5
        R15 = 0.2 * (15#) ^ 0.5
        NIP = UBound(s, 1)
        w(1) = 5# / 9#
        w(2) = 8# / 9#
        w(3) = 5# / 9#
        For i = 1 To 7 Step 3
            v(i) = 5# / 9# * w(1)
            v(i + 1) = 8# / 9# * w(2)
            v(i + 2) = 5# / 9# * w(3)
        Next i
        Select Case element
        Case "line"
     
            Select Case (NIP)
            Case (1)
                s(1, 1) = 0#
                wt(1) = 2#
            Case (2)
                s(1, 1) = root3
                s(2, 1) = -s(1, 1)
                wt(1) = 1#
                wt(2) = 1#
            Case (3)
                s(1, 1) = R15
                s(2, 1) = 0#
                s(3, 1) = -s(1, 1)
                wt = w
            Case (4)
                s(1, 1) = 0.861136311594053
                s(2, 1) = 0.339981043584856
                s(3, 1) = -s(2, 1)
                s(4, 1) = -s(1, 1)
     
                wt(1) = 0.347854845137454
                wt(2) = 0.652145154862546
                wt(3) = wt(2)
                wt(4) = wt(1)
            Case (5)
                s(1, 1) = 0.906179845938664
                s(2, 1) = 0.538469310105683
                s(3, 1) = 0#
                s(4, 1) = -s(2, 1)
                s(5, 1) = -s(1, 1)
                wt(1) = 0.236926885056189
                wt(2) = 0.478628670499366
                wt(3) = 0.568888888888889
                wt(4) = wt(2)
                wt(5) = wt(1)
            Case (6)
                s(1, 1) = 0.932469514203152
                s(2, 1) = 0.661209386466265
                s(3, 1) = 0.238619186083197
                s(4, 1) = -s(3, 1)
                s(5, 1) = -s(2, 1)
                s(6, 1) = -s(1, 1)
                wt(1) = 0.17132449237917
                wt(2) = 0.360761573048139
                wt(3) = 0.467913934572691
                wt(4) = wt(3)
                wt(5) = wt(2)
                wt(6) = wt(1)
            Case Default
                '      print*,"wrong number of integrating points for a line"
            End Select
        End Select
     
    End Sub
     
     
    Sub fmbeam(der2, fun, points, i, ell)
     
     
    ' !
    ' ! this subroutine forms the beam shape functions
    ' ! and their 2nd derivatives in local coordinates
    ' !
     
        Dim xi As Double, xi2 As Double, xi3
     
        xi = points(i, 1)
        xi2 = xi * xi
        xi3 = xi2 * xi
     
        fun(1) = 0.25 * (xi3 - 3# * xi + 2#)
        fun(2) = 0.125 * ell * (xi3 - xi2 - xi + 1#)
        fun(3) = 0.25 * (-xi3 + 3# * xi + 2#)
        fun(4) = 0.125 * ell * (xi3 + xi2 - xi - 1#)
        der2(1) = 1.5 * xi
        der2(2) = 0.25 * ell * (3# * xi - 1#)
        der2(3) = -1.5 * xi
        der2(4) = 0.25 * ell * (3# * xi + 1#)
     
     
    End Sub






    Code : Sélectionner tout - Visualiser dans une fenêtre à part
    1
    2
    3
    4
    5
    6
    7
    8
    9
    10
    11
    12
    13
    14
    15
    16
    17
    18
    19
    20
    21
    22
    23
    24
    25
    26
    27
    28
    29
    30
    31
    32
    33
    34
    35
    36
    37
    38
    39
    40
    41
    42
    43
    44
    45
    46
    47
    48
    49
    50
    51
    52
    53
    54
    55
    56
    57
    58
    59
    60
    61
    62
    63
    64
    65
    66
    67
    68
    69
    70
    71
    72
    73
    74
    75
    76
    77
    78
    79
    80
    81
    82
    83
    84
    85
    86
    87
    88
    89
    90
    91
    92
    93
    94
    95
    96
    97
    98
    99
    100
    101
    102
    103
    104
    105
    106
    107
    108
    109
    110
    111
    112
    113
    114
    115
    116
    117
    118
    119
    120
    121
    122
    123
    124
    125
    126
    127
    128
    129
    130
    131
    132
    133
    134
    135
    136
    137
    138
    139
    140
    141
    142
    143
    144
    145
    146
    147
    148
    149
    150
    151
    152
    153
    154
    155
    156
    157
    158
    159
    160
    161
    162
    163
    164
    165
    166
    167
    168
    169
    170
    171
    172
    173
    174
    175
    176
    177
    178
    179
    180
    181
    182
    183
    184
    185
    186
    187
    188
    189
    190
    191
    192
    193
    194
    195
    196
    197
    198
    199
    200
    201
    202
    203
    204
    205
    206
    207
    208
    209
    210
    211
    212
    213
    214
    215
    216
    217
    218
    219
    220
    221
    222
    223
    224
    225
    226
    227
    228
    229
    230
    231
    232
    233
    234
    235
    236
    237
    238
    239
    240
    241
    242
    243
    244
    245
    246
    247
    248
    249
    250
    251
    252
    253
    254
    255
    256
    257
    258
    259
    260
    261
    262
    263
    264
    265
    266
    267
    268
    269
    270
    271
    272
    273
    274
    275
    276
    277
    278
    279
    280
    281
    282
    283
    284
    285
    286
    287
    288
    289
    290
    291
    292
    293
    294
    295
    296
    297
    298
    299
    300
    301
    302
    303
    304
    305
    306
    307
    308
    309
    310
    311
    312
    313
    314
    315
    316
    317
    318
    319
    320
    321
    322
    323
    324
    325
    326
    327
    328
    329
    330
    331
    332
    333
    334
    335
    336
    337
    338
    339
    340
    341
    342
    343
    344
    345
    346
    347
    348
    349
    350
    351
    352
    353
    354
    355
    356
    357
    358
    359
    360
    361
    362
    363
    364
    365
    366
    367
    368
    369
    370
    371
    372
    373
    374
    375
    376
    377
    378
    379
    380
    381
    382
    383
    384
    385
    386
    387
    388
    389
    390
    391
    392
    393
    394
    395
    396
    397
    398
    399
    400
    401
    402
    403
    404
    405
    406
    407
    408
    409
    410
    411
    412
    413
    414
    415
    416
    417
    418
    419
    420
    421
    422
    423
    424
    425
    426
    427
    428
    429
    430
    431
    432
    433
    434
    435
    436
    437
    438
    439
    440
    441
    442
    443
    444
    445
    446
    447
    448
    449
    450
    451
    452
    453
    454
    455
    456
    457
    458
    459
    460
    461
    462
    463
    464
    465
    466
    467
    468
    469
    470
    471
    472
    473
    474
    475
    476
    477
    478
    479
    480
    481
    482
    483
    484
    485
    486
    487
    488
    489
    490
    491
    492
    493
    494
    495
    496
    497
    498
    499
    500
    501
    502
    503
    504
    505
    506
    507
    508
    509
    510
    511
    512
    513
    514
    515
    516
    517
    518
    519
    520
    521
    522
    523
    524
    525
    526
    527
    528
    529
    530
    531
    532
    533
    534
    535
    536
    537
    538
    539
    540
    541
    542
    543
    544
    545
    546
    547
    548
    549
    550
    551
    552
    553
    554
    555
    556
    557
    558
    559
    560
    561
    562
    563
    564
    565
    566
    567
    568
    569
    570
    571
    572
    573
    574
    575
    576
    577
    578
    579
    580
    581
    582
    583
    584
    585
    586
    587
    588
    589
    590
    591
    592
    593
    594
    595
    596
    597
    598
    599
    600
    601
    602
    603
    604
    605
    606
    607
    608
    609
    610
    611
    612
    613
    614
    615
    616
    617
    618
    619
    620
    621
    622
    623
    624
    625
    626
    627
    628
    629
    630
    631
    632
    633
    634
    635
    636
    637
    638
    639
    640
    641
    642
    643
    644
    645
    646
    647
    648
    649
    650
    651
    652
    653
    654
    655
    656
    657
    658
    659
    660
    661
    662
    663
    664
    665
    666
        ' Programs in this module are based on Fortran programs taken from "Prgramming The Finite Element Method"
        ' by IM Smith ans D.V. Griffiths, published by John Wiley & Sons
     
     
        Sub BoEF(beamA As Array, LoadA As Array, FixedNodeA As Object, NIP As Long, ByRef resa(,) As Double)
     
            Dim nels As Long, neq As Long, nn As Long, nband As Long, nr As Long
            Dim iel As Long, i As Long, j As Long, k As Long, l As Long
            Dim loaded_nodes As Long, fixed_nodes As Long, FixedCols As Long, Sense As Long, Value As Double, Node As Long
            Dim fs As Double, fs0 As Double, fs1 As Double, x As Double, samp_pt As Double, EI0 As Double, EI1 As Double, EI As Double
     
            Dim Loads() As Double
            Dim action As Object
            Dim no() As Long
            Dim n As Long
     
            Const nod As Long = 2, nodof As Long = 2, ndof As Long = 4, ndim As Long = 1
            Const element As String = "line"
     
     
            loaded_nodes = UBound(LoadA)
            nels = UBound(beamA) - 1
            If NIP < 2 Then NIP = 2
            If NIP > 6 Then NIP = 6
     
            nn = nels + 1
     
            Dim nf As Array = Array.CreateInstance(GetType(Long), {nodof, nn}, {1, 1})
            Dim km As Array = Array.CreateInstance(GetType(Double), {ndof, ndof}, {1, 1})
            Dim coord As Array = Array.CreateInstance(GetType(Double), {nod, ndim}, {1, 1})
            Dim g_coord As Array = Array.CreateInstance(GetType(Double), {ndim, nn}, {1, 1})
            Dim eld As Array = Array.CreateInstance(GetType(Double), {ndof}, {1})
            Dim g_num As Array = Array.CreateInstance(GetType(Long), {nod, nels}, {1, 1})
            Dim num As Array = Array.CreateInstance(GetType(Long), {nod}, {1})
            Dim g As Array = Array.CreateInstance(GetType(Long), {ndof}, {1})
            Dim g_g As Array = Array.CreateInstance(GetType(Long), {ndof, nels}, {1, 1})
            Dim mm As Array = Array.CreateInstance(GetType(Double), {ndof, ndof}, {1, 1})
            Dim ftf As Array = Array.CreateInstance(GetType(Double), {ndof, ndof}, {1, 1})
            Dim ell As Array = Array.CreateInstance(GetType(Double), {nels}, {1})
            Dim dtd As Array = Array.CreateInstance(GetType(Double), {ndof, ndof}, {1, 1})
            Dim store_km As Array = Array.CreateInstance(GetType(Double), {ndof, ndof, nels}, {1, 1, 1})
            Dim der2 As Array = Array.CreateInstance(GetType(Double), {ndof}, {1})
            Dim fun As Array = Array.CreateInstance(GetType(Double), {ndof}, {1})
            Dim mom As Array = Array.CreateInstance(GetType(Double), {nn}, {1})
            Dim points As Array = Array.CreateInstance(GetType(Double), {NIP, ndim}, {1, 1})
            Dim weights As Array = Array.CreateInstance(GetType(Double), {NIP}, {1})
            Dim etype As Array = Array.CreateInstance(GetType(Long), {nels}, {1})
     
     
     
            fs0 = beamA(1, 3)
            fs1 = beamA(2, 3)
            EI = beamA(1, 2)
     
            For i = 1 To nels
                ell(i) = beamA(i + 1, 1) - beamA(i, 1)
                etype(i) = 1
            Next i
            nr = 0
     
            For i = 1 To nodof
                For j = 1 To nn
                    nf(i, j) = 1
                Next j
            Next i
     
            If nr > 0 Then
                'read(10,*)(k,nf(:,k),i=1,nr)
            End If
            Call formnf(nf)
            neq = WorksheetMax(nf)
     
            '!--------------loop the elements to find global array sizes--------------------
            nband = 0
     
            For iel = 1 To nels
                Call geometry_2l(iel, ell(iel), coord, num)
                Call num_to_g(num, nf, g, ndof)
     
                For i = 1 To nod
                    g_num(i, iel) = num(i)
                Next i
     
                For j = 1 To nod
                    g_coord(1, num(j)) = coord(j, 1)
                Next j
     
                For i = 1 To ndof
                    g_g(i, iel) = g(i)
                Next i
     
     
                If nband < bandwidth(g) Then nband = bandwidth(g)
            Next iel
     
            Dim kv As Array = Array.CreateInstance(GetType(Double), {neq * (nband + 1)}, {1})
            Loads = Array.CreateInstance(GetType(Double), {neq + 1}, {0})
     
            Call sample(element, points, weights)
     
            x = 0#
     
            For iel = 1 To nels
                km = Array.CreateInstance(GetType(Double), {ndof, ndof}, {1, 1})
                mm = Array.CreateInstance(GetType(Double), {ndof, ndof}, {1, 1})
     
                For i = 1 To ndof
                        g(i) = g_g(i, iel)
                    Next i
     
                    fs0 = beamA(iel, 3)
                    fs1 = beamA(iel + 1, 3)
                    EI0 = beamA(iel, 2)
                    EI1 = beamA(iel + 1, 2)
                    For i = 1 To NIP
                        samp_pt = ell(iel) * 0.5 * (points(i, 1) + 1.0#)
                        fs = samp_pt / (ell(iel)) * (fs1 - fs0) + fs0
                        EI = samp_pt / (ell(iel)) * (EI1 - EI0) + EI0
                        Call fmbeam(der2, fun, points, i, ell(iel))
                        For k = 1 To ndof
                            For l = 1 To ndof
                                ftf(k, l) = fun(k) * fun(l) * weights(i) * 0.5 * ell(iel) * fs
                                dtd(k, l) = der2(k) * der2(l) * weights(i) * 8.0# * EI / (ell(iel) ^ 3)
                                mm(k, l) = mm(k, l) + ftf(k, l)
                                km(k, l) = km(k, l) + dtd(k, l)
     
                            Next l
                        Next k
                    Next i
     
                    For i = 1 To ndof
                        For j = 1 To ndof
                            km(i, j) = km(i, j) + mm(i, j)
                            store_km(i, j, iel) = km(i, j)
                        Next j
                    Next i
     
                    x = beamA(iel, 1)
                Call formkv(kv, km, g, neq, iel)
     
            Next iel
     
            '!-----------------------------read loads---------------------------------------
     
            If loaded_nodes <> 0 Then
                '
                For i = 1 To loaded_nodes
                    k = LoadA(i, 1)
                    For j = 1 To nodof
                        Loads(nf(j, k)) = LoadA(i, j + 1)
                    Next j
                Next i
            End If
     
            'read (10,*)fixed_nodes
            If IsArray(FixedNodeA) = True Then
                fixed_nodes = UBound(FixedNodeA)
                FixedCols = UBound(FixedNodeA, 2)
            Else
                fixed_nodes = 0
            End If
            If fixed_nodes <> 0 Then
                For i = 1 To fixed_nodes
                    Node = FixedNodeA(i, 1)
                    Sense = 0
                    If FixedNodeA(i, 2) = 1 Then Sense = 1
                    If FixedNodeA(i, 3) = 1 Then Sense = Sense + 2
                    If Sense = 1 Or Sense = 3 Then
                        If FixedCols > 3 Then Value = FixedNodeA(i, 4) Else Value = 0
                        n = nf(1, Node)
                        kv(n) = kv(n) + 1.0E+20
                        Loads(n) = kv(n) * Value
                        Sense = Sense - 1
                    End If
                    If Sense = 2 Then
                        If FixedCols > 4 Then Value = FixedNodeA(i, 5) Else Value = 0
                        n = nf(2, Node)
                        kv(n) = kv(n) + 1.0E+20
                        Loads(n) = kv(n) * Value
                    End If
                Next i
     
            End If
     
            '!-----------------------------equation solution -------------------------------
            Call banred(kv, neq)
            Call bacsub(kv, Loads)
     
            '!-----------------------retrieve element end actions---------------------------
            For iel = 1 To nels
                For i = 1 To ndof
                    For j = 1 To ndof
                        km(i, j) = store_km(i, j, iel)
                    Next j
                    g(i) = g_g(i, iel)
                    eld(i) = Loads(g(i))
                Next i
     
     
                action = MATMUL(km, eld)
                resa(iel, 3) = -action(2)
                resa(iel, 4) = action(1)
     
            Next iel
            resa(iel, 3) = -action(4)
            resa(iel, 4) = -action(3)
            For i = 1 To nn
                resa(i, 1) = beamA(i, 1)
                resa(i, 2) = Loads(2 * i - 1)
            Next i
     
     
        End Sub
     
        Sub formnf(nf)
     
            Static TimeofFunction
            Dim Stopwatch1 As Stopwatch = Stopwatch.StartNew()
     
            '  ! reform nf
            Dim i As Long, j As Long, M As Long
            M = 0
            For j = 1 To UBound(nf, 2)
                For i = 1 To UBound(nf, 1)
                    If (nf(i, j) <> 0) Then
                        M = M + 1
                        nf(i, j) = M
                    End If
                Next i
            Next j
     
            Stopwatch1.Stop()
            TimeofFunction = TimeofFunction + Stopwatch1.Elapsed
     
     
        End Sub
     
     
        '!----------------Node to freedom number conversion ----------------------------
        Sub num_to_g(num, nf, g, ndof)
     
            Static TimeofFunction
            Dim Stopwatch1 As Stopwatch = Stopwatch.StartNew()
     
            ' !finds the g vector from num and nf
            Dim i As Long, j As Long, k As Long, nod As Long, nodof As Long
     
            nod = UBound(num, 1)
            nodof = UBound(nf, 1)
            For i = 1 To nod
                k = i * nodof
                For j = 1 To nodof
                    g(k - nodof + j) = nf(j, num(i))
                Next j
            Next i
            Stopwatch1.Stop()
            TimeofFunction = TimeofFunction + Stopwatch1.Elapsed
     
     
        End Sub
     
     
        ' !-------------------------------- Lines  --------------------------------------
     
        Sub geometry_2l(iel, ell, coord, num)
     
            Static TimeofFunction
            Dim Stopwatch1 As Stopwatch = Stopwatch.StartNew()
     
            ' ! node numbers, nodal coordinates and steering vectors for
            ' ! a line of (nonuniform) beam elements
            num(1) = iel
            num(2) = iel + 1
            If (iel = 1) Then
                coord(1, 1) = 0#
                coord(2, 1) = ell
            Else
                coord(1, 1) = coord(2, 1)
                coord(2, 1) = coord(2, 1) + ell
            End If
     
            Stopwatch1.Stop()
            TimeofFunction = TimeofFunction + Stopwatch1.Elapsed
     
        End Sub
     
        Sub rod_km(km, ea, Length)
     
            Static TimeofFunction
            Dim Stopwatch1 As Stopwatch = Stopwatch.StartNew()
     
            Dim Stiff As Double
            Stiff = ea / Length
            km(1, 1) = Stiff
            km(2, 2) = Stiff
            km(1, 2) = -Stiff
            km(2, 1) = -Stiff
     
            Stopwatch1.Stop()
            TimeofFunction = TimeofFunction + Stopwatch1.Elapsed
        End Sub
     
        Sub formkv(bk, km, g, n, iel)
     
            Static TimeofFunction
            Dim Stopwatch1 As Stopwatch = Stopwatch.StartNew()
     
            ' !global stiffness matrix stored as a vector (upper triangle)
            Dim idof As Long, i As Long, j As Long, icd As Long, ival As Long
            idof = UBound(km, 1) - LBound(km, 1) + 1
     
            For i = 1 To idof
                If (g(i) <> 0) Then
                    For j = 1 To idof
                        If (g(j) <> 0) Then
                            icd = g(j) - g(i) + 1
                            If (icd - 1 >= 0) Then
                                ival = n * (icd - 1) + g(i)
                                bk(ival) = bk(ival) + km(i, j)
                            End If
                        End If
                    Next j
     
                End If
     
            Next i
     
            Stopwatch1.Stop()
            TimeofFunction = TimeofFunction + Stopwatch1.Elapsed
     
            ' Return
        End Sub
     
     
        Sub banred(bk, n)
     
            Static TimeofFunction
            Dim Stopwatch1 As Stopwatch = Stopwatch.StartNew()
     
            ' ! gaussian reduction on a vector stored as an upper triangle
            Dim i As Long, il1 As Long, kbl As Long, j As Long, ij As Long, nkb As Long, M As Long, ni As Long, nj As Long, iw As Long
            Dim sum As Double
     
            iw = UBound(bk, 1) / n - 1
            For i = 2 To n
                il1 = i - 1
                kbl = il1 + iw + 1
                If (kbl - n > 0) Then kbl = n
                For j = i To kbl
                    ij = (j - i) * n + i
                    sum = bk(ij)
                    nkb = j - iw
                    If (nkb <= 0) Then nkb = 1
                    If (nkb - il1 <= 0) Then
                        For M = nkb To il1
                            ni = (i - M) * n + M
                            nj = (j - M) * n + M
                            sum = sum - bk(ni) * bk(nj) / bk(M)
                        Next M
                    End If
                    bk(ij) = sum
                Next j
            Next i
     
            Stopwatch1.Stop()
            TimeofFunction = TimeofFunction + Stopwatch1.Elapsed
     
        End Sub
     
        Sub bacsub(bk, Loads)
            Static TimeofFunction
            Dim Stopwatch1 As Stopwatch = Stopwatch.StartNew()
     
            ' ! performs the complete gaussian backsubstitution
            Dim nkb As Long, k As Long, i As Long, jn As Long, jj As Long, i1 As Long, n As Long, iw As Long
            Dim sum As Double
            n = UBound(Loads, 1)
            iw = UBound(bk, 1) / n - 1
            Loads(1) = Loads(1) / bk(1)
            For i = 2 To n
                sum = Loads(i)
                i1 = i - 1
                nkb = i - iw
                If (nkb <= 0) Then nkb = 1
                For k = nkb To i1
                    jn = (i - k) * n + k
                    sum = sum - bk(jn) * Loads(k)
                Next k
                Loads(i) = sum / bk(i)
            Next i
            For jj = 2 To n
                i = n - jj + 1
                sum = 0#
                i1 = i + 1
                nkb = i + iw
                If (nkb - n > 0) Then nkb = n
                For k = i1 To nkb
                    jn = (k - i) * n + i
                    sum = sum + bk(jn) * Loads(k)
                Next k
                Loads(i) = Loads(i) - sum / bk(i)
            Next jj
            ' Return
     
            Stopwatch1.Stop()
            TimeofFunction = TimeofFunction + Stopwatch1.Elapsed
     
        End Sub
     
     
        Function MATMUL(A, b)
            Static TimeofFunction
            Dim Stopwatch1 As Stopwatch = Stopwatch.StartNew()
     
            '      PRODUCT OF TWO MATRICES
            '
            Dim i As Long, j As Long, k As Long, x As Double, l As Long, n As Long
            l = UBound(A, 2)
            n = UBound(A)
            Dim Mat = Array.CreateInstance(GetType(Double), {l}, {1})
            For i = 1 To l
                For j = 1 To n
                    x = 0#
                    For k = 1 To l
                        x = x + A(i, k) * b(k)
                    Next k
                    Mat(i) = x
                Next j
            Next i
            MATMUL = Mat
     
            Stopwatch1.Stop()
            TimeofFunction = TimeofFunction + Stopwatch1.Elapsed
     
        End Function
     
     
        '
        Sub MATRAN(A, IA, b, IB, M, n)
            '
            '      FORMS THE TRANSPOSE OF A MATRIX
            '
            '     REAL A(IA,*),B(IB,*)
            '     DO 1 I = 1,M
            '         DO 1 J = 1,N
            '     A(j, i) = B(i, j)
     
        End Sub
     
        Sub MSMULT(A, IA, c, M, n)
            '
            '      MULTIPLIES A MATRIX BY A SCALAR
            '
            '      REAL A(IA,*)
            '      DO 1 I = 1,M
            '          DO 1 J = 1,N
            '     A(i, j) = A(i, j) * c
        End Sub
     
        Sub MVMULT(M, v, k, l, y)
     
            Static TimeofFunction
            Dim Stopwatch1 As Stopwatch = Stopwatch.StartNew()
     
            '      MULTIPLIES A MATRIX BY A VECTOR
     
     
            Dim i As Long, j As Long, x As Double
            ' ReDim Preserve M(k, l), V(l), Y(k)
            ReDim y(k)
            For i = 1 To k
                x = 0#
                For j = 1 To l
                    x = x + M(i, j) * v(j)
                Next j
                y(i) = x
            Next i
     
            Stopwatch1.Stop()
            TimeofFunction = TimeofFunction + Stopwatch1.Elapsed
     
        End Sub
     
        Function bandwidth(g)
     
            Static TimeofFunction
            Dim Stopwatch1 As Stopwatch = Stopwatch.StartNew()
     
            '! finds the element bandwidth from g
            If WorksheetMin(g) > 0 Then
                bandwidth = WorksheetMax(g) - WorksheetMin(g)
            Else
                bandwidth = 0
            End If
     
            Stopwatch1.Stop()
            TimeofFunction = TimeofFunction + Stopwatch1.Elapsed
     
        End Function
     
     
        Sub sample(element, s, wt)
     
            Static TimeofFunction
            Dim Stopwatch1 As Stopwatch = Stopwatch.StartNew()
            ' ! returns the local coordinates of the integrating points
     
            ' ****  Edited for line element only  *****
     
            Dim NIP As Long, root3 As Double, R15 As Double
            Dim w = Array.CreateInstance(GetType(Double), {3}, {1})
            Dim v = Array.CreateInstance(GetType(Double), {9}, {1})
            Dim b As Double, c As Double, i As Long
     
            root3 = 1.0# / (3.0#) ^ 0.5
            R15 = 0.2 * (15.0#) ^ 0.5
            NIP = UBound(s, 1)
            w(1) = 5.0# / 9.0#
            w(2) = 8.0# / 9.0#
            w(3) = 5.0# / 9.0#
            For i = 1 To 7 Step 3
                v(i) = 5.0# / 9.0# * w(1)
                v(i + 1) = 8.0# / 9.0# * w(2)
                v(i + 2) = 5.0# / 9.0# * w(3)
            Next i
            Select Case element
                Case "line"
     
                    Select Case (NIP)
                        Case (1)
                            s(1, 1) = 0#
                            wt(1) = 2.0#
                        Case (2)
                            s(1, 1) = root3
                            s(2, 1) = -s(1, 1)
                            wt(1) = 1.0#
                            wt(2) = 1.0#
                        Case (3)
                            s(1, 1) = R15
                            s(2, 1) = 0#
                            s(3, 1) = -s(1, 1)
                            wt = w
                        Case (4)
                            s(1, 1) = 0.861136311594053
                            s(2, 1) = 0.339981043584856
                            s(3, 1) = -s(2, 1)
                            s(4, 1) = -s(1, 1)
     
                            wt(1) = 0.347854845137454
                            wt(2) = 0.652145154862546
                            wt(3) = wt(2)
                            wt(4) = wt(1)
                        Case (5)
                            s(1, 1) = 0.906179845938664
                            s(2, 1) = 0.538469310105683
                            s(3, 1) = 0#
                            s(4, 1) = -s(2, 1)
                            s(5, 1) = -s(1, 1)
                            wt(1) = 0.236926885056189
                            wt(2) = 0.478628670499366
                            wt(3) = 0.568888888888889
                            wt(4) = wt(2)
                            wt(5) = wt(1)
                        Case (6)
                            s(1, 1) = 0.932469514203152
                            s(2, 1) = 0.661209386466265
                            s(3, 1) = 0.238619186083197
                            s(4, 1) = -s(3, 1)
                            s(5, 1) = -s(2, 1)
                            s(6, 1) = -s(1, 1)
                            wt(1) = 0.17132449237917
                            wt(2) = 0.360761573048139
                            wt(3) = 0.467913934572691
                            wt(4) = wt(3)
                            wt(5) = wt(2)
                            wt(6) = wt(1)
     
                    End Select
            End Select
            Stopwatch1.Stop()
            TimeofFunction = TimeofFunction + Stopwatch1.Elapsed
        End Sub
     
     
        Sub fmbeam(der2, fun, points, i, ell)
     
            Static TimeofFunction
            Dim Stopwatch1 As Stopwatch = Stopwatch.StartNew()
     
            ' !
            ' ! this subroutine forms the beam shape functions
            ' ! and their 2nd derivatives in local coordinates
            ' !
     
            Dim xi As Double, xi2 As Double, xi3
     
            xi = points(i, 1)
            xi2 = xi * xi
            xi3 = xi2 * xi
     
            fun(1) = 0.25 * (xi3 - 3.0# * xi + 2.0#)
            fun(2) = 0.125 * ell * (xi3 - xi2 - xi + 1.0#)
            fun(3) = 0.25 * (-xi3 + 3.0# * xi + 2.0#)
            fun(4) = 0.125 * ell * (xi3 + xi2 - xi - 1.0#)
            der2(1) = 1.5 * xi
            der2(2) = 0.25 * ell * (3.0# * xi - 1.0#)
            der2(3) = -1.5 * xi
            der2(4) = 0.25 * ell * (3.0# * xi + 1.0#)
     
            Stopwatch1.Stop()
            TimeofFunction = TimeofFunction + Stopwatch1.Elapsed
     
        End Sub
     
     
     
        Function WorksheetMax(tr As Object, Optional InitialRowNumber As Integer = 0, Optional InitialColumnNumber As Integer = 0)
     
            Static TimeofFunction
            Dim Stopwatch1 As Stopwatch = Stopwatch.StartNew()
     
            WorksheetMax = Nothing
     
            For Each element In tr
     
                If IsNumeric(element) Then
     
                    If IsNothing(WorksheetMax) Then
                        WorksheetMax = element
                    Else
                        If element > WorksheetMax Then WorksheetMax = element
                    End If
     
                End If
            Next
     
            Stopwatch1.Stop()
            TimeofFunction = TimeofFunction + Stopwatch1.Elapsed
     
        End Function
     
        Function WorksheetMin(tr As Object, Optional InitialRowNumber As Integer = 0, Optional InitialColumnNumber As Integer = 0)
            Static TimeofFunction
            Dim Stopwatch1 As Stopwatch = Stopwatch.StartNew()
     
            WorksheetMin = Nothing
     
            For Each element In tr
     
                If IsNumeric(element) Then
     
     
                    If IsNothing(WorksheetMin) Then
                        WorksheetMin = element
                    Else
                        If element < WorksheetMin Then WorksheetMin = element
                    End If
     
                End If
     
            Next
     
            Stopwatch1.Stop()
            TimeofFunction = TimeofFunction + Stopwatch1.Elapsed
     
        End Function

  2. #2
    Inactif  

    Homme Profil pro
    Développeur .NET
    Inscrit en
    janvier 2012
    Messages
    4 905
    Détails du profil
    Informations personnelles :
    Sexe : Homme
    Âge : 66
    Localisation : Canada

    Informations professionnelles :
    Activité : Développeur .NET
    Secteur : Finance

    Informations forums :
    Inscription : janvier 2012
    Messages : 4 905
    Points : 10 171
    Points
    10 171
    Billets dans le blog
    36
    Par défaut
    Bonjour,

    Je n'ai pas lu tout ton code. Et il faudrait qu'il soit indenté et balisé.

    à la différence de VBA, en VB.net tous les indices (de tableaux. de boucles, de collections etc) commencent obligatoirement à zéro. Et en conséquence:

    Code : Sélectionner tout - Visualiser dans une fenêtre à part
    Dim nf as long (1 to ndof, 1 to nn)
    ne peut pas être utilisé. Tu as donc fondamentalement deux choix :

    Code : Sélectionner tout - Visualiser dans une fenêtre à part
    Dim nf as long (ndof, nn)
    Cela te permet de continuer à partir de 1 (un) mais ton indice 0 (zéro) reste vide. Cela gaspille un peu de mémoire, mais cela évite de reprendre toutes les numérotations qui sont liées au tableau. Et c'est plus rapide si tu es pressé

    ou bien:

    Code : Sélectionner tout - Visualiser dans une fenêtre à part
    Dim nf as long (ndof -1, nn -1)
    Là tu pars de zéro, mais tu vas (probablement) devoir ajuster tes numérotations. C'est la même chose avec Redim

    P.S. si cela t'intéresse; il existe un paquet de bibliothèques (même des gratuites) toutes faites pour des calculs scientifiques en .net. Et beaucoup sont disponibles en packages NuGet.
    À ma connaissance, le seul personnage qui a été diagnostiqué comme étant allergique au mot effort. c'est Gaston Lagaffe.

    Ô Saint Excel, Grand Dieu de l'Inutile.

    Excel n'a jamais été, n'est pas et ne sera jamais un SGBD, c'est pour cela que Excel s'appelle Excel et ne s'appelle pas Access junior.

  3. #3
    Responsable
    Office & Excel


    Homme Profil pro
    Formateur et développeur chez EXCELLEZ.net
    Inscrit en
    novembre 2003
    Messages
    18 963
    Détails du profil
    Informations personnelles :
    Sexe : Homme
    Âge : 55
    Localisation : Belgique

    Informations professionnelles :
    Activité : Formateur et développeur chez EXCELLEZ.net
    Secteur : Enseignement

    Informations forums :
    Inscription : novembre 2003
    Messages : 18 963
    Points : 55 503
    Points
    55 503
    Billets dans le blog
    128
    Par défaut
    Salut.

    Ca n'a pour moi aucun sens de vouloir traduire du vba en vb.net. Le premier est procédural avec un peu d'objet, le second est POO. Il me semblerait plus sage d'analyser ce que fait le code en VBA, puis de repenser "l'application" en orienté objet et de l'écrire en vb.net. Traduire chaque ligne de VBA à vb.net ne peut amener à mon sens qu'à un code mal foutu.
    "Plus les hommes seront éclairés, plus ils seront libres" (Voltaire)
    ---------------
    Mes billets de blog sur DVP
    Mes remarques et critiques sont purement techniques. Ne les prenez jamais pour des attaques personnelles...
    Pensez à utiliser les tableaux structurés. Ils vous simplifieront la vie, tant en Excel qu'en VBA ==> mon tuto
    Le VBA ne palliera jamais une mauvaise conception de classeur ou un manque de connaissances des outils natifs d'Excel...
    Ce ne sont pas des bonnes pratiques parce que ce sont les miennes, ce sont les miennes parce que ce sont des bonnes pratiques
    VBA pour Excel? Pensez D'ABORD en EXCEL avant de penser en VBA...
    ---------------

  4. #4
    Candidat au Club
    Homme Profil pro
    Développeur .NET
    Inscrit en
    octobre 2019
    Messages
    3
    Détails du profil
    Informations personnelles :
    Sexe : Homme
    Localisation : Espagne

    Informations professionnelles :
    Activité : Développeur .NET
    Secteur : Bâtiment Travaux Publics

    Informations forums :
    Inscription : octobre 2019
    Messages : 3
    Points : 4
    Points
    4
    Par défaut
    Bonjour, et merci pour vôtres réponses, (et pour la balisation du code ).


    clementmarcotte:

    Quelle est le probleme avec "Dim nf As array = Array.CreateInstance(GetType(Long), {ndof,nn}, {1,1})" ??

    Je vais regarder la possibilité de ré-écrire le code avec des base 0:
    Dim nf as long (ndof -1, nn -1)



    Pierre Fauconnier

    En général je serais d'accord, mais pour cette fonction il me semble que les structures sont pas trop différentes ... pourrais tu donner un exemple plus précis a propos de cette ré-structuration?

  5. #5
    Expert éminent sénior Avatar de Pol63
    Homme Profil pro
    .NET / SQL SERVER
    Inscrit en
    avril 2007
    Messages
    14 039
    Détails du profil
    Informations personnelles :
    Sexe : Homme
    Âge : 40
    Localisation : France, Puy de Dôme (Auvergne)

    Informations professionnelles :
    Activité : .NET / SQL SERVER

    Informations forums :
    Inscription : avril 2007
    Messages : 14 039
    Points : 24 836
    Points
    24 836
    Par défaut
    .net est un langage beaucoup plus évolué que vba, et avec une méthodologie différente
    biensur on peut encore écrire du code spaghetti comme il y a 20 ans, ca fonctionnera, mais ce n'est pas fait pour

    en .net on n'utilise plus de tableau de variable, et on ne fait pas de méthode qui demande 50 paramètres, avec la programmation orienté objet tout est beaucoup simple et on peut diviser le nombre de lignes de code nécessaire par 10 ou 20 dans le meilleur des cas.

    par contre ca demande d'apprendre, et de reréfléchir l'application, et malheureusement ceux qui portent des logiciels n'ont parfois ni le temps ni l'envie de faire ca, et font donc un copier coller avec des adaptations
    Cours complets, tutos et autres FAQ ici : C# - VB.NET

  6. #6
    Responsable
    Office & Excel


    Homme Profil pro
    Formateur et développeur chez EXCELLEZ.net
    Inscrit en
    novembre 2003
    Messages
    18 963
    Détails du profil
    Informations personnelles :
    Sexe : Homme
    Âge : 55
    Localisation : Belgique

    Informations professionnelles :
    Activité : Formateur et développeur chez EXCELLEZ.net
    Secteur : Enseignement

    Informations forums :
    Inscription : novembre 2003
    Messages : 18 963
    Points : 55 503
    Points
    55 503
    Billets dans le blog
    128
    Par défaut
    Citation Envoyé par Pol63 Voir le message
    [...]
    Je n'ai pas envie d'entrer dans un débat, mais j'ai quand même envie de dire ceci:
    • en vba non plus, on n'utilise pas 50 paramètres pour une fonction (on peut même créer ses propres classes et ses propres interfaces, si si );
    • en vba, on peut ne pas écrire du code spaghetti;
    • en vba, on peut aussi, en codant proprement, réduire drastiquement ses lignes de code;
    • Il n'y a pas de honte à utiliser des tableaux, surtout lorsque l'on manipule des feuilles Excel (qui permettent notamment de réduire le nombre de paramètres passés à une fonction), et je ne vois pas pourquoi on n'utiliserait pas des tableaux en dot.net;
    • En vba, on peut écrire des codes génériques et donc facilement réutilisables;
    • en VBA, on peut typer ses variables, on peut même s'obliger à les déclarer;
    • en VBA, on peut architecturer son code en "trois tiers";
    • en VBA, on peut gérer les exceptions avec des Try...Catch...Finally (même si ça ne s'appelle pas comme ça);
    • ...


    Bref, on peut écrire du code propre en VBA, comme on peut fabriquer de la merde en vb.net (surtout lorsque l'on essaie de "copier" vba en vb.net).


    Ce n'est vraiment pas sur les exemples cités que l'on peut affirmer une "supériorité" de vb.net par rapport à VBA. Il y a d'un côté un langage OO et de l'autre un langage procéduralEt je préfère un code bien écrit en vba qu'une bouillie infâme en vb.net.
    "Plus les hommes seront éclairés, plus ils seront libres" (Voltaire)
    ---------------
    Mes billets de blog sur DVP
    Mes remarques et critiques sont purement techniques. Ne les prenez jamais pour des attaques personnelles...
    Pensez à utiliser les tableaux structurés. Ils vous simplifieront la vie, tant en Excel qu'en VBA ==> mon tuto
    Le VBA ne palliera jamais une mauvaise conception de classeur ou un manque de connaissances des outils natifs d'Excel...
    Ce ne sont pas des bonnes pratiques parce que ce sont les miennes, ce sont les miennes parce que ce sont des bonnes pratiques
    VBA pour Excel? Pensez D'ABORD en EXCEL avant de penser en VBA...
    ---------------

  7. #7
    Expert éminent sénior Avatar de Pol63
    Homme Profil pro
    .NET / SQL SERVER
    Inscrit en
    avril 2007
    Messages
    14 039
    Détails du profil
    Informations personnelles :
    Sexe : Homme
    Âge : 40
    Localisation : France, Puy de Dôme (Auvergne)

    Informations professionnelles :
    Activité : .NET / SQL SERVER

    Informations forums :
    Inscription : avril 2007
    Messages : 14 039
    Points : 24 836
    Points
    24 836
    Par défaut
    Citation Envoyé par Pierre Fauconnier Voir le message
    Je n'ai pas envie d'entrer dans un débat
    et pourquoi pas, ca permet d'apporter son savoir et ses opinions

    je suis pas expert vba, je suis parti avec une comparaison avec vb6 qui dans mon inconscient ressemble pas mal
    mais tu vas me dire qu'en vb6 aussi on pouvait faire des classes et des interfaces (il parait ^^)

    enfin en vb6 (et ici sur du vba) voir methode (param1 as string, param2 as integer, … param29 as boolean) c'est courant, dès lors qu'on a une classe ca devient plus pratique et plus lisible

    concernant les tableaux, c'est surtout qu'en .net écrire tab(x)(y)(z) c'est ptete un pouillème plus performant mais c'est dommage de passer à côté de myList(x).Parent.Size
    après si tu me dis à nouveau qu'on pouvait réellement faire en ca en vba/6 avec des classes, je suis étonné de n'avoir vu aucun développeur faire ca et que tout le monde faisait des tableaux de tableaux illisibles
    (après vb6 ca remonte à mes débuts et j'en ai pas fait longtemps)

    et je suis pas totalement contre les tableaux, en interne dans une classe
    pour manipuler des octets
    c'est bien, mais les exposer et en empiler je suis pas fan ...

    concernant le nombre de ligne de code, par contre vba doit rester loin de .net (linq, héritage etc...)
    Cours complets, tutos et autres FAQ ici : C# - VB.NET

  8. #8
    Responsable
    Office & Excel


    Homme Profil pro
    Formateur et développeur chez EXCELLEZ.net
    Inscrit en
    novembre 2003
    Messages
    18 963
    Détails du profil
    Informations personnelles :
    Sexe : Homme
    Âge : 55
    Localisation : Belgique

    Informations professionnelles :
    Activité : Formateur et développeur chez EXCELLEZ.net
    Secteur : Enseignement

    Informations forums :
    Inscription : novembre 2003
    Messages : 18 963
    Points : 55 503
    Points
    55 503
    Billets dans le blog
    128
    Par défaut
    Citation Envoyé par Pol63 Voir le message
    et pourquoi pas, ca permet d'apporter son savoir et ses opinions[...]
    Tout à fait, mes réticences relevaient plus de la volonté de ne pas "trop" polluer la discussion existante.


    Citation Envoyé par Pol63 Voir le message
    [...]tu vas me dire qu'en vb6 aussi on pouvait faire des classes et des interfaces (il parait ^^)[...]
    Voir mon tuto sur les classes en VB6/VBA. Je n'y parle pas d'interface, mais c'est très facile à créer en VBA.


    Citation Envoyé par Pol63 Voir le message
    [...]
    enfin en vb6 (et ici sur du vba) voir methode (param1 as string, param2 as integer, … param29 as boolean) c'est courant. dès lors qu'on a une classe ca devient plus pratique et plus lisible[...]
    Je ne vois pas cela très souvent... Mais c'est aussi pourquoi je parle de code bien écrit, qui évite cela en VBA (soit en passant par une structure ou une classe perso,soit en passant un tableau ou une collection en arguments). Cela dit, rien n'empêche d'écrire une fonction vb.net qui reçoit 29 arguments. J'ai déjà vu des constructeurs qui fonctionnaient ainsi, et même si dans certains cas, je peux "admettre", je pense que ce n'est pas dans l'esprit POO. Et pour terminer, si c'est pour ne pas les passer en arguments mais avoir les 29 lignes d'affectation de propriétés, je ne suis pas tout à fait persuadé que c'est vraiment mieux


    Citation Envoyé par Pol63 Voir le message
    [...]
    concernant les tableaux, c'est surtout qu'en .net écrire tab(x)(y)(z) c'est ptete un pouillème plus performant mais c'est dommage de passer à côté de myList(x).Parent.Size[...]
    J'admets, j'ai d'ailleurs un peu tendance à considérer les LIST comme des tableaux avec deux ou trois trucs chouettes en plus, comme je considère les collections VBA comme des tableaux avec deux trois trucs chouettes en plus...


    Citation Envoyé par Pol63 Voir le message
    [...]
    après si tu me dis à nouveau qu'on pouvait réellement faire en ca en vba/6 avec des classes[...]
    qu'on peut faire quoi? Je n'ai pas compris... Au sujet des classes que l'on peut créer en VBA, c'est clair qu'il n'y pas l'héritage (qu'on peut "modéliser" jusqu'à un certain point), ni le polymorphisme, ni même la possibilité de surcharger le constructeur. Donc ok, il y a des trucs gênants, mais rien d'insurmontable. La surcharge n'est pas possible non plus, mais peut être simulée avec les arguments optionnels, qui ont d'ailleurs fait leur apparition en dot.net...



    Citation Envoyé par Pol63 Voir le message
    [...]
    c'est bien, mais les exposer et en empiler je suis pas fan ...[...]
    Je ne suis pas fan non plus des tableaux de tableaux...



    Citation Envoyé par Pol63 Voir le message
    [...]
    concernant le nombre de ligne de code, par contre vba doit rester loin de .net (linq, héritage etc...)
    Je ne nie pas le fait, mais je maintiens qu'un code VBA bien écrit et architecturé permet une réutilisation de code intéressante, et donc moins de lignes de code. Et j'ai déjà vu des héritages avec tellement de trucs overridés que je peux dire qu'en termes de gain (et d'héritage), c'était bof bof.

    Cela étant, je ne dis pas que la POO c'est mal et qu'il faut se cantonner au procédural, et je ne doute pas des avantages de la POO dans une série de situations, loin de là. C'est d'ailleurs le but de ma première intervention dans cette discussion: ne pas transposer du code vba en vb.net, mais le refondre dans l'esprit POO.

    "Plus les hommes seront éclairés, plus ils seront libres" (Voltaire)
    ---------------
    Mes billets de blog sur DVP
    Mes remarques et critiques sont purement techniques. Ne les prenez jamais pour des attaques personnelles...
    Pensez à utiliser les tableaux structurés. Ils vous simplifieront la vie, tant en Excel qu'en VBA ==> mon tuto
    Le VBA ne palliera jamais une mauvaise conception de classeur ou un manque de connaissances des outils natifs d'Excel...
    Ce ne sont pas des bonnes pratiques parce que ce sont les miennes, ce sont les miennes parce que ce sont des bonnes pratiques
    VBA pour Excel? Pensez D'ABORD en EXCEL avant de penser en VBA...
    ---------------

+ Répondre à la discussion
Cette discussion est résolue.

Discussions similaires

  1. Syntaxe VBA DDEpoke sous Access
    Par Mad Ant dans le forum Access
    Réponses: 3
    Dernier message: 13/03/2007, 15h02
  2. Erreur de syntaxe - Access et Ado.net
    Par tibblou dans le forum Accès aux données
    Réponses: 6
    Dernier message: 25/01/2007, 20h26
  3. Questions de syntaxe VBA
    Par antoine46 dans le forum Access
    Réponses: 5
    Dernier message: 18/12/2006, 00h23
  4. equivalent avec vba access de isnull
    Par celiaaa dans le forum Access
    Réponses: 1
    Dernier message: 07/12/2006, 19h01
  5. Lien entre VBA et VB.net ?
    Par RaelRiaK dans le forum Access
    Réponses: 9
    Dernier message: 10/11/2006, 16h09

Partager

Partager
  • Envoyer la discussion sur Viadeo
  • Envoyer la discussion sur Twitter
  • Envoyer la discussion sur Google
  • Envoyer la discussion sur Facebook
  • Envoyer la discussion sur Digg
  • Envoyer la discussion sur Delicious
  • Envoyer la discussion sur MySpace
  • Envoyer la discussion sur Yahoo