Structured matrix construction: fix int type usage

Use ordinal_type to index over grid points and rows, and size_type to
index over matrix entries.

This applies to all backends, but on HIP it also works around a compiler
bug with division-with-remainder.