[deal.II] FESystem in tutorial 31

[Huan Sun]

Hi all,

I was reading tutorial 31 and got confused by the following statements
in output_results() function:

             if (joint_fe.system_to_base_index(i).first.first == 0)
             {
               Assert (joint_fe.system_to_base_index(i).second
                       <
                       local_stokes_dof_indices.size(),
                       ExcInternalError());
               joint_solution(local_joint_dof_indices[i])
                 = 
stokes_solution(local_stokes_dof_indices[joint_fe.system_to_base_index(i).second]);
             }
           else
             {
               Assert (joint_fe.system_to_base_index(i).first.first == 1,
                       ExcInternalError());
               Assert (joint_fe.system_to_base_index(i).second
                       <
                       local_stokes_dof_indices.size(),
                       ExcInternalError());
               joint_solution(local_joint_dof_indices[i])
                 = 
temperature_solution(local_temperature_dof_indices[joint_fe.system_to_base_index(i).second]);
             }

Should the third assert compare joint_fe.system_to_base_index(i).second with 
local_temperature_dof_indices.size() instead of local_stokes_dof_indices.size()?

thanks,

Huan
          
                     

_______________________________________________
dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii