TimSzuHsienWu/cm — comment on commit 10e13388c2699b8a2ad737897f90e28928f1733e
File: src/data_point_routines.f90 line 241
By TimSzuHsienWu:

Hi Chris

The initial idea was to have a gauss_points_type separate from the
data_points_type, where they both use data_projections_type for
projections. Each data_projection_type however is still associated with one
gauss_points_type or data_points_type and one mesh which it is projected
onto.

But, for gauss_points_type especially, we would probably want to project
the same group of gauss points onto multiple meshes, so when creating the
gauss_points_type (or data_points_type), I initialize the array to store
the results that may come from data projection. Since results are stored in
gauss_points_type or data_points_type respectively, it needs to know how
many dimensions (number of projections) it may possibly have.

An alternative is to store the projection results in data_projections_type
(instead of gauss_points_type or data_points_type), but the problem is that
for each projected data point, we need to link it back to each data_point
(not a group of data_points), which I think will require alot of extra
memories for the pointers.

There are probably something I have overlooked, please let me know if I
shall still remove number of data projections from data points create start
routines

Regards
Tim