Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Jump to bottom

BandpassFilter results are empty #20

Open
dan-gies opened this issue Mar 17, 2024 · 1 comment
Open

BandpassFilter results are empty #20

dan-gies opened this issue Mar 17, 2024 · 1 comment

Comments

@dan-gies
Copy link

I'm able to successfully run the BandpassFilter example, but as it stands, the results appear empty.

The data.txt file:

   0.000000000000E+000   0.000000000000E+000
   0.000000000000E+000   0.000000000000E+000
   0.000000000000E+000   0.000000000000E+000

The output of ElmerSolver:

Console log (click to expand) 
ELMER SOLVER (v 9.0) STARTED AT: 2024/03/17 01:05:11
ParCommInit:  Initialize #PEs:            1
MAIN: OMP_NUM_THREADS not set. Using only 1 thread per task.
MAIN:
MAIN: =============================================================
MAIN: ElmerSolver finite element software, Welcome!
MAIN: This program is free software licensed under (L)GPL
MAIN: Copyright 1st April 1995 - , CSC - IT Center for Science Ltd.
MAIN: Webpage http://www.csc.fi/elmer, Email elmeradm@csc.fi
MAIN: Version: 9.0 (Rev: 182610423, Compiled: 2024-03-15)
MAIN:  Running one task without MPI parallelization.
MAIN:  Running with just one thread per task.
MAIN: =============================================================
MAIN:
MAIN:
MAIN: -------------------------------------
MAIN: Reading Model: Filter_Zhai.sif
LoadInputFile: Scanning input file: Filter_Zhai.sif
LoadInputFile: Scanning only size info
LoadInputFile: First time visiting
LoadInputFile: Reading base load of sif file
LoadInputFile: Loading input file: Filter_Zhai.sif
LoadInputFile: Reading base load of sif file
Loading user function library: [EMParam_WR28]...[omega]
Loading user function library: [EMParam_WR28]...[MagnBndLoad]
Loading user function library: [EMParam_WR28]...[betaNull]
Loading user function library: [EMParam_WR28]...[betaNull]
LoadInputFile: Number of BCs: 3
LoadInputFile: Number of Body Forces: 0
LoadInputFile: Number of Initial Conditions: 0
LoadInputFile: Number of Materials: 1
LoadInputFile: Number of Equations: 1
LoadInputFile: Number of Solvers: 4
LoadInputFile: Number of Bodies: 1
ListTagKeywords: Setting weight for keywords!
ListTagKeywords: No parameters width suffix: normalize by area
ListTagKeywords: Setting weight for keywords!
ListTagKeywords: No parameters width suffix: normalize by volume
Loading user function library: [VectorHelmholtz]...[VectorHelmholtzSolver_Init0]
EdgeElementStyle: "Quadratic Approximation" requested without Piola. Setting "Use Piola Transform = True"
Loading user function library: [VectorHelmholtz]...[VectorHelmholtzCalcFields_Init0]
Loading user function library: [SaveData]...[SaveScalars_Init0]
Loading user function library: [ResultOutputSolve]...[ResultOutputSolver_Init0]
Loading user function library: [VectorHelmholtz]...[MagnetoDynamics_Dummy_Init0]
ElmerAsciiMesh: Base mesh name: mesh/Hl_Filter_Zhai_a08
PermuteNodeNumbering: Performing node mapping
MapBodiesAndBCs: Minimum initial body index: 1
MapBodiesAndBCs: Maximum initial body index: 1
MapBodiesAndBCs: Minimum initial boundary index: 1
MapBodiesAndBCs: Maximum initial boundary index: 3
LoadMesh: Elapsed REAL time:     0.0980 (s)
MAIN: -------------------------------------
Loading user function library: [VectorHelmholtz]...[VectorHelmholtzSolver_Init]
Loading user function library: [VectorHelmholtz]...[VectorHelmholtzSolver_bulk]
Loading user function library: [VectorHelmholtz]...[VectorHelmholtzSolver]
OptimizeBandwidth: ---------------------------------------------------------
OptimizeBandwidth: Computing matrix structure for: vector helmholtz equation
OptimizeBandwidth: Initial bandwidth for vector helmholtz equation: 40578
OptimizeBandwidth: Optimized bandwidth for vector helmholtz equation: 2258
OptimizeBandwidth: ---------------------------------------------------------
Loading user function library: [VectorHelmholtz]...[VectorHelmholtzCalcFields_Init]
Loading user function library: [VectorHelmholtz]...[VectorHelmholtzCalcFields_bulk]
Loading user function library: [VectorHelmholtz]...[VectorHelmholtzCalcFields]
OptimizeBandwidth: ---------------------------------------------------------
OptimizeBandwidth: Computing matrix structure for: vector helmholtz post process
OptimizeBandwidth: Initial bandwidth for vector helmholtz post process: 1660
OptimizeBandwidth: Optimized bandwidth for vector helmholtz post process: 153
OptimizeBandwidth: ---------------------------------------------------------
Loading user function library: [SaveData]...[SaveScalars_Init]
Loading user function library: [SaveData]...[SaveScalars_bulk]
Loading user function library: [SaveData]...[SaveScalars]
Loading user function library: [ResultOutputSolve]...[ResultOutputSolver_Init]
Loading user function library: [ResultOutputSolve]...[ResultOutputSolver_bulk]
Loading user function library: [ResultOutputSolve]...[ResultOutputSolver]
Loading user function library: [VectorHelmholtz]...[MagnetoDynamics_Dummy_Init]
Loading user function library: [VectorHelmholtz]...[MagnetoDynamics_Dummy_bulk]
Loading user function library: [VectorHelmholtz]...[MagnetoDynamics_Dummy]
MAIN: Number of timesteps to be saved: 3
MAIN:
MAIN: -------------------------------------
MAIN: Time: 1/3:   1.000E+00
MAIN: -------------------------------------
MAIN:
VectorHelmholtzSolver:
VectorHelmholtzSolver: -------------------------------------------------
VectorHelmholtzSolver: Solving harmonic electromagnetic wave equation!
EnforceDirichletConditions: Dirichlet conditions enforced for dofs: 26516
EnforceDirichletConditions: Dirichlet conditions enforced for dofs: 26516
SolveLinearSystem: Assuming complex valued linear system
CRS_ComplexILUT: ILU(T) (Complex), Starting factorization:
CRS_ComplexILUT:  ILU(T) (Complex), NOF nonzeros:      1577145
CRS_ComplexILUT:  ILU(T) (Complex), filling (%):          100
CRS_ComplexILUT: ILU(T) (Complex), Factorization ready at (s):     1.26
       1 0.1639E+00
       2 0.1111E+00
       3 0.6372E-01
       4 0.2984E-01
       5 0.3528E-01
       6 0.8439E-02
       7 0.5484E-02
       8 0.4529E-02
       9 0.3714E-02
      10 0.8877E-03
      11 0.4062E-03
      12 0.4368E-03
      13 0.8840E-04
      14 0.7870E-05
      15 0.5200E-06
      16 0.1090E-06
      17 0.7460E-07
ComputeChange: Nonlinear iteration did not converge to tolerance
ComputeChange: NS (ITER=1) (NRM,RELC): ( 0.29375131E-03  2.0000000     ) :: vector helmholtz equation
SaveScalars: -----------------------------------------
SaveScalars: Saving scalar values of various kinds
SaveScalars: Found 2 values to save in total
SaveScalars: Saving names of values to file: res/data.txt.names
SaveScalars: Saving values to file: res/data.txt
ResultOutputSolver: -------------------------------------
ResultOutputSolver: Saving with prefix: str
ResultOutputSolver: Creating list for saving - if not present
CreateListForSaving: Field Variables for Saving
CreateListForSaving: Vector Field 1: electric field re
CreateListForSaving: Vector Field 2: electric field im
CreateListForSaving: Vector Field 3: magnetic field strength re
CreateListForSaving: Vector Field 4: magnetic field strength im
CreateListForSaving: Vector Field 5: poynting vector re
CreateListForSaving: Vector Field 6: poynting vector im
CreateListForSaving: Vector Field 7: electric field re e
CreateListForSaving: Vector Field 8: electric field im e
CreateListForSaving: Vector Field 9: magnetic field strength re e
CreateListForSaving: Vector Field 10: magnetic field strength im e
CreateListForSaving: Vector Field 11: poynting vector re e
CreateListForSaving: Vector Field 12: poynting vector im e
ResultOutputSolver: Saving in unstructured VTK XML (.vtu) format
VtuOutputSolver: Saving results in VTK XML format with prefix: str
VtuOutputSolver: Saving number of partitions: 1
VtuOutputSolver: Setting offset for boundary entities: 100
ResultOutputSolver: -------------------------------------
MAIN:
MAIN: -------------------------------------
MAIN: Time: 2/3:   2.000E+00
MAIN: Elapsed time:  7.1 seconds
MAIN: Estimated time left: 10.8 seconds
MAIN: -------------------------------------
MAIN:
VectorHelmholtzSolver:
VectorHelmholtzSolver: -------------------------------------------------
VectorHelmholtzSolver: Solving harmonic electromagnetic wave equation!
EnforceDirichletConditions: Dirichlet conditions enforced for dofs: 26516
EnforceDirichletConditions: Dirichlet conditions enforced for dofs: 26516
SolveLinearSystem: Assuming complex valued linear system
CRS_ComplexILUT: ILU(T) (Complex), Starting factorization:
CRS_ComplexILUT:  ILU(T) (Complex), NOF nonzeros:      1577303
CRS_ComplexILUT:  ILU(T) (Complex), filling (%):          100
CRS_ComplexILUT: ILU(T) (Complex), Factorization ready at (s):     1.27
       1 0.6656E-02
       2 0.2196E-02
       3 0.1157E-02
       4 0.8285E-03
       5 0.3123E-03
       6 0.2498E-03
       7 0.2103E-03
       8 0.4581E-04
       9 0.1587E-04
      10 0.1309E-04
      11 0.8812E-05
      12 0.1269E-04
      13 0.4836E-06
      14 0.3779E-07
ComputeChange: Nonlinear iteration did not converge to tolerance
ComputeChange: NS (ITER=1) (NRM,RELC): ( 0.29869249E-03 0.16680668E-01 ) :: vector helmholtz equation
SaveScalars: -----------------------------------------
SaveScalars: Saving scalar values of various kinds
SaveScalars: Found 2 values to save in total
SaveScalars: Saving values to file: res/data.txt
ResultOutputSolver: -------------------------------------
ResultOutputSolver: Saving in unstructured VTK XML (.vtu) format
VtuOutputSolver: Setting offset for boundary entities: 100
ResultOutputSolver: -------------------------------------
MAIN:
MAIN: -------------------------------------
MAIN: Time: 3/3:   3.000E+00
MAIN: Elapsed time: 12.1 seconds
MAIN: Estimated time left:  5.0 seconds
MAIN: -------------------------------------
MAIN:
VectorHelmholtzSolver:
VectorHelmholtzSolver: -------------------------------------------------
VectorHelmholtzSolver: Solving harmonic electromagnetic wave equation!
EnforceDirichletConditions: Dirichlet conditions enforced for dofs: 26516
EnforceDirichletConditions: Dirichlet conditions enforced for dofs: 26516
SolveLinearSystem: Assuming complex valued linear system
CRS_ComplexILUT: ILU(T) (Complex), Starting factorization:
CRS_ComplexILUT:  ILU(T) (Complex), NOF nonzeros:      1577544
CRS_ComplexILUT:  ILU(T) (Complex), filling (%):          100
CRS_ComplexILUT: ILU(T) (Complex), Factorization ready at (s):     1.28
       1 0.6928E-02
       2 0.3007E-02
       3 0.1199E-02
       4 0.3534E-02
       5 0.3793E-03
       6 0.2850E-03
       7 0.2699E-03
       8 0.6368E-04
       9 0.3000E-04
      10 0.1094E-04
      11 0.6679E-05
      12 0.8440E-05
      13 0.4205E-06
      14 0.4692E-07
ComputeChange: Nonlinear iteration did not converge to tolerance
ComputeChange: NS (ITER=1) (NRM,RELC): ( 0.30399058E-03 0.17581677E-01 ) :: vector helmholtz equation
SaveScalars: -----------------------------------------
SaveScalars: Saving scalar values of various kinds
SaveScalars: Found 2 values to save in total
SaveScalars: Saving values to file: res/data.txt
ResultOutputSolver: -------------------------------------
ResultOutputSolver: Saving in unstructured VTK XML (.vtu) format
VtuOutputSolver: Setting offset for boundary entities: 100
ResultOutputSolver: -------------------------------------
MAIN: Reporting unused list entries for sif improvement!
MAIN: If you do not want these lines undefine > DEVEL_LISTUSAGE < !
   Unused keywords:
       Simulation            bdf order
       Constants             gravity
       Constants             stefan boltzmann
       Constants             boltzmann constant
       Constants             unit charge
       Material 1            relative permeability
       Solver 1              steady state convergence tolerance
       Solver 1              linear system abort not converged
       Solver 1              idrs parameter
       Solver 2              steady state convergence tolerance
       Solver 2              nonlinear system max iterations
       Solver 2              linear system iterative method
       Solver 2              linear system max iterations
       Solver 2              linear system convergence tolerance
       Solver 2              bicgstabl polynomial degree
       Solver 2              linear system ilut tolerance
       Solver 2              linear system abort not converged
MAIN: *** Elmer Solver: ALL DONE ***
ReleaseMeshFaceTables: Releasing number of faces: 12099
MAIN: The end
SOLVER TOTAL TIME(CPU,REAL):        16.40       17.14
ELMER SOLVER FINISHED AT: 2024/03/17 01:05:28

In Paraview, all fields appear as 0.00:

image
@gbassi94
Copy link

Hi @dan-gies, I am definitely not an expert but removing/commenting the line Exec Solver = never under Solver 2 in file Filter_Zhai.sif should do the job. As far as I understand Solver 2 computes the actual fields so it should be run.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Assignees
No one assigned
Labels
None yet
Projects
None yet
Milestone
No milestone
Development

No branches or pull requests
2 participants
@gbassi94 @dan-gies