in reply to
Re: What will scientific computing in Perl 6 look like?
in thread What will scientific computing in Perl 6 look like?
By the way complex numbers are natively supported in Perl 6.
Can you construct arbitrary vector spaces (e.g., the quadratics) or are the complex numbers a special case?
UPDATE: I think abstract algebra starts at the undergrad level. Also, wow, perl6 is like a whole new language.
Re^3: What will scientific computing in Perl 6 look like? by moritz (Cardinal) on Jul 14, 2008 at 12:22 UTC 
# suppose you inherit from List to implement your vector:
class MyVector is List {
sub infix:<+>(MyVector @self, MyVector @other) is export {
return @self »+« @other;
}
# adding a scalar:
sub infix:<+>(MyVector @self, Num $other) is export {
return @self »+ $other
}
}
No more hassle of overloading literal constants in the source, since multi method dispatch takes care of it all.  [reply] [d/l] 

 [reply] 

Well, Perl 6's type system has difficulties with respect to numeric types, because they don't form a nice hierarchy. Luckily we have roles, and thus can have a type of all (scalar) numeric types without pressing them into a hierarchy.
Implementing the numeric types correctly is a bit tricky, but the fact remains that the user (aka Perl 6 programmer) can simply add (or override) operators for each pair of types. If one of them is always a userdefined one you can guarantee no clashes with existing semantics.
So I don't see how Perl 6 isn't a solution wrt to adding new numeric types to the language  would you care to elaborate?
 [reply] 


 

Don't you need to use the multi keyword on those subs?
 [reply] [d/l] 

It's better style, but if the prelude defines a proto infix:<+> (and I'm quite sure it does) it's redundant. Good catch anyway.
 [reply] [d/l] 

Re^3: What will scientific computing in Perl 6 look like? by blazar (Canon) on Jul 20, 2008 at 10:12 UTC 
Can you construct arbitrary vector spaces (e.g., the quadratics) or are the complex numbers a special case?
I personally believe that this is slightly nonsense, since vector spaces are "simply" vector spaces, while the complex numbers can be described in a variety of different algebraic structures: they're a field, a two dimensional real algebra, a one dimensional complex algebra, etc. Now, I'm sure that binary field operations will be supported for complex numbers: thus whether arbitrary vector spaces will be supported or not, it won't be a matter of "special case."
 [reply] [d/l] 

It isn't nonsense at all. Subspaces, fields, and vector spaces are all interrelated. The complex numbers are a vector space and the quadratics are a subspace of them. I was hoping that instead of just complex numbers, you could use the special type for all sorts of vector math, tensors, quadratics, values with units attached. Why make complex numbers the special case? Why not have two tuples or ntuples that know about coefficients so you can do all kinds of interesting algebra in the same way you can use complex numbers.
All you have to do to make this really really flexible is allow users to change the value of that sqrt(1) coefficient(s) and change the length of the vector/sum/value. Tada.
UPDATE: Yes, well, I believe the quadratics are a 2dimensional vector space and a field of reals. It happens to be a subspace of the complexes, so maybe it doesn't matter, but I fail to see why there couldn't be an interface to use the complex number system for more.
 [reply] 

I personally believe it is still a nonsense, wrt the "special case" (your literal words) bit: do not misunderstand me! I am too hoping that to quote you verbatim we "could use the special type for all sorts of vector math, tensors, quadratics, values with units attached." (And more!) Nevertheless, complex numbers are something more: i.e. a field, and that cannot be told of general ndimensional vector spaces: given the confidence you talk about these topics, you certainly know that the only three real algebras (with division) are the reals themselves, the complex numbers and quaternions. (Then you have the more esoteric Cayley numbers, if you're willing to give up on associativity.) Or else, are you taking, say, two ndimensional vectors and multiply them together?!? Certainly, you can have an inner product, but that's an entirely different beast, and not a generalization of complex multiplication. Thus, definitely, not as a "special case."
 [reply] [d/l] 

