IEEE_UNDERFLOW_FLAG IEEE_DENORMAL

[ad10]

Dear all,

I am facing a strange issue here. Most of time to check my formulation I am using write(7,*)'......',..... but when I use this in Ubuntu 20.04 for same code it gives following errors.

Code: [Select]

  *ERROR* Residual norm is NaN or Inf
 --> ERRORS OCCURRED: For details see file: Osolid.txt

Note: The following floating-point exceptions are signalling: IEEE_INVALID_FLAG IEEE_OVERFLOW_FLAG IEEE_UNDERFLOW_FLAG IEEE_DENORMAL

As I have seen some post on forum I tried to debug using the GDB. I got following,

Code: [Select]

  *ERROR* Residual norm is NaN or Inf
 --> ERRORS OCCURRED: For details see file: Osolid.txt

Note: The following floating-point exceptions are signalling: IEEE_UNDERFLOW_FLAG IEEE_DENORMAL
[Inferior 1 (process 177227) exited normally]

I also tried to debug using valgrind, I got following

Code: [Select]

==177410== Memcheck, a memory error detector
==177410== Copyright (C) 2002-2017, and GNU GPL'd, by Julian Seward et al.
==177410== Using Valgrind-3.15.0-608cb11914-20190413 and LibVEX; rerun with -h for copyright info
==177410== Command: ./feap
==177410==
--177410-- Valgrind options:
--177410--    -v
--177410--    --leak-check=full
--177410--    --show-reachable=yes
--177410-- Contents of /proc/version:
--177410--   Linux version 5.11.0-37-generic (buildd@lcy01-amd64-021) (gcc (Ubuntu 9.3.0-17ubuntu1~20.04) 9.3.0, GNU ld (GNU Binutils for Ubuntu) 2.34) #41~20.04.2-Ubuntu SMP Fri Sep 24 09:06:38 UTC 2021
--177410--
--177410-- Arch and hwcaps: AMD64, LittleEndian, amd64-cx16-lzcnt-rdtscp-sse3-ssse3-avx-avx2-bmi-f16c-rdrand
--177410-- Page sizes: currently 4096, max supported 4096
--177410-- Valgrind library directory: /usr/lib/x86_64-linux-gnu/valgrind
--177410-- Reading syms from /home/ajinkya/Desktop/Viscous/user2/feap
--177410-- Reading syms from /usr/lib/x86_64-linux-gnu/ld-2.31.so
--177410--   Considering /usr/lib/x86_64-linux-gnu/ld-2.31.so ..
--177410--   .. CRC mismatch (computed 975d0390 wanted 30bd717f)
--177410--   Considering /lib/x86_64-linux-gnu/ld-2.31.so ..
--177410--   .. CRC mismatch (computed 975d0390 wanted 30bd717f)
--177410--   Considering /usr/lib/debug/lib/x86_64-linux-gnu/ld-2.31.so ..
--177410--   .. CRC is valid
--177410-- Reading syms from /usr/lib/x86_64-linux-gnu/valgrind/memcheck-amd64-linux
--177410--    object doesn't have a symbol table
--177410--    object doesn't have a dynamic symbol table
--177410-- Scheduler: using generic scheduler lock implementation.
--177410-- Reading suppressions file: /usr/lib/x86_64-linux-gnu/valgrind/default.supp
==177410== embedded gdbserver: reading from /tmp/vgdb-pipe-from-vgdb-to-177410-by-ajinkya-on-???
==177410== embedded gdbserver: writing to   /tmp/vgdb-pipe-to-vgdb-from-177410-by-ajinkya-on-???
==177410== embedded gdbserver: shared mem   /tmp/vgdb-pipe-shared-mem-vgdb-177410-by-ajinkya-on-???
==177410==
==177410== TO CONTROL THIS PROCESS USING vgdb (which you probably
==177410== don't want to do, unless you know exactly what you're doing,
==177410== or are doing some strange experiment):
==177410==   /usr/lib/x86_64-linux-gnu/valgrind/../../bin/vgdb --pid=177410 ...command...
==177410==
==177410== TO DEBUG THIS PROCESS USING GDB: start GDB like this
==177410==   /path/to/gdb ./feap
==177410== and then give GDB the following command
==177410==   target remote | /usr/lib/x86_64-linux-gnu/valgrind/../../bin/vgdb --pid=177410
==177410== --pid is optional if only one valgrind process is running
==177410==
--177410-- REDIR: 0x4022e10 (ld-linux-x86-64.so.2:strlen) redirected to 0x580c9ce2 (???)
--177410-- REDIR: 0x4022be0 (ld-linux-x86-64.so.2:index) redirected to 0x580c9cfc (???)
--177410-- Reading syms from /usr/lib/x86_64-linux-gnu/valgrind/vgpreload_core-amd64-linux.so
--177410--    object doesn't have a symbol table
--177410-- Reading syms from /usr/lib/x86_64-linux-gnu/valgrind/vgpreload_memcheck-amd64-linux.so
--177410--    object doesn't have a symbol table
==177410== WARNING: new redirection conflicts with existing -- ignoring it
--177410--     old: 0x04022e10 (strlen              ) R-> (0000.0) 0x580c9ce2 ???
--177410--     new: 0x04022e10 (strlen              ) R-> (2007.0) 0x0483f060 strlen
--177410-- REDIR: 0x401f5f0 (ld-linux-x86-64.so.2:strcmp) redirected to 0x483ffd0 (strcmp)
--177410-- REDIR: 0x4023370 (ld-linux-x86-64.so.2:mempcpy) redirected to 0x4843a20 (mempcpy)
--177410-- Reading syms from /usr/lib/x86_64-linux-gnu/libX11.so.6.3.0
--177410--    object doesn't have a symbol table
--177410-- Reading syms from /usr/lib/x86_64-linux-gnu/libgfortran.so.5.0.0
--177410--    object doesn't have a symbol table
--177410-- Reading syms from /usr/lib/x86_64-linux-gnu/libm-2.31.so
--177410--   Considering /usr/lib/x86_64-linux-gnu/libm-2.31.so ..
--177410--   .. CRC mismatch (computed fcb42c76 wanted f6c95789)
--177410--   Considering /lib/x86_64-linux-gnu/libm-2.31.so ..
--177410--   .. CRC mismatch (computed fcb42c76 wanted f6c95789)
--177410--   Considering /usr/lib/debug/lib/x86_64-linux-gnu/libm-2.31.so ..
--177410--   .. CRC is valid
--177410-- Reading syms from /usr/lib/x86_64-linux-gnu/libjpeg.so.8.2.2
--177410--    object doesn't have a symbol table
--177410-- Reading syms from /usr/lib/x86_64-linux-gnu/libgcc_s.so.1
--177410--    object doesn't have a symbol table
--177410-- Reading syms from /usr/lib/x86_64-linux-gnu/libc-2.31.so
--177410--   Considering /usr/lib/x86_64-linux-gnu/libc-2.31.so ..
--177410--   .. CRC mismatch (computed 86b78530 wanted e380f01c)
--177410--   Considering /lib/x86_64-linux-gnu/libc-2.31.so ..
--177410--   .. CRC mismatch (computed 86b78530 wanted e380f01c)
--177410--   Considering /usr/lib/debug/lib/x86_64-linux-gnu/libc-2.31.so ..
--177410--   .. CRC is valid
--177410-- Reading syms from /usr/lib/x86_64-linux-gnu/libxcb.so.1.1.0
--177410--    object doesn't have a symbol table
--177410-- Reading syms from /usr/lib/x86_64-linux-gnu/libdl-2.31.so
--177410--   Considering /usr/lib/x86_64-linux-gnu/libdl-2.31.so ..
--177410--   .. CRC mismatch (computed b5d487d3 wanted 6ef97e7c)
--177410--   Considering /lib/x86_64-linux-gnu/libdl-2.31.so ..
--177410--   .. CRC mismatch (computed b5d487d3 wanted 6ef97e7c)
--177410--   Considering /usr/lib/debug/lib/x86_64-linux-gnu/libdl-2.31.so ..
--177410--   .. CRC is valid
--177410-- Reading syms from /usr/lib/x86_64-linux-gnu/libquadmath.so.0.0.0
--177410--    object doesn't have a symbol table
--177410-- Reading syms from /usr/lib/x86_64-linux-gnu/libXau.so.6.0.0
--177410--   Considering /usr/lib/x86_64-linux-gnu/libXau.so.6.0.0 ..
--177410--   .. CRC mismatch (computed 7de9b6ad wanted e8a17129)
--177410--   Considering /lib/x86_64-linux-gnu/libXau.so.6.0.0 ..
--177410--   .. CRC mismatch (computed 7de9b6ad wanted e8a17129)
--177410--    object doesn't have a symbol table
--177410-- Reading syms from /usr/lib/x86_64-linux-gnu/libXdmcp.so.6.0.0
--177410--    object doesn't have a symbol table
--177410-- Reading syms from /usr/lib/x86_64-linux-gnu/libbsd.so.0.10.0
--177410--    object doesn't have a symbol table
--177410-- REDIR: 0x4ee7600 (libc.so.6:memmove) redirected to 0x48311d0 (_vgnU_ifunc_wrapper)
--177410-- REDIR: 0x4ee6900 (libc.so.6:strncpy) redirected to 0x48311d0 (_vgnU_ifunc_wrapper)
--177410-- REDIR: 0x4ee7930 (libc.so.6:strcasecmp) redirected to 0x48311d0 (_vgnU_ifunc_wrapper)
--177410-- REDIR: 0x4ee6220 (libc.so.6:strcat) redirected to 0x48311d0 (_vgnU_ifunc_wrapper)
--177410-- REDIR: 0x4ee6960 (libc.so.6:rindex) redirected to 0x48311d0 (_vgnU_ifunc_wrapper)
--177410-- REDIR: 0x4ee8dd0 (libc.so.6:rawmemchr) redirected to 0x48311d0 (_vgnU_ifunc_wrapper)
--177410-- REDIR: 0x4f03e60 (libc.so.6:wmemchr) redirected to 0x48311d0 (_vgnU_ifunc_wrapper)
--177410-- REDIR: 0x4f039a0 (libc.so.6:wcscmp) redirected to 0x48311d0 (_vgnU_ifunc_wrapper)
--177410-- REDIR: 0x4ee7760 (libc.so.6:mempcpy) redirected to 0x48311d0 (_vgnU_ifunc_wrapper)
--177410-- REDIR: 0x4ee7590 (libc.so.6:bcmp) redirected to 0x48311d0 (_vgnU_ifunc_wrapper)
--177410-- REDIR: 0x4ee6890 (libc.so.6:strncmp) redirected to 0x48311d0 (_vgnU_ifunc_wrapper)
--177410-- REDIR: 0x4ee62d0 (libc.so.6:strcmp) redirected to 0x48311d0 (_vgnU_ifunc_wrapper)
--177410-- REDIR: 0x4ee76c0 (libc.so.6:memset) redirected to 0x48311d0 (_vgnU_ifunc_wrapper)
--177410-- REDIR: 0x4f03960 (libc.so.6:wcschr) redirected to 0x48311d0 (_vgnU_ifunc_wrapper)
--177410-- REDIR: 0x4ee67f0 (libc.so.6:strnlen) redirected to 0x48311d0 (_vgnU_ifunc_wrapper)
--177410-- REDIR: 0x4ee63b0 (libc.so.6:strcspn) redirected to 0x48311d0 (_vgnU_ifunc_wrapper)
--177410-- REDIR: 0x4ee7980 (libc.so.6:strncasecmp) redirected to 0x48311d0 (_vgnU_ifunc_wrapper)
--177410-- REDIR: 0x4ee6350 (libc.so.6:strcpy) redirected to 0x48311d0 (_vgnU_ifunc_wrapper)
--177410-- REDIR: 0x4ee7ad0 (libc.so.6:memcpy@@GLIBC_2.14) redirected to 0x48311d0 (_vgnU_ifunc_wrapper)
--177410-- REDIR: 0x4f050d0 (libc.so.6:wcsnlen) redirected to 0x48311d0 (_vgnU_ifunc_wrapper)
--177410-- REDIR: 0x4f039e0 (libc.so.6:wcscpy) redirected to 0x48311d0 (_vgnU_ifunc_wrapper)
--177410-- REDIR: 0x4ee69a0 (libc.so.6:strpbrk) redirected to 0x48311d0 (_vgnU_ifunc_wrapper)
--177410-- REDIR: 0x4ee6280 (libc.so.6:index) redirected to 0x48311d0 (_vgnU_ifunc_wrapper)
--177410-- REDIR: 0x4ee67b0 (libc.so.6:strlen) redirected to 0x48311d0 (_vgnU_ifunc_wrapper)
--177410-- REDIR: 0x4eefd20 (libc.so.6:memrchr) redirected to 0x48311d0 (_vgnU_ifunc_wrapper)
--177410-- REDIR: 0x4ee79d0 (libc.so.6:strcasecmp_l) redirected to 0x48311d0 (_vgnU_ifunc_wrapper)
--177410-- REDIR: 0x4ee7550 (libc.so.6:memchr) redirected to 0x48311d0 (_vgnU_ifunc_wrapper)
--177410-- REDIR: 0x4f03ab0 (libc.so.6:wcslen) redirected to 0x48311d0 (_vgnU_ifunc_wrapper)
--177410-- REDIR: 0x4ee6c60 (libc.so.6:strspn) redirected to 0x48311d0 (_vgnU_ifunc_wrapper)
--177410-- REDIR: 0x4ee78d0 (libc.so.6:stpncpy) redirected to 0x48311d0 (_vgnU_ifunc_wrapper)
--177410-- REDIR: 0x4ee7870 (libc.so.6:stpcpy) redirected to 0x48311d0 (_vgnU_ifunc_wrapper)
--177410-- REDIR: 0x4ee8e10 (libc.so.6:strchrnul) redirected to 0x48311d0 (_vgnU_ifunc_wrapper)
--177410-- REDIR: 0x4ee7a20 (libc.so.6:strncasecmp_l) redirected to 0x48311d0 (_vgnU_ifunc_wrapper)
--177410-- REDIR: 0x4fcf490 (libc.so.6:__strrchr_avx2) redirected to 0x483ea10 (rindex)
--177410-- REDIR: 0x4ee2c90 (libc.so.6:calloc) redirected to 0x483dce0 (calloc)
--177410-- REDIR: 0x4ee1260 (libc.so.6:malloc) redirected to 0x483b780 (malloc)
--177410-- REDIR: 0x4fd2670 (libc.so.6:__memcpy_avx_unaligned_erms) redirected to 0x48429f0 (memmove)
--177410-- REDIR: 0x4fcf660 (libc.so.6:__strlen_avx2) redirected to 0x483ef40 (strlen)
--177410-- REDIR: 0x4fcf800 (libc.so.6:__strnlen_avx2) redirected to 0x483eee0 (strnlen)
--177410-- REDIR: 0x4ee1850 (libc.so.6:free) redirected to 0x483c9d0 (free)
--177410-- REDIR: 0x4fcd6b0 (libc.so.6:__strncasecmp_avx) redirected to 0x483f910 (strncasecmp)


    F I N I T E   E L E M E N T   A N A L Y S I S   P R O G R A M

--177410-- REDIR: 0x4fd2af0 (libc.so.6:__memset_avx2_unaligned_erms) redirected to 0x48428e0 (memset)
           FEAP (C) Regents of the University of California
                         All Rights Reserved.
                       VERSION: Release 8.6.1e     
                          DATE: 10 October 2020     

         Files are set as:   Status  Filename

           Input   (read ) : Exists  Isolid.txt                     
           Output  (write) : Exists  Osolid.txt                     
           Restart (read ) : New     Rsolid.txt                     
           Restart (write) : New     Rsolid.txt                     
           Plots   (write) : New     Psolid.txt                     

         Caution, existing write files will be overwritten.

         Are filenames correct?( y or n; r = redefine all; s = stop) : y
--177410-- REDIR: 0x4fcb4c0 (libc.so.6:__memchr_avx2) redirected to 0x4840050 (memchr)

         R U N N I N G    F E A P    P R O B L E M    N O W

          --> Post errors at forum: feap.berkeley.edu/forum
                 or send e-mail to: feap@berkeley.edu

--177410-- REDIR: 0x4ee7120 (libc.so.6:__GI_strstr) redirected to 0x4843ca0 (__strstr_sse2)
--177410-- REDIR: 0x4fcbc50 (libc.so.6:__memcmp_avx2_movbe) redirected to 0x48421e0 (bcmp)
--177410-- REDIR: 0x4fcab60 (libc.so.6:__strcmp_avx2) redirected to 0x483fed0 (strcmp)
--177410-- REDIR: 0x4fcf2a0 (libc.so.6:__strchrnul_avx2) redirected to 0x4843540 (strchrnul)
--177410-- REDIR: 0x4fcafa0 (libc.so.6:__strncmp_avx2) redirected to 0x483f670 (strncmp)

     T w o   D i m e n s i o n a l   S o l i d   E l e m e n t

 2D Quadratic element with linear shape functions
--177410-- REDIR: 0x4fd2650 (libc.so.6:__mempcpy_avx_unaligned_erms) redirected to 0x4843660 (mempcpy)
  E          :    200.00000000000000     
  nu         :   0.29999999999999999     
 Thickness   :    1.0000000000000000     
 xi:   -0.57735026918962573     
 xi:    0.57735026918962573     
 xi:    0.57735026918962573     
 xi:   -0.57735026918962573     
==177410== Conditional jump or move depends on uninitialised value(s)
==177410==    at 0x4BFA2D1: ??? (in /usr/lib/x86_64-linux-gnu/libgfortran.so.5.0.0)
==177410==    by 0x4BFBC46: ??? (in /usr/lib/x86_64-linux-gnu/libgfortran.so.5.0.0)
==177410==    by 0x4BF2410: ??? (in /usr/lib/x86_64-linux-gnu/libgfortran.so.5.0.0)
==177410==    by 0x4BF27E3: ??? (in /usr/lib/x86_64-linux-gnu/libgfortran.so.5.0.0)
==177410==    by 0x389B63: p_dof_norm_ (p_dof_norm.f:31)
==177410==    by 0x35839A: pmacr1_ (pmacr1.f:567)
==177410==    by 0x17C462: pmacr_ (pmacr.f:578)
==177410==    by 0x10E97B: pcontr_ (pcontr.f:1307)
==177410==    by 0x10B912: feap (feap86.f:166)
==177410==    by 0x10B912: main (feap86.f:172)
==177410==
==177410== Conditional jump or move depends on uninitialised value(s)
==177410==    at 0x4BF8989: ??? (in /usr/lib/x86_64-linux-gnu/libgfortran.so.5.0.0)
==177410==    by 0x4BFBC46: ??? (in /usr/lib/x86_64-linux-gnu/libgfortran.so.5.0.0)
==177410==    by 0x4BF2410: ??? (in /usr/lib/x86_64-linux-gnu/libgfortran.so.5.0.0)
==177410==    by 0x4BF27E3: ??? (in /usr/lib/x86_64-linux-gnu/libgfortran.so.5.0.0)
==177410==    by 0x389B63: p_dof_norm_ (p_dof_norm.f:31)
==177410==    by 0x35839A: pmacr1_ (pmacr1.f:567)
==177410==    by 0x17C462: pmacr_ (pmacr.f:578)
==177410==    by 0x10E97B: pcontr_ (pcontr.f:1307)
==177410==    by 0x10B912: feap (feap86.f:166)
==177410==    by 0x10B912: main (feap86.f:172)
==177410==
==177410== Conditional jump or move depends on uninitialised value(s)
==177410==    at 0x4BF8B6B: ??? (in /usr/lib/x86_64-linux-gnu/libgfortran.so.5.0.0)
==177410==    by 0x4BFBC46: ??? (in /usr/lib/x86_64-linux-gnu/libgfortran.so.5.0.0)
==177410==    by 0x4BF2410: ??? (in /usr/lib/x86_64-linux-gnu/libgfortran.so.5.0.0)
==177410==    by 0x4BF27E3: ??? (in /usr/lib/x86_64-linux-gnu/libgfortran.so.5.0.0)
==177410==    by 0x389B63: p_dof_norm_ (p_dof_norm.f:31)
==177410==    by 0x35839A: pmacr1_ (pmacr1.f:567)
==177410==    by 0x17C462: pmacr_ (pmacr.f:578)
==177410==    by 0x10E97B: pcontr_ (pcontr.f:1307)
==177410==    by 0x10B912: feap (feap86.f:166)
==177410==    by 0x10B912: main (feap86.f:172)
==177410==
==177410== Conditional jump or move depends on uninitialised value(s)
==177410==    at 0x1AAA9E: naninfck_ (naninfck.f:46)
==177410==    by 0x4773DD: p_res_norm_ (p_res_norm.f:49)
==177410==    by 0x3548BB: pmacr1_ (pmacr1.f:569)
==177410==    by 0x17C462: pmacr_ (pmacr.f:578)
==177410==    by 0x10E97B: pcontr_ (pcontr.f:1307)
==177410==    by 0x10B912: feap (feap86.f:166)
==177410==    by 0x10B912: main (feap86.f:172)
==177410==
  *ERROR* Residual norm is NaN or Inf
 --> ERRORS OCCURRED: For details see file: Osolid.txt

==177410==
==177410== HEAP SUMMARY:
==177410==     in use at exit: 0 bytes in 0 blocks
==177410==   total heap usage: 4,706 allocs, 4,706 frees, 3,728,185 bytes allocated
==177410==
==177410== All heap blocks were freed -- no leaks are possible
==177410==
==177410== Use --track-origins=yes to see where uninitialised values come from
==177410== ERROR SUMMARY: 13 errors from 4 contexts (suppressed: 0 from 0)
==177410==
==177410== 1 errors in context 1 of 4:
==177410== Conditional jump or move depends on uninitialised value(s)
==177410==    at 0x1AAA9E: naninfck_ (naninfck.f:46)
==177410==    by 0x4773DD: p_res_norm_ (p_res_norm.f:49)
==177410==    by 0x3548BB: pmacr1_ (pmacr1.f:569)
==177410==    by 0x17C462: pmacr_ (pmacr.f:578)
==177410==    by 0x10E97B: pcontr_ (pcontr.f:1307)
==177410==    by 0x10B912: feap (feap86.f:166)
==177410==    by 0x10B912: main (feap86.f:172)
==177410==
==177410==
==177410== 4 errors in context 2 of 4:
==177410== Conditional jump or move depends on uninitialised value(s)
==177410==    at 0x4BF8B6B: ??? (in /usr/lib/x86_64-linux-gnu/libgfortran.so.5.0.0)
==177410==    by 0x4BFBC46: ??? (in /usr/lib/x86_64-linux-gnu/libgfortran.so.5.0.0)
==177410==    by 0x4BF2410: ??? (in /usr/lib/x86_64-linux-gnu/libgfortran.so.5.0.0)
==177410==    by 0x4BF27E3: ??? (in /usr/lib/x86_64-linux-gnu/libgfortran.so.5.0.0)
==177410==    by 0x389B63: p_dof_norm_ (p_dof_norm.f:31)
==177410==    by 0x35839A: pmacr1_ (pmacr1.f:567)
==177410==    by 0x17C462: pmacr_ (pmacr.f:578)
==177410==    by 0x10E97B: pcontr_ (pcontr.f:1307)
==177410==    by 0x10B912: feap (feap86.f:166)
==177410==    by 0x10B912: main (feap86.f:172)
==177410==
==177410==
==177410== 4 errors in context 3 of 4:
==177410== Conditional jump or move depends on uninitialised value(s)
==177410==    at 0x4BF8989: ??? (in /usr/lib/x86_64-linux-gnu/libgfortran.so.5.0.0)
==177410==    by 0x4BFBC46: ??? (in /usr/lib/x86_64-linux-gnu/libgfortran.so.5.0.0)
==177410==    by 0x4BF2410: ??? (in /usr/lib/x86_64-linux-gnu/libgfortran.so.5.0.0)
==177410==    by 0x4BF27E3: ??? (in /usr/lib/x86_64-linux-gnu/libgfortran.so.5.0.0)
==177410==    by 0x389B63: p_dof_norm_ (p_dof_norm.f:31)
==177410==    by 0x35839A: pmacr1_ (pmacr1.f:567)
==177410==    by 0x17C462: pmacr_ (pmacr.f:578)
==177410==    by 0x10E97B: pcontr_ (pcontr.f:1307)
==177410==    by 0x10B912: feap (feap86.f:166)
==177410==    by 0x10B912: main (feap86.f:172)
==177410==
==177410==
==177410== 4 errors in context 4 of 4:
==177410== Conditional jump or move depends on uninitialised value(s)
==177410==    at 0x4BFA2D1: ??? (in /usr/lib/x86_64-linux-gnu/libgfortran.so.5.0.0)
==177410==    by 0x4BFBC46: ??? (in /usr/lib/x86_64-linux-gnu/libgfortran.so.5.0.0)
==177410==    by 0x4BF2410: ??? (in /usr/lib/x86_64-linux-gnu/libgfortran.so.5.0.0)
==177410==    by 0x4BF27E3: ??? (in /usr/lib/x86_64-linux-gnu/libgfortran.so.5.0.0)
==177410==    by 0x389B63: p_dof_norm_ (p_dof_norm.f:31)
==177410==    by 0x35839A: pmacr1_ (pmacr1.f:567)
==177410==    by 0x17C462: pmacr_ (pmacr.f:578)
==177410==    by 0x10E97B: pcontr_ (pcontr.f:1307)
==177410==    by 0x10B912: feap (feap86.f:166)
==177410==    by 0x10B912: main (feap86.f:172)
==177410==
==177410== ERROR SUMMARY: 13 errors from 4 contexts (suppressed: 0 from 0)


I am not able to figure out what making the program ill. Please your suggestion will help me alot.

[Prof. S. Govindjee]

The underflow/denormal warning can be ignored, that is not your main issue.

For the debugging, try write(*,*).

The main issue is that your code is producing NaN and Inf.  There can be lots of reasons for this, but most likely it is a coding error in your element.  You should first confirm that FEAP works with its own built in elements, so as to confirm that your input file is correct.  Then you can move on to debugging your element.

There is also way to have the program halt when it computes the NaN/Inf for the first time to make it easier to figure out where the problem is coming from, but I forget the details.  If you google it you should be able to learn how and then you can use gdb to exactly locate the problem.

[ad10]

Prof. Thanks a lot for replying back.

The user subroutine works perfectly fine (without any NaN and Inf), when i don't use any write/print statements during loop.

I get error when I use them whether it write(7,*) or any in code. I am not able to find reason behind that.

[Prof. S. Govindjee]

That is an indication that there is a memory overwrite error in your code.  When you have the write statement the compiler is producing different object code that is exposing the memory error.  These kind of errors are really hard to track down.   If you are lucky you can figure out where in your element you are writing out of bounds, just by proofreading the code.  If not, you will need to back track from the NaN/Inf residual to figure out what is being over written, you can then use gdb and the watch feature to identify the line of code that is writing out of bounds.

[ad10]

Thanks a lot professor for response.

I do figure out the location from my code starts deviating (from where Nan is coming).

[ad10]

Dear all,

I have a one issue here. I am trying to learn the GDB to debug the FEAP. but's only working on feap86.f not on the other corresponding files. I am interested to know the flow of program so I see where it has went wrong. Can you please suggest me how to see the actual path followed by program through gdb?

Code: [Select]

Copyright (C) 2021 Free Software Foundation, Inc.
License GPLv3+: GNU GPL version 3 or later <http://gnu.org/licenses/gpl.html>
This is free software: you are free to change and redistribute it.
There is NO WARRANTY, to the extent permitted by law.
Type "show copying" and "show warranty" for details.
This GDB was configured as "x86_64-linux-gnu".
Type "show configuration" for configuration details.
For bug reporting instructions, please see:
<https://www.gnu.org/software/gdb/bugs/>.
Find the GDB manual and other documentation resources online at:
    <http://www.gnu.org/software/gdb/documentation/>.

For help, type "help".
Type "apropos word" to search for commands related to "word"...
Reading symbols from ./feap...
(gdb) start
Temporary breakpoint 1 at 0x384a: file feap86.f, line 116.
Starting program: /home/ajinkya/Desktop/Solid/feap

Temporary breakpoint 1, feap () at feap86.f:116
116       versn(1) = 'Release 8.6.1e'
(gdb)
(gdb)
(gdb) n
131       ior = 15
(gdb) n
116       versn(1) = 'Release 8.6.1e'
(gdb) n
131       ior = 15
(gdb) n
127       iop = 11
(gdb) n
129       ird = 13
(gdb) n
134       icl = 18
(gdb) n
149       prompt = .true. ! Prompt for graphics inputs when defalt = .false.
(gdb) n
116       versn(1) = 'Release 8.6.1e'
(gdb) n
117       versn(2) = '10 October 2020'
(gdb) n
149       prompt = .true. ! Prompt for graphics inputs when defalt = .false.
(gdb) n
163       call pstart()
(gdb) n
117       versn(2) = '10 October 2020'
(gdb) n
122       ipr = 2             ! ipr = 2 for double length real to integers
(gdb) n
126       ilg = 10
(gdb) n
138       coflg = .true.  ! Parse all input as expressions   (slower mode)
(gdb) n
153       pscolr = .true. ! PostScript outputs are in color
(gdb) n
158       hlplev = 0      ! Basic
(gdb) n
163       call pstart()
(gdb) n


    F I N I T E   E L E M E N T   A N A L Y S I S   P R O G R A M

           FEAP (C) Regents of the University of California
                         All Rights Reserved.
                       VERSION: Release 8.6.1e     
                          DATE: 10 October 2020     

         Files are set as:   Status  Filename

           Input   (read ) : Exists  Isolid.txt                     
           Output  (write) : Exists  Osolid.txt                     
           Restart (read ) : New     Rsolid.txt                     
           Restart (write) : New     Rsolid.txt                     
           Plots   (write) : New     Psolid.txt                     

         Caution, existing write files will be overwritten.

         Are filenames correct?( y or n; r = redefine all; s = stop) : y

         R U N N I N G    F E A P    P R O B L E M    N O W

          --> Post errors at forum: feap.berkeley.edu/forum
                 or send e-mail to: feap@berkeley.edu

166       call pcontr()
(gdb)

     T w o   D i m e n s i o n a l   S o l i d   E l e m e n t

 2D Quadratic element with linear shape functions
  E          :    200.00000000000000     
  nu         :   0.29999999999999999     
 Thickness   :    1.0000000000000000     
 Input a solution command: Enter "help" for list of commands.
 Enter "exit" to end with restart save, "quit" to end without restart save.
   List 10  Command  1>
   List 10  Command  1>
   List 10  Command  1>
   List 10  Command  1>
   List 10  Command  1>
   List 10  Command  1>



[Prof. S. Govindjee]

A useful thing to do is to put a break point in a routine that you are interested to know how FEAP gets there.  Then you can issue the command WHERE

Code: [Select]

(gdb) break pmacr1
(gdb) run
Then when FEAP stops and you are back to the GDB prompt try

Code: [Select]

(gdb) where

I suggest looking for some online GDB tutorials, looking for a local class at your university, or trying to find someone that can show you how to use gdb.  It is enormously helpful, so it is worth your time to learn it.

[ad10]

Thanks for quick reply.

Prof. I am facing issue here. My formulated element should work each and every time without errors but why it is happening like sometimes its works and sometimes not. you can see the below output it gives me correct values. But for inbuilt FEAP functions it gives correct each time. Is it a programming error?


I have followed the same steps given in user manual. For your reference I have also attached my element.

Same goes to stress plot if I add the Plot,stre,1 in Input it gives error. Any comments on this will be really helpful.

Code: [Select]

ajinkya@ajinkya-HP-Pavilion-Gaming-Laptop-15-dk0xxx:~/Desktop/Solid$ ./feap


    F I N I T E   E L E M E N T   A N A L Y S I S   P R O G R A M

           FEAP (C) Regents of the University of California
                         All Rights Reserved.
                       VERSION: Release 8.6.1e     
                          DATE: 10 October 2020     

         Files are set as:   Status  Filename

           Input   (read ) : Exists  Isolid.txt                     
           Output  (write) : Exists  Osolid.txt                     
           Restart (read ) : New     Rsolid.txt                     
           Restart (write) : New     Rsolid.txt                     
           Plots   (write) : New     Psolid.txt                     

         Caution, existing write files will be overwritten.

         Are filenames correct?( y or n; r = redefine all; s = stop) : y

         R U N N I N G    F E A P    P R O B L E M    N O W

          --> Post errors at forum: feap.berkeley.edu/forum
                 or send e-mail to: feap@berkeley.edu


     T w o   D i m e n s i o n a l   S o l i d   E l e m e n t

 2D Quadratic element with linear shape functions
  E          :    200.00000000000000     
  nu         :   0.29999999999999999     
 Thickness   :    1.0000000000000000     
 Input a solution command: Enter "help" for list of commands.
 Enter "exit" to end with restart save, "quit" to end without restart save.
   List  8  Command  1> ^C
ajinkya@ajinkya-HP-Pavilion-Gaming-Laptop-15-dk0xxx:~/Desktop/Solid$ ./feap


    F I N I T E   E L E M E N T   A N A L Y S I S   P R O G R A M

           FEAP (C) Regents of the University of California
                         All Rights Reserved.
                       VERSION: Release 8.6.1e     
                          DATE: 10 October 2020     

         Files are set as:   Status  Filename

           Input   (read ) : Exists  Isolid.txt                     
           Output  (write) : Exists  Osolid.txt                     
           Restart (read ) : New     Rsolid.txt                     
           Restart (write) : New     Rsolid.txt                     
           Plots   (write) : New     Psolid.txt                     

         Caution, existing write files will be overwritten.

         Are filenames correct?( y or n; r = redefine all; s = stop) : y

         R U N N I N G    F E A P    P R O B L E M    N O W

          --> Post errors at forum: feap.berkeley.edu/forum
                 or send e-mail to: feap@berkeley.edu


     T w o   D i m e n s i o n a l   S o l i d   E l e m e n t

 2D Quadratic element with linear shape functions
  E          :    200.00000000000000     
  nu         :   0.29999999999999999     
 Thickness   :    1.0000000000000000     
  *ERROR* Residual norm is NaN or Inf
 --> ERRORS OCCURRED: For details see file: Osolid.txt

Note: The following floating-point exceptions are signalling: IEEE_OVERFLOW_FLAG IEEE_UNDERFLOW_FLAG IEEE_DENORMAL


[Prof. R.L. Taylor]

I get a compile error that you have too many arguments on the call to pzero for Cmat - at line 388

This could cause the erratic behavior.

[Prof. R.L. Taylor]

there is a second error.  The routine "pinput" is a logical function used as

setval = pinput(.......

[ad10]

Thanks professor. I have changed those. But I think there are some more issues. I will be careful for these errors!