@@ -275,41 +275,54 @@ int main(int argc, char** argv) {
275275 for (const auto & block_params : interface_block_params) {
276276 std::vector<int > IDs;
277277 const std::string& block_name = block_params.first ;
278+ std::cout << " [dbg] Processing block: " << block_name << " \n " ; flush_now ();
278279 interface.get_block_node_IDs (block_name, IDs);
279280 if (IDs.size () < 2 ) throw std::runtime_error (" IDs too small for " + block_name);
280281 int num_inlet_nodes = IDs[0 ];
281282 int num_outlet_nodes = IDs[1 + num_inlet_nodes * 2 ];
283+ std::cout << " [dbg] " << block_name << " : inlets=" << num_inlet_nodes
284+ << " , outlets=" << num_outlet_nodes << " \n " ; flush_now ();
282285 if (block_name == " outlet_aorta" ) {
283- if ((num_inlet_nodes != 1 ) || (num_outlet_nodes != 0 )) {
284- throw std::runtime_error (" Wrong number of inlets/outlets for outlet_aorta" );
285- }
286+ if ((num_inlet_nodes != 1 ) || (num_outlet_nodes != 0 )) {
287+ throw std::runtime_error (" Wrong number of inlets/outlets for outlet_aorta" );
288+ }
286289 } else {
287- if ((num_inlet_nodes != 0 ) || (num_outlet_nodes != 1 )) {
288- throw std::runtime_error (" Wrong number of inlets/outlets for " + block_name);
289- }
290+ if ((num_inlet_nodes != 0 ) || (num_outlet_nodes != 1 )) {
291+ throw std::runtime_error (" Wrong number of inlets/outlets for " + block_name);
292+ }
290293 }
291294 }
292295 // Outlet from heart block
293296 std::vector<int > IDs;
294297 std::string block_name = " J_heart_outlet" ;
298+ std::cout << " [dbg] Processing heart outlet block\n " ; flush_now ();
295299 interface.get_block_node_IDs (block_name, IDs);
296300 int num_inlet_nodes = IDs[0 ];
297301 int num_outlet_nodes = IDs[1 + num_inlet_nodes * 2 ];
302+ std::cout << " [dbg] Heart outlet: inlets=" << num_inlet_nodes
303+ << " , outlets=" << num_outlet_nodes << " \n " ; flush_now ();
298304 if ((num_inlet_nodes != 1 ) && (num_outlet_nodes != 1 )) {
299305 throw std::runtime_error (" Wrong number of inlets/outlets for J_heart_outlet" );
300306 }
301307 int aortic_inlet_flow_id = IDs[1 ];
302308 int aortic_inlet_pressure_id = IDs[2 ];
309+ std::cout << " [dbg] Aortic IDs: flow=" << aortic_inlet_flow_id
310+ << " , pressure=" << aortic_inlet_pressure_id << " \n " ; flush_now ();
303311 // Outlet from coronary
304312 block_name = " BC_lca1" ;
313+ std::cout << " [dbg] Processing coronary block\n " ; flush_now ();
305314 interface.get_block_node_IDs (block_name, IDs);
306315 num_inlet_nodes = IDs[0 ];
307316 num_outlet_nodes = IDs[1 + num_inlet_nodes * 2 ];
317+ std::cout << " [dbg] Coronary: inlets=" << num_inlet_nodes
318+ << " , outlets=" << num_outlet_nodes << " \n " ; flush_now ();
308319 if ((num_inlet_nodes != 1 ) && (num_outlet_nodes != 1 )) {
309320 throw std::runtime_error (" Wrong number of inlets/outlets for BC_lca1" );
310321 }
311322 int bc_lca1_outlet_flow_id = IDs[4 ];
312323 int bc_lca1_outlet_pressure_id = IDs[5 ];
324+ std::cout << " [dbg] Coronary IDs: flow=" << bc_lca1_outlet_flow_id
325+ << " , pressure=" << bc_lca1_outlet_pressure_id << " \n " ; flush_now ();
313326 std::cout << " [ok ] block node IDs\n " ; flush_now ();
314327
315328 // Update block parameters with current flow from 3D solver
0 commit comments