Problem:
F(t,y,dy)=dy+y
- Noniterative
| cmake_minimum_required(VERSION 2.8) | |
| project( mmpde6_harmonic Fortran ) | |
| find_program(ifort "ifort") | |
| # message(${ifort}) | |
| set(CMAKE_Fortran_COMPILER ${ifort}) | |
| # message(${comp}) |
| #include <stdio.h> | |
| #include <stdint.h> | |
| #include <endian.h> | |
| typedef union | |
| { | |
| __float128 val; | |
| struct | |
| { |
| (setq require-final-newline t) | |
| (add-hook 'before-save-hook 'delete-trailing-whitespace) |
| program test_function_registry | |
| use class_func_registry | |
| use flu_module | |
| use class_platonic | |
| use platonic_from_lua_module | |
| use flu_get_module | |
| use class_function | |
| use class_boundary | |
| use boundary_box_module | |
| use class_boundary_dirichlet |
| pkgname=camlimages | |
| pkgver=4.0.1 | |
| pkgrel=1 | |
| pkgdesc="Objective Caml image processing library" | |
| arch=(i686 x86_64) | |
| url="http://gallium.inria.fr/camlimages" | |
| license=('GPL2') | |
| depends=('libpng' 'libtiff' 'libxpm' 'libungif' 'freetype2' 'libjpeg') | |
| makedepends=('ocaml' 'omake') | |
| md5sums=('e2e1c7360f6a13fa1c5bbc51e7d05444') |
| default: script_embed | |
| ./script_embed | |
| script_embed: script_embed.f90 script_data.o | |
| gfortran -o $@ $^ | |
| script_data.o: script_data.s script.lua | |
| gcc -c -o $@ $< |
| default: script_embed | |
| ./script_embed | |
| script_embed: script_embed.f90 script_data.o | |
| $(FC) -o $@ $^ | |
| script_data.o: script_data.s script.lua | |
| gcc -c -o $@ $< | |
| clean: |
| $MaxExtraPrecision = 1000; | |
| der[f_, k_, rk_, h_] := | |
| NDSolve`FiniteDifferenceDerivative[Derivative[k], h Range[0, 100], | |
| f /@ (h Range[0, 100]), DifferenceOrder -> rk]; | |
| abserr[f_, k_, rk_] := | |
| Abs[der[f, k, rk, #][[2]] - | |
| Derivative[k][f][#]] & /@ (2^-Range[1, 30]); | |
| relerr[f_, k_, rk_] := | |
| abserr[f, k, rk]/Abs@(Derivative[k][f] /@ (2^-Range[1, 30])); | |
| graph[{f1_, f2_}, k_, rk_] := |