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_] := |