Commit 07304b3f authored by m.tavelli's avatar m.tavelli

Bug fix, seems to do some time steps

parent 34bc6724
......@@ -98,13 +98,14 @@ void FOCCZ4::FOCCZ4Solver_ADERDG::boundaryValues(const double* const x,const dou
std::memset(fluxOut , 0, nVar * sizeof(double));
std::copy_n(stateIn,nVar,stateOut);
for(int dd=0; dd<nDim; dd++) F[dd] = Fs[dd];
std::copy_n(fluxIn,nVar,fluxOut);
/*for(int dd=0; dd<nDim; dd++) F[dd] = Fs[dd];
for(int i=0; i < basisSize; i++) { // i == time
const double weight = kernels::legendre::weights[order][i];
const double xi = kernels::legendre::nodes[order][i];
double ti = t + xi * dt;
// initialdata_(x, &ti, Qgp,&md,&cms,&order);
initialdata_(x_3, &ti, Qgp);
flux(Qgp, F);
......@@ -112,7 +113,7 @@ void FOCCZ4::FOCCZ4Solver_ADERDG::boundaryValues(const double* const x,const dou
stateOut[m] += weight * Qgp[m];
//fluxOut[m] += weight * Fs[direction][m];
}
}
}*/
}
exahype::solvers::Solver::RefinementControl FOCCZ4::FOCCZ4Solver_ADERDG::refinementCriterion(const double* const luh,const tarch::la::Vector<DIMENSIONS,double>& cellCentre,const tarch::la::Vector<DIMENSIONS,double>& cellSize,double t,const int level) {
......
......@@ -25,6 +25,8 @@ RECURSIVE SUBROUTINE InitParameters(STRLEN,PARSETUP)
EQN%CCZ4GLMepsP = 1.0 ! 5.
EQN%CCZ4GLMepsD = 1.0 ! 0.1
!
EQN%CCZ4itau = 0.0
EQN%CCZ4k1 = 0.0
EQN%CCZ4k2 = 0.0
EQN%CCZ4k3 = 0.0
......@@ -68,7 +70,7 @@ RECURSIVE SUBROUTINE InitialData(xGP, tGp, Q)
!
r = SQRT( xGP(1)**2 + xGP(2)**2 )
!
!V0(60) = 0.001*EXP(-0.5*( (xGP(1)-2.0)**2+xGP(2)**2+0*xGP(3)**2)/1.0**2 ) + 0.001*EXP(-0.5*( (xGP(1)+2.0)**2+xGP(2)**2+0*xGP(3)**2)/1.0**2 )
V0(60) = 0.001*EXP(-0.5*( (xGP(1)-2.0)**2+xGP(2)**2+0*xGP(3)**2)/1.0**2 ) + 0.001*EXP(-0.5*( (xGP(1)+2.0)**2+xGP(2)**2+0*xGP(3)**2)/1.0**2 )
IF( r>5.0 .AND. r<10. ) THEN
V0(61) = -0.2*xGP(2)*( 10.0 - r )/5.0
V0(62) = +0.2*xGP(1)*( 10.0 - r )/5.0
......@@ -79,7 +81,7 @@ RECURSIVE SUBROUTINE InitialData(xGP, tGp, Q)
V0(61:62) = 0.0
ENDIF
V0(63) = 0.0
!V0(64) = 0.5*0.001*EXP(-0.5*( (xGP(1)-2.0)**2+xGP(2)**2+0*xGP(3)**2)/1.0**2 ) + 0.5*0.001*EXP(-0.5*( (xGP(1)+2.0)**2+xGP(2)**2+0*xGP(3)**2)/1.0**2 )
V0(64) = 0.5*0.001*EXP(-0.5*( (xGP(1)-2.0)**2+xGP(2)**2+0*xGP(3)**2)/1.0**2 ) + 0.5*0.001*EXP(-0.5*( (xGP(1)+2.0)**2+xGP(2)**2+0*xGP(3)**2)/1.0**2 )
!
CASE DEFAULT
PRINT *, 'NOT IMPLEMENTED'
......
......@@ -2605,10 +2605,10 @@ RECURSIVE SUBROUTINE PDEFusedSrcNCP(Src_BgradQ,Q,gradQin)
!print *, 'g_cov',g_cov
!print *, 'DD',DD
!print *, 'PP',PP
print *, 'mytmp1',mytmp1
print *, 'mytmp2',mytmp2
print *, 'Christoffel',Christoffel
pause
!print *, 'mytmp1',mytmp1
!print *, 'mytmp2',mytmp2
!print *, 'Christoffel',Christoffel
!pause
DO l = 1, 3
DO j = 1, 3
......@@ -2772,13 +2772,13 @@ RECURSIVE SUBROUTINE PDEFusedSrcNCP(Src_BgradQ,Q,gradQin)
!print *, 'dtraceK=', dtraceK
!print *, 'Theta=', Theta
!print *, 'Christoffel',Christoffel
print *, 'Christoffel',Christoffel
!print *, 'Christoffel',Christoffel
!print *, 'Gtilde',Gtilde
!print *, 'Christoffel_tilde',Christoffel_tilde
!print *, 'dChristoffelsrc',dChristoffelsrc
!print *, 'dChristoffel_tildesrc',dChristoffel_tildesrc
!print *, 'Z',Z
pause
!pause
traceB = BB(1,1) + BB(2,2) + BB(3,3)
dtphi = beta(1)*PP(1) + beta(2)*PP(2) + beta(3)*PP(3) + 1./3.*alpha*traceK - 1./3.*traceB
......
# MPI Logfilter.
# black or white list entry
# -1 means all ranks
# default entries
debug tarch -1 black
debug peano -1 black
info tarch -1 black
info tarch::logging::CommandLineLogger -1 white
info peano -1 black
info peano::utils::UserInterface -1 white
info peano::performanceanalysis -1 white
info exahype -1 white
# on first node
info mpibalancing -1 black
# too many messages
info peano::parallel::SendReceiveBufferAbstractImplementation::releaseSentMessages -1 black
info sharedmemoryoracles -1 black
info Z4Solver -1 black
info exahype::solvers::Plotter -1 black
info exahype::runner::Runner::runAsWorker -1 black
info exahype::mappings::LoadBalancing -1 white
info exahype::runners::Runner::runAsWorker -1 black
info exahype::runners::Runner::init -1 black
info exahype::runners::Runner::createRepository -1 white
info exahype::mappings::LoadBalancing::endIteration -1 black
......@@ -105,7 +105,7 @@
"type": "user::defined",
"name": "TecplotWriter",
"time": 0.0,
"repeat": 0.1,
"repeat": 1.0,
"output": "./output/tecplot",
"variables": 96
}
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment