[BUG]:

[ruf10]

How to reproduce the bug

When I run my code I have the following errors and sometimes my code got ended earlier

Minimal Example (Python)

kn_form = form(kn*dx)
        kn_scalar = fem.assemble.assemble_scalar(kn_form)

        for jj in range(k_on, t_num):
            t = t+dt
            count = count + 1
            kn_scalar = 0.5*(kn_scalar+ abs(kn_scalar))
            nu_t = kn_scalar*mu*tau*(func_wall_n_dist/LL)**2 #no-slip BCs
            F = (1./dt)*inner(u,v)*dx + eps*p*q*dx + b(2*un-unMinus1,u,v)*dx +2.*nu*a_sym(u,v)*dx
            F += - c(p,v)*dx + c(q,u)*dx + nu_t*a_sym(u,v)*dx
            F -= (1./dt)*(inner(un,v))*dx

            unPlus1, pnPlus = solve_u(F, t)
            #update solutions
            unMinus1.x.array[:] = un.x.array[:]
            un.x.array[:] = unPlus1.x.array[:]
            u_Sym = grad_sym(unPlus1)
            frobenius_norm_squared_u_Sym = frobenius_norm_squared(u_Sym)
            M_form = form(nu_t*frobenius_norm_squared_u_Sym*dx)
            # Assemble integral over the domain
            M = fem.assemble.assemble_scalar(M_form)

            epsilon_turb = 1/vol * M
            knPlus1_scalar  = (kn_scalar+ epsilon_turb*dt)/(1 + np.sqrt(2.0)*dt/(2.0*tau))
            #update solutions, onlyt take the postive part of kn
            kn_scalar = knPlus1_scalar

Output (Python)

INFO:root:running build_ext
INFO:root:building 'libffcx_forms_aa98e3af6fcc9a609fc690a72c688b994578c0c8' extension
INFO:root:clang -Wno-unused-result -Wsign-compare -Wunreachable-code -DNDEBUG -fwrapv -O2 -Wall -fPIC -O2 -isystem /Users/ruifang/opt/anaconda3/envs/fenicsx-env/include -fPIC -O2 -isystem /Users/ruifang/opt/anaconda3/envs/fenicsx-env/include -I/Users/ruifang/opt/anaconda3/envs/fenicsx-env/lib/python3.10/site-packages/ffcx/codegeneration -I/Users/ruifang/opt/anaconda3/envs/fenicsx-env/include/python3.10 -c libffcx_forms_aa98e3af6fcc9a609fc690a72c688b994578c0c8.c -o ./libffcx_forms_aa98e3af6fcc9a609fc690a72c688b994578c0c8.o -O2 -g0
INFO:root:clang -bundle -undefined dynamic_lookup -Wl,-rpath,/Users/ruifang/opt/anaconda3/envs/fenicsx-env/lib -L/Users/ruifang/opt/anaconda3/envs/fenicsx-env/lib -Wl,-rpath,/Users/ruifang/opt/anaconda3/envs/fenicsx-env/lib -L/Users/ruifang/opt/anaconda3/envs/fenicsx-env/lib ./libffcx_forms_aa98e3af6fcc9a609fc690a72c688b994578c0c8.o -o ./libffcx_forms_aa98e3af6fcc9a609fc690a72c688b994578c0c8.cpython-310-darwin.so
dyld[1773]: Assertion failed: (this->magic == kMagic), function contains, file Loader.cpp, line 144.

zsh: abort      python -u ct6_eps0_tau1.py > june2_log6.txt

Version

0.6.0

DOLFINx git commit

No response

Installation

No response

Additional information

No response

[jorgensd]

The code provided is not sufficient to reproduce the issue. Please post a minimal reproducible example

[garth-wells]

No follow up on request for more information, so closing.