MPCと変位境界条件の両方が与えられた節点の計算結果がおかしい

MPC自由度消去の変形をした後に変位境界条件処理をしなければいけないところ、順番が逆になっているので、変位境界を与えられた節点の変位がおかしい。接触がある場合（solve_lineq_contact）はソルバ内部で自由度消去を行っているので、どう順番を変えるか。

下記、発生したときの報告

!EQUATIONの入力節点に!BOUNDARYで変位境界条件を与えると、!BOUNDARYが無視される動作が不具合かもしれないのでご連絡いたします。 下記サンプルデータでは、節点1の解析結果Z変位が!BOUNDARYを無視して0になります。 【.msh】 !HEADER HECMW_MSH File generated by REVOCAP !NODE 1, 0.00000000e+00, 0.00000000e+00, 1.00000000e+00 2, 1.00000000e+00, 0.00000000e+00, 1.00000000e+00 3, 1.00000000e+00, 1.00000000e+00, 1.00000000e+00 4, 0.00000000e+00, 1.00000000e+00, 1.00000000e+00 5, 1.00000000e+00, 0.00000000e+00, 0.00000000e+00 6, 0.00000000e+00, 0.00000000e+00, 0.00000000e+00 7, 0.00000000e+00, 1.00000000e+00, 0.00000000e+00 8, 1.00000000e+00, 1.00000000e+00, 0.00000000e+00 !ELEMENT, TYPE=361, EGRP=Pconcrete 7,7,4,1,6,8,3,2,5 !MATERIAL, NAME=concrete, ITEM=1 !ITEM=1, SUBITEM=1 0.000000e+00 !SECTION, TYPE=SOLID, EGRP=Pconcrete, MATERIAL=concrete !AMPLITUDE, NAME=step, VALUE=RELATIVE 0.000000e+00,0.000000e+00 1.000000e+00,1.000000e+06 1.000000e+00,1.000000e+07 !NGROUP, NGRP=zdn1 1 !NGROUP, NGRP=BND_fix_1 5,6,7,8 1 !NGROUP, NGRP=BND_fix_2 1,2,5,6 3,4,7,8 101 !NGROUP, NGRP=BND_fix_3 5,6,7,8 !EQUATION 3,0 1, 3, -1.0 2, 3, 0.5 3, 3, 0.5 !END

【.cnt】 !VERSION 3 !SOLUTION,TYPE=DYNAMIC !DYNAMIC, TYPE=NONLINEAR 1, 1 0.000000e+00, 2000000.0, 2, 1000000.0, 0 5.000000e-01, 2.500000e-01 1, 1, 0.000000e+00, 0.000100e+00 1, 7, 1 0, 0, 0, 0, 0, 0 !SOLVER,METHOD=MUMPS 20000, 2, 10, 10 1.000000e-06, 1.000000e+00, 0.000000e+00 !STEP, CONVERG=1.0E-3, ITMAX=10 !BOUNDARY, GRPID=1 BND_fix_1, 1, 1, 0.000000e+00 BND_fix_2, 2, 2, 0.000000e+00 BND_fix_3, 3, 3, 0.000000e+00 !BOUNDARY, GRPID=101, AMP=step zdn1, 3, 3, 1.0 !MATERIAL, NAME=concrete !ELASTIC, TYPE=ISOTROPIC, INFINITE=1 28545.45455e+6, 0.2 !DENSITY 2400.0 !WRITE,VISUAL, FREQUENCY=1 #!WRITE,RESULT, FREQUENCY=1 !OUTPUT_VIS DISP, ON REACTION, ON NSTRAIN, OFF NSTRESS, OFF NMISES, OFF ESTRAIN, ON ESTRESS, ON !VISUAL, method=PSR !surface_num = 1 !surface !display_method = 1 !surface_style = 1 !initial_style = 1 !output_type = COMPLETE_AVS !x_resolution = 800 !y_resolution = 600 !num_of_lights = 1 !position_of_lights = 100.000, 100.000, 50.0000 !up_direction = 0.00000, 0.00000, 1.00000 !ambient_coef = 0.3 !diffuse_coef = 0.7 !specular_coef = 0.5 !color_mapping_bar_on = 1 !scale_marking_on = 1 !num_of_scale = 3 !font_size = 1.5 !font_color = 1.00000, 1.00000, 1.00000 !background_color = 0.00000, 0.00000, 0.00000 !END