 <OpenGeoSysProjectDiff base_file="cube_1e0.prj">

   <add pos="after" sel=".//linear_solvers/linear_solver">

     <linear_solver>

       <name>general_linear_solver</name>

       <eigen>

         <solver_type>PardisoLU</solver_type>

       </eigen>

     </linear_solver>

   </add>

   <remove sel=".//linear_solvers/linear_solver[1]" />

 </OpenGeoSysProjectDiff>