Step 23: Renumber the nodes to reduce bandwidth
Use my utility program OrbNumber,
which runs in text console mode and requires very few decisions from the user.
It will read in your current .feg file, change all the nodes numbers
(but nothing else), and write a modified .feg file.
Sample output is below (with some lines omitted in the middle).
Note that the bandwidth was reduced from 1195 to 44.
When NeoKinema is run, this can speed up the calculation by as much as a
factor of (1195/44)2 = 738,
and also can reduce memory usage (for the banded linear system) by a factor of
(1195/44) = 27.
Initializing the
clock;
Attempting to read existing .feg file on UNIT 1...
File name missing or blank - please enter file name
UNIT 1?
Iberia1.feg
TITLE OF FINITE ELEMENT GRID =
Iberia1: Iberian margin, with faults
(SINCE OPTION BRIEF=.TRUE., GRID WILL NOT BE ECHOED HERE. BE CAREFUL!!!)
SINCE FAULT OFFSETS ARE ALL ZERO, INPUT PARAMETER BYERLY WILL HAVE NO EFFECT
....DONE.
6.00 SECONDS HAVE ELAPSED SINCE LAST CHECK.
SQUARE IS CHECKING FOR CORRECT TOPOLOGY...
IN SQUARE: ARE ALL NODES CONNECTED?
IN SQUARE: AVERAGING NODE COORDINATES
IN SQUARE: CHECKING FOR NEGATIVE AREAS
IN SQUARE: CREATING LIST OF BOUNDARY NODES
IN SQUARE: ORDERING LIST OF BOUNDARY NODES
....DONE.
1.00 SECONDS HAVE ELAPSED SINCE LAST CHECK.
CALLING NUMBER...
CREATING THE GRAPH MATRIX...
....DONE.
0.00 SECONDS HAVE ELAPSED SINCE LAST CHECK.
DETERMINING THE ORIGINAL BANDWIDTH...
....DONE.
0.00 SECONDS HAVE ELAPSED SINCE LAST CHECK.
NUMBERING FROM INITIAL NODE = 1
ATTEMPT NUMBER 1 BEGINS WITH NODE 1.
WHEN THIS NODE IS #1, THE BANDWIDTH IS 44.
NUMBERING FROM INITIAL NODE = 2
ATTEMPT NUMBER 2 BEGINS WITH NODE 2.
WHEN THIS NODE IS #1, THE BANDWIDTH IS 49.
NUMBERING FROM INITIAL NODE = 6
ATTEMPT NUMBER 3 BEGINS WITH NODE 6.
WHEN THIS NODE IS #1, THE BANDWIDTH IS 52.
NUMBERING FROM INITIAL NODE = 14
ATTEMPT NUMBER 4 BEGINS WITH NODE 14.
WHEN THIS NODE IS #1, THE BANDWIDTH IS 56.
NUMBERING FROM INITIAL NODE = 26
ATTEMPT NUMBER 5 BEGINS WITH NODE 26.
WHEN THIS NODE IS #1, THE BANDWIDTH IS 109.
NUMBERING FROM INITIAL NODE = 40
ATTEMPT NUMBER 6 BEGINS WITH NODE 40.
WHEN THIS NODE IS #1, THE BANDWIDTH IS 58.
NUMBERING FROM INITIAL NODE = 60
ATTEMPT NUMBER 7 BEGINS WITH NODE 60.
WHEN THIS NODE IS #1, THE BANDWIDTH IS 57.
NUMBERING FROM INITIAL NODE = 84
ATTEMPT NUMBER 8 BEGINS WITH NODE 84.
WHEN THIS NODE IS #1, THE BANDWIDTH IS 57.
=-=-=-=-=-=- (lines omitted) -=-=-=-=-=-=-=-=
NUMBERING FROM INITIAL NODE = 39
ATTEMPT NUMBER 144 BEGINS WITH NODE 39.
WHEN THIS NODE IS #1, THE BANDWIDTH IS 48.
NUMBERING FROM INITIAL NODE = 25
ATTEMPT NUMBER 145 BEGINS WITH NODE 25.
WHEN THIS NODE IS #1, THE BANDWIDTH IS 47.
NUMBERING FROM INITIAL NODE = 13
ATTEMPT NUMBER 146 BEGINS WITH NODE 13.
WHEN THIS NODE IS #1, THE BANDWIDTH IS 47.
NUMBERING FROM INITIAL NODE = 5
ATTEMPT NUMBER 147 BEGINS WITH NODE 5.
WHEN THIS NODE IS #1, THE BANDWIDTH IS 84.
NUMBERING FROM INITIAL NODE = 1
OLD NODES 494 AND 538 HAVE NEW NUMBERS 495 AND 539
WHOSE DIFFERENCE OF 44 DEFINES THE NEW BANDWIDTH.
OLD NODES 523 AND 566 HAVE NEW NUMBERS 523 AND 567
WHOSE DIFFERENCE OF 44 DEFINES THE NEW BANDWIDTH.
OLD NODES 534 AND 577 HAVE NEW NUMBERS 535 AND 579
WHOSE DIFFERENCE OF 44 DEFINES THE NEW BANDWIDTH.
OLD NODES 535 AND 580 HAVE NEW NUMBERS 537 AND 581
WHOSE DIFFERENCE OF 44 DEFINES THE NEW BANDWIDTH.
OLD NODES 579 AND 623 HAVE NEW NUMBERS 580 AND 624
WHOSE DIFFERENCE OF 44 DEFINES THE NEW BANDWIDTH.
OLD NODES 619 AND 663 HAVE NEW NUMBERS 620 AND 664
WHOSE DIFFERENCE OF 44 DEFINES THE NEW BANDWIDTH.
ATTEMPT NUMBER 1 BEGINS WITH NODE 1.
WHEN THIS NODE IS #1, THE BANDWIDTH IS 44.
....DONE.
68.00 SECONDS HAVE ELAPSED SINCE LAST CHECK.
BANDWIDTH BEFORE RENUMBERING = 1195
BANDWIDTH AFTER RENUMBERING = 44
RENUMBERING COMPLETED.
==========================================
Attempting to write output .feg file on UNIT 2...
File name missing or blank - please enter file name
UNIT 2?
Iberia2.feg
|
You can open the newly-renumbered .feg file in OrbWin, if you
like, to confirm that it looks just the same as before.
Only the node-numbers which appear in the lower-right of the screen (during
certain commands) will have changed.
Now, return to your NeoKinema-parameter file (introduced in Step #8; sample file is parameters_for_NeoKinema.nki.txt)
and open it with a plain-ASCII text editor such as NotePad or EditPad
Pro.
In line #20, replace the name of the .feg file with the filename of this
newly-renumbered .feg file from OrbNumber.
Note that, when OrbNumber is finished, it will ALSO provide you with
a new file listing the numbers (and locations) of the nodes around the grid
boundary.
Name this file something like yourFEGridName.per (where .per
stands for “perimeter”).
Save this (plain-ASCII) file, which will be very useful soon, in Step 26 of this Guide.