diff --git a/src/language/algorithms/HeatDiamondAlgorithm.cpp b/src/language/algorithms/HeatDiamondAlgorithm.cpp index 33d6783e59be016e532441e3b6669e56c6a094c5..32cdef7ab8046e31133c0c6087bb5963fa7ea12c 100644 --- a/src/language/algorithms/HeatDiamondAlgorithm.cpp +++ b/src/language/algorithms/HeatDiamondAlgorithm.cpp @@ -573,18 +573,8 @@ HeatDiamondScheme<Dimension>::HeatDiamondScheme( } Vector<double> T{number_of_dof}; - T = 0; - { - parallel_for( - mesh->numberOfCells(), PUGS_LAMBDA(CellId cell_id) { T[cell_dof_number[cell_id]] = Tj[cell_id]; }); - parallel_for( - mesh->numberOfFaces(), PUGS_LAMBDA(FaceId face_id) { - if (primal_face_is_on_boundary[face_id]) { - T[face_dof_number[face_id]] = Tl[face_id]; - } - }); - T = 1; - } + T = 1; + BiCGStab{b, A, T, 1000, 1e-9}; CellValue<double> Temperature{mesh->connectivity()};