[Gmsh] Magnetodynamics with cohomology conditions

[Matti Pellikka]

Hello,

Such modification reduces the number of tunnels in the air region by one,
and therefore reduces the number of required "cohomology basis functions"
by one. The solution is to remove all the references to the non-existent
second cohomology basis function in the .pro-file, and check that the
remaining one refers to the correct physical group of the mesh. Here's a
.geo and a .pro-file that should work.

-- 
Matti Pellikka, Researcher, Electromagnetics
Tampere University of Technology, Finland



On 31 July 2013 09:17, Maik Franz <maikfranz25 at yahoo.com> wrote:

> I have big problems with the "Magnetodynamics with cohomology
> conditions"-Solver. If I modify the geometry of the tube e.g. if I create a
> massive cylinder insted of the tube, I recieve the following error. Please
> can you help me?? I'm trying to solve this problem for several days without
> success...
>
>
> Skindepth Al 0.0129949
> Skindepth Cu 0.00918881
> Info    : Selected Resolution 'MagDynTOComplex'
> Info    : Loading Geometric data 'indheat3d_therm5.msh'
> Info    : System 'A' : Complex, Frequency = 50 Hz
> P r e - P r o c e s s i n g . . .
> Info    : Treatment Formulation 'MagDynTO'
> Info    :   Generate ExtendedGroup '_BF_Entity_17' (NodesOf)
> Info    :   Generate ExtendedGroup '_BF_Entity_18' (EdgesOf)
> Info    :   Generate ExtendedGroup '_CO_Entity_21' (EdgesOfTreeIn)
> Info    :   Generate ExtendedGroup '_BF_Entity_19' (GroupsOfEdgesOf)
> Info    :   Generate ExtendedGroup '_BF_Entity_20' (GroupsOfEdgesOf)
> Info    : (Tue Jul 30 13:06:03 2013, CPU = 0.795605s, Mem = 21Mb)
> E n d   P r e - P r o c e s s i n g
> P r o c e s s i n g . . .
> Info    : InitSolution[A]
> Info    : Generate[A]
> Warning : Dt not implemented, using DtDof instead
> Info    : (Tue Jul 30 13:06:04 2013, CPU = 1.68481s, Mem = 49Mb)
> Info    : Solve[A]
> Info    : N: 13992 - preonly lu mumps
> [0]PETSC ERROR: --------------------- Error Message
> ----------------------------
> --------
> [0]PETSC ERROR: Error in external library!
> [0]PETSC ERROR: Error reported by MUMPS in numerical factorization phase:
> INFO(1
> )=-10, INFO(2)=0
> !
> [0]PETSC ERROR:
> ----------------------------------------------------------------
> --------
> [0]PETSC ERROR: Petsc Release Version 3.1.0, Patch 8, Thu Mar 17 13:37:48
> CDT 20
> 11
> [0]PETSC ERROR: See docs/changes/index.html for recent updates.
> [0]PETSC ERROR: See docs/faq.html for hints about trouble shooting.
> [0]PETSC ERROR: See docs/index.html for manual pages.
> [0]PETSC ERROR:
> ----------------------------------------------------------------
> --------
> [0]PETSC ERROR: indheat3d_therm5.pro on a win64_com named TURE-OFFICE by
> Ture Tu
> e Jul 30 13:06:03 2013
> [0]PETSC ERROR: Libraries linked from
> /home/geuzaine/src/petsc-3.1-p8/win64_comp
> lex_mumps_seq/lib
> [0]PETSC ERROR: Configure run at Thu Apr  7 15:27:36 2011
> [0]PETSC ERROR: Configure options --CC=/usr/bin/x86_64-w64-mingw32-gcc.exe
> --CXX
> =/usr/bin/x86_64-w64-mingw32-g++.exe
> --FC=/usr/bin/x86_64-w64-mingw32-gfortran.e
> xe --with-blas-lapack-dir=/usr/local/lib --with-debugging=0
> --with-clanguage=cxx
>  --with-shared=0 --with-x=0 --with-mpi=0
> --with-mumps-dir=/home/geuzaine/src/mum
> ps_seq --with-scalar-type=complex --useThreads=0
> [0]PETSC ERROR:
> ----------------------------------------------------------------
> --------
> [0]PETSC ERROR: MatFactorNumeric_MUMPS() line 517 in
> src/mat/impls/aij/mpi/mumps
> /mumps.c
> [0]PETSC ERROR: MatLUFactorNumeric() line 2587 in
> src/mat/interface/matrix.c
> [0]PETSC ERROR: PCSetUp_LU() line 158 in src/ksp/pc/impls/factor/lu/lu.c
> [0]PETSC ERROR: PCSetUp() line 795 in src/ksp/pc/interface/precon.c
> [0]PETSC ERROR: KSPSetUp() line 237 in src/ksp/ksp/interface/itfunc.c
> [0]PETSC ERROR: KSPSolve() line 353 in src/ksp/ksp/interface/itfunc.c
> [0]PETSC ERROR: User provided function() line 67 in
> "unknowndirectory/"/home/geu
> zaine/src/getdp/Legacy/LinAlg_PETSC.cpp
> Error   : PETSc error 76
>
>
>
> _______________________________________________
> gmsh mailing list
> gmsh at geuz.org
> http://www.geuz.org/mailman/listinfo/gmsh
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://www.geuz.org/pipermail/gmsh/attachments/20130731/506aebc7/attachment.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: indheat_massive.geo
Type: application/octet-stream
Size: 5161 bytes
Desc: not available
URL: <http://www.geuz.org/pipermail/gmsh/attachments/20130731/506aebc7/attachment.geo>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: indheat_massive.pro
Type: application/octet-stream
Size: 15920 bytes
Desc: not available
URL: <http://www.geuz.org/pipermail/gmsh/attachments/20130731/506aebc7/attachment.pro>