diff options
| -rw-r--r-- | src/ChangeLog | 4 | ||||
| -rw-r--r-- | src/algebra/Makefile.in | 3 | ||||
| -rw-r--r-- | src/algebra/Makefile.pamphlet | 3 | ||||
| -rw-r--r-- | src/algebra/catdef.spad.pamphlet | 5 | ||||
| -rw-r--r-- | src/algebra/exposed.lsp.pamphlet | 1 | ||||
| -rw-r--r-- | src/algebra/vector.spad.pamphlet | 52 |
6 files changed, 64 insertions, 4 deletions
diff --git a/src/ChangeLog b/src/ChangeLog index fe260b86..bcab5bca 100644 --- a/src/ChangeLog +++ b/src/ChangeLog @@ -1,5 +1,9 @@ 2010-06-30 Gabriel Dos Reis <gdr@cs.tamu.edu> + * algebra/vector.spad.pamphlet (LinearElement): New. + +2010-06-30 Gabriel Dos Reis <gdr@cs.tamu.edu> + * interp/c-util.boot (getRepresentation): Tidy. * interp/define.boot (checkRepresentation): Don't compile the base domain here yet. diff --git a/src/algebra/Makefile.in b/src/algebra/Makefile.in index ca09293a..5d04110f 100644 --- a/src/algebra/Makefile.in +++ b/src/algebra/Makefile.in @@ -612,6 +612,7 @@ $(OUT)/DIRPROD.$(FASLEXT): $(OUT)/DIRPCAT.$(FASLEXT) $(OUT)/DIRPCAT.$(FASLEXT): $(OUT)/VSPACE.$(FASLEXT) $(OUT)/IVECTOR.$(FASLEXT) $(OUT)/MATRIX.$(FASLEXT): $(OUT)/MATCAT.$(FASLEXT) $(OUT)/BTAGG.$(FASLEXT): $(OUT)/BOOLE.$(FASLEXT) +$(OUT)/LINELT.$(FASLEXT): $(OUT)/VSPACE.$(FASLEXT) $(OUT)/FM.$(FASLEXT) axiom_algebra_layer_10 = \ RESULT BFUNCT BPADIC ANY \ @@ -640,7 +641,7 @@ axiom_algebra_layer_10 = \ FUNDESC XPBWPOLY SMATCAT SMATCAT- \ RMATRIX RMATCAT RMATCAT- DIRPROD \ DIRPCAT DIRPCAT- IVECTOR MATRIX \ - MATCAT MATCAT- IIARRAY2 + MATCAT MATCAT- IIARRAY2 LINELT axiom_algebra_layer_10_nrlibs = \ diff --git a/src/algebra/Makefile.pamphlet b/src/algebra/Makefile.pamphlet index b51a660e..2a600ccb 100644 --- a/src/algebra/Makefile.pamphlet +++ b/src/algebra/Makefile.pamphlet @@ -619,6 +619,7 @@ $(OUT)/DIRPROD.$(FASLEXT): $(OUT)/DIRPCAT.$(FASLEXT) $(OUT)/DIRPCAT.$(FASLEXT): $(OUT)/VSPACE.$(FASLEXT) $(OUT)/IVECTOR.$(FASLEXT) $(OUT)/MATRIX.$(FASLEXT): $(OUT)/MATCAT.$(FASLEXT) $(OUT)/BTAGG.$(FASLEXT): $(OUT)/BOOLE.$(FASLEXT) +$(OUT)/LINELT.$(FASLEXT): $(OUT)/VSPACE.$(FASLEXT) $(OUT)/FM.$(FASLEXT) axiom_algebra_layer_10 = \ RESULT BFUNCT BPADIC ANY \ @@ -647,7 +648,7 @@ axiom_algebra_layer_10 = \ FUNDESC XPBWPOLY SMATCAT SMATCAT- \ RMATRIX RMATCAT RMATCAT- DIRPROD \ DIRPCAT DIRPCAT- IVECTOR MATRIX \ - MATCAT MATCAT- IIARRAY2 + MATCAT MATCAT- IIARRAY2 LINELT axiom_algebra_layer_10_nrlibs = \ diff --git a/src/algebra/catdef.spad.pamphlet b/src/algebra/catdef.spad.pamphlet index 104c4f10..1c22c456 100644 --- a/src/algebra/catdef.spad.pamphlet +++ b/src/algebra/catdef.spad.pamphlet @@ -295,12 +295,15 @@ OrderedType(): Category == BasicType with ++ a model of the category \spadtype{OrderedType}. )abbrev domain ORDSTRCT OrderedStructure OrderedStructure(T: Type,f: (T,T) -> Boolean): Public == Private where - Public == Join(OrderedType,HomotopicTo T) + Public == Join(OrderedType,HomotopicTo T) with + if T has CoercibleTo OutputForm then CoercibleTo OutputForm Private == add Rep == T coerce(x: %): T == rep x coerce(y: T): % == per y x < y == f(rep x,rep y) + if T has CoercibleTo OutputForm then + coerce(x: %): OutputForm == rep(x)::OutputForm @ diff --git a/src/algebra/exposed.lsp.pamphlet b/src/algebra/exposed.lsp.pamphlet index 66f6278f..ba666e3a 100644 --- a/src/algebra/exposed.lsp.pamphlet +++ b/src/algebra/exposed.lsp.pamphlet @@ -231,6 +231,7 @@ (|LetAst| . LETAST) (|Library| . LIB) (|LieSquareMatrix| . LSQM) + (|LinearElement| . LINELT) (|LinearOrdinaryDifferentialOperator| . LODO) (|LinearSystemMatrixPackage| . LSMP) (|LinearSystemMatrixPackage1| . LSMP1) diff --git a/src/algebra/vector.spad.pamphlet b/src/algebra/vector.spad.pamphlet index 08ec4b60..05474370 100644 --- a/src/algebra/vector.spad.pamphlet +++ b/src/algebra/vector.spad.pamphlet @@ -447,11 +447,60 @@ DirectProductFunctions2(dim, A, B): Exports == Implementation where reduce(f, v, b) == reduce(f, v::VA, b) @ + +\section{Vector space of finite dimension given by a basis} + +<<domain LINELT LinearElement>>= +)abbrev domain LINELT LinearElement +++ Author: Gabriel Dos Reis +++ Date Created: June 30, 2010 +++ Date Last Modified: June 30, 2010 +++ Description: +++ A simple data structure for elements that form a +++ vector space of finite dimension over a given field, +++ with a given symbolic basis. +LinearElement(K,B): Public == Private where + K: Field + B: List Symbol + macro Basis == OrderedVariableList B + Public == Join(VectorSpace K,CoercibleFrom Basis) with + linearElement: (List K,List Basis) -> % + ++ \spad{linearElement([x1,..,xn],[b1,..,bn]) constructs + ++ a linear element with coordinates \spad{[x1,..,xn]} with + ++ respect to the basis elements \spad{b1},...\spad{bn}. + coordinates: % -> Vector K + ++ \spad{coordinates x} returns the coordinates of the linear + ++ element with respect to the basis \spad{B}. + Private == FreeModule(K,Basis) add + coerce(b: Basis): % == + per monomial(1$K,b) + + dimension() == + size()$Basis::CardinalNumber + + linearElement(cs,bs) == + null cs or null bs => per 0$Rep + reduce(_+@(%,%)->%, + [per monomial(c,b) for c in cs for b in bs]) + + coordinates x == + n := #B + v: Vector K := new(n,0$K) + ts := terms rep x + for i in 1..n repeat + t := first ts + lookup first t = i => + v.i := second t + ts := rest ts + v +@ + + \section{License} <<license>>= --Copyright (c) 1991-2002, The Numerical ALgorithms Group Ltd. --All rights reserved. ---Copyright (C) 2007-2009, Gabriel Dos Reis. +--Copyright (C) 2007-2010, Gabriel Dos Reis. --All rights reserved. -- --Redistribution and use in source and binary forms, with or without @@ -492,6 +541,7 @@ DirectProductFunctions2(dim, A, B): Exports == Implementation where <<category DIRPCAT DirectProductCategory>> <<domain DIRPROD DirectProduct>> <<package DIRPROD2 DirectProductFunctions2>> +<<domain LINELT LinearElement>> @ \eject |