Commit 531c0ee3 authored by Shuang Chen's avatar Shuang Chen
Browse files

update

parent 83ad83e9
......@@ -38,12 +38,12 @@ void MatrixTranslatorGeneral<ODESystemTag::FirstOrderImplicitQuasilinear>::
_time_disc.getWeightedOldX(tmp, x_prev);
// debugging
tmp.write(
"C:\\george\\PhD\\UFZ\\task\\89\\opti_BHEprocess\\debug\\weighted_old_x.txt");
"C:\\george\\PhD\\UFZ\\task\\89\\opti_BHEprocess\\debug\\weighted_pre_x.txt");
// rhs = M * weighted_old_x + b
LinAlg::matMultAdd(M, tmp, b, rhs);
// debugging
rhs.write(
"C:\\george\\PhD\\UFZ\\task\\89\\opti_BHEprocess\\debug\\global_B.txt");
"C:\\george\\PhD\\UFZ\\task\\89\\opti_BHEprocess\\debug\\global_B_before_applyBC.txt");
NumLib::GlobalVectorProvider::provider.releaseVector(tmp);
}
......
......@@ -93,7 +93,7 @@ NonlinearSolverStatus NonlinearSolver<NonlinearSolverTag::Picard>::solve(
sys.getA(A);
// debugging
A.write(
"C:\\george\\PhD\\UFZ\\task\\89\\opti_BHEprocess\\debug\\global_A.txt");
"C:\\george\\PhD\\UFZ\\task\\89\\opti_BHEprocess\\debug\\global_A_before_applyBC.txt");
sys.getRhs(*x_prev[process_id], rhs);
INFO("[time] Assembly took {:g} s.", time_assembly.elapsed());
......@@ -107,6 +107,10 @@ NonlinearSolverStatus NonlinearSolver<NonlinearSolverTag::Picard>::solve(
sys.applyKnownSolutionsPicard(A, rhs, *x_new[process_id]);
time_dirichlet += timer_dirichlet.elapsed();
INFO("[time] Applying Dirichlet BCs took {:g} s.", time_dirichlet);
A.write(
"C:\\george\\PhD\\UFZ\\task\\89\\opti_BHEprocess\\debug\\global_A_after_applyBC.txt");
rhs.write(
"C:\\george\\PhD\\UFZ\\task\\89\\opti_BHEprocess\\debug\\global_B_after_applyBC.txt");
if (!sys.isLinear() && _convergence_criterion->hasResidualCheck()) {
GlobalVector res;
......
Supports Markdown
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