I would like to seek clarification regarding mapping of cell (scalar) node to vector node while filling up the matrix.
As mentioned in the documentation, the type_ls_map does the needful job.
What does equation_ls_map%l(i,j,k) refers to?
The vector or matrix node corresponding to a particular scalar node has been described based on the basis of above definition (ia=l*equation_stencil%size + 1, as in add_cell_diffusion_term_centered_o2).
May I know the significance of adding 1.