Skip to content
[Jenkins]Jenkins
log in

Valgrind Error Details

back to process details
back to processes overview
Executable/home1/jenkins/workspace/ProjectX_Weekly/PX_ARC/relwithdebinfo_gnu/buildnode/reynolds/platforms/relwithdebinfo_gnu/bin/ProjectXExec
PID26196
PPID16382
Unique ID0xf77
Text3,568 bytes in 1 blocks are still reachable in loss record 3,958 of 4,065
Auxiliary
Suppression{
   <insert_a_suppression_name_here>
   Memcheck:Leak
   match-leak-kinds: reachable
   fun:malloc
   fun:PXAllocate
   fun:PXPrecalcElemMasterPolynomial
   fun:PXGetMasterPolynomialElem
   fun:PXElementGalerkinPrimalResidual
   fun:PXElementGalerkinResidual
   fun:PXGalerkinResidual
   fun:PXCalculateResidual
   fun:PXSolveNonlinearPrimalNewton
   fun:PXSolveNonlinearProblem
   fun:PXSolveUnsteadyStep
   fun:PXSolveUnsteadyProblem
   fun:PXSolvePrimal
   fun:PXPSequencing
   fun:main
}

Stacktrace

Object/usr/lib/valgrind/vgpreload_memcheck-amd64-linux.so
Functionmalloc
File/Line
CodeSource code not available

Object/home1/jenkins/workspace/ProjectX_Weekly/PX_ARC/relwithdebinfo_gnu/buildnode/reynolds/platforms/relwithdebinfo_gnu/bin/ProjectXExec
FunctionPXAllocate
File/Line/home1/jenkins/workspace/ProjectX_Weekly/PX_ARC/relwithdebinfo_gnu/buildnode/reynolds/src/Fundamentals/PXMemory.c:99
Code
089 
090 
091   tot = n*size;
092   
093   if (likely(tot)){
094     if (unlikely(n<0)){
095       printf("Error, requesting allocation of negative memory size.\n"); fflush(stdout);
096       *po = NULL;
097       return PXError(PX_MEMORY_ERROR);
098     }
099     else if ( unlikely(( _po = (char *)malloc(tot)) == NULL) ){
100       *po = NULL;
101       return PXError(PX_MEMORY_ERROR);
102     }
103   }
104   *po = (void *) _po;

Object/home1/jenkins/workspace/ProjectX_Weekly/PX_ARC/relwithdebinfo_gnu/buildnode/reynolds/platforms/relwithdebinfo_gnu/bin/ProjectXExec
FunctionPXPrecalcElemMasterPolynomial
File/Line/home1/jenkins/workspace/ProjectX_Weekly/PX_ARC/relwithdebinfo_gnu/buildnode/reynolds/src/Reference/PXMasterPolynomial.c:137
Code
127   int porder;
128   enum PXE_Shape shape;
129   enum PXE_BASISTYPE basisType;
130   PX_REAL * RESTRICT phi = NULL, * RESTRICT gphi = NULL;
131 
132   PXOrder2nbf(order, &nbf);
133   PXOrder2porder(order, &porder);
134   PXGetOrderBasisType(order, &basisType);
135   PXOrder2Shape(order, &shape);
136 
137   PXErrorReturn( PXAllocate( sizeof(PX_MasterPolynomial)+(Dim+1)*nbf*nquad, sizeof(PX_REAL), (void **)&(*masterPoly) ) ); 
138   phi = (PX_REAL *)  ((*masterPoly) + 1);
139   gphi = phi + nbf*nquad;
140 
141   (**masterPoly).phi = phi;
142   (**masterPoly).order = order;

Object/home1/jenkins/workspace/ProjectX_Weekly/PX_ARC/relwithdebinfo_gnu/buildnode/reynolds/platforms/relwithdebinfo_gnu/bin/ProjectXExec
FunctionPXGetMasterPolynomialElem
File/Line/home1/jenkins/workspace/ProjectX_Weekly/PX_ARC/relwithdebinfo_gnu/buildnode/reynolds/src/Reference/PXMasterPolynomial.c:329
Code
319 }
320 
321 int
322 PXGetMasterPolynomialElem(PX_PolynomialPrecalcStruct * PolyPrecalcStruct, int quad_order, enum PXE_Shape shape, enum PXE_QuadratureRule quadRule, enum PXE_SolutionOrder order, PX_MasterPolynomial const ** MasterPolyData){
323   if( unlikely(PolyPrecalcStruct->MasterPolynomialElemLookup[shape][quadRule][order] == NULL) ){
324     /* make sure quadrature rules exist */
325     if( unlikely(PolyPrecalcStruct->QuadElemLookup[shape][quadRule] == NULL) ){
326       PXComputeQuadElemPrecalc(PolyPrecalcStruct->Dim, quad_order, shape, &(PolyPrecalcStruct->QuadElemLookup[shape][quadRule]));
327     }
328 
329     PXErrorReturn( PXPrecalcElemMasterPolynomial(PolyPrecalcStruct->Dim, order, PolyPrecalcStruct->QuadElemLookup[shape][quadRule]->nquad, PolyPrecalcStruct->QuadElemLookup[shape][quadRule]->xquad, &PolyPrecalcStruct->MasterPolynomialElemLookup[shape][quadRule][order]) );
330   }
331 
332   if(MasterPolyData != NULL)
333     *MasterPolyData = PolyPrecalcStruct->MasterPolynomialElemLookup[shape][quadRule][order];
334 

Object/home1/jenkins/workspace/ProjectX_Weekly/PX_ARC/relwithdebinfo_gnu/buildnode/reynolds/platforms/relwithdebinfo_gnu/bin/ProjectXExec
FunctionPXElementGalerkinPrimalResidual
File/Line/home1/jenkins/workspace/ProjectX_Weekly/PX_ARC/relwithdebinfo_gnu/buildnode/reynolds/src/PXResidual.c:2041
Code
02031     PXErrorReturn( PXQuadReference(shape, quad_order, &nquad, &quad_rule, NULL, NULL) );
02032 
02033     PXErrorReturn( PXGetMasterPolynomialElemWithQuad(PolynomialPrecalcStruct_GLOBAL, quad_order, shape, quad_rule, order, &QuadElemData, &MasterPolyData) );
02034     phiStart = MasterPolyData->phi;
02035     gphiStart = MasterPolyData->phi + nbf*nquad;
02036 
02037     PXErrorReturn( PXGetMasterPolynomialElem(PolynomialPrecalcStruct_GLOBAL, quad_order, shape, quad_rule, orderQ, &MasterPolyData) );
02038     phiQStart = MasterPolyData->phi;
02039     gphiQStart = MasterPolyData->phi + nbfQ*nquad;
02040 
02041     PXErrorReturn( PXGetMasterPolynomialElem(PolynomialPrecalcStruct_GLOBAL, quad_order, shape, quad_rule, orderH, &MasterPolyData) );
02042     phiHStart = MasterPolyData->phi;
02043 
02044     /*-----------------------*/
02045     /* Get quadrature points */
02046     /*-----------------------*/

Object/home1/jenkins/workspace/ProjectX_Weekly/PX_ARC/relwithdebinfo_gnu/buildnode/reynolds/platforms/relwithdebinfo_gnu/bin/ProjectXExec
FunctionPXElementGalerkinResidual
File/Line/home1/jenkins/workspace/ProjectX_Weekly/PX_ARC/relwithdebinfo_gnu/buildnode/reynolds/src/PXResidual.c:2730
Code
02720   if (AdjointFlag == PXE_True || TangentFlag == PXE_True){
02721     PXErrorReturn( PXAllocate( Rank*Rank, sizeof(PX_REAL), (void **)&(pR_U) ) );
02722     for (k = 0; k < Rank*Rank; k++) pR_U[k] = 0.0;
02723   }
02724   else
02725     pR_U = R_U; /* Set pointer to R_U to R_U */
02726 
02727   /*************************************************************/
02728 
02729   /* Evaluate Primal Residual and Jacobian */
02730   PXErrorReturn( PXElementGalerkinPrimalResidual(pxa, ResJobParams, ResElemData, egrp, elem, order, CurrentTime, U, R, pR_U) );
02731 
02732   /*************************************************************/
02733 
02734   /*--------------------------*/
02735   /* Adjoint/Tangent Specific */

Object/home1/jenkins/workspace/ProjectX_Weekly/PX_ARC/relwithdebinfo_gnu/buildnode/reynolds/platforms/relwithdebinfo_gnu/bin/ProjectXExec
FunctionPXGalerkinResidual
File/Line/home1/jenkins/workspace/ProjectX_Weekly/PX_ARC/relwithdebinfo_gnu/buildnode/reynolds/src/PXResidual.c:2928
Code
02918       /*-------------------*/
02919       /* Evaluate Residual */
02920       /*-------------------*/
02921 
02922       /* Get pointer into Jacobian Matrix */
02923       if (DGSolverStruct != NULL) {
02924         PXErrorReturn( PXDGSolver_GetMatPointer(DGSolverStruct, pg, egrp, elem, -1, &JM) );
02925       }
02926 
02927       /* Evaluate Residual for an element */
02928       PXErrorReturn( PXElementGalerkinResidual(pxa, ResJobParams, &ResElemData, egrp, elem, order, CurrentTime,
02929                                                State->value[egrp][elem], Adj, Tan,
02930                                                R, AdjR, TanR, JM) );
02931 
02932       /* Insert values into jacobian matrix */
02933       if (DGSolverStruct != NULL) {

Object/home1/jenkins/workspace/ProjectX_Weekly/PX_ARC/relwithdebinfo_gnu/buildnode/reynolds/platforms/relwithdebinfo_gnu/bin/ProjectXExec
FunctionPXCalculateResidual
File/Line/home1/jenkins/workspace/ProjectX_Weekly/PX_ARC/relwithdebinfo_gnu/buildnode/reynolds/src/PXResidual.c:16319
Code
16309   //PXErrorReturn( PXPrecalcElemThings(pxa, State->order) );
16310 #if PREALLOC == 1
16311   PXErrorReturn( PXPreAllocateThings(pxa, &ResJobParams, State->order, -1, -1, calcJacobianFlag) );
16312   //printf("prealloc_Galerkin_Flag = %d\n",prealloc_Galerkin_Flag);
16313 #endif
16314 
16315   /* Element residuals */
16316 #ifdef PX_TIMING
16317   c2 = PXclock();
16318 #endif
16319   PXErrorReturn( PXGalerkinResidual(pxa, &ResJobParams, State, AdjState, TanState,
16320                                     Res, AdjRes, TanRes, DGSolverStruct));
16321 #ifdef PX_TIMING
16322   c3 = PXclock();
16323 #endif
16324 

Object/home1/jenkins/workspace/ProjectX_Weekly/PX_ARC/relwithdebinfo_gnu/buildnode/reynolds/platforms/relwithdebinfo_gnu/bin/ProjectXExec
FunctionPXSolveNonlinearPrimalNewton
File/Line/home1/jenkins/workspace/ProjectX_Weekly/PX_ARC/relwithdebinfo_gnu/buildnode/reynolds/src/PXSolver.c:6898
Code
06888     /*--------------------------*/
06889     /* Compute Spatial Residual */
06890     /*--------------------------*/
06891 
06892     if ( unlikely(PingFlag == PXE_PingLocal) ) {
06893       PXErrorReturn(PXPingResidualFast(pxa, State));
06894     }else if( unlikely(PingFlag == PXE_PingGlobal) ){
06895       PXErrorReturn(PXPingResidualFull(pxa, State));
06896     }
06897 
06898     terr =  PXCalculateResidual(pxa, mode, State, NULL, NULL, Residual, NULL, NULL,
06899                                 DGSolverStruct);
06900     /* Check that the flow state does not need to be rewound */
06901     PXErrorReturn( PXCheckRewind( pxa, terr, UpdateState, State, ValidState,
06902                                   NULL, NULL, &CFL, &CFLValid,
06903                                   &RewindFlag ) );

Object/home1/jenkins/workspace/ProjectX_Weekly/PX_ARC/relwithdebinfo_gnu/buildnode/reynolds/platforms/relwithdebinfo_gnu/bin/ProjectXExec
FunctionPXSolveNonlinearProblem
File/Line/home1/jenkins/workspace/ProjectX_Weekly/PX_ARC/relwithdebinfo_gnu/buildnode/reynolds/src/PXSolver.c:8862
Code
08852     /* Set number of adjoint outputs */
08853     nSensParam = pxa->SensitivityStruct->nSensParam;
08854     /* Allocate adjoint residual norm array */
08855     PXErrorReturn( PXAllocate( nSensParam, sizeof(PX_REAL),
08856                                (void **)&(TanResidualNorm) ) );
08857   }
08858 
08859   switch ( NonlinearSolver ){
08860   case PXE_Newton:
08861     if ( FlowFlag == PXE_True ){
08862       ierr = PXSolveNonlinearPrimalNewton( pxa, MaxIter, unsteadyWeight,
08863                                            mode, CFLStart,
08864                                            &ResidualTol, &ResidualNorm,
08865                                            State0, Source, State, Residual, &smResIter);
08866       if ( ierr !=PX_NO_ERROR && ierr !=PX_NOT_CONVERGED ) PXErrorReturn(ierr);
08867     }

Object/home1/jenkins/workspace/ProjectX_Weekly/PX_ARC/relwithdebinfo_gnu/buildnode/reynolds/platforms/relwithdebinfo_gnu/bin/ProjectXExec
FunctionPXSolveUnsteadyStep
File/Line/home1/jenkins/workspace/ProjectX_Weekly/PX_ARC/relwithdebinfo_gnu/buildnode/reynolds/src/PXSolver.c:9320
Code
09310                                                            TanSpatialRes ) );
09311 
09312         /* Set flag to converged for this stage */
09313         ConvergedFlag = PXE_True;
09314       }
09315       else
09316       {
09317         /*----------------------------------------------*/
09318         /* Solve the problem using the nonlinear solver */
09319         /*----------------------------------------------*/
09320         PXErrorReturn(
09321             PXSolveNonlinearProblem( pxa, unsteadyWeight, mode,
09322                                      NULL, Source, State, StageRes[iStage],
09323                                      NULL, NULL, NULL, NULL, /* No Adjoint */
09324                                      NULL, TanSource, TanState, TanRes,
09325                                      &ResidualNorm, NULL, NULL,

Object/home1/jenkins/workspace/ProjectX_Weekly/PX_ARC/relwithdebinfo_gnu/buildnode/reynolds/platforms/relwithdebinfo_gnu/bin/ProjectXExec
FunctionPXSolveUnsteadyProblem
File/Line/home1/jenkins/workspace/ProjectX_Weekly/PX_ARC/relwithdebinfo_gnu/buildnode/reynolds/src/PXSolver.c:9855
Code
09845 
09846 
09847     }
09848     else{
09849 
09850       StageState[0] = State;
09851 
09852     }
09853 
09854     /* Call time step solver */
09855     PXErrorReturn( PXSolveUnsteadyStep(pxa, currentIndex, UnsteadyStruct,
09856                                        UnsteadyFlag, PXE_Mode_Solve,
09857                                        StageState,
09858                                        (PrevState!=NULL)?PrevState[0]:NULL,
09859                                            PrevState,
09860                                            StageResidual, Source, StateSource,

Object/home1/jenkins/workspace/ProjectX_Weekly/PX_ARC/relwithdebinfo_gnu/buildnode/reynolds/platforms/relwithdebinfo_gnu/bin/ProjectXExec
FunctionPXSolvePrimal
File/Line/home1/jenkins/workspace/ProjectX_Weekly/PX_ARC/relwithdebinfo_gnu/buildnode/reynolds/src/PXRunSolver.c:386
Code
0376   PXPrintf("\n***************\n");
0377   PXPrintf(  " Primal Solver \n");
0378   PXPrintf(  "***************\n");
0379 
0380   if ( Discretization == PXE_Discretization_HDG ) {
0381     /* HDG discretization */
0382     PXErrorReturn( PXHDGSolve(pxa, &ConvergedFlag) );
0383   }
0384   else {
0385     /* DG discretization */
0386     PXErrorReturn( PXSolveUnsteadyProblem(pxa, &ConvergedFlag) );
0387   }
0388 
0389   /*-------------------*/
0390   /* Write Output File */
0391   /*-------------------*/

Object/home1/jenkins/workspace/ProjectX_Weekly/PX_ARC/relwithdebinfo_gnu/buildnode/reynolds/platforms/relwithdebinfo_gnu/bin/ProjectXExec
FunctionPXPSequencing
File/Line/home1/jenkins/workspace/ProjectX_Weekly/PX_ARC/relwithdebinfo_gnu/buildnode/reynolds/src/PXRunSolver.c:681
Code
0671 
0672     /*----------------------*/
0673     /* Solve Primal Problem */
0674     /*----------------------*/
0675 
0676 #ifdef PX_TIMING
0677     PXBarrier();
0678     c0 = PXclock();
0679 #endif
0680 
0681     PXErrorReturn( PXSolvePrimal(pxa, AdaptFlag, aiter, &ConvergedFlag, &iSolOrder) );
0682 
0683 
0684 #ifdef PX_TIMING
0685     PXBarrier();
0686     c1 = PXclock();

Object/home1/jenkins/workspace/ProjectX_Weekly/PX_ARC/relwithdebinfo_gnu/buildnode/reynolds/platforms/relwithdebinfo_gnu/bin/ProjectXExec
Functionmain
File/Line/home1/jenkins/workspace/ProjectX_Weekly/PX_ARC/relwithdebinfo_gnu/buildnode/reynolds/src/PXRunSolver.c:999
Code
0989     PXErrorReturn( PXInitializeAdaptIter( pxa, aiter ) );
0990     //pg = pxa->pg;
0991 
0992 #ifdef PX_TIMING
0993     PXBarrier();
0994     c2=PXclock();
0995     PXPrintf("Time : (InitializeAdaptIter) = %f\n", c2 - c1 );
0996 #endif
0997 
0998     /* Perform PSequencing */
0999     PXErrorReturn( PXPSequencing( pxa, AdaptFlag, aiter, &ConvergedFlag ) );
1000 
1001 #ifdef PX_TIMING
1002     PXBarrier();
1003     c1=PXclock();
1004     PXPrintf("Time : (PXPSequencing) = %f\n",  c1 - c2);