Dellnitz's talk

1. http://www-math.uni-paderborn.de/~agdellnitz/index.html
2. GAIO/Mangen for LCS; NTG/DMOC(Discrete Mechanics and Optimal Control) for optimization
3. Kuhn-Tucker condition and Houman's problem; relation to Pareto's equilibrium?
4. Application to molecular dynamics, find invariant manifold and invraiant measures.

量子场论，弦理论与数学

人们通常认为，近代科学与以前的科学的区分别是近代科学有实验。这种看法是值得商榷的。著名物理学家杨振宁教授和著名哲学家海德格尔认为近代科学的最根本的特征是数学和实验的结合，自然科学的定律用抽象的数学形式表达，从而达到前所未有的深度和广度。作为近代科学标志的两大发明，万有引力和微积分都是由牛顿创造的。在牛顿以后的科学发展中也反复印证了这一点。
近代科学史上许多有伟大贡献的自然科学家也是数学家。这种状况一直延续到20世纪20年代。此后形式化的数学一度占据数学的中心，数学在很长一段时间淡化了和其他科学，尤其是理论物理的联系。从20世纪20年代，量子场论开始出现并逐步成为理论物理的中心。到20世纪70年代中数学和量子场论才开始建立起密切的联系。从80年代以来，获得菲尔兹奖的数学家中其工作和量子场论或弦论有直接联系的占一半。
对称性和量子化：支配物理和数学的两个基本原则
也许我们要问：为什麽量子场论和弦论会和数学有密切的关系？一个答案是，它们被相同的原则所支配。其中最重要的原则是：对称性和量子化。
什麽是对称性？从一些建筑设计，巴赫的音乐和粒子物理中的CPT破缺（杨振宁和李政道的诺贝尔奖工作）我们体验到各种离散对称性。伽罗瓦是第一个系统研究离散对称性并用于解决高次多项式方程不可解的问题的。对于自然界连续对称性似更重要。例如我们有：
。从伽里略的相对性原理导出牛顿第一定律，
。从洛伦茨对称性导出狭义相对论，
。从坐标变换不变性和局域洛伦茨不变性导出广义相对论，
。经魏耳等人的努力，电动力学可以表述为阿贝尔规范场，即具有局域变换不变性，规范群是阿贝尔群。非阿贝尔规范场，即杨－Mills场，是粒子物理的基础，也具有局域变换不变性，规范群是非阿贝尔群
这里我们也许可以用两个原理来表述对称性的重要作用：
爱因斯坦原理：物理世界的规律应该和我们的表述无关。
杨振宁原理：对称性支配相互作用。
上述原理在几何中也是基本的。几何量，如长度，面积，体积等也是和描述他们的方式无关。这一点充分反映在以下理论中：
嘉当和陈省身：活动标架法。
在70年代中杨振宁意识到规范场和陈省身先生研究的联络是一回事，似就是局域对称性在物理和几何两个领域的各自实现。
以下我们解释一下什麽是量子化。
量子化原理：微观世界的描述不能用决定性的方式来描述，他们是几率式的。事件的几率全体组成Hilbert 空间。动力学变量实现为Hilbert空间上的算子。
玻尔相容性原理：我们对于世界的每一种描述是不完备的，但是他们是相容，自洽的。
测不准原理是玻尔相容性原理的具体实现。
我们知道，量子力学已成为了解微观世界的基本工具。在量子力学发明后不久，人们把它用到电动力学的研究上。这时我们必须引入场的概念。经典的麦克斯韦方程是线性方程。它的解就是无穷多个波的叠加。其量子化乃是将无穷多个谐振子放在一起而无相互作用。当人们作计算时发现有许多无穷大。一直到1948年，量子电动力学才在引入重正化以后有了有限的定义并和实验吻合的极好。在1954年杨振宁－Mills 将规范场推广到非阿贝尔群。其量子化经许多人的努力得到实现。人们发现量子规范场理论是唯一具有渐进自由性质的量子场论。物理学家对于围扰场论用费曼图给出了定义。到1974年物理学家建立了基本粒子的标准模型。从此物质场基本被标准模型所描述。在此过程中杨先生的“对称性支配相互作用”起了重要作用。拉氏量中的相互作用往往被对称性的考虑所决定。人们也试图在此框架下将引力量子化，没有成功。实际上，引力场是不可重整的。
为什麽要研究超弦理论？
由上我们也许可以得到一点启示，即相互作用的统一实际上是对称性的统一。从20世纪70年代起，人们又发现了超对称。它是一种将对易和反对易关系非平凡的合在一起的代数结构。将这种代数局域化我们得到局域超对称。在此类变换下不变的就是所谓超引力。在超引力中我们所知道的4种相互作用合在一起。所以我们说在经典的意义下超引力把4种相互作用统一起来了。超引力的量子理论就是超弦理论。这就是为什麽我们认为超弦理论中包涵了量子引力。
弦理论把粒子不再看成一个点，而是看成一根弦。弦的运动扫出一条曲面，弦的振动给出粒子。当粒子碰撞时，他们不在某个特定的点碰撞，因而免去场论中令人头疼的无穷大问题。到了1985年人们发现共有5种协调的超弦理论。他们都在10维时空中运动。在我们将其中6维空间紧致化以后，我们可以得到通常的4维规范场论。从保持部分超对称的考虑，紧致化的6维空间必须是卡拉比－丘成桐空间。弦理论里自然包涵引力子，超引力是超弦理论的低能极限. 在1985年人们面临的问题是，在5种超弦理论中，哪一种是描述自然的？超弦理论如何和实验建立联系？
在1995－1998的第二次超弦革命中，上述问题取得了突破。人们发现了对偶性，即不同理论在其适当的范围内可以相互等价。其中最让人惊奇的是一些强相互作用的理论和某些弱相互作用的理论等价。这就为人们研究强相互作用开辟了道路。人们最初在超引力方程中找到了孤立子解，p-膜，后来在超弦中发现了在某些超对称变换下不变的超对称态，D－膜。由于保持某些超对称，他们的量子性质与相互作用强度无关。因而人们可以得到一些强耦合下的信息。人们发现上述5种超弦理论是等价的。他们都是M理论的极限，M理论在低能下的极限就是11维的超引力。
上面所及的量子场论只是在微扰的情况下有意义。这相当于在很小的尺度下经典近似是非常好的近似。反过来，当尺度变大，相互作用变强，上述理论失效。在粒子物理里，人们猜测当尺度变大，相互作用变强，从而无法把夸克分开。这就是著名的夸克幽禁猜测。这是标准模型中的核心问题之一。弦论前几年的发展为我们建立夸克幽禁开辟了一条全新的道路。
实际上，前几年超弦理论的第二次革命使我们可以系统的处理非围扰的量子场论。
在超导，超流等研究中，最困难的是处理强耦合的系统。超弦理论因为具有较高的超对称，目前还无法直接应用到超导，超流等系统中。
也许人们会认为，量子引力只在Planck尺度以下（10^{-33}cm）才起作用，这个尺度目前和我们没有多大关系。弦论前几年的进展从第一原理导出黑洞熵的公式。这对于超弦理论是强有力的实验支持。
另外，弦理论和数学有极其密切的关系。数学为弦理论提供了很多理想实验并得到许多令人惊奇的结果。
量子场论和弦论的数学基础
从70年代以来，数学和场论及弦理论发生了密切的关系。70年代中杨振宁先生的关于规范场和微分几何关系的工作,70年代末指标定理和反常的关系等起了很重要的作用.
在代数的研究中，人们发现无穷维李代数如Kac-Moody代数及其表示理论为共形场论及围扰弦理论建立了基础。而由特征标的对偶性质也可建立其它量子场论的对偶性质。Borcherds将顶点算子数学化和应用到理解例外有限群使他荣获菲尔兹奖。
80年代，在低维拓扑的研究中有若干重大突破。有些数学事实很难被理解。例如Donaldson（菲尔兹奖获得者）理论给出4维时空有无穷多种微分结构。这些结果被Witten在量子场论的框架下得到自然的解释。Donaldson不变量即是某种N＝2超对称Yang-Mills场的相关函数。后来从对偶性考虑，Seiberg-Witten引入新的不变量，
使这一理论得到极大的简化。这一对偶性对于研究弦理论中的对偶性有启发性，是引发第二次弦理论革命的重要线索。
还有许多和量子场论有关的工作,例如纽结多项式,模空间的相交理论,椭圆上同调,镜对称等.这些工作大都是考虑场论的经典解并考虑附近的量子修正得到.数学家们抛开物理背景直接从有限维构造这些理论.
我们对这种状况显然不能满意.到目前为止量子场论还没有建立起数学基础.量子场论的考虑可以提供猜测,但无法提供证明.我们希望这种状况能够改变.在量子场论的框架下直接考虑数学问题,使很多问题的理解变得直接明了.
如Witten最近在一些文章中所强调的，有两个问题是非常基本的。一个是量子Yang-Mills规范场的有限性，这可从渐进自由看出。但是目前数学上还没有证明。另一个是Yang-Mills场的质量界猜测，这和夸克幽禁有极其密切的联系。这也是Clay研究所提出的7个千年僖数学难题。目前这问题最有希望的解答是通过和弦理论的对偶得到。Maldacena前几年猜测具有极大超对称的以SU（N）为规范群的场论和某些以1/N为耦合常数的弦论对偶.这种规范场/引力对偶近两年拓展到N=2, 1 的超对称Yang-Mills场论. 夸克幽禁问题很可能在不远的将来得到解决.
Witten 建议数学家在作4维的量子场论的问题之前作2维和3维的场论.对于2维Sigma模型,质量下界对于特定情形建立起来.我们应设法拓展到广泛的情形并得到一些几何上的应用.对于3维场论他建议在Chern-Simons项前增加Yang-Mills项.这种场论的质量也应当是有下界的.
弦理论的对偶性为数学提出许多深刻的问题. 例如Sen指出弦理论的某些对偶蕴涵某些模空间上调和形式的关系.从物理学家的角度考虑，Seiberg-Witten－Donaldson的对偶性可从弦论的对偶性解释。Seiberg-Witten－Donaldson的等价性是富有挑战性的问题。也许我们需要建立某种无穷维的微积分，在这里BRST算子相当于无穷维的微分算子。Seiberg-Witten的工作相当于对于有超对称的特别的Yang-Mills场建立了夸克幽禁。
这些问题的实质性进展无疑将量子场论,弦论变为数学的一章.这是我们期待以久的.由于数学和物理长期的隔阂,在国外将两者真正结合起来作的也是凤毛麟角.这对于我们来说是个很好的机会. 我们希望中国的科学家能在此过程中继续作出贡献.

MatPy - Matrix package for Python

MatPy - Matrix package for Python

MatPy is a Python package for numerical computation and plotting with a MatLab-like interface. It currently consists of wrappers around the Numeric, Gnuplot and SpecialFuncs packages. It provides an alternative interface to NumPy that is somewhat more convenient for matrix and vector computation. Eventually both will be based directly on the same low level routines. We are also looking for the possibility of interface to Octave.

Why another package of numerical computation?

• We want Python's capability of object oriented programming and its rich data structures.

• We want Matlab/Octave style expression that is convenient for matrix computations.

• The interface of the Numeric package is more optimized for multiarray than linear algebra.

The goal of this project is to provide a package with an interface that combines the elegance of Python for OO programming and convenience of Matlab/Octave for numerical computation.

The documentations are in html, ps and pdf, generated from LaTeX source.

This is an open source project. The Project Page below has many useful links, including how to join the developement.

You can get the releases from the Download page. If you want to try the latest and greatest grab the CVS tree.

Any comments, suggestion, questions or codes are welcome. Just drop me an email. But if you're going to hack a lot, the most efficient way is to join the development team so that you can check in to the CVS tree yourself. See the SourceForge Site Documentation for how to do it.

Scilab Documentation and Support Pages

Scilab Documentation and Support Pages

Documentation and support

• On-line help. Each man is a HTML file generated from XML file. In the Scilab hierarchy they are in the man directory:
  • /man/eng/: all man files in English
  • /man/fr/: all man files with a number of them translated into French
  
  You can download below the HTML files for the stable version:
  • English version: zip, chm or tar.gz "eng" directory
  • Partially French version: zip, chm or tar.gz "fr" directory

• Next to the on line help of the software, we encourage you (specially if you are novice) to read our webpage dedicated to free books or reports to make the first move with Scilab or to consult the Scilab step-to-step web site. Scilab team express one's thanks to Jacques-Deric Rouault from NAMC laboratory for his work - CNRS UMR 8620, bat 446 Université Paris-sud, F91400 Orsay.

• Some documentation not yet updated are still available and can be useful.

last update: February 23 2006

Simula+

Le projet Simula+ réunit les laboratoires L.A.M.A.V. (EA 4015) de l'Université de Valenciennes et du Hainaut Cambrésis (UVHC) et L.P.M.M. (UMR 7554) de l'université, l'ENSAM et l'ENIM de Metz.

Objectifs Machines de calcul Rencontres LPMM/LAMAV Participants au projet Développements en cours Structure du code, aide et téléchargement Publications ayant utilisé Simula+ Présentation de Simula+ : Université de Lille 1, 10 Mars 2006 Contact : simula+@ensam.fr

Objectifs :

Ce projet consiste à développer un code orienté objet (langage de programmation C++) commun aux deux laboratoires en vue d'objectifs identifiés, à savoir :

- Permettre des simulations numériques de types éléments finis 2D/3D sur maillages non structurés, l'étude d'estimateurs à posteriori et la mise en oeuvre de techniques de raffinement de maillages, dans le cadre de problèmes anisotropes;

- Obtenir des bibliothèques performantes permettant la résolution de problèmes algébriques standards (par exemple la résolution de très grands systèmes linéaires, la recherche de valeurs propres,...), sans aucune « boite noire » sur laquelle devoir s'appuyer;

- Developper un code pour modeliser le comportement multiphysique des materiaux intelligents (alliages a memoire de forme, materiaux piezoelectriques). La loi de comportement pour ces materiaux est obtenue par une approche d'homogeneisation, ou "micro-macro", qui utilise une transition d'echelles permettant de prendre en consideration l'influence des parametres microstructuraux sur le comportement global;

- Concevoir un code qui pourra être vectorisé et/ou parallélisé en fonction des besoins futurs;

- Développer une interface graphique compatible avec tous les OS.

En outre, de par son orientation objet, Simula+ doit permettre aux futurs arrivants dans les laboratoires (étudiants, thésards, ingénieurs, enseignants chercheurs,...) de comprendre très rapidement sa structure globale et de pouvoir y travailler rapidement. Pour ce faire, sont organisées des formations au C++ pour les doctorants, et le code comporte une aide pour chacune des classes développées. Une "charte du projet" régit les modalités de développement et d'utilisation du code, afin d'assurer une bonne communication entre tous les participants, la construction d'un outil numérique de qualité, et de permettre une valorisation satisfaisante du travail accompli.

Willem Ottevanger

Willem Ottevanger
(e-mail: w.ottevanger AT math.utwente.nl)

Report:
Discontinuous Finite Element Modeling of River Hydraulics and Morphology with Application to the Paraná River

Code:
Compiling and running the code
Adding a test case
Documentation of the code using doxygen.

Presentations:
Final presentation (15.15.2005)
Build up of the code (16.9.2005)
Jonge Honden Seminar (21.7.2005)

IFISS Software Package

IFISS Software Package
This is an open-source MATLAB software package associated with the book
Finite Elements and Fast Iterative Solvers with applications in incompressible fluid dynamics
by Howard C. Elman, David J. Silvester, and Andrew J. Wathen, published by Oxford University Press.
See the publisher's web pages (UK Edition) / (US Edition) for additional information about the book.

The IFISS software can be used to generate typical linear systems arising from finite element discretizations of four important PDE applications: diffusion, convection-diffusion, Stokes flow and Navier-Stokes flow problems. It has built-in multigrid and Krylov subspace solvers and includes a variety of appropriate preconditioning strategies for each problem. We have used early versions of the software to support technical workshops we have given in the last decade on fast solvers for incompressible flow problems.

Key features include

• implementation of a variety of mixed finite element approximation methods
• automatic calculation of stabilization parameters where appropriate
• a posteriori error estimation
• a range of preconditioned Krylov subspace solvers (including MINRES and BICGSTAB(ell))
• a built-in geometric multigrid solver/preconditioner
• an interface to the algebraic multigrid solver of FEMLAB (see http://www.comsol.com/ )
• useful visualization tools.

IFISS has been tested under Matlab Versions 5.3 to 7.0 and can be run under Windows, Unix and Mac architectures. Documentation is available at ifiss_guide.pdf.

FEniCS

FEniCS is free software for the Automation of Computational Mathematical Modeling (ACMM).

Projects

FEniCS is organized as a collection of sub projects/components. Currently hosted projects are

• DOLFIN, the C++ interface of FEniCS,
• FErari, optimizes the evaluation of variational forms,
• FFC, automates the evaluation of variational forms,
• FIAT, automates the generation of finite elements,
• Ko, automates the simulation of mechanical systems,
• Puffin, a light-weight version of FEniCS for MATLAB/Octave.

To take part in the development of FEniCS, subscribe to one of the mailing lists.