sort_and_sum Subroutine

public subroutine sort_and_sum(nz, j, val)

Uniques and sorts values according to j-ascending ordering Same j-values are summed up in val

Arguments

Type IntentOptional Attributes Name
integer, intent(inout) :: nz

On input: Dimension of j and val On output: Number of unique

integer, intent(inout), dimension(nz) :: j

Integer list, on output uniqued and sorted in ascending order

real(kind=FP), intent(inout), dimension(nz) :: val

Values, on output sorted and summed according to new j-ordering