Skip to content
Snippets Groups Projects
Commit 36a46619 authored by Christophe Geuzaine's avatar Christophe Geuzaine
Browse files 

Adding python solver interface from Carl Osterwisch <osterwischc@asme.org>.

The interface comes with a example python "solver" that converts an Abacus
result file (in CSV format) into Gmsh post-processing views. Very cool ;-)
Thanks Carl!
parent a75e4fc5
No related branches found
No related tags found
No related merge requests found
Hide whitespace changes
Inline Side-by-side
Showing
with 1165 additions and 9 deletions
TODO
+ 5
− 1
View file @ 36a46619
$Id: TODO,v 1.78 2005-01-14 22:53:20 geuzaine Exp $
$Id: TODO,v 1.79 2005-01-20 01:25:49 geuzaine Exp $
********************************************************************
......@@ -6,6 +6,10 @@ add GUI for the Layer stuff in Extrude commands
********************************************************************
add vertex arrays for isolines in post-pro
********************************************************************
Labels:
- add dynamic variables? E.g., if a string contains %string.string,
replace it dynamically in Draw_Text2D/3D by calling
......
doc/CREDITS
+ 3
− 2
View file @ 36a46619
$Id: CREDITS,v 1.24 2005-01-01 19:35:39 geuzaine Exp $
$Id: CREDITS,v 1.25 2005-01-20 01:25:49 geuzaine Exp $
Gmsh is copyright (C) 1997-2005
......@@ -95,4 +95,5 @@ Jose Miguel Pasini <jmp84 at cornell.edu>, Philippe Lussou <plussou at
necs.fr>, Jacques Kools <JKools at veeco.com>, Bayram Yenikaya
<yenikaya at math.umn.edu>, Peter Hornby <p.hornby at arrc.csiro.au>,
Krishna Mohan Gundu <gkmohan at gmail.com>, Christopher Stott
<C.Stott@surrey.ac.uk>, Timmy Schumacher <Tim.Schumacher@colorado.edu>.
<C.Stott@surrey.ac.uk>, Timmy Schumacher <Tim.Schumacher@colorado.edu>,
Carl Osterwisch <osterwischc@asme.org>.
doc/texinfo/gmsh.texi
+ 8
− 6
View file @ 36a46619
\input texinfo.tex @c -*-texinfo-*-
@c $Id: gmsh.texi,v 1.169 2005-01-18 16:13:35 geuzaine Exp $
@c $Id: gmsh.texi,v 1.170 2005-01-20 01:25:49 geuzaine Exp $
@c
@c Copyright (C) 1997-2005 C. Geuzaine, J.-F. Remacle
@c
......@@ -404,11 +404,13 @@ evaluated at the center of its circumscribed circle/sphere.
@node Solver, Post-processing, Mesh, Overview
@section Solver: external solver interface
External solvers can be interfaced with Gmsh through Unix sockets, which
permits to easily launch external computations and to collect and exploit
the simulation results within Gmsh's post-processing module. The default
solver interfaced with Gmsh is GetDP
(@uref{http://www.geuz.org/getdp/}).
External solvers can be interfaced with Gmsh through Unix or TCP/IP sockets,
which permits to easily launch external computations and to collect and
exploit the simulation results within Gmsh's post-processing module. The
default solver interfaced with Gmsh is GetDP
(@uref{http://www.geuz.org/getdp/}). Small examples of solvers written in C,
C++, Perl and Python are available in the source distribution (in the
@file{utils/solvers/} directory).
@c -------------------------------------------------------------------------
@c Post-processing: scalar, vector and tensor field visualization
......
utils/solvers/python/Gmsh.py 0 → 100644
+ 210
− 0
View file @ 36a46619
#!/usr/bin/env python
"""Gmsh.py - Python module for interacting with the Gmsh pre/post processor.
Module functions:
getopt() -- Parse command line options
send() -- Send a string to the Gmsh socket if open
Included classes:
Client -- A class for maintaining socket communication with Gmsh.
Mesh -- A class for storing nodes and elements.
Typical use:
import sys
import Gmsh
opt = Gmsh.getopt(sys.argv[1:])
if opt.has_key('--socket'):
c = Gmsh.Client(opt['--socket'])
Gmsh.send('Testing testing', Gmsh.INFO)
Carl Osterwisch <osterwischc@asme.org>, January 2005
"""
import socket
import struct
import time
import os
import sys
START = 1
STOP = 2
INFO = 10
WARNING = 11
ERROR = 12
PROGRESS = 13
VIEW = 20
OPTION = 100
OPTION_1 = OPTION + 0
OPTION_2 = OPTION + 1
OPTION_3 = OPTION + 2
OPTION_4 = OPTION + 3
OPTION_5 = OPTION + 4
# Module variable used to store the open communication socket with Gmsh
sock = None
def getopt(opts):
"""Parse the given list (usually sys.argv) into a dict.
This is a very crude function which scans for options of the form
['-option=value'] and returns a dict of the form {'-option':'value'}.
"""
options = {}
for o in opts:
s = o.split('=', 1)
if len(s) < 2:
if s[0].startswith('-'):
# Option with no value given.
options[s[0].lower()] = None
else:
# Value with no option given.
options[None] = s[0]
else:
# Assign the given value to the option.
options[s[0].lower()] = s[1]
return options
def send(message, type=INFO):
"""Send information to Gmsh using the open socket.
Keyword arguments:
message -- the string message to send
type -- the type of message being sent (default INFO)
As of Gmsh version 1.56, type should be one of:
START STOP INFO WARNING ERROR PROGRESS VIEW OPTION
OPTION_1 OPTION_2 OPTION_3 OPTION_4 OPTION_5
If no socket has been opened (using a Gmsh.Client instance) then print the
string sys.stderr.
"""
if sock:
sock.sendall(struct.pack('II', type, len(message)))
sock.sendall(message)
elif type is not PROGRESS:
print >>sys.stderr, message
class Client:
"""Initiate and maintain connection with Gmsh socket.
Instance this class to open and initiate socket communication with Gmsh.
Communication will be automatically and gracefully terminated when the
instance is destroyed. Only one instance should exist at a time.
"""
def __init__(self, filename):
"""Open and initiate connection to Gmsh socket.
Keyword arguments:
filename -- the name of the Unix socket given by Gmsh
Assigns the new socket to the module variable sock and sends the
process id to Gmsh to initiate the communication.
"""
self.initialized = 0
global sock
if not sock:
sock = socket.socket(socket.AF_UNIX, socket.SOCK_STREAM)
time.sleep(0.1)
sock.connect(filename)
send(str(os.getpid()), START)
self.initialized = 1
else:
send('Gmsh.sock already exists as'+sock.getsockname(), ERROR)
def __del__(self):
"""Gracefully close communications upon destruction.
Sends the STOP signal to Gmsh and closes the socket.
"""
global sock
if self.initialized:
send('Done.', STOP)
sock.close()
sock = None
class Mesh:
"""This is a class for storing nodes and elements.
Members:
nodes -- A dict of the form { nodeID: [ xcoord, ycoord, zcoord] }
elements -- A dict of the form { elemID: (type, [tags], [nodeIDs]) }
Methods:
read(file) -- Parse a Gmsh version 1.0 or 2.0 mesh file
write(file) -- Output a Gmsh version 2.0 mesh file
"""
def __init__(self):
self.nodes = {}
self.elements = {}
def read(self, mshfile):
"""Read a Gmsh .msh file.
Reads Gmsh format 1.0 and 2.0 mesh files, storing the nodes and
elements in the appropriate dicts.
"""
readmode = 0
send('Reading %s'%mshfile.name)
for line in mshfile:
line = line.strip()
if line.startswith('$'):
if line == '$NOD' or line == '$Nodes':
readmode = 1
elif line == '$ELM':
readmode = 2
elif line == '$Elements':
readmode = 3
else:
readmode = 0
elif readmode:
columns = line.split()
if readmode == 1 and len(columns) == 4:
# Version 1.0 or 2.0 Nodes
try:
self.nodes[int(columns[0])] = map(float, columns[1:])
except ValueError:
send('Node format error: '+line, ERROR)
readmode = 0
elif readmode > 1 and len(columns) > 5:
# Version 1.0 or 2.0 Elements
try:
columns = map(int, columns)
except ValueError:
send('Element format error: '+line, ERROR)
readmode = 0
else:
(id, type) = columns[0:2]
if readmode == 2:
# Version 1.0 Elements
tags = columns[2:4]
nodes = columns[5:]
else:
# Version 2.0 Elements
ntags = columns[2]
tags = columns[3:3+ntags]
nodes = columns[3+ntags:]
self.elements[id] = (type, tags, nodes)
send(' %d Nodes'%len(self.nodes))
send(' %d Elements'%len(self.elements))
def write(self, file):
"""Dump the mesh out to a Gmsh 2.0 msh file."""
print >>file, '$MeshFormat\n2.0 0 8\n$EndMeshFormat'
print >>file, '$Nodes\n%d'%len(self.nodes)
for (id, coord) in self.nodes.items():
print >>file, id, ' '.join(map(str, coord))
print >>file, '$EndNodes'
print >>file, '$Elements\n%d'%len(self.elements)
for (id, elem) in self.elements.items():
(type, tags, nodes) = elem
print >>file, id, type, len(tags),
print >>file, ' '.join(map(str, tags)),
print >>file, ' '.join(map(str, nodes))
print >>file, '$EndElements'
if __name__ == '__main__':
print __doc__
utils/solvers/python/postcsv.opt 0 → 100644
+ 20
− 0
View file @ 36a46619
// These options define 'postcsv.py' as the second solver in Gmsh's
// solver module, under the name 'PostCSV'.
Solver.Name1 = "PostCSV";
Solver.Help1 = "A filter to read node or element
indexed, comma-separated-value results into Gmsh.
Carl Osterwisch <osterwischc@asme.org>
January 2005";
Solver.Executable1 = "./postcsv.py";
Solver.Extension1 = ".csv";
Solver.MeshName1 = "";
Solver.MeshCommand1 = "--mesh=%s";
Solver.SocketCommand1 = "--socket=%s";
Solver.NameCommand1 = "--results=%s";
Solver.OptionCommand1 = "--options";
Solver.FirstOption1 = "Option";
Solver.FirstButton1 = "Import";
Solver.FirstButtonCommand1 = "%s";
Solver.ClientServer1 = 1;
Solver.MergeViews1 = 1;
Solver.PopupMessages1 = 1;
utils/solvers/python/postcsv.py 0 → 100755
+ 306
− 0
View file @ 36a46619
#!/usr/bin/env python
"""Convert csv finite element data into Gmsh pos format.
Usage:
postcsv.py --mesh=gmshfile.msh --results=csvfile.csv [--average]
If --average is specified, the data for each node ID and each element ID is
averaged, otherwise repeated IDs are assumed to be new time steps.
csvfile.csv is in the following format:
# comments start with #
# time step 1 nodal scalar and vector data
node, heading1, heading2
1, 132.5, .0032 .0012 .0011
2, 100.2, .0030 .0004 -.0120
3, 98.3, .0021 .0032 .0002
# time step 1 elemental scalar data
element, heading3
1, 32.63
2, 30.34
3, 21.32
# time step 2 nodal scalar and vector data
node, heading1, heading2
1, 232.5, .0022 .0022 .0031
2, 200.2, .0020 .0014 .0110
3, 68.3, .0011 .0022 .0012
Carl Osterwisch <osterwischc@asme.org>, January 2005
"""
import Gmsh
import os, sys, re
class PostMesh(Gmsh.Mesh):
"""Store and convert FE results into Gmsh pos format.
Extends the Gmsh.Mesh class by adding the following attributes:
Members:
dataTypes -- dict to convert data length into type (scalar, vector, ...)
elemTypes -- dict to convert mesh element type into pos element type
posTypes -- list in pos defined order of pos object types
nodalData -- dict of nodal results
elementData -- dict of element centroid results
Methods:
csvparse(csvfile): Read a comma-separated-value file of nodal results
makeviews(file): Write all the data out to a Gmsh pos file
saveview(posfile): Called by makeviews to output the data for one view
"""
# dataTypes is used to translate the length of a vector into its
# corresponding Gmsh type.
dataTypes = {1:'S', 3:'V', 9:'T'}
# elemTypes is used to convert from Gmsh mesh element type to Gmsh
# post-processing object type.
elemTypes = {1:'L', 2:'T', 3:'Q', 4:'S', 5:'H', 6:'I', 7:'Y'}
# postTypes lists the valid Gmsh post-processing objects (dataType-elemType
# combinations) in the order they are defined in the pos output file
# format.
postTypes = ['SP', 'VP', 'TP',
'SL', 'VL', 'TL',
'ST', 'VT', 'TT',
'SQ', 'VQ', 'TQ',
'SS', 'VS', 'TS',
'SH', 'VH', 'TH',
'SI', 'VI', 'TI',
'SY', 'VY', 'TY',
'T2', 'T2C', 'T3', 'T3C']
def __init__(self):
Gmsh.Mesh.__init__(self)
self.nodalData = {}
#{ viewName: { nodeID: [ [time1data], [time2data] ] } }
self.elementData = {}
#{ viewName: { elemID: [ [time1data], [time2data] ] } }
# Delimiter between columns in the data file
self.separator = '[,;]'
def csvparse(self, csvfile):
"""Parse ascii file of results into elementData and nodeData.
The file should be comma separated columns with matching column
headings at the top. The first column is assumed to be the node or
element number. Each unique heading will generate a new view within
Gmsh. Repeated headings will generate time steps for that view.
"""
csv = re.compile('\s*' + self.separator + '\s*')
elem = re.compile('el', re.I)
datatype = self.nodalData
Gmsh.send('Parsing results from %s'%csvfile.name)
headings = []
for line in csvfile:
columns = csv.split(line.strip())
if len(columns) > 1:
# There are at least two columns.
linedata = []
try:
id = int(columns[0])
linedata = [map(float, column.split()) for column in columns[1:]]
except ValueError:
# Assume this is a heading.
headings = columns
if elem.match(headings[0]):
datatype = self.elementData
else:
datatype = self.nodalData
else:
for viewName, viewData in zip(headings[1:], linedata):
if not datatype.has_key(viewName):
datatype[viewName] = {}
if not datatype[viewName].has_key(id):
datatype[viewName][id] = [ viewData ]
else:
datatype[viewName][id].append(viewData)
for (viewName, nodes) in self.nodalData.items():
Gmsh.send(' %s %d Nodal Results'%(viewName, len(nodes)))
for (viewName, elems) in self.elementData.items():
Gmsh.send(' %s %d Elemental Results'%(viewName, len(elems)))
def saveview(self, posfile, viewName=None, viewData=[]):
"""Store the given viewData to the posfile.
If viewName is None it will output a header line to the file.
"""
if not viewName:
print >>posfile, '$PostFormat\n1.2 0 8\n$EndPostFormat'
elif len(viewData):
Gmsh.send(' View '+viewName)
nts = len(viewData.values()[0][0][1][0])
print >>posfile, '$View\n%s %d'%(viewName, nts)
# Print section lengths.
for type in self.postTypes:
if viewData.has_key(type):
print >>posfile, len(viewData[type]),
Gmsh.send(' %s %d'%(type, len(viewData[type])))
else:
print >>posfile, 0,
# Print dummy time-step values.
print >>posfile, '\n', ' '.join([str(i) for i in range(nts)])
# Print data.
for type in self.postTypes:
if viewData.has_key(type):
for object in viewData[type]:
for x in range(3):
# Print the coordinates in the appropriate order.
print >>posfile, ' '.join([str(i[x]) for i in object[0]])
for timestep in range(nts):
for nodalResult in object[1]:
# Print the results.
try:
print >>posfile, ' '.join([str(i)
for i in nodalResult[timestep]])
except IndexError:
Gmsh.send(' Insufficient time steps for all locations.',
Gmsh.ERROR)
Gmsh.send(' Consider using --average if this is one time step.',
Gmsh.ERROR)
print >>posfile, '$EndView'
def makeviews(self, file):
"""Generate post-processing views for the results.
1. Loop through all defined elements and generate corresponding
post-processor objects for elements with nodalData.
2. Loop through all elementData and generate corresponding
post-processor objects.
"""
self.saveview(file) # This just prints the header information
if len(self.nodalData):
Gmsh.send('Saving nodal results to '+file.name)
for viewName in self.nodalData.keys():
viewData = {}
for elem in self.elements.values():
# Iterate over all defined elements.
coord = []
result = []
for nid in elem[2]:
# Gather node coordinates and results for this element.
try:
coord.append(self.nodes[nid])
except KeyError:
Gmsh.send(' Mesh is missing node %d'%nid,
Gmsh.ERROR)
try:
result.append(self.nodalData[viewName][nid])
except KeyError:
pass
if len(result) == len(coord) and \
self.elemTypes.has_key(elem[0]):
# There was data available at every node for this
# element and this element has a corresponding
# post-processing object.
if self.dataTypes.has_key(len(result[0][0])):
postType = self.dataTypes[len(result[0][0])] + \
self.elemTypes[elem[0]]
if not viewData.has_key(postType):
viewData[postType] = []
viewData[postType].append( (coord, result) )
else:
Gmsh.send(" No data type for "+repr(result[0]),
Gmsh.ERROR)
self.saveview(file, viewName, viewData)
if len(self.elementData):
Gmsh.send('Saving elemental results to %s'%file.name)
for (viewName, data) in self.elementData.items():
viewData = {}
for (eid, edata) in data.items():
try:
elem = self.elements[eid]
except KeyError:
Gmsh.send(' Mesh is missing element %d'%eid,
Gmsh.ERROR)
else:
coord = []
for nid in elem[2]:
# Gather node coordinates and results for this
# element.
try:
coord.append(self.nodes[nid])
except KeyError:
Gmsh.send(' Mesh is missing node %d'%nid,
Gmsh.ERROR)
if len(elem[2]) == len(coord) and \
self.elemTypes.has_key(elem[0]):
# There was data available at every node for this
# element and this element has a corresponding
# post-processing object.
result = [edata]*len(coord)
if self.dataTypes.has_key(len(result[0][0])):
postType = (self.dataTypes[len(result[0][0])]
+ self.elemTypes[elem[0]])
if not viewData.has_key(postType):
viewData[postType] = []
viewData[postType].append( (coord, result) )
else:
Gmsh.send(" No data type for "+repr(result[0]),
Gmsh.ERROR)
self.saveview(file, viewName, viewData)
###############################################################################
# Main program begins here
###############################################################################
cmdline = Gmsh.getopt(sys.argv[1:])
if cmdline.has_key('--socket'):
connection = Gmsh.Client(cmdline['--socket'])
if cmdline.has_key('--help') or cmdline.has_key('-h'):
print __doc__
elif cmdline.has_key('--options'):
Gmsh.send('--noaverage', Gmsh.OPTION)
Gmsh.send('--average', Gmsh.OPTION)
else:
if not cmdline.has_key('--mesh'):
Gmsh.send("No mesh specified.", Gmsh.ERROR)
sys.exit(1)
if not cmdline.has_key('--results'):
Gmsh.send("No results specified.", Gmsh.ERROR)
sys.exit(1)
try:
mshfile = open(cmdline['--mesh'])
csvfile = open(cmdline['--results'])
except:
Gmsh.send(str(sys.exc_info()[1]), Gmsh.ERROR)
else:
pmesh = PostMesh()
pmesh.read(mshfile)
mshfile.close()
pmesh.csvparse(csvfile)
csvfile.close()
if cmdline.has_key('--average'):
Gmsh.send('Averaging data to one time step.')
for result in [ pmesh.nodalData, pmesh.elementData ]:
for view in result.keys():
for (id, data) in result[view].items():
nts = len(data)
averaged = []
for col in range(min(map(len, data))):
sum = 0.0
for d in data:
sum += d[col]
averaged.append(sum/nts)
result[view][id] = [averaged]
try:
posfile = open(os.path.splitext(csvfile.name)[0]+'.pos', 'w')
except:
Gmsh.send(str(sys.exc_info()[1]), Gmsh.ERROR)
else:
pmesh.makeviews(posfile)
posfile.close()
Gmsh.send(posfile.name, Gmsh.VIEW)
utils/solvers/python/testcsv.csv 0 → 100644
+ 200
− 0
View file @ 36a46619
1
ABAQUS VERSION 6.4-4 DATE 18-Jan-2005 TIME 17:42:16 PAGE 1
AAAAAA BBBBBBBBB AAAAAA QQQQQQQQ U U SSSSSSSS
A A B B A A Q Q U U S
A A B B A A Q Q U U S
A A B B A A Q Q U U S
AAAAAAAAAA BBBBBBBBB AAAAAAAAAA Q Q U U SSSSSSSS
A A B B A A Q Q Q U U S
A A B B A A Q Q Q U U S
A A B B A A Q Q Q U U S
A A BBBBBBBBB A A QQQQQQQQ UUUUUUUU SSSSSSSS
This program has been developed by
ABAQUS, Inc.
1080 Main Street
Pawtucket, R.I. 02860
P R O B L E M S I Z E
NUMBER OF ELEMENTS IS 218
NUMBER OF NODES IS 131
NUMBER OF NODES DEFINED BY THE USER 131
TOTAL NUMBER OF VARIABLES IN THE MODEL 262
(DEGREES OF FREEDOM PLUS ANY LAGRANGE MULTIPLIER VARIABLES)
END OF USER INPUT PROCESSING
S T E P 1 S T A T I C A N A L Y S I S
INCREMENT 1 SUMMARY
TIME INCREMENT COMPLETED 1.00 , FRACTION OF STEP COMPLETED 1.00
STEP TIME COMPLETED 1.00 , TOTAL TIME COMPLETED 1.00
E L E M E N T O U T P U T
THE FOLLOWING TABLE IS PRINTED FOR ALL ELEMENTS WITH TYPE CPE3 AVERAGED AT THE NODES
node, temperature, displacement, stress
2, 69.86, 8.4194E-04 0.000 0, -31.56 13.43 0 13.43 -329.1 0 0 0 -251.5
3, 39.96, 3.0836E-03 0.000 0, -1.491 -3.133 0 -3.133 -44.68 0 0 0 -97.55
4, 11.16, 1.9388E-03 0.000 0, 11.99 -10.21 0 -10.21 -116.0 0 0 0 -64.08
5, 13.47, 0.000 0.000 0, 37.40 1.292 0 1.292 -73.75 0 0 0 -49.52
6, 67.60, 0.000 -6.6056E-04 0, -21.23 6.614 0 6.614 -8.562 0 0 0 -148.1
7, 69.83, 8.3356E-04 -6.7009E-05 0, -21.90 33.08 0 33.08 -337.4 0 0 0 -251.5
8, 69.76, 8.0950E-04 -1.4135E-04 0, -23.03 62.61 0 62.61 -345.2 0 0 0 -254.9
9, 69.62, 7.9671E-04 -2.0956E-04 0, -40.32 70.09 0 70.09 -285.4 0 0 0 -240.2
10, 69.48, 7.5652E-04 -2.7698E-04 0, -54.81 86.73 0 86.73 -249.9 0 0 0 -234.0
11, 69.28, 7.0955E-04 -3.3576E-04 0, -63.35 94.38 0 94.38 -210.6 0 0 0 -224.6
12, 69.07, 6.6124E-04 -3.9676E-04 0, -74.58 95.83 0 95.83 -172.0 0 0 0 -216.2
13, 68.82, 6.0588E-04 -4.4704E-04 0, -74.75 84.32 0 84.32 -129.5 0 0 0 -202.5
14, 68.59, 5.3983E-04 -4.9435E-04 0, -72.84 72.02 0 72.02 -97.01 0 0 0 -191.8
15, 68.35, 4.7004E-04 -5.3725E-04 0, -70.10 59.18 0 59.18 -70.98 0 0 0 -182.3
16, 68.14, 3.9512E-04 -5.7221E-04 0, -59.27 44.19 0 44.19 -43.37 0 0 0 -170.3
17, 67.96, 3.2004E-04 -6.0584E-04 0, -53.38 33.76 0 33.76 -26.85 0 0 0 -164.6
18, 67.80, 2.4153E-04 -6.2869E-04 0, -45.07 24.73 0 24.73 -17.33 0 0 0 -157.7
19, 67.67, 1.6298E-04 -6.4554E-04 0, -35.73 16.37 0 16.37 -11.90 0 0 0 -151.8
20, 67.61, 8.2122E-05 -6.5673E-04 0, -28.73 6.691 0 6.691 -8.362 0 0 0 -148.7
21, 66.81, 9.8446E-04 0.000 0, -29.35 6.504 0 6.504 -289.6 0 0 0 -235.8
22, 63.73, 1.1361E-03 0.000 0, -40.65 5.472 0 5.472 -257.8 0 0 0 -224.4
23, 60.65, 1.2821E-03 0.000 0, -49.06 1.671 0 1.671 -208.1 0 0 0 -204.7
24, 57.54, 1.4246E-03 0.000 0, -50.58 -0.6161 0 -0.6161 -173.6 0 0 0 -187.8
25, 54.50, 1.5756E-03 0.000 0, -47.33 -1.646 0 -1.646 -152.7 0 0 0 -174.7
26, 51.55, 1.7358E-03 0.000 0, -42.15 -3.541 0 -3.541 -133.9 0 0 0 -161.3
27, 48.72, 1.9085E-03 0.000 0, -34.37 -4.128 0 -4.128 -117.9 0 0 0 -148.1
28, 46.10, 2.0976E-03 0.000 0, -28.24 -4.892 0 -4.892 -105.2 0 0 0 -137.1
29, 43.78, 2.3081E-03 0.000 0, -19.19 -5.296 0 -5.296 -89.82 0 0 0 -123.8
30, 41.81, 2.5322E-03 0.000 0, -13.39 -3.509 0 -3.509 -76.96 0 0 0 -114.7
31, 40.45, 2.7922E-03 0.000 0, -5.308 -3.571 0 -3.571 -63.00 0 0 0 -105.7
32, 39.12, 3.0185E-03 2.2115E-04 0, -2.576 -5.276 0 -5.276 -47.55 0 0 0 -96.55
33, 36.19, 2.8256E-03 4.1759E-04 0, -2.044 -9.147 0 -9.147 -65.46 0 0 0 -96.16
34, 30.51, 2.5157E-03 4.8788E-04 0, -0.7721 -14.20 0 -14.20 -90.33 0 0 0 -92.70
35, 22.25, 2.1581E-03 3.8111E-04 0, 2.853 -11.24 0 -11.24 -114.0 0 0 0 -80.06
36, 11.69, 1.3632E-03 0.000 0, 20.52 -7.606 0 -7.606 -92.20 0 0 0 -56.02
37, 12.84, 7.0866E-04 0.000 0, 37.23 -2.644 0 -2.644 -79.34 0 0 0 -47.26
38, 26.51, 0.000 1.1136E-04 0, 18.80 2.425 0 2.425 -73.57 0 0 0 -71.25
39, 36.98, 0.000 5.4140E-05 0, -1.709 3.663 0 3.663 -69.27 0 0 0 -96.91
40, 45.62, 0.000 -8.1302E-05 0, -18.69 6.841 0 6.841 -61.46 0 0 0 -118.2
41, 52.82, 0.000 -2.4234E-04 0, -31.92 9.526 0 9.526 -41.54 0 0 0 -131.9
42, 58.77, 0.000 -4.0527E-04 0, -41.24 9.057 0 9.057 -26.89 0 0 0 -143.7
43, 63.66, 0.000 -5.4837E-04 0, -33.13 7.239 0 7.239 -10.12 0 0 0 -147.6
44, 39.11, 1.4706E-03 2.1356E-04 0, -2.215 -16.24 0 -16.24 -107.1 0 0 0 -115.4
45, 51.64, 1.3526E-03 9.5846E-05 0, -22.39 -6.105 0 -6.105 -141.9 0 0 0 -157.7
46, 47.75, 5.5026E-04 -1.9718E-05 0, -12.23 17.50 0 17.50 -88.57 0 0 0 -130.3
47, 60.44, 1.2091E-03 -7.5841E-06 0, -44.20 8.877 0 8.877 -194.6 0 0 0 -197.9
48, 58.57, 2.2758E-04 -3.4992E-04 0, -35.73 26.28 0 26.28 -39.02 0 0 0 -146.2
49, 53.82, 8.5472E-04 -1.6590E-05 0, -14.35 18.67 0 18.67 -126.4 0 0 0 -154.2
50, 59.09, 1.0037E-03 -4.4003E-05 0, -31.44 22.41 0 22.41 -168.6 0 0 0 -184.0
51, 42.40, 2.0545E-03 1.6521E-04 0, -13.94 -17.69 0 -17.69 -98.98 0 0 0 -122.4
52, 29.42, 7.8391E-04 1.7522E-04 0, 11.99 -3.246 0 -3.246 -86.77 0 0 0 -85.40
53, 53.34, 1.5467E-03 5.7672E-05 0, -40.47 -7.413 0 -7.413 -148.0 0 0 0 -168.6
54, 61.16, 7.3601E-04 -1.7843E-04 0, -36.58 46.11 0 46.11 -139.8 0 0 0 -182.0
55, 64.37, 1.0845E-03 -3.5633E-05 0, -44.93 17.37 0 17.37 -250.6 0 0 0 -223.5
56, 58.10, 4.9318E-04 -2.3073E-04 0, -28.35 39.79 0 39.79 -86.64 0 0 0 -156.9
57, 33.97, 2.1313E-03 2.8372E-04 0, -2.154 -17.92 0 -17.92 -90.18 0 0 0 -98.58
58, 43.42, 1.0689E-03 1.4500E-04 0, -3.284 -3.999 0 -3.999 -112.8 0 0 0 -126.0
59, 62.17, 1.1621E-04 -4.9096E-04 0, -38.48 14.31 0 14.31 -19.57 0 0 0 -148.0
60, 63.27, 9.8677E-04 -9.1671E-05 0, -41.88 39.07 0 39.07 -213.4 0 0 0 -209.8
61, 47.33, 1.8777E-03 9.3115E-05 0, -28.91 -12.50 0 -12.50 -117.6 0 0 0 -143.9
62, 24.09, 1.4647E-03 2.1595E-04 0, 9.910 -12.59 0 -12.59 -94.44 0 0 0 -74.41
63, 43.83, 1.6560E-03 1.7798E-04 0, -12.67 -17.58 0 -17.58 -113.9 0 0 0 -129.5
64, 39.96, 2.5284E-03 1.6466E-04 0, -7.909 -12.34 0 -12.34 -77.77 0 0 0 -109.8
65, 63.48, 8.1289E-04 -1.8230E-04 0, -43.51 54.68 0 54.68 -168.6 0 0 0 -197.0
66, 55.78, 1.2850E-03 4.3255E-05 0, -32.66 -1.098 0 -1.098 -159.4 0 0 0 -174.1
67, 48.37, 1.1958E-03 1.1577E-04 0, -10.63 -4.579 0 -4.579 -128.9 0 0 0 -143.2
68, 36.16, 4.7265E-04 1.0862E-04 0, 0.7402 4.679 0 4.679 -82.59 0 0 0 -100.9
69, 56.96, 1.4146E-03 2.3575E-05 0, -46.39 -1.231 0 -1.231 -171.6 0 0 0 -185.1
70, 61.98, 5.6338E-04 -2.8586E-04 0, -38.12 51.00 0 51.00 -104.0 0 0 0 -172.5
71, 65.16, 9.8242E-04 -9.3900E-05 0, -37.02 43.60 0 43.60 -237.8 0 0 0 -219.4
72, 53.17, 1.0946E-03 4.9133E-05 0, -17.30 6.233 0 6.233 -141.0 0 0 0 -158.6
73, 48.36, 2.4764E-04 -1.1023E-04 0, -21.28 13.88 0 13.88 -63.61 0 0 0 -126.9
74, 42.44, 2.3173E-03 1.0627E-04 0, -14.60 -11.03 0 -11.03 -89.29 0 0 0 -120.0
75, 35.01, 1.1293E-03 2.0386E-04 0, 4.220 -9.489 0 -9.489 -99.26 0 0 0 -102.5
76, 64.86, 6.2305E-04 -3.2075E-04 0, -54.70 67.14 0 67.14 -112.6 0 0 0 -187.0
77, 50.21, 1.7736E-03 5.1510E-05 0, -37.30 -7.163 0 -7.163 -127.9 0 0 0 -155.1
78, 47.98, 1.6505E-03 1.2818E-04 0, -23.30 -13.91 0 -13.91 -125.9 0 0 0 -146.0
79, 60.89, 1.2466E-03 -5.0573E-06 0, -48.24 6.321 0 6.321 -208.8 0 0 0 -205.2
80, 38.36, 1.8932E-03 2.2949E-04 0, -5.402 -20.85 0 -20.85 -99.62 0 0 0 -112.0
81, 48.10, 8.1422E-04 5.0337E-05 0, -7.133 11.99 0 11.99 -109.7 0 0 0 -135.8
82, 60.72, 3.3586E-04 -3.6942E-04 0, -38.51 35.77 0 35.77 -53.49 0 0 0 -154.7
83, 39.97, 7.4362E-04 1.2307E-04 0, 0.1859 3.982 0 3.982 -94.25 0 0 0 -112.7
84, 65.54, 8.4336E-04 -1.9118E-04 0, -43.71 68.36 0 68.36 -201.5 0 0 0 -211.5
85, 32.16, 1.6376E-03 2.5023E-04 0, 2.514 -18.72 0 -18.72 -96.50 0 0 0 -95.80
86, 57.42, 6.8273E-04 -1.3380E-04 0, -22.90 35.16 0 35.16 -119.5 0 0 0 -163.4
87, 22.02, 4.5206E-04 1.1849E-04 0, 27.57 -0.5780 0 -0.5780 -78.37 0 0 0 -61.70
88, 27.06, 1.9435E-03 3.0431E-04 0, 2.174 -15.48 0 -15.48 -98.32 0 0 0 -85.54
89, 57.29, 9.1126E-04 -4.6903E-05 0, -24.34 23.77 0 23.77 -147.5 0 0 0 -171.9
90, 20.35, 1.0270E-03 1.4077E-04 0, 23.95 -6.210 0 -6.210 -86.24 0 0 0 -62.78
91, 53.78, 3.7698E-04 -1.8593E-04 0, -25.62 27.49 0 27.49 -69.66 0 0 0 -142.6
92, 53.39, 1.6126E-03 3.0206E-05 0, -45.13 -3.757 0 -3.757 -146.9 0 0 0 -170.3
93, 65.62, 7.0106E-04 -2.9006E-04 0, -55.91 74.27 0 74.27 -144.5 0 0 0 -198.0
94, 38.27, 2.2964E-03 2.2631E-04 0, -6.964 -17.67 0 -17.67 -86.84 0 0 0 -108.7
95, 66.38, 8.8448E-04 -1.6366E-04 0, -42.42 66.73 0 66.73 -239.6 0 0 0 -224.3
96, 59.54, 1.1434E-03 -1.0713E-05 0, -37.34 14.42 0 14.42 -183.4 0 0 0 -191.1
97, 45.13, 2.1104E-03 7.5866E-05 0, -23.85 -10.18 0 -10.18 -103.1 0 0 0 -132.7
98, 63.31, 2.2765E-04 -4.8975E-04 0, -41.71 25.96 0 25.96 -22.51 0 0 0 -152.7
99, 52.78, 6.2969E-04 -7.5140E-05 0, -15.72 26.76 0 26.76 -103.1 0 0 0 -147.0
100, 43.88, 1.3641E-03 1.6792E-04 0, -7.276 -12.54 0 -12.54 -119.0 0 0 0 -130.8
101, 47.33, 1.4363E-03 1.4140E-04 0, -14.20 -12.38 0 -12.38 -127.2 0 0 0 -141.6
102, 56.84, 1.3404E-03 3.2032E-05 0, -40.35 -0.8910 0 -0.8910 -168.0 0 0 0 -181.3
103, 63.68, 4.2738E-04 -4.0621E-04 0, -49.22 48.06 0 48.06 -64.35 0 0 0 -168.2
104, 50.09, 1.6750E-03 8.8554E-05 0, -33.07 -11.61 0 -11.61 -131.7 0 0 0 -154.9
105, 18.82, 1.7143E-03 1.8591E-04 0, 8.110 -9.748 0 -9.748 -106.3 0 0 0 -69.59
106, 66.56, 9.9793E-04 -3.3011E-05 0, -31.92 15.61 0 15.61 -288.6 0 0 0 -235.6
107, 40.03, 2.3416E-04 3.0424E-05 0, -8.431 6.413 0 6.413 -72.11 0 0 0 -108.2
108, 66.84, 9.5633E-04 -8.8428E-05 0, -30.30 42.13 0 42.13 -290.6 0 0 0 -236.8
109, 55.88, 1.1496E-03 2.8041E-05 0, -27.72 7.216 0 7.216 -160.8 0 0 0 -174.2
110, 47.46, 1.9630E-03 4.6760E-05 0, -30.85 -6.555 0 -6.555 -112.4 0 0 0 -142.7
111, 27.87, 1.1867E-03 2.0507E-04 0, 10.06 -11.70 0 -11.70 -92.79 0 0 0 -83.81
112, 65.42, 5.3018E-04 -3.9857E-04 0, -58.41 62.36 0 62.36 -84.20 0 0 0 -180.3
113, 55.60, 1.3108E-04 -3.0201E-04 0, -35.26 16.36 0 16.36 -41.02 0 0 0 -139.8
114, 39.83, 2.8067E-03 1.2332E-04 0, -3.517 -5.677 0 -5.677 -57.45 0 0 0 -101.9
115, 52.72, 1.4654E-03 7.8496E-05 0, -32.42 -8.561 0 -8.561 -145.4 0 0 0 -163.9
116, 42.43, 4.4960E-04 3.6423E-05 0, -7.868 10.95 0 10.95 -79.02 0 0 0 -115.2
117, 66.50, 7.5078E-04 -2.7189E-04 0, -54.98 83.04 0 83.04 -176.9 0 0 0 -209.3
118, 41.24, 2.5562E-03 7.5637E-05 0, -10.34 -5.624 0 -5.624 -75.38 0 0 0 -112.5
119, 60.45, 8.7062E-04 -1.1044E-04 0, -35.49 34.88 0 34.88 -164.6 0 0 0 -187.2
120, 61.95, 9.4380E-04 -9.8688E-05 0, -38.04 36.30 0 36.30 -183.3 0 0 0 -196.4
121, 35.48, 2.4940E-03 3.2925E-04 0, -2.175 -13.00 0 -13.00 -77.42 0 0 0 -98.28
122, 49.70, 9.8927E-04 6.7429E-05 0, -10.44 5.383 0 5.383 -126.0 0 0 0 -144.8
123, 63.08, 1.1594E-03 -1.3116E-05 0, -43.33 7.414 0 7.414 -234.5 0 0 0 -216.0
124, 56.01, 1.4914E-03 1.3061E-05 0, -49.84 -1.514 0 -1.514 -164.5 0 0 0 -182.1
125, 30.47, 3.0653E-04 1.2008E-04 0, 8.389 2.497 0 2.497 -76.10 0 0 0 -86.24
126, 58.42, 1.3757E-03 8.2717E-06 0, -49.91 0.3047 0 0.3047 -182.9 0 0 0 -192.5
127, 64.90, 3.0955E-04 -5.0667E-04 0, -50.18 35.88 0 35.88 -32.63 0 0 0 -161.6
128, 62.35, 1.0849E-03 -4.8287E-05 0, -40.94 22.78 0 22.78 -211.7 0 0 0 -206.7
129, 40.75, 2.2941E-03 1.6746E-04 0, -10.97 -16.67 0 -16.67 -89.00 0 0 0 -115.6
130, 38.08, 2.7138E-03 2.5436E-04 0, -3.250 -10.01 0 -10.01 -64.39 0 0 0 -100.3
131, 45.32, 1.8808E-03 1.3680E-04 0, -21.81 -16.15 0 -16.15 -112.0 0 0 0 -135.3
132, 63.26, 1.1019E-03 -4.0239E-05 0, -42.03 20.67 0 20.67 -226.9 0 0 0 -213.6
utils/solvers/python/testcsv.geo 0 → 100644
+ 18
− 0
View file @ 36a46619
Point(1) = {0, 0, 0, 0.1};
Point(2) = {1, 0, 0, 0.1};
Point(3) = {4, 0, 0, 0.5};
Point(4) = {4, 4, 0, 1.0};
Point(5) = {0, 4, 0, 1.0};
Point(6) = {0, 1, 0, 0.1};
Circle(1) = {2,1,6};
Line(2) = {2,3};
Line(3) = {3,4};
Line(4) = {4,5};
Line(5) = {5,6};
Line Loop(6) = {1,-5,-4,-3,-2};
Plane Surface(7) = {6};
Physical Line(101) = {5};
Physical Line(102) = {2};
Physical Line(112) = {4};
Physical Line(200) = {1};
Physical Surface(11) = {-7};
utils/solvers/python/testcsv.msh 0 → 100644
+ 395
− 0
View file @ 36a46619
$MeshFormat
2 0 8
$EndMeshFormat
$Nodes
131
2 1 0 0
3 4 0 0
4 4 4 0
5 0 4 0
6 0 1 0
7 0.9945218953682734 0.104528463267653 0
8 0.9781476007338058 0.2079116908177585 0
9 0.951056516295154 0.3090169943749462 0
10 0.9135454576426018 0.4067366430757982 0
11 0.8660254037844408 0.4999999999999963 0
12 0.8090169943749513 0.5877852522924678 0
13 0.7431448254774 0.6691306063588518 0
14 0.6691306063588659 0.7431448254773874 0
15 0.5877852522924825 0.8090169943749406 0
16 0.5000000000000111 0.8660254037844323 0
17 0.4067366430758131 0.9135454576425951 0
18 0.3090169943749584 0.95105651629515 0
19 0.2079116908177668 0.978147600733804 0
20 0.1045284632676574 0.994521895368273 0
21 1.107660966482431 0 0
22 1.230764829752582 0 0
23 1.371542876294749 0 0
24 1.532522866946499 0 0
25 1.71659942352474 0 0
26 1.927090777472968 0 0
27 2.167782956522458 0 0
28 2.443096045924229 0 0
29 2.757877025314089 0 0
30 3.117761880505327 0 0
31 3.529324854187504 0 0
32 4 0.5949066135025001 0
33 4 1.278145619573775 0
34 4 2.063201578232027 0
35 4 2.964633608956098 0
36 2.666666666666685 4 0
37 1.33333333333337 4 0
38 0 3.065647789360017 0
39 0 2.393288956030129 0
40 0 1.909232044233197 0
41 0 1.560915832167684 0
42 0 1.310258795174613 0
43 0 1.129841025580155 0
44 2.072632706027311 1.69598688973258 0
45 1.613773926086385 0.8500227436870333 0
46 0.8294669788088875 1.658857826312972 0
47 1.32513596375101 0.3476839123937975 0
48 0.3521126881546655 1.277363555436711 0
49 1.12466853437021 1.153418043518669 0
50 1.192696518113869 0.7306415166144991 0
51 2.576053729862485 0.9053498230115412 0
52 1.340126667186485 2.734096405958066 0
53 1.726422310283826 0.4059954321843577 0
54 0.9353736948832166 0.8715848956617313 0
55 1.18615207058655 0.1857246235325533 0
56 0.7052561075612429 1.175916841150846 0
57 3.128587449887757 1.766432616500146 0
58 1.502973959682794 1.639173253302476 0
59 0.1759133514351944 1.169780059600397 0
60 1.126589503584839 0.4922655384184934 0
61 2.196140328827049 0.5068889305674 0
62 2.612345788958327 2.921450948015479 0
63 2.112581290158194 1.167800481700336 0
64 3.237383018383969 0.7234638462140814 0
65 0.9840394726232626 0.6906803311775085 0
66 1.468485632956789 0.6195916609164185 0
67 1.543585212840636 1.22337304019121 0
68 0.7847572529703151 2.373911411047601 0
69 1.536612492066048 0.2546061344524544 0
70 0.7563814323183446 0.9674850965584817 0
71 1.106950139898796 0.3511105542101679 0
72 1.35889190895776 0.9986524025805176 0
73 0.3959273532525635 1.738325897497147 0
74 2.844593724249942 0.5138633427199154 0
75 1.754062226045712 2.193549370361175 0
76 0.7938495145246145 0.8047643924259087 0
77 2.001875040491023 0.2863548721871211 0
78 1.970700830475661 0.7993756598722095 0
79 1.34411212289718 0.1825137374421063 0
80 2.591371751876443 1.47181123899766 0
81 1.145266577772124 1.501068978824429 0
82 0.4941440667927327 1.148854074736235 0
83 1.149429186474111 2.031817360579769 0
84 0.9939713247765112 0.5426636213837798 0
85 2.53934709873342 2.170371417857803 0
86 0.9170333339922072 1.089415598156251 0
87 0.8299773866185869 3.344838685871387 0
88 3.289379674794993 2.516276442461847 0
89 1.137453345295398 0.9199637046208775 0
90 1.882033197768151 3.350776213150626 0
91 0.5747040540405435 1.429957302128625 0
92 1.775313987913787 0.1942086133302341 0
93 0.8643434939808741 0.7017184961964265 0
94 3.072389280032639 1.182610643612223 0
95 1.025398445390052 0.425064443790196 0
96 1.291599018349766 0.5379400141021424 0
97 2.500950996588131 0.3796311220360604 0
98 0.329249276924329 1.100214750608649 0
99 0.894324716202368 1.353915876710289 0
100 1.813941250117164 1.40844809571916 0
101 1.794878469096937 1.088460304995382 0
102 1.486769023536358 0.4427525772291758 0
103 0.5900512883039151 0.9886452832900361 0
104 1.926700328468347 0.5373153173291814 0
105 3.279545588397917 3.324720918759852 0
106 1.111952245445099 0.09914654476876072 0
107 0.387323186063824 2.191864376335631 0
108 1.078521106254799 0.230024464668231 0
109 1.35920433836455 0.7749771722215344 0
110 2.265435202084209 0.2409669473795033 0
111 2.013801144774503 2.711101094307226 0
112 0.6907699614193001 0.8590816528292089 0
113 0.2143922726422754 1.42492042583774 0
114 3.598069662403972 0.4358901037284767 0
115 1.680368656509361 0.6158390637855682 0
116 0.7133557167561021 1.999543410730382 0
117 0.9072102561528353 0.6052161494598445 0
118 3.178430153760137 0.3236365929559723 0
119 1.063899507509137 0.7888937937290703 0
120 1.1075300813509 0.6309016457330627 0
121 3.544094384876648 1.445283348829069 0
122 1.313711830501423 1.292934084685645 0
123 1.253184469019452 0.1023441938573933 0
124 1.616253868583924 0.119775562455931 0
125 0.5355241979412432 2.779799497471038 0
126 1.478800932191749 0.1185154109768451 0
127 0.4286494767521648 1.012418382992747 0
128 1.208825921620229 0.4088813796481784 0
129 2.923716641100324 0.8451929254086457 0
130 3.629773810408995 0.9401286651846473 0
131 2.278190977419342 0.7701736465013063 0
132 1.209925918352208 0.2871708528743117 0
$EndNodes
$Elements
255
1 1 2 101 5 5 38
2 1 2 101 5 38 39
3 1 2 101 5 39 40
4 1 2 101 5 40 41
5 1 2 101 5 41 42
6 1 2 101 5 42 43
7 1 2 101 5 43 6
8 1 2 102 2 2 21
9 1 2 102 2 21 22
10 1 2 102 2 22 23
11 1 2 102 2 23 24
12 1 2 102 2 24 25
13 1 2 102 2 25 26
14 1 2 102 2 26 27
15 1 2 102 2 27 28
16 1 2 102 2 28 29
17 1 2 102 2 29 30
18 1 2 102 2 30 31
19 1 2 102 2 31 3
20 1 2 112 4 4 36
21 1 2 112 4 36 37
22 1 2 112 4 37 5
23 1 2 200 1 2 7
24 1 2 200 1 7 8
25 1 2 200 1 8 9
26 1 2 200 1 9 10
27 1 2 200 1 10 11
28 1 2 200 1 11 12
29 1 2 200 1 12 13
30 1 2 200 1 13 14
31 1 2 200 1 14 15
32 1 2 200 1 15 16
33 1 2 200 1 16 17
34 1 2 200 1 17 18
35 1 2 200 1 18 19
36 1 2 200 1 19 20
37 1 2 200 1 20 6
38 2 2 11 7 75 58 44
39 2 2 11 7 80 63 51
40 2 2 11 7 59 43 20
41 2 2 11 7 72 45 67
42 2 2 11 7 101 67 45
43 2 2 11 7 80 44 63
44 2 2 11 7 59 20 19
45 2 2 11 7 85 75 44
46 2 2 11 7 131 97 51
47 2 2 11 7 83 46 81
48 2 2 11 7 83 81 58
49 2 2 11 7 98 59 19
50 2 2 11 7 83 58 75
51 2 2 11 7 43 6 20
52 2 2 11 7 100 44 58
53 2 2 11 7 97 74 51
54 2 2 11 7 79 69 47
55 2 2 11 7 93 65 54
56 2 2 11 7 89 86 54
57 2 2 11 7 111 75 85
58 2 2 11 7 95 9 71
59 2 2 11 7 89 49 86
60 2 2 11 7 113 91 73
61 2 2 11 7 113 48 91
62 2 2 11 7 108 71 9
63 2 2 11 7 84 11 10
64 2 2 11 7 102 47 69
65 2 2 11 7 118 74 29
66 2 2 11 7 108 7 106
67 2 2 11 7 108 8 7
68 2 2 11 7 91 82 56
69 2 2 11 7 88 57 34
70 2 2 11 7 118 29 30
71 2 2 11 7 131 61 97
72 2 2 11 7 91 48 82
73 2 2 11 7 86 70 54
74 2 2 11 7 95 84 10
75 2 2 11 7 127 16 103
76 2 2 11 7 127 17 16
77 2 2 11 7 76 14 13
78 2 2 11 7 121 34 57
79 2 2 11 7 98 19 18
80 2 2 11 7 93 54 76
81 2 2 11 7 109 50 96
82 2 2 11 7 113 42 59
83 2 2 11 7 90 87 52
84 2 2 11 7 99 49 81
85 2 2 11 7 83 75 52
86 2 2 11 7 76 54 70
87 2 2 11 7 73 40 41
88 2 2 11 7 115 78 45
89 2 2 11 7 85 44 80
90 2 2 11 7 59 42 43
91 2 2 11 7 77 26 27
92 2 2 11 7 119 54 65
93 2 2 11 7 131 51 63
94 2 2 11 7 112 103 15
95 2 2 11 7 120 96 50
96 2 2 11 7 83 52 68
97 2 2 11 7 109 89 50
98 2 2 11 7 90 37 87
99 2 2 11 7 109 72 89
100 2 2 11 7 115 104 78
101 2 2 11 7 99 86 49
102 2 2 11 7 123 106 22
103 2 2 11 7 113 73 41
104 2 2 11 7 94 80 51
105 2 2 11 7 92 53 69
106 2 2 11 7 109 96 66
107 2 2 11 7 97 29 74
108 2 2 11 7 90 62 36
109 2 2 11 7 106 7 2
110 2 2 11 7 93 76 13
111 2 2 11 7 125 68 52
112 2 2 11 7 106 2 21
113 2 2 11 7 122 49 72
114 2 2 11 7 103 16 15
115 2 2 11 7 108 9 8
116 2 2 11 7 95 10 9
117 2 2 11 7 122 81 49
118 2 2 11 7 106 21 22
119 2 2 11 7 101 45 78
120 2 2 11 7 105 36 62
121 2 2 11 7 120 60 96
122 2 2 11 7 113 59 48
123 2 2 11 7 100 63 44
124 2 2 11 7 112 14 76
125 2 2 11 7 87 37 5
126 2 2 11 7 85 80 57
127 2 2 11 7 132 55 79
128 2 2 11 7 115 102 53
129 2 2 11 7 125 52 87
130 2 2 11 7 88 85 57
131 2 2 11 7 110 77 27
132 2 2 11 7 92 26 77
133 2 2 11 7 124 92 69
134 2 2 11 7 89 72 49
135 2 2 11 7 102 69 53
136 2 2 11 7 125 39 107
137 2 2 11 7 122 72 67
138 2 2 11 7 100 58 67
139 2 2 11 7 129 94 51
140 2 2 11 7 109 66 45
141 2 2 11 7 123 79 55
142 2 2 11 7 99 91 56
143 2 2 11 7 122 67 58
144 2 2 11 7 115 66 102
145 2 2 11 7 86 56 70
146 2 2 11 7 111 85 62
147 2 2 11 7 104 53 77
148 2 2 11 7 116 46 83
149 2 2 11 7 97 28 29
150 2 2 11 7 87 5 38
151 2 2 11 7 116 73 46
152 2 2 11 7 125 107 68
153 2 2 11 7 92 77 53
154 2 2 11 7 91 46 73
155 2 2 11 7 115 45 66
156 2 2 11 7 107 40 73
157 2 2 11 7 132 79 47
158 2 2 11 7 99 56 86
159 2 2 11 7 93 13 12
160 2 2 11 7 109 45 72
161 2 2 11 7 130 94 64
162 2 2 11 7 103 70 56
163 2 2 11 7 131 63 78
164 2 2 11 7 127 98 18
165 2 2 11 7 117 11 84
166 2 2 11 7 130 121 94
167 2 2 11 7 123 55 106
168 2 2 11 7 92 25 26
169 2 2 11 7 88 62 85
170 2 2 11 7 95 71 60
171 2 2 11 7 122 58 81
172 2 2 11 7 112 70 103
173 2 2 11 7 102 96 47
174 2 2 11 7 114 3 32
175 2 2 11 7 123 23 79
176 2 2 11 7 101 78 63
177 2 2 11 7 114 31 3
178 2 2 11 7 128 60 71
179 2 2 11 7 101 100 67
180 2 2 11 7 125 38 39
181 2 2 11 7 126 79 23
182 2 2 11 7 129 51 74
183 2 2 11 7 112 76 70
184 2 2 11 7 121 33 34
185 2 2 11 7 118 30 31
186 2 2 11 7 103 56 82
187 2 2 11 7 88 34 35
188 2 2 11 7 105 35 4
189 2 2 11 7 126 69 79
190 2 2 11 7 127 18 17
191 2 2 11 7 105 62 88
192 2 2 11 7 98 48 59
193 2 2 11 7 104 61 78
194 2 2 11 7 118 31 114
195 2 2 11 7 115 53 104
196 2 2 11 7 94 57 80
197 2 2 11 7 120 119 65
198 2 2 11 7 111 90 52
199 2 2 11 7 108 106 55
200 2 2 11 7 101 63 100
201 2 2 11 7 117 12 11
202 2 2 11 7 131 78 61
203 2 2 11 7 123 22 23
204 2 2 11 7 111 52 75
205 2 2 11 7 120 50 119
206 2 2 11 7 107 39 40
207 2 2 11 7 104 77 61
208 2 2 11 7 110 27 28
209 2 2 11 7 102 66 96
210 2 2 11 7 99 46 91
211 2 2 11 7 113 41 42
212 2 2 11 7 110 28 97
213 2 2 11 7 112 15 14
214 2 2 11 7 105 88 35
215 2 2 11 7 95 60 84
216 2 2 11 7 117 93 12
217 2 2 11 7 99 81 46
218 2 2 11 7 116 107 73
219 2 2 11 7 110 61 77
220 2 2 11 7 111 62 90
221 2 2 11 7 119 89 54
222 2 2 11 7 132 71 108
223 2 2 11 7 132 108 55
224 2 2 11 7 118 64 74
225 2 2 11 7 98 82 48
226 2 2 11 7 124 25 92
227 2 2 11 7 120 65 84
228 2 2 11 7 129 74 64
229 2 2 11 7 126 23 24
230 2 2 11 7 128 96 60
231 2 2 11 7 120 84 60
232 2 2 11 7 119 50 89
233 2 2 11 7 125 87 38
234 2 2 11 7 128 47 96
235 2 2 11 7 110 97 61
236 2 2 11 7 117 84 65
237 2 2 11 7 116 83 68
238 2 2 11 7 121 57 94
239 2 2 11 7 124 24 25
240 2 2 11 7 130 64 114
241 2 2 11 7 127 103 82
242 2 2 11 7 129 64 94
243 2 2 11 7 117 65 93
244 2 2 11 7 130 114 32
245 2 2 11 7 132 128 71
246 2 2 11 7 118 114 64
247 2 2 11 7 126 24 124
248 2 2 11 7 130 33 121
249 2 2 11 7 116 68 107
250 2 2 11 7 126 124 69
251 2 2 11 7 132 47 128
252 2 2 11 7 127 82 98
253 2 2 11 7 90 36 37
254 2 2 11 7 105 4 36
255 2 2 11 7 130 32 33
$EndElements
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment