Hi srvaldez.
I made up a fractions udt with operators.
It should keep values, it only uses integers.
I kept your matrix type, but had to tweak the multiply slightly to suit type fraction.
Code: Select all
Type fraction
As String n,d 'numerator,denominator
Declare Constructor
Declare Constructor(As String)
Declare Operator Let(As String)
Declare Operator Cast() As String
Declare Function Val() As Double
#define f Type<fraction>
End Type
Sub Remove(Text As String,Char As String)
Var index = 0,asci=Asc(char)
For i As Integer = 0 To Len(Text) - 1
If Text[i] <> asci Then Text[index] = Text[i] : index =index+ 1
Next : Text = Left(Text,index)
End Sub
Function hcf(a As Longint,b As Longint) As Longint
a=abs(a):b=abs(b)
If b>a Then Swap b,a
Dim As Longint c
While b
c = a
a = b
b = c Mod b
Wend
Return a
End Function
Constructor fraction
End Constructor
Constructor fraction (s As String)
remove(s," ") 'keep tight,optional if tight anyway
Var i=Instr(s,"/")
n=Mid(s,1,i-1)
d=Mid(s,i+1)
If Instr(s,"/")=0 Then n=s:d="1"
End Constructor
Operator fraction.let(s As String)
Var i=Instr(s,"/")
n=Mid(s,1,i-1)
d=Mid(s,i+1)
End Operator
Operator fraction.cast() As String
if valint(n)<0 and valint(d)<0 then n=ltrim(n,"-"):d=ltrim(d,"-")
If Instr(d,"-") Then d=Ltrim(d,"-"):n="-"+n
Return Rtrim(n+"/"+d,"/1")
End Operator
Function fraction.val() As Double
Return Valint(n)/Valint(d)
End Function
Operator +(a As fraction,b As fraction) As fraction
Var p=(Valint(a.d)*Valint(b.d))
Var n1=Valint(a.n)*Valint(b.d)
Var n2=Valint(b.n)*Valint(a.d)
Var n3=n1+n2
Var h= hcf(n3,p)
If h>1 Then
n3/=h:p/=h
End If
Return Str(n3)+"/"+Str(p)
End Operator
Operator -(a As fraction,b As fraction) As fraction
Var p=(Valint(a.d)*Valint(b.d))
Var n1=Valint(a.n)*Valint(b.d)
Var n2=Valint(b.n)*Valint(a.d)
Var n3=n1-n2
Var h= hcf(n3,p)
If h>1 Then
n3/=h:p/=h
End If
Return Str(n3)+"/"+Str(p)
End Operator
Operator *(a As fraction,b As fraction) As fraction
Var n2=Valint(a.n)*Valint(b.n)
Var d2=Valint(a.d)*Valint(b.d)
Var h= hcf(Abs(n2),Abs(d2))
n2/=h:d2/=h
Return Str(n2)+"/"+Str(d2)
End Operator
Operator /(a As fraction,b As fraction) As fraction
Swap b.n,b.d
Return a*b
End Operator
'============= srvaldez ==============
Type Matrix
Dim As fraction m( Any , Any )
Declare Constructor ( )
Declare Constructor ( Byval x As Uinteger , Byval y As Uinteger )
End Type
Constructor Matrix ( )
End Constructor
Constructor Matrix ( Byval x As Uinteger ,Byval y As Uinteger )
Redim this.m( x - 1 , y - 1 )
End Constructor
Operator *(m1 As matrix,m2 As matrix) As matrix
Dim rows As Integer=Ubound(m1.m,1)
Dim columns As Integer=Ubound(m2.m,2)
If Ubound(m1.m,2)<>Ubound(m2.m,1) Then
Print "Can't do multiply"
Exit Operator
End If
Dim As matrix ans=matrix(rows+1,columns+1)
Dim rxc As fraction
For r As Integer=0 To rows
For c As Integer=0 To columns
rxc=f("0")
For k As Integer = 0 To Ubound(m1.m,2)
rxc=rxc+(m1.m(r,k)*m2.m(k,c))
Next k
ans.m(r,c)=rxc
Next c
Next r
Operator= ans
End Operator
'some garbage matrices for demonstration
Dim As Matrix a = Matrix(4 , 4)
a.m(0 , 0) = f("1") : a.m(0 , 1) = f("1") : a.m(0 , 2) = f("1") : a.m(0 , 3) = f("1")
a.m(1 , 0) = f("2") : a.m(1 , 1) = f("4") : a.m(1 , 2) = f("8") : a.m(1 , 3) = f("16")
a.m(2 , 0) = f("3") : a.m(2 , 1) = f("9") : a.m(2 , 2) = f("27"): a.m(2 , 3) = f("81")
a.m(3 , 0) = f("4") : a.m(3 , 1) = f("16") : a.m(3 , 2) = f("64") : a.m(3 , 3) = f("256")
Dim As Matrix b = Matrix( 4 , 4 )
b.m(0 , 0) = f("4") : b.m(0 , 1) = f("-3") : b.m(0 , 2) = f("4/3") : b.m (0, 3) = f("-1/4")
b.m(1 , 0) = f("-13/3") : b.m(1 , 1) = f("19/4") : b.m(1 , 2) = f("-7/3") : b.m (1, 3) = f("11/24")
b.m(2 , 0) = f("3/2") : b.m(2 , 1) = f("-2") : b.m(2 , 2) = f("7/6") : b.m (2, 3) = f("-1/4")
b.m(3 , 0) = f("-1/6") : b.m(3 , 1) = f("1/4") : b.m(3 , 2) = f("-1/6") : b.m (3, 3) = f("1/24")
Dim As Matrix c = a * a * b
Print c.m(0, 0), c.m(0, 1), c.m(0, 2), c.m(0, 3)
Print c.m(1, 0), c.m(1, 1), c.m(1, 2), c.m(1, 3)
Print c.m(2, 0), c.m(2, 1), c.m(2, 2), c.m(2, 3)
Print c.m(3, 0), c.m(3, 1), c.m(3, 2), c.m(3, 3)
?"=============================================="
For i As Long=1 To 4
c = c * c * b
Print c.m(0, 0), c.m(0, 1), c.m(0, 2), c.m(0, 3)
Print c.m(1, 0), c.m(1, 1), c.m(1, 2), c.m(1, 3)
Print c.m(2, 0), c.m(2, 1), c.m(2, 2), c.m(2, 3)
Print c.m(3, 0), c.m(3, 1), c.m(3, 2), c.m(3, 3)
?"=============================================="
Next
Sleep
seems to work for me
