http://didgeridoo.une.edu.au/womwiki/ 2026-05-20T01:41:40+1000 http://didgeridoo.une.edu.au/womwiki/ http://didgeridoo.une.edu.au/womwiki/lib/tpl/vector/images/favicon.ico text/html 2018-07-02T12:37:16+1000 Anonymous (anonymous@undisclosed.example.com) http://didgeridoo.une.edu.au/womwiki/doku.php?id=fortran:basicsplines&rev=1530499036&do=diff Click tab to download ! BSPLINES ! Purpose: evaluate B-spline basis functions ! linear, quadratic or cubic ! equidistant (default) or at user-specified knots ! Input: via command line options (-h will list them) or interactively ! program will prompt for esential options not given on command line; ! ages and knot positions can be given in separate files ! Choose ! degree: command line option -d1, -d2 or -d3 for linear, quadratic or cubic… text/html 2018-08-01T12:03:37+1000 Anonymous (anonymous@undisclosed.example.com) http://didgeridoo.une.edu.au/womwiki/doku.php?id=fortran:fortran&rev=1533089017&do=diff Karin's FORTRAN gems This page gathers some small programs or routines, with limited documentation. If you like it, cite it! B-splines Simple program to generate B-splines functions for given range or given values of covariable, with choice of equi-distant splines or user-defined knots. text/html 2018-07-02T12:37:16+1000 Anonymous (anonymous@undisclosed.example.com) http://didgeridoo.une.edu.au/womwiki/doku.php?id=fortran:frontgeninv&rev=1530499036&do=diff ! calculate generalised inverse of real symmetric, halfstored matrix ! input file: non-zero elements of matrix, one line per element ! -> row no., column no., element (space separated) ! program will ask for file name & no. of rows/columns ! karin meyer program ginverse implicit none real(8), dimension(:), allocatable :: avec real(8) :: xx, det, zero=1.d-9 character(len=30) :: fname … text/html 2018-07-02T12:37:16+1000 Anonymous (anonymous@undisclosed.example.com) http://didgeridoo.une.edu.au/womwiki/doku.php?id=fortran:geninverse&rev=1530499036&do=diff ! ================================================ subroutine dkmxhf( a, det, zero, nrank, n, iopt) ! ================================================ implicit none integer, intent(in) :: n ! no. of rows/columns logical, intent(in) :: iopt ! .true./.false. : print warnings real(8), dimension(n*(n+1)/2), intent(inout) :: a ! matrix; upper triangle, stored row-wise real(8), intent(in) :: zero ! operat… text/html 2018-07-02T12:37:16+1000 Anonymous (anonymous@undisclosed.example.com) http://didgeridoo.une.edu.au/womwiki/doku.php?id=fortran:harville&rev=1530499036&do=diff module inv_harville !© karin meyer integer, save :: nrczhz, maxrec integer, dimension(:),allocatable, save :: ifirst, ivcol, inext real(8), dimension(:),allocatable, save :: dia, zhz integer, dimension(:), allocatable, save :: ivec real(8), dimension(:), allocatable, save :: xvec contains ! ========================================================== subroutine harville( nrow, ngibbs, nburn, iseed1, iseed2 )… text/html 2018-07-02T12:37:16+1000 Anonymous (anonymous@undisclosed.example.com) http://didgeridoo.une.edu.au/womwiki/doku.php?id=fortran:harvilletest&rev=1530499036&do=diff !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! ! example matrix ! 2.000 1.000 1.000 0.000 0.000 ! 1.000 3.000 0.000 1.000 1.000 ! 1.000 0.000 3.000 0.000 1.000 ! 0.000 1.000 0.000 1.000 1.000 ! 0.000 1.000 1.000 1.000 2.000 ! exact inverse ! 1.000 -0.500 -0.500 0.000 0.500 ! -0.500 0.750 0.250 -0.500 -0.250 ! -0.500 0.250 0.750 0.500 -0.750 ! 0.000 -0.500 0.500 3.000 -1.500 ! … text/html 2018-07-02T12:37:16+1000 Anonymous (anonymous@undisclosed.example.com) http://didgeridoo.une.edu.au/womwiki/doku.php?id=fortran:hashfront&rev=1530499036&do=diff program xxx use hash_tables implicit none integer, parameter :: nex = 10 integer, dimension(nex) :: my_ids =(/ 772, 9999, 263643, 789, 1281, 66661, & & 55, 12132, 789, 111111 /) integer :: i, idorig, idnew, ipos ! define how many lists mlists = 1 ! allocate arrays call hash_setup ! recode in running order do i = 1, nex idorig = my_ids(i) call hash_it(idorig, 1, idnew, ipos, .false.) ! call hash_it(idorig, ivar, id… text/html 2018-07-02T12:37:16+1000 Anonymous (anonymous@undisclosed.example.com) http://didgeridoo.une.edu.au/womwiki/doku.php?id=fortran:hashtabs&rev=1530499036&do=diff MODULE hash_tables ! mlevels and step define the hash function; set to prime numbers as required integer, parameter :: step = 123 integer, parameter :: mlevels = 99139 ! examples of prime numbers ! small(ish) : 50021, 99139 ! medium : 333667, 524827 ! large : 1299877, 2000003, ! huge : 4000079, 6026179, 8034709 integer, save :: mlists integer, dim…