ngsolve 6.2.2506.post216.dev0__cp314-cp314-manylinux_2_28_x86_64.whl

ngsolve-6.2.2506.post216.dev0.data/data/share/ngsolve/d10_DGdoubleglazing.pde

@@ -0,0 +1,50 @@
+geometry = doubleglazing.in2d
+mesh = doubleglazing.vol
+
+####################### DISCONTINUOUS GALERKIN DEMO #############
+
+#### the double glazing problem: 
+#### taken from Howard Elman, David Silvester, Andy Wathen: 
+#### Finite Elements and Fast Iterative Solvers 
+#### with applications in incompressible fluid dynamics
+#### 3.1.4 Example
+
+define coefficient lam 
+(0.005),(0.005),
+
+#for inflow
+define coefficient dirichletcoef
+0, 1, 
+
+define coefficient calpha 
+10,10,
+
+#convection velocity
+define coefficient b1 
+(2*y*(1-x*x)),(2*y*(1-x*x)),
+
+#convection velocity
+define coefficient b2
+(-2*x*(1-y*y)),(-2*x*(1-y*y)),
+
+
+define fespace vdisc -type=l2ho -order=4 -dgjumps 
+define gridfunction udisc -fespace=vdisc
+
+define bilinearform adisc -fespace=vdisc
+DGIP_innfac_laplace lam calpha
+DGIP_bndfac_laplace lam calpha -definedon=[1,2]
+DG_innfac_convection b1 b2
+DG_bndfac_convection b1 b2
+convection b1 b2
+laplace lam
+
+define linearform fdisc -fespace=vdisc
+DGIP_bndfac_dir lam dirichletcoef calpha -definedon=2
+DG_bndfac_convdir b1 b2 dirichletcoef -definedon=2
+
+define preconditioner c -type=direct -bilinearform=adisc
+
+numproc bvp npdisc -bilinearform=adisc -linearform=fdisc -gridfunction=udisc -preconditioner=c
+
+numproc visualization npvis -scalarfunction=udisc -subdivision=4 -deformationscale=1

ngsolve-6.2.2506.post216.dev0.data/data/share/ngsolve/d11_chip_nitsche.pde

@@ -0,0 +1,40 @@
+#modified version of the d2_chip.pde example. here boundary conditions 
+#are incorporated in a weak sense with nitsche penalization
+geometry = chip.in2d
+mesh = chip.vol
+
+define coefficient lam
+1, 1000, 10,
+
+
+define coefficient coef_source
+0, 0, 1,
+
+define fespace v -type=h1ho -order=3 #-dirichlet=[1]
+define fespace verr -type=l2 -order=0
+
+define gridfunction u -fespace=v -nested
+define gridfunction err -fespace=verr
+
+define constant alpha = 10
+
+define bilinearform a -fespace=v -symmetric
+laplace lam
+nitsche lam alpha -definedon=1
+
+define linearform f -fespace=v
+source coef_source
+#nitsche lam dir alpha -definedon=1
+
+define preconditioner c -type=multigrid -bilinearform=a -smoothingsteps=3
+
+numproc bvp np1 -bilinearform=a -linearform=f -gridfunction=u -preconditioner=c -maxsteps=1000
+
+numproc drawflux np2 -bilinearform=a -solution=u -applyd -label=flux 
+
+numproc zzerrorestimator np3 -bilinearform=a -linearform=f -solution=u -error=err -minlevel=1
+numproc markelements np4 -error=err -minlevel=1 -factor=0.5 
+
+
+
+numproc visualization npv1 -scalarfunction=u -subdivision=2 -nolineartexture

ngsolve-6.2.2506.post216.dev0.data/data/share/ngsolve/d1_square.pde

@@ -0,0 +1,43 @@
+#
+# solve the Poisson equation -Delta u = f
+#
+# with boundary conditions
+#      u = 0  on Gamma1
+#  du/dn = 1  on Gamma2
+
+
+# load geometry
+geometry = square.in2d
+
+# and mesh
+mesh = square.vol
+
+
+# define a finite element space
+# Dirichlet boundary is Gamma_1 
+# play around with -order=...
+fespace v -type=h1ho -order=3 -dirichlet=[1]
+
+# the solution field
+gridfunction u -fespace=v -nested
+
+linearform f -fespace=v
+source x*sin(pi*y)
+# neumann 1 --definedon=[2]    # Neumann on Gamma_2
+
+# the bilinear-form 
+bilinearform a -fespace=v -symmetric -eliminate_internal
+laplace 1
+
+
+# preconditioner c -type=direct -bilinearform=a
+# preconditioner c -type=local -bilinearform=a 
+# preconditioner c -type=multigrid -bilinearform=a -smoother=block
+preconditioner c -type=bddc -bilinearform=a
+
+
+numproc bvp np1 -bilinearform=a -linearform=f -gridfunction=u -preconditioner=c -maxsteps=1000
+
+numproc drawflux np2 -bilinearform=a -solution=u -label=flux -applyd
+
+numproc visualization npv1 -scalarfunction=u -subdivision=2 -nolineartexture

ngsolve-6.2.2506.post216.dev0.data/data/share/ngsolve/d2_chip.pde

@@ -0,0 +1,35 @@
+geometry = chip.in2d
+mesh = chip.vol
+
+# one coefficient per sub-domain
+coefficient lam
+1, 1000, 10,
+
+# source in sub-domain 3
+coefficient coef_source
+0, 0, 1,
+
+fespace v -type=h1ho -order=3 -dirichlet=[1]
+fespace verr -type=l2 -order=0
+
+gridfunction u -fespace=v -nested
+gridfunction err -fespace=verr
+
+bilinearform a -fespace=v -symmetric
+laplace lam
+
+linearform f -fespace=v
+source coef_source
+
+preconditioner c -type=multigrid -bilinearform=a -smoothingsteps=3
+
+numproc bvp np1 -bilinearform=a -linearform=f -gridfunction=u -preconditioner=c -maxsteps=1000
+
+numproc drawflux np2 -bilinearform=a -solution=u -applyd -label=flux 
+
+numproc zzerrorestimator np3 -bilinearform=a -solution=u -error=err -minlevel=1
+numproc markelements np4 -error=err -minlevel=1 -factor=0.5 
+
+numproc writefile wf -variables=[mesh.levels,fes.v.ndof,ZZerrest.np3.err]  -filename=d2_error.out
+
+numproc visualization npv1 -scalarfunction=u -subdivision=2 -nolineartexture

ngsolve-6.2.2506.post216.dev0.data/data/share/ngsolve/d3_helmholtz.pde

@@ -0,0 +1,22 @@
+mesh = piezo2d40round4.vol.gz
+
+define constant k = 1.8
+
+define fespace v -type=h1ho -order=3 -complex -dirichlet=[1,2]
+
+define gridfunction u -fespace=v -nested
+
+define bilinearform a -fespace=v -symmetric
+laplace 1
+mass (-k*k)
+robin (-k*I) --definedon=4
+
+define linearform f -fespace=v
+
+define coefficient coef_dirichlet
+1, -1, 0, 0,
+numproc setvalues npsv -coefficient=coef_dirichlet -gridfunction=u -boundary
+
+numproc bvp np1 -bilinearform=a -linearform=f -gridfunction=u -solver=direct
+
+numproc visualization npv1 -scalarfunction=u -subdivision=2 -nolineartexture

ngsolve-6.2.2506.post216.dev0.data/data/share/ngsolve/d4_cube.pde

@@ -0,0 +1,46 @@
+geometry = cube.geo
+mesh = cube.vol
+
+constant heapsize = 10000000
+
+constant pen = 1e7
+
+coefficient lam
+1, 
+
+coefficient rho
+1,
+
+coefficient penalty
+(pen), 0,
+
+coefficient dirich_bc
+(pen*x*y), 0, 
+
+fespace v -order=4 -type=h1ho
+fespace vp -order=3 -dim=3 -type=h1ho
+
+gridfunction u -fespace=v -nested
+gridfunction p -fespace=vp
+
+linearform f -fespace=v
+source x*x*x*x
+neumann dirich_bc
+
+bilinearform a -fespace=v -symmetric -linearform=f -eliminate_internal
+laplace lam
+mass rho
+robin penalty
+
+
+
+# preconditioner c -type=direct -bilinearform=a
+# preconditioner c -type=local -bilinearform=a
+preconditioner c -type=multigrid -bilinearform=a -smoothingsteps=1 -smoother=block -notest -blocktype=9
+# preconditioner c -type=amg -bilinearform=a -coefe=lam -notiming -test
+
+numproc bvp np1 -bilinearform=a -linearform=f -gridfunction=u -preconditioner=c -maxsteps=200 -noprint -prec=1e-8
+
+numproc calcflux np2 -bilinearform=a -solution=u -flux=p -applyd
+        
+numproc visualization npv1 -scalarfunction=u -subdivision=2 -nolineartexture

ngsolve-6.2.2506.post216.dev0.data/data/share/ngsolve/d5_beam.pde

@@ -0,0 +1,74 @@
+#
+# PDE example file for solving a linear elasticity model
+#
+
+
+# load prepared geometry and mesh file:
+#
+geometry = beam.geo
+mesh = beam.vol
+
+
+define constant heapsize = 100000000
+
+
+# the elastic modulus for each sub-domain (here just one)
+define coefficient E
+2.1E11,
+
+# the Poisson ratio
+define coefficient nu
+0.2,
+
+
+# Displacement boundary conditions (clumped) are set by penalty
+# Give a large value for each clamped piece of the boundary 
+# The index corresponds to the value set in 'Edit Boundary Conditons'
+
+define coefficient penalty
+1e20, 0, 
+
+
+# A volume force in z-direction
+define coefficient coef_force_z
+7e4, 
+
+# A surface force in z-direction for each piece of the boundary
+define coefficient coef_surface_force_z
+0, 1e5, 
+
+
+define fespace v -type=h1ho -dim=3 -order=4 -dirichlet=[1]
+define fespace vp -type=h1ho -dim=6  -order=3
+
+define gridfunction u -fespace=v
+define gridfunction stress -fespace=vp
+
+# generate load vector, volume force in z-direction (-comp=3)
+define linearform f -fespace=v
+# source coef_force_z -comp=1
+neumann coef_surface_force_z -comp=3
+
+# define system matrix. robin adds penalty terms to the x,y, and z-components
+define bilinearform a -fespace=v -symmetric  -eliminate_internal -linearform=f
+elasticity E nu
+
+
+# use either a direct factorization, or a multigrid preconditioner.
+# for problems smaller than 10000 nodes, the direct one is usually faster,
+# otherwise the multigrid is faster (and needs much less memory)
+ 
+# define preconditioner c -type=direct -bilinearform=a
+define preconditioner c -type=multigrid -bilinearform=a  # -inverse=sparsecholesky
+
+#solve the problem by calling the cg-iteration
+numproc bvp np1 -bilinearform=a -linearform=f -gridfunction=u -preconditioner=c -maxsteps=200
+
+
+# postprocessing computes the stresses:
+numproc calcflux np2 -bilinearform=a -solution=u -flux=stress -applyd
+
+
+
+
+numproc visualization npv1 -vectorfunction=u -subdivision=2 -nolineartexture -deformationscale=100 -scalarfunction=stress.1

ngsolve-6.2.2506.post216.dev0.data/data/share/ngsolve/d6_shaft.pde

@@ -0,0 +1,73 @@
+geometry = shaft.geo
+mesh = shaft.vol
+
+define constant geometryorder = 3
+
+define coefficient E
+1e6, 
+
+define coefficient nu
+0, 
+
+define coefficient transmissionbc
+0, 0, 0, 0, 0.01, 0, 0, 0, 0,
+
+define coefficient bearingbc1
+1, 0, 0, 0, 0, 0, 0,
+
+define coefficient bearingbc2
+0, 1, 0, 0, 0, 0, 0,
+
+define coefficient bearingbc3
+0, 0, 1, 0, 0, 0, 0, 
+
+define coefficient surfaceloadz
+0, 0, 0, 1, 0, 0, 0, 0, 0,
+
+
+define fespace v -type=h1ho -dim=3 -order=3 
+define fespace vp -type=h1ho -dim=6 -order=2
+define fespace verr -type=l2 -order=0
+
+define gridfunction u -fespace=v -nested #-addcoef
+define gridfunction p -fespace=vp
+define gridfunction error -fespace=verr
+
+define bilinearform a -fespace=v -symmetric
+elasticity E nu
+robin transmissionbc -comp=1
+robin transmissionbc -comp=2
+robin transmissionbc -comp=3
+
+define linearform f -fespace=v
+neumann surfaceloadz -comp=3
+
+define linearform cnsty1 -fespace=v
+neumann bearingbc1 -comp=2
+define linearform cnstz1 -fespace=v
+neumann bearingbc1 -comp=3
+
+define linearform cnsty2 -fespace=v
+neumann bearingbc2 -comp=2
+define linearform cnstz2 -fespace=v
+neumann bearingbc2 -comp=3
+
+define linearform cnsty3 -fespace=v
+neumann bearingbc3 -comp=2
+define linearform cnstz3 -fespace=v
+neumann bearingbc3 -comp=3
+
+
+
+define preconditioner c -type=multigrid -bilinearform=a -smoothingsteps=2
+# define preconditioner c -type=local -bilinearform=a
+# define preconditioner c -type=direct -bilinearform=a
+
+
+numproc constrainedbvp np1 -bilinearform=a -linearform=f -gridfunction=u -preconditioner=c -maxsteps=50 -constraints=[cnsty1,cnstz1,cnsty2,cnstz2,cnsty3,cnstz3] 
+
+numproc calcflux np2 -bilinearform=a -solution=u -flux=p -applyd
+numproc drawflux np2a -bilinearform=a -solution=u -label=stress -applyd
+
+numproc zzerrorestimator np3 -bilinearform=a -linearform=f -solution=u -error=error -flux=p
+numproc markelements np3a -error=error -factor=0.5 -minlevel=1 

ngsolve-6.2.2506.post216.dev0.data/data/share/ngsolve/d7_coil.pde

@@ -0,0 +1,50 @@
+geometry = coil.geo
+mesh = coil.vol
+
+
+define constant geometryorder = 4
+# define constant heapsize = 50000000
+
+define constant mu0 = 1.257e-7
+
+# nu = 1/mu with mu permeability
+define coefficient nu
+(1/mu0), (1/mu0), 
+
+# artificial conductivity for regularization
+define coefficient sigma 1e-8*nu
+
+# 1000 Ampere-turns
+define coefficient r (sqrt(x*x+y*y))
+define coefficient cs 
+( 1000/16 * (y/r,-x/r,0))
+(0, 0, 0),
+
+define fespace v -type=hcurlho -order=4 -nograds -dirichlet=[1]
+
+define gridfunction u -fespace=v  -novisual
+
+define linearform f -fespace=v
+sourceedge cs  -definedon=1
+
+
+define bilinearform a -fespace=v -spd -nonsym_storage -linearform=f -eliminate_internal 
+curlcurledge nu 
+massedge sigma
+
+
+define bilinearform acurl -fespace=v -symmetric -nonassemble
+curlcurledge nu 
+
+
+# define preconditioner c -type=multigrid -bilinearform=a  -smoother=block
+define preconditioner c -type=bddc -bilinearform=a
+
+
+numproc bvp np1 -bilinearform=a -linearform=f -gridfunction=u -preconditioner=c  -prec=1.e-9
+
+numproc drawflux np5 -bilinearform=acurl -solution=u  -label=flux
+
+
+
+numproc visualization npv1 -scalarfunction=flux:0 -clipsolution=scalar -subdivision=2 -clipvec=[0,1,0] -nolineartexture

ngsolve-6.2.2506.post216.dev0.data/data/share/ngsolve/d8_coilshield.pde

@@ -0,0 +1,49 @@
+geometry = coilshield.geo
+mesh = coilshield.vol
+
+define constant geometryorder = 4
+define constant refinep = 1 
+define constant heapsize = 200000000
+
+define coefficient nu material
+iron 1e-4
+default 1
+
+define coefficient kappa material
+default 1e-6
+
+define coefficient cs1
+0, ( y ), 0, 0, 0, 
+define coefficient cs2
+0, ( -x ), 0, 0, 0, 
+define coefficient cs3
+0, 0, 0, 0, 0, 
+
+define coefficient penalty
+0,0,0,0,0,0,0,0,
+#1e+6, 0, 0, 0, 0, 0
+
+
+define fespace v -type=hcurlho -order=5  -nograds
+
+define gridfunction u -fespace=v
+
+define linearform f -fespace=v -noprintelvec
+sourceedge cs1 cs2 cs3  -definedon=2
+
+
+define bilinearform a -fespace=v -symmetric  -eliminate_internal  -linearform=f
+curlcurledge nu 
+massedge kappa  -order=2
+robinedge penalty
+	
+define bilinearform acurl -fespace=v -symmetric  -nonassemble
+curlcurledge nu
+
+
+define preconditioner c -type=multigrid -bilinearform=a -cylce=1 -smoother=block -coarsetype=direct -coarsesmoothingsteps=5 -notest
+
+
+numproc bvp np1 -bilinearform=a -linearform=f -gridfunction=u -preconditioner=c -maxsteps=400 -prec=1.e-8
+
+numproc drawflux np3 -bilinearform=acurl -solution=u -label=flux

ngsolve-6.2.2506.post216.dev0.data/data/share/ngsolve/d9_hybridDG.pde

@@ -0,0 +1,72 @@
+########################################################################################
+# solve poisson equation
+#    - eps Delta u = f in Omega,  u=g on Gamma0
+########################################################################################
+
+
+geometry = square.in2d
+mesh = square.vol
+
+# geometry = cube.geo
+# mesh = cube.vol
+
+
+define constant heapsize = 10000000
+
+
+# ho-fespace: compound space for (u, lambda) ######
+define fespace v -type=HDG -order=3 -dirichlet=[1,2]
+
+
+define gridfunction u -fespace=v
+
+
+
+
+## boundary terms ########################################
+
+define coefficient coef_dirichlet
+0, (y * (1-y)), 0, 0,  
+
+numproc setvalues npsv -coefficient=coef_dirichlet -gridfunction=u -boundary -comp=2
+
+
+## some coefficients #####################################
+define coefficient one
+1, 
+
+define coefficient alpha
+2,
+
+
+define coefficient b
+(50,10),
+
+define coefficient cf 
+1,
+
+
+define linearform f -fespace=v
+# source cf -comp=1
+
+
+define bilinearform a -fespace=v -eliminate_internal -linearform=f -printelmat
+HDG_laplace one alpha
+HDG_convection b
+
+
+define preconditioner c -type=direct -bilinearform=a # -inverse=pardiso
+
+numproc bvp np1 -bilinearform=a -linearform=f -gridfunction=u -preconditioner=c
+
+
+
+
+# visualization
+
+## automatic visualizaiton for 2D
+numproc visualization vis1 -scalarfunction=u -subdivision=2
+
+## automatic visualization for 3D
+# numproc visualization vis2 -scalarfunction=u -subdivision=1 -clipvec=[0,0,-1] -clipsolution=scalar
+

ngsolve-6.2.2506.post216.dev0.data/data/share/ngsolve/doubleglazing.in2d

@@ -0,0 +1,27 @@
+splinecurves2dv2
+2
+points
+#point number
+#	x-coordinate
+#		y-coordinate
+#			local refinement
+1	-1	 1 
+2	 1	 1	-maxh=0.01
+3	 1	-1	-maxh=0.01
+4	-1	-1
+
+segments
+#domain left of the segment
+#	domain right of the segment
+#		number of points the segment consists of
+#			first point of the segment
+#				second point of the segment
+#					boundary number
+#						local refinement	
+0	1	2	1	2	-maxh=0.25 -bc=1
+0	1	2	2	3	-maxh=0.1 -bc=2
+0	1	2	3	4	-maxh=0.25 -bc=1
+0	1	2	4	1	-maxh=0.25 -bc=1
+
+materials
+1	inner	-maxh=0.5