Loading src/operators/diffusion/nonlinear-diffusion-operators/tnlFiniteVolumeNonlinearOperator_impl.h +1 −1 Original line number Diff line number Diff line Loading @@ -294,7 +294,7 @@ updateLinearSystem( const RealType& time, operatorQ.getValue(mesh, index, coordinates, u, time, 1, 0, 0 ); const RealType fCoef = - tau * operatorQ.getValue(mesh, index, coordinates, u, time ) * mesh.getHySquareInverse() / operatorQ.getValue(mesh, index, coordinates, u, time, 0, 1, 0 ); const RealType gCoef = - tau * operatorQ.getValue(mesh, index, coordinates, u, time * mesh.getHzSquareInverse() / const RealType gCoef = - tau * operatorQ.getValue(mesh, index, coordinates, u, time ) * mesh.getHzSquareInverse() / operatorQ.getValue(mesh, index, coordinates, u, time, 0, 0, 1 ); matrixRow.setElement( 0, mesh.template getCellNextToCell< 0,0,-1 >( index ), aCoef ); matrixRow.setElement( 1, mesh.template getCellNextToCell< 0,-1,0 >( index ), bCoef ); Loading Loading
src/operators/diffusion/nonlinear-diffusion-operators/tnlFiniteVolumeNonlinearOperator_impl.h +1 −1 Original line number Diff line number Diff line Loading @@ -294,7 +294,7 @@ updateLinearSystem( const RealType& time, operatorQ.getValue(mesh, index, coordinates, u, time, 1, 0, 0 ); const RealType fCoef = - tau * operatorQ.getValue(mesh, index, coordinates, u, time ) * mesh.getHySquareInverse() / operatorQ.getValue(mesh, index, coordinates, u, time, 0, 1, 0 ); const RealType gCoef = - tau * operatorQ.getValue(mesh, index, coordinates, u, time * mesh.getHzSquareInverse() / const RealType gCoef = - tau * operatorQ.getValue(mesh, index, coordinates, u, time ) * mesh.getHzSquareInverse() / operatorQ.getValue(mesh, index, coordinates, u, time, 0, 0, 1 ); matrixRow.setElement( 0, mesh.template getCellNextToCell< 0,0,-1 >( index ), aCoef ); matrixRow.setElement( 1, mesh.template getCellNextToCell< 0,-1,0 >( index ), bCoef ); Loading