| 414 | } // initializePatchHierarchy |
| 415 | |
| 416 | void |
| 417 | PenaltyIBMethod::initializeLevelData(Pointer<BasePatchHierarchy<NDIM>> hierarchy, |
| 418 | int level_number, |
| 419 | double init_data_time, |
| 420 | bool can_be_refined, |
| 421 | bool initial_time, |
| 422 | Pointer<BasePatchLevel<NDIM>> old_level, |
| 423 | bool allocate_data) |
| 424 | { |
| 425 | IBMethod::initializeLevelData( |
| 426 | hierarchy, level_number, init_data_time, can_be_refined, initial_time, old_level, allocate_data); |
| 427 | |
| 428 | if (initial_time && d_l_data_manager->levelContainsLagrangianData(level_number)) |
| 429 | { |
| 430 | // Initialize Mass and Spring constant data. |
| 431 | // Position and Velocity will be copied later. |
| 432 | Pointer<LData> M_data = d_l_data_manager->createLData("M", |
| 433 | level_number, |
| 434 | 1, |
| 435 | /*manage_data*/ true); |
| 436 | Pointer<LData> K_data = d_l_data_manager->createLData("K", |
| 437 | level_number, |
| 438 | 1, |
| 439 | /*manage_data*/ true); |
| 440 | static const int global_index_offset = 0; |
| 441 | static const int local_index_offset = 0; |
| 442 | d_l_initializer->initializeMassDataOnPatchLevel(global_index_offset, |
| 443 | local_index_offset, |
| 444 | M_data, |
| 445 | K_data, |
| 446 | hierarchy, |
| 447 | level_number, |
| 448 | init_data_time, |
| 449 | can_be_refined, |
| 450 | initial_time, |
| 451 | d_l_data_manager); |
| 452 | |
| 453 | if (d_silo_writer) |
| 454 | { |
| 455 | d_silo_writer->registerVariableData("M", M_data, level_number); |
| 456 | } |
| 457 | } |
| 458 | return; |
| 459 | } |
| 460 | |
| 461 | void |
| 462 | PenaltyIBMethod::putToDatabase(Pointer<Database> db) |
nothing calls this directly
no test coverage detected