Re: Types with physical dimension

[mab@dst17.wdl.loral.com (Mark A Biggar)]

Another possibility is to use a discriminated type to keep track of dimensions
like so:

generic
	type Base is digits <>;
package Units is

type Unit(Length_DIM: Integer; Mass_DIM: Integer; Time_DIM: Integer) is
	private;

function "*"(Left, Right: Unit) return Unit;
function "*"(Left: Unit; Right: Base) return Unit;
function "+"(Left, Right: Unit) return Unit;
function "-"(X: Unit) return Unit;

-- etc;

subtype Length is Unit(1,0,0);
subtype Area is Unit(2,0,0);
subtype Mass is Unit(0,1,0);
subtype Time is Unit(0,0,1);
suntype Frequency is Unit(0,0,-1);
subtype Energy is Unit(2,1,-2);

-- etc;

private

type Unit(Length_DIM: Integer; Mass_DIM: Integer; Time_DIM: Integer) is
	record
		Value: Base;
	end record;
end Units;

package body Units is

procedure Check_DIM(Left, Right: Unit) is
begin
	if Left.Length_DIM \= Right.Length_DIM or else
	   Left.Mass_DIM \= Right.Mass_DIM or else
	   Left.Time_DIM \= Right.Time_DIM then
		raise CONSTRAINT_ERROR;
end Check_DIM;

function "*"(Left, Right: Unit) return Unit is
begin
	return Unit'(Left.Length_DIM+Right.Length_DIM,
		     Left.Mass_DIM+Right.Mass DIM,
		     Left.Time_DIM+Right.Time_DIM,
		     Left.Value * Right.Value);
end "*";

	

function "+"(Left, RIght: Unit) return Unit is
begin
	Check_DIM(Keft, Right);
	return UniL'(Left.Length_DIM,
		     Left.Mass_DIM,
		     Left.Time_DIM,
		     Left.Value + Right.Value);
end "+";

function "-"(X: Unit) is
begin
        return Unit'(X.Length_DIM,
                     X.Mass_DIM,
                     X.Time_DIM,
                     -X.Value);
end "-";

-- etc;

end Units;

Yes I know that this requires runtime dimension checking, but you only
have to use it until you are through debugging your code then
you can replace the body with one that turns off all the checking.

--
Mark Biggar
mab@wdl.loral.com