DG/HDG Splitting iTutorial

ddrake
Topic Author
Offline
Premium Member

7 years 2 months ago - 7 years 1 month ago #589 by ddrake

DG/HDG Splitting iTutorial was created by ddrake

Hi,

I found a small bug in the DG/HDG Splitting iTutorial (Unit 3.4). In input block 7 there is the function:

Code:

def SolveConvectionSteps(gfuvec, res, tau, steps):
    bfmixed.Apply (gfuvec, gfuL2.vec) 
    VL2.SolveM(CoefficientFunction(1), gfuL2.vec)
    conv_applied = gfuL2.vec.CreateVector()
    for i in range(steps):
        convL2.Apply(gfuL2.vec,conv_applied)
        VL2.SolveM(CoefficientFunction(1), conv_applied)
        gfuL2.vec.data -= tau/steps * conv_applied
        #Redraw()    
    bfmixedT.Apply (gfuL2.vec, res)

The two calls to SolveM give:

TypeError: SolveM(): incompatible function arguments. The following argument types are supported:
1. (self: ngsolve.comp.FESpace, vec: ngsolve.la.BaseVector, rho: ngsolve.fem.CoefficientFunction=None) -> None

The error can be fixed by reversing the order of the arguments to SolveM or providing keywords like this:

Code:

VL2.SolveM(rho=CoefficientFunction(1), vec=conv_applied)

Edit: The same error is present in block 9 of the DG-method for convection (Unit 3.3)

Best,
Dow

Last edit: 7 years 1 month ago by ddrake.

ddrake
Topic Author
Offline
Premium Member

7 years 1 month ago #603 by ddrake

Replied by ddrake on topic DG/HDG Splitting iTutorial

These issues are fixed in the current iTutorials:
DG-method for convection (Unit 3.3) and DG/HDG Splitting iTutorial (Unit 3.4)

Thanks!
Dow

Time to create page: 0.120 seconds