<p class="MsoNormal">
        </p><p class="MsoNormal">
                Dear
Mr./Ms.，
        </p>
        <p class="MsoNormal">
                <br>
        </p>
        <p class="MsoNormal" style="text-indent:24.0pt;">
                I have used
‘DMCreateMatrix’
to create a matrix <b>K</b>, and also the ‘DMCreateGlobalVector’
to create two vectors <b>U</b> (to be solved)
and <b>F </b>(right-hand side), i.e. <b>KU</b>=<b>F</b>.
Now, I want to add some complex constraints to this system through lagrangian
multiplier method, and the constraint matrix is <b>Q</b>. The <b>KU</b>=<b>F </b>transforms to
        </p>
        <p class="MsoNormal" style="text-align:center;text-indent:24pt;">
                <img src="cid:74f9957c$2$18003a469e0$Coremail$gaoyahe$buaa.edu.cn">
        </p>
        <p class="MsoNormal" align="center" style="text-align:center;">
                <b></b>
        </p>
        <p class="MsoNormal">
                   How
to create <b>K</b><b>_λ</b>, and how to effectively copy values<b> K</b> and <b>Q</b> to <b>K</b><b>_λ</b>? Does the newly created <b>K</b><b>_λ</b> and <b>F</b><b>_λ</b>still have an advantage
of DMDA? Or do you have any other good suggestions for this kind of problem? 
        </p>
        <p class="MsoNormal">
                <br>
        </p>
        <p class="MsoNormal">
                Thank you very
much!
        </p>
        <p class="MsoNormal">
                <br>
        </p>
        <p class="MsoNormal">
                Best
regards,
        </p>
        <p class="MsoNormal">
                A
PETSc user
        </p>