Skip to content

Commit 3b6fdd9

Browse files
dalerobertsdaleroberts
committed
Memory and performance improvements
- Matrix class: packed/symmetric storage, attach-to-mmap, calloc allocator - LM solver + Anderson acceleration CLI - Static build fixes (glibc 2.34+ pthread) - Matrix test suite
1 parent 216a1d3 commit 3b6fdd9

15 files changed

Lines changed: 3541 additions & 573 deletions

‎dynadjust/CMakeLists.txt‎

Lines changed: 55 additions & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -195,7 +195,7 @@ else()
195195
endif()
196196
message(STATUS "Found XercesC (${XercesC_INCLUDE_DIRS} ${XercesC_LIBRARIES})")
197197

198-
# Static builds need ICU for Xerces transcoding (glibc iconv can't dlopen gconv modules in static binaries)
198+
# Static builds need ICU for Xerces transcoding (glibc iconv can't dlopen gconv modules in static binaries).
199199
if(BUILD_STATIC AND UNIX AND NOT APPLE)
200200
find_package(ICU COMPONENTS uc data)
201201
if(ICU_FOUND)
@@ -1132,6 +1132,53 @@ if (BUILD_TESTING)
11321132
add_test (NAME test-urban-phased-network COMMAND $<TARGET_FILE:${DNADIFF_TARGET}> urban.phased.adj urban.phased.adj.expected --skip-to-marker "M Station 1" -t 0.001)
11331133
add_test (NAME test-urban-thread-network COMMAND $<TARGET_FILE:${DNADIFF_TARGET}> urban_mt.phased-mt.adj urban_mt.phased-mt.adj.expected --skip-to-marker "M Station 1" -t 0.01 -v)
11341134

1135+
# ........................................................................
1136+
# LM / Trust-Region solver tests
1137+
# ........................................................................
1138+
1139+
# LM-1. GNSS simultaneous with LM — should match GN results
1140+
add_test (NAME import-gnss-lm COMMAND $<TARGET_FILE:${DNAIMPORT_TARGET}> -n gnss_lm gnss-network.stn gnss-network.msr -r GDA2020 --flag-unused-stations --quiet)
1141+
add_test (NAME geoid-gnss-lm COMMAND $<TARGET_FILE:${DNAGEOID_TARGET}> gnss_lm -g gnss-network-geoid.gsb --export-dna-geo)
1142+
add_test (NAME reftran-gnss-lm COMMAND $<TARGET_FILE:${DNAREFTRAN_TARGET}> gnss_lm --export-dna --export-xml)
1143+
add_test (NAME adjust-gnss-lm COMMAND $<TARGET_FILE:${DNAADJUST_TARGET}> gnss_lm --output-adj-msr --scale-normals-to-unity --lm-enabled)
1144+
add_test (NAME test-gnss-lm-vs-expected COMMAND $<TARGET_FILE:${DNADIFF_TARGET}> gnss_lm.simult.adj gnss.simult.adj.expected --skip-headers 52 -t 0.001)
1145+
1146+
# LM-2. Urban phased with LM — should match GN results
1147+
add_test (NAME import-urban-lm COMMAND $<TARGET_FILE:${DNAIMPORT_TARGET}> -n urban_lm urban-network.stn urban-network.msr --flag-unused-stations --quiet)
1148+
add_test (NAME geoid-urban-lm COMMAND $<TARGET_FILE:${DNAGEOID_TARGET}> urban_lm -g urban-network-geoid.gsb --export-dna-geo)
1149+
add_test (NAME segment-urban-lm COMMAND $<TARGET_FILE:${DNASEGMENT_TARGET}> urban_lm --min 50 --max 150)
1150+
add_test (NAME adjust-urban-lm COMMAND $<TARGET_FILE:${DNAADJUST_TARGET}> urban_lm --output-adj-msr --phased --lm-enabled)
1151+
add_test (NAME test-urban-lm-vs-expected COMMAND $<TARGET_FILE:${DNADIFF_TARGET}> urban_lm.phased.adj urban.phased.adj.expected --skip-to-marker "M Station 1" -t 0.001)
1152+
1153+
# LM-3. Urban multi-threaded phased with LM — should match GN results
1154+
add_test (NAME import-urban-mt-lm COMMAND $<TARGET_FILE:${DNAIMPORT_TARGET}> -n urban_mt_lm urban-network.stn urban-network.msr --quiet)
1155+
add_test (NAME reftran-urban-mt-lm COMMAND $<TARGET_FILE:${DNAREFTRAN_TARGET}> urban_mt_lm -r gda2020)
1156+
add_test (NAME geoid-urban-mt-lm COMMAND $<TARGET_FILE:${DNAGEOID_TARGET}> urban_mt_lm -g urban-network-geoid.gsb)
1157+
add_test (NAME segment-urban-mt-lm COMMAND $<TARGET_FILE:${DNASEGMENT_TARGET}> urban_mt_lm --min 50 --max 85)
1158+
add_test (NAME adjust-urban-mt-lm COMMAND $<TARGET_FILE:${DNAADJUST_TARGET}> urban_mt_lm --output-adj-msr --multi --lm-enabled)
1159+
add_test (NAME test-urban-mt-lm-vs-expected COMMAND $<TARGET_FILE:${DNADIFF_TARGET}> urban_mt_lm.phased-mt.adj urban_mt.phased-mt.adj.expected --skip-to-marker "M Station 1" -t 0.01)
1160+
1161+
# LM-4. DSG distance-only simultaneous — multi-iteration stress test
1162+
add_test (NAME import-dsg-lm COMMAND $<TARGET_FILE:${DNAIMPORT_TARGET}> -n dsg_lm dsg.stn dsg.msr --flag-unused-stations --quiet)
1163+
add_test (NAME adjust-dsg-lm COMMAND $<TARGET_FILE:${DNAADJUST_TARGET}> dsg_lm --output-adj-msr --max-iterations 50 --lm-enabled)
1164+
1165+
# LM-5. DSG distance-only simultaneous — GN baseline for comparison
1166+
add_test (NAME import-dsg-gn COMMAND $<TARGET_FILE:${DNAIMPORT_TARGET}> -n dsg_gn dsg.stn dsg.msr --flag-unused-stations --quiet)
1167+
add_test (NAME adjust-dsg-gn COMMAND $<TARGET_FILE:${DNAADJUST_TARGET}> dsg_gn --output-adj-msr --max-iterations 50)
1168+
add_test (NAME test-dsg-lm-vs-gn COMMAND $<TARGET_FILE:${DNADIFF_TARGET}> dsg_lm.simult.adj dsg_gn.simult.adj --skip-headers 52 -t 0.01)
1169+
1170+
# LM-6. DSG distance-only phased with LM
1171+
add_test (NAME import-dsg-ph-lm COMMAND $<TARGET_FILE:${DNAIMPORT_TARGET}> -n dsg_ph_lm dsg.stn dsg.msr --flag-unused-stations --quiet)
1172+
add_test (NAME segment-dsg-ph-lm COMMAND $<TARGET_FILE:${DNASEGMENT_TARGET}> dsg_ph_lm --min 10 --max 30)
1173+
add_test (NAME adjust-dsg-ph-lm COMMAND $<TARGET_FILE:${DNAADJUST_TARGET}> dsg_ph_lm --output-adj-msr --phased --max-iterations 50 --lm-enabled)
1174+
1175+
# LM-7. GNSS simultaneous with large initial lambda (robustness check)
1176+
add_test (NAME import-gnss-lm-large-lambda COMMAND $<TARGET_FILE:${DNAIMPORT_TARGET}> -n gnss_lm_ll gnss-network.stn gnss-network.msr -r GDA2020 --flag-unused-stations --quiet)
1177+
add_test (NAME geoid-gnss-lm-large-lambda COMMAND $<TARGET_FILE:${DNAGEOID_TARGET}> gnss_lm_ll -g gnss-network-geoid.gsb --export-dna-geo)
1178+
add_test (NAME reftran-gnss-lm-large-lambda COMMAND $<TARGET_FILE:${DNAREFTRAN_TARGET}> gnss_lm_ll --export-dna --export-xml)
1179+
add_test (NAME adjust-gnss-lm-large-lambda COMMAND $<TARGET_FILE:${DNAADJUST_TARGET}> gnss_lm_ll --output-adj-msr --scale-normals-to-unity --lm-enabled --lm-lambda-init 10.0)
1180+
add_test (NAME test-gnss-lm-large-lambda COMMAND $<TARGET_FILE:${DNADIFF_TARGET}> gnss_lm_ll.simult.adj gnss.simult.adj.expected --skip-headers 52 -t 0.001)
1181+
11351182
# 8. Source tag preservation in --export-xml (issue #317)
11361183
# Import XML data with <Source> tags and verify they are preserved in export
11371184
add_test (NAME import-source-test COMMAND $<TARGET_FILE:${DNAIMPORT_TARGET}> -n srctest source-test-stn.xml source-test-msr.xml --export-xml -r GDA2020)
@@ -1998,6 +2045,13 @@ if (BUILD_TESTING)
19982045
set_tests_properties(test-urban-phased-network PROPERTIES DEPENDS adjust-urban-network)
19992046
set_tests_properties(test-urban-thread-network PROPERTIES DEPENDS adjust-urban-network-thread-01)
20002047

2048+
# LM test dependencies
2049+
set_tests_properties(test-gnss-lm-vs-expected PROPERTIES DEPENDS adjust-gnss-lm)
2050+
set_tests_properties(test-urban-lm-vs-expected PROPERTIES DEPENDS adjust-urban-lm)
2051+
set_tests_properties(test-urban-mt-lm-vs-expected PROPERTIES DEPENDS adjust-urban-mt-lm)
2052+
set_tests_properties(test-dsg-lm-vs-gn PROPERTIES DEPENDS "adjust-dsg-lm;adjust-dsg-gn")
2053+
set_tests_properties(test-gnss-lm-large-lambda PROPERTIES DEPENDS adjust-gnss-lm-large-lambda)
2054+
20012055
set_tests_properties(ref-itrf-pmm-06 PROPERTIES DEPENDS ref-itrf-pmm-05)
20022056
#set_tests_properties(ref-itrf-pmm-07 PROPERTIES DEPENDS ref-itrf-pmm-06)
20032057

‎dynadjust/dynadjust/dnaadjust/dnaadjust-multi.cpp‎

Lines changed: 11 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -91,6 +91,12 @@ bool prepareAdjustmentExceptionThrown(const std::vector<std::exception_ptr>& pre
9191
//
9292
void dna_adjust::AdjustPhasedMultiThread()
9393
{
94+
if (projectSettings_.a.lm_enabled)
95+
{
96+
AdjustPhasedMultiThreadLM();
97+
return;
98+
}
99+
94100
initialiseIteration();
95101

96102
std::string corr_msg;
@@ -295,7 +301,7 @@ void dna_adjust::SolveMT(bool COMPUTE_INVERSE, const UINT32& block)
295301
{
296302
// Compute inverse of normals (aposteriori variance matrix)
297303
// (AT * V-1 * A)-1
298-
FormInverseVarianceMatrix(&(v_normalsR_.at(block)), false);
304+
FormInverseVarianceMatrix(&(v_normalsR_.at(block)), false, true);
299305
}
300306

301307
// compute weighted "measured minus computed"
@@ -304,7 +310,10 @@ void dna_adjust::SolveMT(bool COMPUTE_INVERSE, const UINT32& block)
304310

305311
// Solve corrections from normal equations
306312
v_correctionsR_.at(block).redim(v_designR_.at(block).columns(), 1);
307-
v_correctionsR_.at(block).multiply(v_normalsR_.at(block), "N", At_Vinv_m, "N");
313+
if (v_normalsR_.at(block).is_symmetric())
314+
v_correctionsR_.at(block).multiply_sym(v_normalsR_.at(block), At_Vinv_m);
315+
else
316+
v_correctionsR_.at(block).multiply(v_normalsR_.at(block), "N", At_Vinv_m, "N");
308317

309318
// debug output?
310319
if (projectSettings_.g.verbose > 3)

‎dynadjust/dynadjust/dnaadjust/dnaadjust-stage.cpp‎

Lines changed: 79 additions & 39 deletions
Original file line numberDiff line numberDiff line change
@@ -173,7 +173,7 @@ void dna_adjust::DeserialiseBlockFromMappedFile(const UINT32& block, const int f
173173
break;
174174
case sf_normals_r:
175175
addr = normalsR_map_.GetBlockRegionAddr(block);
176-
v_normalsR_.at(block).ReadMappedFileRegion(addr);
176+
v_normalsR_.at(block).AttachMappedFileRegion(addr);
177177
break;
178178
case sf_atvinv:
179179
v_AtVinv_.at(block).allocate();
@@ -183,47 +183,47 @@ void dna_adjust::DeserialiseBlockFromMappedFile(const UINT32& block, const int f
183183
break;
184184
case sf_meas_minus_comp:
185185
addr = measMinusComp_map_.GetBlockRegionAddr(block);
186-
v_measMinusComp_.at(block).ReadMappedFileRegion(addr);
186+
v_measMinusComp_.at(block).AttachMappedFileRegion(addr);
187187
break;
188188
case sf_estimated_stns:
189189
addr = estimatedStations_map_.GetBlockRegionAddr(block);
190-
v_estimatedStations_.at(block).ReadMappedFileRegion(addr);
190+
v_estimatedStations_.at(block).AttachMappedFileRegion(addr);
191191
break;
192192
case sf_original_stns:
193193
addr = originalStations_map_.GetBlockRegionAddr(block);
194-
v_originalStations_.at(block).ReadMappedFileRegion(addr);
194+
v_originalStations_.at(block).AttachMappedFileRegion(addr);
195195
break;
196196
case sf_rigorous_stns:
197197
addr = rigorousStations_map_.GetBlockRegionAddr(block);
198-
v_rigorousStations_.at(block).ReadMappedFileRegion(addr);
198+
v_rigorousStations_.at(block).AttachMappedFileRegion(addr);
199199
break;
200200
case sf_junction_vars:
201201
addr = junctionVariances_map_.GetBlockRegionAddr(block);
202-
v_junctionVariances_.at(block).ReadMappedFileRegion(addr);
202+
v_junctionVariances_.at(block).AttachMappedFileRegion(addr);
203203
break;
204204
case sf_junction_vars_f:
205205
addr = junctionVariancesFwd_map_.GetBlockRegionAddr(block);
206-
v_junctionVariancesFwd_.at(block).ReadMappedFileRegion(addr);
206+
v_junctionVariancesFwd_.at(block).AttachMappedFileRegion(addr);
207207
break;
208208
case sf_junction_ests_f:
209209
addr = junctionEstimatesFwd_map_.GetBlockRegionAddr(block);
210-
v_junctionEstimatesFwd_.at(block).ReadMappedFileRegion(addr);
210+
v_junctionEstimatesFwd_.at(block).AttachMappedFileRegion(addr);
211211
break;
212212
case sf_junction_ests_r:
213213
addr = junctionEstimatesRev_map_.GetBlockRegionAddr(block);
214-
v_junctionEstimatesRev_.at(block).ReadMappedFileRegion(addr);
214+
v_junctionEstimatesRev_.at(block).AttachMappedFileRegion(addr);
215215
break;
216216
case sf_rigorous_vars:
217217
addr = rigorousVariances_map_.GetBlockRegionAddr(block);
218-
v_rigorousVariances_.at(block).ReadMappedFileRegion(addr);
218+
v_rigorousVariances_.at(block).AttachMappedFileRegion(addr);
219219
break;
220220
case sf_prec_adj_msrs:
221221
addr = precAdjMsrs_map_.GetBlockRegionAddr(block);
222-
v_precAdjMsrsFull_.at(block).ReadMappedFileRegion(addr);
222+
v_precAdjMsrsFull_.at(block).AttachMappedFileRegion(addr);
223223
break;
224224
case sf_corrections:
225225
addr = corrections_map_.GetBlockRegionAddr(block);
226-
v_corrections_.at(block).ReadMappedFileRegion(addr);
226+
v_corrections_.at(block).AttachMappedFileRegion(addr);
227227

228228
if (v_blockMeta_.at(block)._blockLast)
229229
v_correctionsR_.at(block).allocate();
@@ -240,8 +240,8 @@ void dna_adjust::SerialiseBlockToMappedFile(const UINT32& block, const int file_
240240
{
241241
// serialise all. That is, call this function again, but with
242242
// arguments for all files
243-
SerialiseBlockToMappedFile(block, 16,
244-
sf_normals, sf_normals_r, sf_atvinv, sf_design, sf_meas_minus_comp,
243+
SerialiseBlockToMappedFile(block, 12,
244+
sf_normals_r, sf_meas_minus_comp,
245245
sf_estimated_stns, sf_original_stns, sf_rigorous_stns,
246246
sf_junction_vars, sf_junction_vars_f, sf_junction_ests_f, sf_junction_ests_r,
247247
sf_rigorous_vars, sf_prec_adj_msrs, sf_corrections);
@@ -257,16 +257,10 @@ void dna_adjust::SerialiseBlockToMappedFile(const UINT32& block, const int file_
257257
{
258258
switch (va_arg(vlist, int))
259259
{
260-
case sf_normals:
261-
break;
262260
case sf_normals_r:
263261
addr = normalsR_map_.GetBlockRegionAddr(block);
264262
v_normalsR_.at(block).WriteMappedFileRegion(addr);
265263
break;
266-
case sf_atvinv:
267-
break;
268-
case sf_design:
269-
break;
270264
case sf_meas_minus_comp:
271265
addr = measMinusComp_map_.GetBlockRegionAddr(block);
272266
v_measMinusComp_.at(block).WriteMappedFileRegion(addr);
@@ -395,7 +389,11 @@ void dna_adjust::OpenStageFileStreams(const int file_count, ...)
395389
}
396390

397391
std::stringstream ss;
398-
ss << projectSettings_.g.output_folder << FOLDER_SLASH << projectSettings_.g.network_name << "-";
392+
if (projectSettings_.a.stage_path.empty())
393+
ss << projectSettings_.g.output_folder;
394+
else
395+
ss << projectSettings_.a.stage_path;
396+
ss << FOLDER_SLASH << projectSettings_.g.network_name << "-";
399397
std::string filePath(ss.str());
400398

401399
v_stageFileStreams_.clear();
@@ -851,6 +849,46 @@ void dna_adjust::OffloadBlockToMappedFile(const UINT32& block)
851849

852850
// Unload block matrix data from memory
853851
UnloadBlock(block);
852+
853+
// Advise kernel that mapped pages for this block are no longer needed,
854+
// freeing page cache for upcoming blocks
855+
AdviseBlockDontNeed(block);
856+
}
857+
858+
859+
void dna_adjust::AdviseBlockDontNeed(const UINT32& block)
860+
{
861+
using advice = boost::interprocess::mapped_region::advice_types;
862+
normalsR_map_.AdviseRegion(block, advice::advice_dontneed);
863+
measMinusComp_map_.AdviseRegion(block, advice::advice_dontneed);
864+
estimatedStations_map_.AdviseRegion(block, advice::advice_dontneed);
865+
originalStations_map_.AdviseRegion(block, advice::advice_dontneed);
866+
rigorousStations_map_.AdviseRegion(block, advice::advice_dontneed);
867+
junctionVariances_map_.AdviseRegion(block, advice::advice_dontneed);
868+
junctionVariancesFwd_map_.AdviseRegion(block, advice::advice_dontneed);
869+
junctionEstimatesFwd_map_.AdviseRegion(block, advice::advice_dontneed);
870+
junctionEstimatesRev_map_.AdviseRegion(block, advice::advice_dontneed);
871+
rigorousVariances_map_.AdviseRegion(block, advice::advice_dontneed);
872+
precAdjMsrs_map_.AdviseRegion(block, advice::advice_dontneed);
873+
corrections_map_.AdviseRegion(block, advice::advice_dontneed);
874+
}
875+
876+
877+
void dna_adjust::AdviseBlockWillNeed(const UINT32& block)
878+
{
879+
using advice = boost::interprocess::mapped_region::advice_types;
880+
normalsR_map_.AdviseRegion(block, advice::advice_willneed);
881+
measMinusComp_map_.AdviseRegion(block, advice::advice_willneed);
882+
estimatedStations_map_.AdviseRegion(block, advice::advice_willneed);
883+
originalStations_map_.AdviseRegion(block, advice::advice_willneed);
884+
rigorousStations_map_.AdviseRegion(block, advice::advice_willneed);
885+
junctionVariances_map_.AdviseRegion(block, advice::advice_willneed);
886+
junctionVariancesFwd_map_.AdviseRegion(block, advice::advice_willneed);
887+
junctionEstimatesFwd_map_.AdviseRegion(block, advice::advice_willneed);
888+
junctionEstimatesRev_map_.AdviseRegion(block, advice::advice_willneed);
889+
rigorousVariances_map_.AdviseRegion(block, advice::advice_willneed);
890+
precAdjMsrs_map_.AdviseRegion(block, advice::advice_willneed);
891+
corrections_map_.AdviseRegion(block, advice::advice_willneed);
854892
}
855893

856894

@@ -896,7 +934,7 @@ void dna_adjust::UnloadBlock(const UINT32& block, const int file_count, ...)
896934
if (file_count == 0)
897935
{
898936
// deserialise all
899-
UnloadBlock(block, 16,
937+
UnloadBlock(block, 15,
900938
sf_normals, sf_normals_r, sf_atvinv, sf_design, sf_meas_minus_comp,
901939
sf_estimated_stns, sf_original_stns, sf_rigorous_stns,
902940
sf_junction_vars, sf_junction_vars_f, sf_junction_ests_f, sf_junction_ests_r,
@@ -907,58 +945,60 @@ void dna_adjust::UnloadBlock(const UINT32& block, const int file_count, ...)
907945
va_list vlist;
908946
va_start(vlist, file_count);
909947

910-
// Unload block matrix data from memory
948+
// Unload block matrix data from memory by freeing buffers.
949+
// Use deallocate() instead of explicit destructor calls to keep
950+
// the matrix objects in a valid state for later reuse.
911951
for (UINT16 file(0); file<file_count; ++file)
912952
{
913953
switch (va_arg(vlist, int))
914954
{
915955
case sf_normals:
916-
v_normals_.at(block).~matrix_2d();
956+
v_normals_.at(block).deallocate();
917957
break;
918958
case sf_normals_r:
919-
v_normalsR_.at(block).~matrix_2d();
959+
v_normalsR_.at(block).deallocate();
920960
break;
921961
case sf_atvinv:
922-
v_AtVinv_.at(block).~matrix_2d();
962+
v_AtVinv_.at(block).deallocate();
923963
break;
924964
case sf_design:
925-
v_design_.at(block).~matrix_2d();
965+
v_design_.at(block).deallocate();
926966
break;
927967
case sf_meas_minus_comp:
928-
v_measMinusComp_.at(block).~matrix_2d();
968+
v_measMinusComp_.at(block).deallocate();
929969
break;
930970
case sf_estimated_stns:
931-
v_estimatedStations_.at(block).~matrix_2d();
971+
v_estimatedStations_.at(block).deallocate();
932972
break;
933973
case sf_original_stns:
934-
v_originalStations_.at(block).~matrix_2d();
974+
v_originalStations_.at(block).deallocate();
935975
break;
936976
case sf_rigorous_stns:
937-
v_rigorousStations_.at(block).~matrix_2d();
977+
v_rigorousStations_.at(block).deallocate();
938978
break;
939979
case sf_junction_vars:
940-
v_junctionVariances_.at(block).~matrix_2d();
980+
v_junctionVariances_.at(block).deallocate();
941981
break;
942982
case sf_junction_vars_f:
943-
v_junctionVariancesFwd_.at(block).~matrix_2d();
983+
v_junctionVariancesFwd_.at(block).deallocate();
944984
break;
945985
case sf_junction_ests_f:
946-
v_junctionEstimatesFwd_.at(block).~matrix_2d();
986+
v_junctionEstimatesFwd_.at(block).deallocate();
947987
break;
948988
case sf_junction_ests_r:
949-
v_junctionEstimatesRev_.at(block).~matrix_2d();
989+
v_junctionEstimatesRev_.at(block).deallocate();
950990
break;
951991
case sf_rigorous_vars:
952-
v_rigorousVariances_.at(block).~matrix_2d();
992+
v_rigorousVariances_.at(block).deallocate();
953993
break;
954994
case sf_prec_adj_msrs:
955-
v_precAdjMsrsFull_.at(block).~matrix_2d();
995+
v_precAdjMsrsFull_.at(block).deallocate();
956996
break;
957997
case sf_corrections:
958-
v_corrections_.at(block).~matrix_2d();
998+
v_corrections_.at(block).deallocate();
959999

9601000
if (v_blockMeta_.at(block)._blockLast)
961-
v_correctionsR_.at(block).~matrix_2d();
1001+
v_correctionsR_.at(block).deallocate();
9621002
break;
9631003
}
9641004
}

0 commit comments

Comments
 (0)