At the core of Catlab is a system for defining generalized algebraic structures, such as categories and monoidal categories, and then creating instances of these structures in Julia code. The objects, morphisms, and even higher-order morphisms can also be represented as typed symbolic expressions, as in a computer algebra system. These expressions can be manipulated abstractly or transformed into more concrete representations, such as [wiring diagrams](@ref wiring_diagrams) or [Julia functions](@ref programs).
The basic elements of this system are:
-
Generalized algebraic theories (GATs), defined using the
@theorymacro. Categories and other typed (multisorted) algebraic structures can be defined as GATs. Alternatively, the@signaturemacro can be used when only the signature (not the axioms) of the GAT are to be specified. -
Instances, or concrete implementations, of theories, defined using the
@instancemacro. -
Syntax systems for theories, defined using the
@syntaxmacro. These are type-safe expression trees constructed using ordinary Julia functions.
We'll explain each of these elements in greater detail in the following sections. From the programming perspective, theories can be thought of as interfaces and bear some resemblance to type classes in languages like Haskell. Both instances and syntax systems can then be thought of as implementations of the interface.
Generalized algebraic theories (GATs) are the natural logical system in which to define categories and related algebraic structures. GATs generalize the typed (multisorted) algebraic theories of universal algebra by incorporating a fragment of dependent type theory; they are perhaps the simplest dependently typed logics.
Catlab implements a version of the GAT formalism on top of Julia's type system,
taking advantage of Julia macros to provide a pleasant syntax. GATs are defined
using the @theory macro.
For example, the theory of categories could be defined by:
using Catlab
import Catlab.Theories: Ob, Hom, ObExpr, HomExpr, dom, codom, compose, ⋅, id
@theory Category{Ob,Hom} begin
@op begin
(→) := Hom
(⋅) := compose
end
Ob::TYPE
Hom(dom::Ob, codom::Ob)::TYPE
id(A::Ob)::(A → A)
compose(f::(A → B), g::(B → C))::(A → C) ⊣ (A::Ob, B::Ob, C::Ob)
(f ⋅ g) ⋅ h == f ⋅ (g ⋅ h) ⊣ (A::Ob, B::Ob, C::Ob, D::Ob,
f::(A → B), g::(B → C), h::(C → D))
f ⋅ id(B) == f ⊣ (A::Ob, B::Ob, f::(A → B))
id(A) ⋅ f == f ⊣ (A::Ob, B::Ob, f::(A → B))
end
nothing # hide
The code is simplified only slightly from the official Catlab definition of
Category. The theory has two type constructors, Ob (object) and Hom
(morphism). The type Hom is a dependent type, depending on two objects, named
dom (domain) and codom (codomain). The theory has two term
constructors, id (identity) and compose (composition).
Notice how the return types of the term constructors depend on the argument
values. For example, the term id(A) has type Hom(A,A). The term constructor
compose also uses context variables, listed to the right of the ⊣
symbol. These context variables can also be defined after a where clause,
but the left hand side must be surrounded by parentheses. This allows us to
write compose(f,g), instead of the more verbose compose(A,B,C,f,g) (for
discussion, see Cartmell, 1986, Sec 10: Informal syntax).
Notice the @op call where we can create method aliases that can then be used
throughout the rest of the theory and outside of definition. We can either use
this block notation, or a single line notation such as @op (⋅) := compose to
define a single alias. Here we utilize this functionality by replacing the Hom
and compose methods with their equivalent Unicode characters, → and ⋅
respectively. These aliases are also automatically available to definitions that
inherit a theory that already has the alias defined.
!!! note
In general, a GAT consists of a *signature*, defining the types and terms of
the theory, and a set of *axioms*, the equational laws satisfied by models
of the theory. The theory of categories, for example, has axioms of
unitality and associativity. At present, Catlab supports the specification
of both signatures and the axioms, but is not currently utilizing the axiom
definitions in any way, reflecting its status as a programming library, not
a proof assistant. It is the programmer's responsibility to ensure any
declared instances of an algebraic structure satisfy its axioms.
- Cartmell, 1986: Generalized algebraic theories and contextual categories, DOI:10.1016/0168-0072(86)90053-9
- Cartmell, 1978, PhD thesis: Generalized algebraic theories and contextual categories
- Pitts, 1995: Categorical logic, Sec 6: Dependent types
A theory can have one or more instances, or instantiations by ordinary Julia types and functions. This feature builds on Julia's support for generic functions with multiple dispatch.
Instances are declared using the @instance macro. In an instance of a
theory, each theory type is mapped to a Julia type and each term is mapped
to a Julia method of the same name. For example, the category of matrices could
be defined as an instance of the theory Category defined above:
using LinearAlgebra: I
struct MatrixDomain
eltype::Type
dim::Int
end
@instance Category{MatrixDomain, Matrix} begin
dom(M::Matrix) = MatrixDomain(eltype(M), size(M,1))
codom(M::Matrix) = MatrixDomain(eltype(M), size(M,2))
id(m::MatrixDomain) = Matrix{m.eltype}(I, m.dim, m.dim)
compose(M::Matrix, N::Matrix) = M*N
end
A = Matrix{Float64}([0 1; 1 0])
id(dom(A))
In this instance, the theory type Ob is mapped to the custom Julia type
MatrixDomain. The latter type has two fields, a Julia type eltype
representing a field dim representing the dimensionality
Hom is mapped to the standard Julia type Matrix.
Theories can also be instantiated as systems of symbolic expressions, using
the @syntax macro. The symbolic expressions are expression trees, as
commonly used in computer algebra systems. They are similar to Julia's Expr
type but they are instead subtyped from Catlab's GATExpr type and they
have a more refined type hierarchy.
A single theory can have different syntax systems, treating different terms as primitive or performing different simplication or normalization procedures. Catlab tries to make it easy to define new syntax systems. Many of the theories included with Catlab have default syntax systems, but the user is encouraged to define their own to suit their needs.
To get started, you can always call the @syntax macro with an empty body.
Below, we subtype from Catlab's abstract types ObExpr and HomExpr to enable
LaTeX pretty-printing and other convenient features, but this is not required.
@syntax CategoryExprs{ObExpr, HomExpr} Category begin
end
A, B, C, D = [ Ob(CategoryExprs.Ob, X) for X in [:A, :B, :C, :D] ]
f, g, h = Hom(:f, A, B), Hom(:g, B, C), Hom(:h, C, D)
compose(compose(f,g),h)
The resulting symbolic expressions perform no simplification. For example, the associativity law is not satisfied:
compose(compose(f,g),h) == compose(f,compose(g,h))
Thus, unlike instances of a theory, syntactic expressions are not expected to obey all the axioms of the theory.
However, the user may supply logic in the body of the @syntax macro to enforce
the axioms or perform other kinds of simplification. Below, we use the
associate function provided by Catlab to convert the binary
expressions representing composition into strict=true tells Catlab to check that the domain and codomain
objects are strictly equal and throw an error if they are not.
@syntax SimplifyingCategoryExprs{ObExpr, HomExpr} Category begin
compose(f::Hom, g::Hom) = associate(new(f,g; strict=true))
end
A, B, C, D = [ Ob(SimplifyingCategoryExprs.Ob, X) for X in [:A, :B, :C, :D] ]
f, g, h = Hom(:f, A, B), Hom(:g, B, C), Hom(:h, C, D)
compose(compose(f,g),h)
Now the associativity law is satisfied:
compose(compose(f,g),h) == compose(f,compose(g,h))
In some algebraic structures, there is a choice as to which operations should be
considered primitive and which should be derived. For example, in a cartesian
monoidal category,
the copy operation
In Catlab, the recommended way to deal with such situations is to define all
the operations in the theory and then allow particular syntax systems to
determine which operations, if any, will be derived from others. In the case of
the cartesian monoidal category, we could define a signature CartesianCategory
by inheriting from the builtin theory SymmetricMonoidalCategory.
using Catlab
import Catlab.Theories: Ob, Hom, ObExpr, HomExpr, SymmetricMonoidalCategory,
dom, codom, compose, id, otimes, munit, braid
@signature CartesianCategory{Ob,Hom} <: SymmetricMonoidalCategory{Ob,Hom} begin
mcopy(A::Ob)::(A → (A ⊗ A))
delete(A::Ob)::(A → munit())
pair(f::(A → B), g::(A → C))::(A → (B ⊗ C)) ⊣ (A::Ob, B::Ob, C::Ob)
proj1(A::Ob, B::Ob)::((A ⊗ B) → A)
proj2(A::Ob, B::Ob)::((A ⊗ B) → B)
end
nothing # hide
We could then define the copying operation in terms of the pairing.
@syntax CartesianCategoryExprsV1{ObExpr,HomExpr} CartesianCategory begin
mcopy(A::Ob) = pair(id(A), id(A))
end
A = Ob(CartesianCategoryExprsV1.Ob, :A)
mcopy(A)
Alternatively, we could define the pairing and projections in terms of the copying and deleting operations.
@syntax CartesianCategoryExprsV2{ObExpr,HomExpr} CartesianCategory begin
pair(f::Hom, g::Hom) = compose(mcopy(dom(f)), otimes(f,g))
proj1(A::Ob, B::Ob) = otimes(id(A), delete(B))
proj2(A::Ob, B::Ob) = otimes(delete(A), id(B))
end
A, B, C = [ Ob(CartesianCategoryExprsV2.Ob, X) for X in [:A, :B, :C] ]
f, g = Hom(:f, A, B), Hom(:g, A, C)
pair(f, g)
Modules = [
GAT,
Syntax,
Rewrite,
Present,
]
Private = false