Hi,
I'm trying to run a parallel simulation with the inflowWithPressureAndTangentialVelocityGiven boundary condition. I recive PETSc error, so I run gdb and it address me to the Boundary Conditions
classes (insBC.C, applyBoundaryConditions.C, ...).
NOTE1: if I substitute the bc with inflowWithVelocityGiven the code works in
parallel.
NOTE2: for both boundary conditions it runs fine in serial.
Is this actually a bug or some conditions are not implemented for parallel?
I report here the output from the debugger.
Regards,
Alessandro
Cgins
...
Parameters::updateTurbulenceModels: turbulenceModel=0
interpolateAndApplyBoundaryConditions (start) steps=-1 t=0: number of array
ID's has increased to 2276
interpolateAndApplyBoundaryConditions (after interpolate): number of array
ID's has increased to 2300
Program received signal SIGSEGV, Segmentation fault.
0x00002aaab10f9385 in
Array_Domain_Type::Error_Checking_For_Index_Operators(Internal_Index* const*)
const ()
from /home/orchini/SW/A++P++-0.8.0/P++/install/lib/libPpp.so
(gdb) backtrace
#0 0x00002aaab10f9385 in
Array_Domain_Type::Error_Checking_For_Index_Operators(Internal_Index* const*)
const ()
from /home/orchini/SW/A++P++-0.8.0/P++/install/lib/libPpp.so
#1 0x00002aaab11093b1 in doubleArray::operator()(Internal_Index const&, Internal_Index
const&, Internal_Index const&, Internal_Index const&) const ()
from /home/orchini/SW/A++P++-0.8.0/P++/install/lib/libPpp.so
#2 0x00002aaaad56a902 in MappedGridOperators::applyBCtangentialComponent(doubleMappedGridFunction&, int, int, Index const&, BCTypes::BCNames const&, int const&, double const&, doubleSerialArray
const&, doubleSerialArray const&, doubleMappedGridFunction const&, double const&, intSerialArray const&, intSerialArray const&, intArray const&, BoundaryConditionParameters const&,
MappedGridOperators::BoundaryConditionOption, int cons) () from /home/orchini/SW/v25/Overture.v24g/lib/libOverture.so
#3 0x00002aaaac55a37b in MappedGridOperators::applyBoundaryCondition(doubleMappedGridFunction&, Index const&, BCTypes::BCNames const&, int const&, double const&, doubleSerialArray const&,
doubleSerialArray const&, doubleSerialArray* (*) [3], doubleMappedGridFunction const&, double const&, BoundaryConditionParameters const&, MappedGridOperators::BoundaryConditionOption, int const&) ()
from /home/orchini/SW/v25/Overture.v24g/lib/libOverture.so
#4 0x00002aaaac55ad06 in MappedGridOperators::applyBoundaryCondition(doubleMappedGridFunction&, Index const&, BCTypes::BCNames const&, int const&, double const&, double const&,
BoundaryConditionParameters const&, int const&) () from /home/orchini/SW/v25/Overture.v24g/lib/libOverture.so
#5 0x00002aaaac9b9e11 in doubleMappedGridFunction::applyBoundaryCondition(Index const&, BCTypes::BCNames const&, int const&, double const&, double const&, BoundaryConditionParameters const&, int
const&) () from /home/orchini/SW/v25/Overture.v24g/lib/libOverture.so
#6 0x00002aaaaabe9456 in Cgins::applyBoundaryConditions (this=0x10624c0,
t=@0x10660b8, u=..., gridVelocity=..., grid=@0x7ffffffed8ac,
option=@0x7ffffffee02c,
puOld=0x0, pGridVelocityOld=0x0, dt=@0x7ffffffed908) at src/insBC.C:2348
#7 0x00002aaaabe691e4 in DomainSolver::applyBoundaryConditions
(this=0x10624c0, cgf=..., option=@0x7ffffffee02c, grid_=-1, puOld=0x0,
dt=@0x7ffffffee668)
at src/applyBoundaryConditions.C:69
#8 0x00002aaaabe606c0 in DomainSolver::interpolateAndApplyBoundaryConditions
(this=0x10624c0, cgf=..., uOld=0x0, dt=@0x7ffffffee668) at src/advance.C:777
#9 0x00002aaaabddff1a in DomainSolver::initializeSolution (this=0x10624c0) at
src/setupGridFunctions.C:370
#10 0x00002aaaaacfb84d in Cgins::initializeSolution (this=0x10624c0) at
src/setupGridFunctions.C:86
#11 0x00002aaaabec4836 in DomainSolver::setup (this=0x10624c0, time=@0x1d77e80)
at src/DomainSolver.C:554
#12 0x00002aaaabcfe083 in DomainSolver::setParametersInteractively
(this=0x10624c0, runSetupOnExit=true) at src/setup.C:2545
#13 0x0000000000418d17 in main (argc=3, argv=0x7fffffffe698) at
src/cginsMain.C:148
