ngsolve 6.2.2501__cp312-cp312-macosx_10_15_universal2.whl → 6.2.2501.post37.dev1__cp312-cp312-macosx_10_15_universal2.whl

ngsolve-6.2.2501.data/data/share/ngsolve/d10_DGdoubleglazing.pde

@@ -1,50 +0,0 @@
-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.2501.data/data/share/ngsolve/d11_chip_nitsche.pde

@@ -1,40 +0,0 @@
-#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.2501.data/data/share/ngsolve/d1_square.pde

@@ -1,43 +0,0 @@
-#
-# 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.2501.data/data/share/ngsolve/d2_chip.pde

@@ -1,35 +0,0 @@
-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.2501.data/data/share/ngsolve/d3_helmholtz.pde

@@ -1,22 +0,0 @@
-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.2501.data/data/share/ngsolve/d4_cube.pde

@@ -1,46 +0,0 @@
-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.2501.data/data/share/ngsolve/d5_beam.pde

@@ -1,74 +0,0 @@
-#
-# 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.2501.data/data/share/ngsolve/d6_shaft.pde

@@ -1,73 +0,0 @@
-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.2501.data/data/share/ngsolve/d7_coil.pde

@@ -1,50 +0,0 @@
-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.2501.data/data/share/ngsolve/d8_coilshield.pde

@@ -1,49 +0,0 @@
-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.2501.data/data/share/ngsolve/d9_hybridDG.pde

@@ -1,72 +0,0 @@
-########################################################################################
-# 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.2501.data/data/share/ngsolve/doubleglazing.in2d

@@ -1,27 +0,0 @@
-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