CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Programming & Development

Problem with run continue of a custom solver

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   January 23, 2013, 09:08
Default Problem with run continue of a custom solver
  #1
Senior Member
 
Hisham's Avatar
 
Hisham Elsafti
Join Date: Apr 2011
Location: Braunschweig, Germany
Posts: 257
Blog Entries: 10
Rep Power: 16
Hisham is on a distinguished road
Dear Foamers,

I am working with a custom solver. Every time I try to continue a run from a previous time step (not the initial), the GAMG (and other, also increased tolerance) solver fails (serial and parallel). The failure occurs at the statement:

Code:
initialResidual = DEqn.solve().initialResidual();
I checked all fields involved in the equation and they are all read properly. I also turned on the debug switch for GAMG but that didn't provide me with useful info. The error is:

Code:
   Normalisation factor = 800.459
DICPCG:  Solving for coarsestLevelCorr, Initial residual = 1, Final residual = 0.21207, No Iterations 1
Post-smoothing scaling factors: 1.03931 1.11069 1.09697 1.38127 1.15639
GAMG:  Solving for Dx, Initial residual = 0.864838, Final residual = 0.22546, No Iterations 0
GAMG:  Solving for Dx, Initial residual = 0.864838, Final residual = 0.22546, No Iterations 1
   Normalisation factor = 2.02601e+247
#0  Foam::error::printStack(Foam::Ostream&) in "/home/hisham/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#1  Foam::sigFpe::sigHandler(int) in "/home/hisham/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#2   in "/lib/x86_64-linux-gnu/libc.so.6"
#3  Foam::PCG::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const in "/home/hisham/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#4  Foam::GAMGSolver::solveCoarsestLevel(Foam::Field<double>&, Foam::Field<double> const&) const in "/home/hisham/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#5  Foam::GAMGSolver::Vcycle(Foam::PtrList<Foam::lduMatrix::smoother> const&, Foam::Field<double>&, Foam::Field<double> const&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::PtrList<Foam::Field<double> >&, Foam::PtrList<Foam::Field<double> >&, unsigned char) const in "/home/hisham/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#6  Foam::GAMGSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const in "/home/hisham/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#7  
 in "/home/hisham/OpenFOAM/hisham-2.1.x/platforms/linux64GccDPOpt/bin/geotechFoam3"
#8  
 in "/home/hisham/OpenFOAM/hisham-2.1.x/platforms/linux64GccDPOpt/bin/geotechFoam3"
#9  
 in "/home/hisham/OpenFOAM/hisham-2.1.x/platforms/linux64GccDPOpt/bin/geotechFoam3"
#10  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#11  
 in "/home/hisham/OpenFOAM/hisham-2.1.x/platforms/linux64GccDPOpt/bin/geotechFoam3"
The solver on the other hand continues the run if the command is issued with valgrind. Valgrind complaints are:


Code:
  
Normalisation factor = 800.459
==7950== Conditional jump or move depends on uninitialised value(s)
==7950==    at 0x6957954: Foam::lduMatrix::solverPerformance::checkConvergence(double, double) (in /home/hisham/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libOpenFOAM.so)
==7950==    by 0x6973AAC: Foam::GAMGSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const (in /home/hisham/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libOpenFOAM.so)
==7950==    by 0x482A60: Foam::fvMatrix<Foam::Vector<double> >::solve(Foam::dictionary const&) (in /home/hisham/OpenFOAM/hisham-2.1.x/platforms/linux64GccDPOpt/bin/geotechFoam3)
==7950==    by 0x48351F: Foam::fvMatrix<Foam::Vector<double> >::solve() (in /home/hisham/OpenFOAM/hisham-2.1.x/platforms/linux64GccDPOpt/bin/geotechFoam3)
==7950==    by 0x42FA16: main (in /home/hisham/OpenFOAM/hisham-2.1.x/platforms/linux64GccDPOpt/bin/geotechFoam3)
==7950== 
==7950== Conditional jump or move depends on uninitialised value(s)
==7950==    at 0x6973AAF: Foam::GAMGSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const (in /home/hisham/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libOpenFOAM.so)
==7950==    by 0x482A60: Foam::fvMatrix<Foam::Vector<double> >::solve(Foam::dictionary const&) (in /home/hisham/OpenFOAM/hisham-2.1.x/platforms/linux64GccDPOpt/bin/geotechFoam3)
==7950==    by 0x48351F: Foam::fvMatrix<Foam::Vector<double> >::solve() (in /home/hisham/OpenFOAM/hisham-2.1.x/platforms/linux64GccDPOpt/bin/geotechFoam3)
==7950==    by 0x42FA16: main (in /home/hisham/OpenFOAM/hisham-2.1.x/platforms/linux64GccDPOpt/bin/geotechFoam3)
==7950== 
==7950== Conditional jump or move depends on uninitialised value(s)
==7950==    at 0x6957954: Foam::lduMatrix::solverPerformance::checkConvergence(double, double) (in /home/hisham/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libOpenFOAM.so)
==7950==    by 0x6960E64: Foam::PCG::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const (in /home/hisham/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libOpenFOAM.so)
==7950==    by 0x6970E9A: Foam::GAMGSolver::solveCoarsestLevel(Foam::Field<double>&, Foam::Field<double> const&) const (in /home/hisham/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libOpenFOAM.so)
==7950==    by 0x69726D9: Foam::GAMGSolver::Vcycle(Foam::PtrList<Foam::lduMatrix::smoother> const&, Foam::Field<double>&, Foam::Field<double> const&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::PtrList<Foam::Field<double> >&, Foam::PtrList<Foam::Field<double> >&, unsigned char) const (in /home/hisham/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libOpenFOAM.so)
==7950==    by 0x6973B62: Foam::GAMGSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const (in /home/hisham/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libOpenFOAM.so)
==7950==    by 0x482A60: Foam::fvMatrix<Foam::Vector<double> >::solve(Foam::dictionary const&) (in /home/hisham/OpenFOAM/hisham-2.1.x/platforms/linux64GccDPOpt/bin/geotechFoam3)
==7950==    by 0x48351F: Foam::fvMatrix<Foam::Vector<double> >::solve() (in /home/hisham/OpenFOAM/hisham-2.1.x/platforms/linux64GccDPOpt/bin/geotechFoam3)
==7950==    by 0x42FA16: main (in /home/hisham/OpenFOAM/hisham-2.1.x/platforms/linux64GccDPOpt/bin/geotechFoam3)
==7950== 
==7950== Conditional jump or move depends on uninitialised value(s)
==7950==    at 0x6960E67: Foam::PCG::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const (in /home/hisham/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libOpenFOAM.so)
==7950==    by 0x6970E9A: Foam::GAMGSolver::solveCoarsestLevel(Foam::Field<double>&, Foam::Field<double> const&) const (in /home/hisham/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libOpenFOAM.so)
==7950==    by 0x69726D9: Foam::GAMGSolver::Vcycle(Foam::PtrList<Foam::lduMatrix::smoother> const&, Foam::Field<double>&, Foam::Field<double> const&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::PtrList<Foam::Field<double> >&, Foam::PtrList<Foam::Field<double> >&, unsigned char) const (in /home/hisham/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libOpenFOAM.so)
==7950==    by 0x6973B62: Foam::GAMGSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const (in /home/hisham/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libOpenFOAM.so)
==7950==    by 0x482A60: Foam::fvMatrix<Foam::Vector<double> >::solve(Foam::dictionary const&) (in /home/hisham/OpenFOAM/hisham-2.1.x/platforms/linux64GccDPOpt/bin/geotechFoam3)
==7950==    by 0x48351F: Foam::fvMatrix<Foam::Vector<double> >::solve() (in /home/hisham/OpenFOAM/hisham-2.1.x/platforms/linux64GccDPOpt/bin/geotechFoam3)
==7950==    by 0x42FA16: main (in /home/hisham/OpenFOAM/hisham-2.1.x/platforms/linux64GccDPOpt/bin/geotechFoam3)
==7950== 
==7950== Conditional jump or move depends on uninitialised value(s)
==7950==    at 0x6961055: Foam::PCG::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const (in /home/hisham/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libOpenFOAM.so)
==7950==    by 0x6970E9A: Foam::GAMGSolver::solveCoarsestLevel(Foam::Field<double>&, Foam::Field<double> const&) const (in /home/hisham/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libOpenFOAM.so)
==7950==    by 0x69726D9: Foam::GAMGSolver::Vcycle(Foam::PtrList<Foam::lduMatrix::smoother> const&, Foam::Field<double>&, Foam::Field<double> const&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::PtrList<Foam::Field<double> >&, Foam::PtrList<Foam::Field<double> >&, unsigned char) const (in /home/hisham/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libOpenFOAM.so)
==7950==    by 0x6973B62: Foam::GAMGSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const (in /home/hisham/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libOpenFOAM.so)
==7950==    by 0x482A60: Foam::fvMatrix<Foam::Vector<double> >::solve(Foam::dictionary const&) (in /home/hisham/OpenFOAM/hisham-2.1.x/platforms/linux64GccDPOpt/bin/geotechFoam3)
==7950==    by 0x48351F: Foam::fvMatrix<Foam::Vector<double> >::solve() (in /home/hisham/OpenFOAM/hisham-2.1.x/platforms/linux64GccDPOpt/bin/geotechFoam3)
==7950==    by 0x42FA16: main (in /home/hisham/OpenFOAM/hisham-2.1.x/platforms/linux64GccDPOpt/bin/geotechFoam3)
==7950== 
==7950== Conditional jump or move depends on uninitialised value(s)
==7950==    at 0x6957954: Foam::lduMatrix::solverPerformance::checkConvergence(double, double) (in /home/hisham/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libOpenFOAM.so)
==7950==    by 0x69613EF: Foam::PCG::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const (in /home/hisham/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libOpenFOAM.so)
==7950==    by 0x6970E9A: Foam::GAMGSolver::solveCoarsestLevel(Foam::Field<double>&, Foam::Field<double> const&) const (in /home/hisham/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libOpenFOAM.so)
==7950==    by 0x69726D9: Foam::GAMGSolver::Vcycle(Foam::PtrList<Foam::lduMatrix::smoother> const&, Foam::Field<double>&, Foam::Field<double> const&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::PtrList<Foam::Field<double> >&, Foam::PtrList<Foam::Field<double> >&, unsigned char) const (in /home/hisham/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libOpenFOAM.so)
==7950==    by 0x6973B62: Foam::GAMGSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const (in /home/hisham/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libOpenFOAM.so)
==7950==    by 0x482A60: Foam::fvMatrix<Foam::Vector<double> >::solve(Foam::dictionary const&) (in /home/hisham/OpenFOAM/hisham-2.1.x/platforms/linux64GccDPOpt/bin/geotechFoam3)
==7950==    by 0x48351F: Foam::fvMatrix<Foam::Vector<double> >::solve() (in /home/hisham/OpenFOAM/hisham-2.1.x/platforms/linux64GccDPOpt/bin/geotechFoam3)
==7950==    by 0x42FA16: main (in /home/hisham/OpenFOAM/hisham-2.1.x/platforms/linux64GccDPOpt/bin/geotechFoam3)
==7950== 
==7950== Conditional jump or move depends on uninitialised value(s)
==7950==    at 0x69613F2: Foam::PCG::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const (in /home/hisham/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libOpenFOAM.so)
==7950==    by 0x6970E9A: Foam::GAMGSolver::solveCoarsestLevel(Foam::Field<double>&, Foam::Field<double> const&) const (in /home/hisham/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libOpenFOAM.so)
==7950==    by 0x69726D9: Foam::GAMGSolver::Vcycle(Foam::PtrList<Foam::lduMatrix::smoother> const&, Foam::Field<double>&, Foam::Field<double> const&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::PtrList<Foam::Field<double> >&, Foam::PtrList<Foam::Field<double> >&, unsigned char) const (in /home/hisham/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libOpenFOAM.so)
==7950==    by 0x6973B62: Foam::GAMGSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const (in /home/hisham/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libOpenFOAM.so)
==7950==    by 0x482A60: Foam::fvMatrix<Foam::Vector<double> >::solve(Foam::dictionary const&) (in /home/hisham/OpenFOAM/hisham-2.1.x/platforms/linux64GccDPOpt/bin/geotechFoam3)
==7950==    by 0x48351F: Foam::fvMatrix<Foam::Vector<double> >::solve() (in /home/hisham/OpenFOAM/hisham-2.1.x/platforms/linux64GccDPOpt/bin/geotechFoam3)
==7950==    by 0x42FA16: main (in /home/hisham/OpenFOAM/hisham-2.1.x/platforms/linux64GccDPOpt/bin/geotechFoam3)
==7950== 
==7950== Conditional jump or move depends on uninitialised value(s)
==7950==    at 0x6957A09: Foam::lduMatrix::solverPerformance::print() const (in /home/hisham/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libOpenFOAM.so)
==7950==    by 0x6970D51: Foam::GAMGSolver::solveCoarsestLevel(Foam::Field<double>&, Foam::Field<double> const&) const (in /home/hisham/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libOpenFOAM.so)
==7950==    by 0x69726D9: Foam::GAMGSolver::Vcycle(Foam::PtrList<Foam::lduMatrix::smoother> const&, Foam::Field<double>&, Foam::Field<double> const&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::PtrList<Foam::Field<double> >&, Foam::PtrList<Foam::Field<double> >&, unsigned char) const (in /home/hisham/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libOpenFOAM.so)
==7950==    by 0x6973B62: Foam::GAMGSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const (in /home/hisham/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libOpenFOAM.so)
==7950==    by 0x482A60: Foam::fvMatrix<Foam::Vector<double> >::solve(Foam::dictionary const&) (in /home/hisham/OpenFOAM/hisham-2.1.x/platforms/linux64GccDPOpt/bin/geotechFoam3)
==7950==    by 0x48351F: Foam::fvMatrix<Foam::Vector<double> >::solve() (in /home/hisham/OpenFOAM/hisham-2.1.x/platforms/linux64GccDPOpt/bin/geotechFoam3)
==7950==    by 0x42FA16: main (in /home/hisham/OpenFOAM/hisham-2.1.x/platforms/linux64GccDPOpt/bin/geotechFoam3)
==7950== 
==7950== Use of uninitialised value of size 8
==7950==    at 0x788D9E7: __printf_fp (printf_fp.c:828)
==7950==    by 0x788B35F: vfprintf (vfprintf.c:2007)
==7950==    by 0x78B1441: vsnprintf (vsnprintf.c:120)
==7950==    by 0x70AF29D: ??? (in /usr/lib/x86_64-linux-gnu/libstdc++.so.6.0.16)
==7950==    by 0x70B4399: std::ostreambuf_iterator<char, std::char_traits<char> > std::num_put<char, std::ostreambuf_iterator<char, std::char_traits<char> > >::_M_insert_float<double>(std::ostreambuf_iterator<char, std::char_traits<char> >, std::ios_base&, char, char, double) const (in /usr/lib/x86_64-linux-gnu/libstdc++.so.6.0.16)
==7950==    by 0x70B47D8: std::num_put<char, std::ostreambuf_iterator<char, std::char_traits<char> > >::do_put(std::ostreambuf_iterator<char, std::char_traits<char> >, std::ios_base&, char, double) const (in /usr/lib/x86_64-linux-gnu/libstdc++.so.6.0.16)
==7950==    by 0x70C104D: std::ostream& std::ostream::_M_insert<double>(double) (in /usr/lib/x86_64-linux-gnu/libstdc++.so.6.0.16)
==7950==    by 0x68CD00C: Foam::OSstream::write(double) (in /home/hisham/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libOpenFOAM.so)
==7950==    by 0x68791BC: Foam::operator<<(Foam::Ostream&, double) (in /home/hisham/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libOpenFOAM.so)
==7950==    by 0x6957A4C: Foam::lduMatrix::solverPerformance::print() const (in /home/hisham/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libOpenFOAM.so)
==7950==    by 0x6970D51: Foam::GAMGSolver::solveCoarsestLevel(Foam::Field<double>&, Foam::Field<double> const&) const (in /home/hisham/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libOpenFOAM.so)
==7950==    by 0x69726D9: Foam::GAMGSolver::Vcycle(Foam::PtrList<Foam::lduMatrix::smoother> const&, Foam::Field<double>&, Foam::Field<double> const&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::PtrList<Foam::Field<double> >&, Foam::PtrList<Foam::Field<double> >&, unsigned char) const (in /home/hisham/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libOpenFOAM.so)
==7950== 
==7950== Use of uninitialised value of size 8
==7950==    at 0x7876458: __mpn_lshift (lshift.S:60)
==7950==    by 0x788D9EB: __printf_fp (printf_fp.c:828)
==7950==    by 0x788B35F: vfprintf (vfprintf.c:2007)
==7950==    by 0x78B1441: vsnprintf (vsnprintf.c:120)
==7950==    by 0x70AF29D: ??? (in /usr/lib/x86_64-linux-gnu/libstdc++.so.6.0.16)
==7950==    by 0x70B4399: std::ostreambuf_iterator<char, std::char_traits<char> > std::num_put<char, std::ostreambuf_iterator<char, std::char_traits<char> > >::_M_insert_float<double>(std::ostreambuf_iterator<char, std::char_traits<char> >, std::ios_base&, char, char, double) const (in /usr/lib/x86_64-linux-gnu/libstdc++.so.6.0.16)
==7950==    by 0x70B47D8: std::num_put<char, std::ostreambuf_iterator<char, std::char_traits<char> > >::do_put(std::ostreambuf_iterator<char, std::char_traits<char> >, std::ios_base&, char, double) const (in /usr/lib/x86_64-linux-gnu/libstdc++.so.6.0.16)
==7950==    by 0x70C104D: std::ostream& std::ostream::_M_insert<double>(double) (in /usr/lib/x86_64-linux-gnu/libstdc++.so.6.0.16)
==7950==    by 0x68CD00C: Foam::OSstream::write(double) (in /home/hisham/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libOpenFOAM.so)
==7950==    by 0x68791BC: Foam::operator<<(Foam::Ostream&, double) (in /home/hisham/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libOpenFOAM.so)
==7950==    by 0x6957A4C: Foam::lduMatrix::solverPerformance::print() const (in /home/hisham/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libOpenFOAM.so)
==7950==    by 0x6970D51: Foam::GAMGSolver::solveCoarsestLevel(Foam::Field<double>&, Foam::Field<double> const&) const (in /home/hisham/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libOpenFOAM.so)
==7950== 
==7950== Use of uninitialised value of size 8
==7950==    at 0x787645B: __mpn_lshift (lshift.S:61)
==7950==    by 0x788D9EB: __printf_fp (printf_fp.c:828)
==7950==    by 0x788B35F: vfprintf (vfprintf.c:2007)
==7950==    by 0x78B1441: vsnprintf (vsnprintf.c:120)
==7950==    by 0x70AF29D: ??? (in /usr/lib/x86_64-linux-gnu/libstdc++.so.6.0.16)
==7950==    by 0x70B4399: std::ostreambuf_iterator<char, std::char_traits<char> > std::num_put<char, std::ostreambuf_iterator<char, std::char_traits<char> > >::_M_insert_float<double>(std::ostreambuf_iterator<char, std::char_traits<char> >, std::ios_base&, char, char, double) const (in /usr/lib/x86_64-linux-gnu/libstdc++.so.6.0.16)
==7950==    by 0x70B47D8: std::num_put<char, std::ostreambuf_iterator<char, std::char_traits<char> > >::do_put(std::ostreambuf_iterator<char, std::char_traits<char> >, std::ios_base&, char, double) const (in /usr/lib/x86_64-linux-gnu/libstdc++.so.6.0.16)
==7950==    by 0x70C104D: std::ostream& std::ostream::_M_insert<double>(double) (in /usr/lib/x86_64-linux-gnu/libstdc++.so.6.0.16)
....
So it seems that some Field of size 8 is uninitialised however the run continues. So does anyone have a clue to what may cause this bug?

Best regards,
Hisham
Hisham is offline   Reply With Quote

Old   January 23, 2013, 10:36
Default
  #2
Senior Member
 
Hisham's Avatar
 
Hisham Elsafti
Join Date: Apr 2011
Location: Braunschweig, Germany
Posts: 257
Blog Entries: 10
Rep Power: 16
Hisham is on a distinguished road
I have rerun valgrind using the option --track-origins=yes
and I have a better idea now where the error comes from

Still if any has an advise ... I am all ears/eyes!
Thanks
Hisham is offline   Reply With Quote

Reply

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
new Solver won't run parallel Chris Lucas OpenFOAM 4 January 10, 2012 11:30
Problem with Phoenics solver Alan Phoenics 1 April 6, 2009 17:02
compressible two phase flow in CFX4.4 youngan CFX 0 July 2, 2003 00:32
Mashing Problem? ji CFX 2 April 4, 2002 06:03
CFX 5.5 Roued CFX 1 October 2, 2001 17:49


All times are GMT -4. The time now is 07:54.