笑死：Welcome to Skip Thompson's Homepage
Welcome to Skip Thompson's Homepage
有一位退休的姓“汤姆森”的老教授，他老人家的个人主页名字如上，直接翻译过来就是“欢迎跳过（略过，忽略）汤姆森的主页”。
我发现这个名字还是挺不错的。
I am an Emeritus Professor of Mathematics and Statistics at Radford University in Radford, Virginia. Please click on the following links for various information, downloads, and the like.
|
这里的东西挺不错的。
|
VODE_F90 Support Page
G.D Byrne
Departments of Applied Mathematics and Chemical & Environmental Engineering
Illinois Institute of Technology
Chicago, IL 60616
gbyrne@wi.rr.com
S. Thompson
Mathematics & Statistics Department
Radford University
Radford, VA 24142
thompson@radford.edu
This page contains the source code and other downloadable material for the VODE_F90 ordinary differential equation (ode) solver. VODE_F90 is a Fortran 90 version of the well-known VODE ode solver by Brown, Byrne, and Hindmarsh. Please contact one of the authors if you have any questions or suggestions about the information on this page.
To date VODE_F90 has been tested with several commericial F90 and F95 compilers. However, we are particularly interested in learning of any difficulties you encounter using your F90 or F95 compiler(s).
Please click on either of the following links for more information.
- VODE_F90 Source Code
- Demo Programs
- VODE_F90 Documentation
- Kyle Niemeyer's OpenMP compatible version of dvode_f90_m.f90
- Choosing Error Tolerances
Click here for information about DVODE_F90's cousins, CVODE and PVODE.
Last updated: 08-27-13
Overview of the Demonstration Programs
G.D. Byrne
Departments of Applied Mathematics and
Chemical & Environmental Engineering
Illinois Institute of Technology
Chicago, IL 60616
gdbyrne847@yahoo.com,
S. Thompson
Department of Mathematics and Statistics
Radford University
Radford, VA 24142,
thompson@radford.edu
Overview
The demonstration programs illustrate the use of the Fortran 90 ordinary differential equation (ODE) VODE_F90 which is an extension of the well-known f77 VODE solver. VODE_F90 is applicable to initial value problems of the form
|
Several extensions were incorporated in VODE_F90. In particular, it contains an option for root finding ("event location"); and it contains an option to accommodate sparse Jacobians. It has also an option to impose bounds on the solution y(t) . The demonstration programs illustrate the use of these options as well as others.
Most of the demonstration programs have an accompanying MATLAB script for plotting the solution generated by the program. To run the demonstration programs you will need to first obtain the solver file vode_double_m.f90 (or vode_single_m.f90) from the solver source code directory. If you encounter any difficulties with the demonstration programs or have questions regarding them, please contact the authors. Please note that the demo programs are still being refined; so they may change at any time.
The programs reside in directories each containing subdirectories for different problems.Following is a brief summary of what the various demonstration programs do and what each is intended to illustrate. If you click on the link for any demo program, you can download it and other related files.
Miscellaneous Demos |
- DemoHarmonic: The program demoharmonic.f90 solves a problem involving two coupled harmonic oscillators from the book Solving ODEs with MATLAB by Shampine, Gladwell, and Thompson. Root finding is used to construct Poincare maps for the system. The script demoharmonic.m generates a MATLAB plot of the Poincare maps.
- DemoLaser: The program demolaser.f90 solves a model for a ruby laser from the paper "Stiff ODE Solvers: A Review of Current and Coming Attractions" by G.D. Byrne and A.C. Hindmarsh, J. Comput. Phys., Vol.70, 1 (1987). The problem is stiff initially and is oscillatory later. The script demolaser.m generates MATLAB plots of the solution components.
- DemoFran: In the paper "Bifurcations for Tori and Phase Locking in a Dissipative System of Differential Equations" by V. Franceschini, (Physica 6D, 1983, pp. 285-304), the author considers properties of attractors for the solution. The shapes of the attractors change as a critical system parameter changes. The program demofran.f90 uses root finding to reproduce several of the attractors. The script demofran.m generates MATLAB plots of the attractors. Note that since many thousands of points for several attractors are located using root finding, this program takes a few minutes to run.
- DemoKnee: This program solves a problem discussed in the paper "Non-negative Solutions of ODEs" by Shampine, Thompson, Kierzenka, and Byrne. The problem is a well-known model in which one of the components becomes negative as it rounds the bend in a "knee" in the problem. It causes difficulty for several ODE solvers. The program demoknee.f90 enforces non-negativity of the solution to prevent this and to ensure that the numerical solution follows the correct solution branch after the bend. The script demoknee.m generates a MATLAB plot of the solution.
- Demosf: In the paper "Efficient Automatic Integration of Ordinary Differential Equations with Discontinuities," by D. Ellison, Math. Comp. Simulation, Volume XXIII, 1981, pp. 190-196, the motion of a block sliding across a surface is considered. The program demosf.f90 uses root finding to track the motion of the block by locating the times at which the block sticks and unsticks. The script demosf.m generates MATLAB plots of the solution components.
MOL Demos |
- DemoHumCol: This problem solves a model of a humidification column from the book Dynamic Modeling of Transport Process Systems by C.A. Silebi and W.E. Schiesser. Bill Schiesser kindly allowed us to translate his f77 DSS/2 spatial differentiation routines to Fortran 90 and include the translated version for use with this program. The program demohum21.f90 recomputes the sparsity pattern after it changes in this method of lines model and solves the desired problem. The script demohum21.m generates MATLAB plots of selected solution components. DSS/2 subroutines are used to generate line plots of the solution, eigenvalues of the Jacobian matrix, and Jacobians maps. We note that DSS/2 is now a part of the programs in the above book and that the original f77 version is now freely available here . In addition, the same site contains information on an entire library of models amenable to solution in a manner similar to this problem. You can do no better than spend time checking out the library of available models to get a taste of applied modelling of PDEs at its best.
- DemoPolymer: The program demop11.f90 solves a model of a convective moving polymer sheet communicated to us by W.E. Schiesser for use in the paper "Non-negative Solutions of ODEs" by Shampine, Thompson, Kierzenka, and Byrne. Near the beginning of the integration the temperatures oscillate outside the physically meaningful region. This program imposes bounds on the solution components to prevent this. The script demop11.m generates a MATLAB plot of the solution which demonstrates that the bounds are satisfied.
- DemoFlow: The program demofow.f90 solves a fluid flow problem. It is required to obtain the numerical steady state solution for Euler's equations starting from an initial guess. A second program demoflow2.f90 solves the same problem using a different set of boundary conditions. The scripts demoflow.m and demoflow2.m generate MATLAB plots of the solution.
Pitagora Demos |
The programs in this directory solve problems from the well-known Pitagora test suite . The authors kindly permitted us to include their pdf descriptions of the problems solved by these programs. Refer to these descriptions for more details about the problems.
- DemoE5: Problem E5 consists of a stiffsystem with four ODEs. The program demoe5.f90 solves the problem. The script demoe5.m generates MATLAB plots of the solution components.
- DemoEmep: Problem EMEP consists of a stiff system with 66 ODEs. The program demoemep.f90 solves the problem. The script demoemep.m generates MATLAB plots of selected solution components.
- DemoHires: Problem HIRES consists of a stiff system with 8 ODEs. The program demohires.f90 solves the problem on the interval [0,6] . The script demohires.m generates MATLAB plots of selected solution components on the interval [0,6] . The program demohires2.f90 solves the problem on the original integration interval [0,321.8122] .
- DemoMedakzo: Problem MEDAKZO is a stiff MOL system consisting of 400 ODEs. The program demomedakzo.f90 solves this problem. The script demomedakzo.m generates MATLAB plots of the solution surfaces.
- DemoOrego: Problem OREGO is the Oregonator system which consists of three stiff equations and whose components have very sharp interfaces. The program demoorego.f90 solves this problem. The script demoorego.m generates MATLAB plots of the solution components.
- DemoRobertson: The program demorobertson.f90 solves the well-known Robertson chemical kinetics problem. Non-negativity of the solution components is enforced.
Toronto Demos |
The programs in this directory solve problems from the well-known Toronto (DETEST) test suites.
- DemoStiff: The program demostiff.f90 solves the problems in the stiff test suite.
- DemoNonstiff: The program demononstiff.f90 solves the problems in the nonstiff test suite.
ODEPACK Demos |
The programs in this directory consist of demo programs for the ODEPACK family of integrators. These programs are included primarily for people familiar with their use in conjunction with the ODEPACK solvers and who would like to compare the results with those for the original solvers.
- DemoFoodWeb: The program demofw.f90 solves one of the f77 demo programs for LSODPK. It is a large, sparse problem.
- DemoDiurnal: The program demodiurn.f90 solves another of the f77 demo programs for LSODPK. It is a large, sparse problem.
- DemoOdepackAR: The two programs in this directory, demoar1.f90 and demoar2.f90 solve two ODEPACK demo problems for LSODAR.
- DemoOdepackDB: The program demodb.f90 solves two of the ODEPACK demos for VODE, the f77 counterpart of VODE_F90.
- DemoOdepackSP: The program demosp.f90 solves three of the ODEPACK demo problems for LSODES.
Single Precision Demo Programs
Several of the programs described above have single precision counterparts. Not all of the double precision demos were converted to single precision because we typically use double precision on computers with short single precision word lengths. Following is a list of the available single precision demos.
- DemoRobertson
- DemoHumCol
- DemoPolymer
- DemoDiurnal
- DemoHarmonic
- DemoFran
- DemoKnee
- DemoOdepackAR
- DemoOdepackDB
- DemoOdepackSP
- DemoToronto
- Distribution Example 1
- Distribution Example 2
笑死：Welcome to Skip Thompson's Homepage相关推荐
- 20多年前微软曾计划收购，任天堂嘲讽道：“笑死我了”
前不久,Xbox主机迎来20周年纪念.在这次纪念活动中,微软大方晒出了一封当年公司副总裁写给任天堂高管的一封信,内容是说服对方同意被收购. 20多年前,当微软准备发布第一代Xbox时,公司对他们缺少重 ...
- 赵本山说日本车。（笑死不偿命）
赵本山说日本车.(笑死不偿命) 赵本山上场说:" 车到山前必有路,可惜丰田刹不住! 丰田就是刹的住,吓得车主尿一裤! 开上丰田雅力士,重症病房把伤治! 买了丰田凯美瑞.太平间里安心睡! 三轮 ...
- 昨晚，B站崩了！看了网友们的评论，我差点笑死...
昨晚,刚洗好澡,突然被群友@出来,说聊聊容灾的话题... 仔细一看,原来是B站挂了... 在很多群里,看到了各种404.500等报错截图: 看看程序员们的反应是这样的: 有把股票联系起来的 有把过去技 ...
- 笑死了，打错一个字，竟然...
淘宝上,让人快活的打字错误: -有大妈吗?(大码) -亲,客服最大的27岁 -你能活到付款吗?(货到付款) -亲,我尽量 -你们有尸体店吗?(实体店) -亲,淘宝不让卖那个-- -你什么时候发火啊?( ...
- php微信段子,年度挤进前十名的微信段子，笑死了
年度挤进前十名的微信段子,笑死了 选送:向国敏 第十名 昨天中午有个男同事外出,没把手机带走,他老婆不停地打电话来.午睡的女同事被吵烦了,拿过手机大吼:"我们在睡觉,你烦 ...
- java 搞笑翻译_超级搞笑的英文翻译图文版！笑死我了
原标题:超级搞笑的英文翻译图文版!笑死我了 I never change whatever U change 电流不随电压改变而改变,当时就把这句话死死的记住了. 2.If you think you ...
- 不同部门员工吃饭时聊些什么，程序员那一段笑死我了…
不同部门员工吃饭时聊些什么,程序员那一段笑死我了- 分类 程序员笑话 身为程序员,大家和同事吃饭的时候喜欢聊天吗?大家又在聊些什么呢?有没有想过,你们聊天的内容已经暴露了自己的职业哦. (一) 甲:你 ...
- 请领导过目文件怎么说_领导说“辛苦了”怎么回才显情商高？网友回复笑死
工作中遇到领导说"辛苦了",该怎么回复? 一般人会说"没事,应该的",但也有人"不一般",非要别出心裁凸显"情商高",还 ...
- 唐人码输入法破解版，笑死我了
单位一老同志需要用唐人码,可是TMD正版只能用5次,还要用软盘,不会安装,我说到网上给下个破解版吧, 找了半天,居然没有免费的,到时有几个叫卖硬盘版,宣称不要软盘,一机任意安装多次,还有一个地址可以下 ...
最新文章
- 扩展jquery实现客户端表格的分页、排序
- [CQOI2010]扑克牌
- OpenMP并行化实例----Mandelbrot集合并行化计算
- asp.net高亮显示搜索关键字的两个方法
- react gps坐标转换_手持GPS的三参数计算方法
- 干货 | 携程数据库发布系统演进之路
- jQuery面试题-区别mouseover和mouseenter的不同之处(看了也许对你有好处)
- MVC 中 Razor 无限分类的展示
- ipython安装成功却无法运营_pyspider显示安装成功但仍无法运行
- 比尔盖茨为什么能成为世界首富？
- Android面试，BroadCastReceiver的两种注册方式的异同
- 视差滚动教程、实现方法、源代码
- windowns10安装httpd
- 水利知识库AI助手+公务员资料【产品经理】
- matlab群速度色散系数,第三章-群速度色散.ppt
- 13位知名科技公司CEO首份工作揭秘
- （2022杭电多校三）1002-Boss Rush（状压DP+二分）
- CyberDuck是一个很棒的 FTP SFTP图形化工具
- 11.3 树的遍历：LDR，LRD，VLR 相关代码
- java代码控制电脑发出警报_java – 如何使用警报构建器对话框显示消息