import numpy as np
#import matplotlib
import matplotlib.pyplot as plt
import numeric
import helper
from helper import Timer
import fem
import assembly
#def rho(x): return 0#math.exp(-0.5/4**2*((x[0]-4)**2+(x[1]-12)**2))*1e-10-math.exp(-0.5/4**2*((x[0]-20)**2+(x[1]-12)**2))*1e-10
#def eps(x): return 1#return if ((x[0] > x1 ) and (x[0] < x2) and (x[1] > y1) and (x[1] < y2))#1
with Timer('importing solution'):
nodes_to_coordinates = helper.read_matrix_from_file(
'data/nodes_to_coordinates.txt')
local_to_global = helper.read_matrix_from_file('data/local_to_global.txt')
boundary_table = helper.read_matrix_from_file('data/boundary_table.txt')
#print nodes_to_coordinates
u_h = helper.read_matrix_from_file('data/u_h.txt')
with Timer('Calculation 1'):
energy = 0
# phi1 = (1-xx) * (1-yy)
# phi2 = xx * (1-yy)
# phi3 = xx * yy
# phi4 = (1-xx) * yy
import math
import preparation
import assembly
import solve
import plots
import helper
from helper import Timer
ELEMENTS_PER_LINE = 12
NODES_PER_ELEMENT = 4
STEP_SIZE = 0.2
points = helper.read_matrix_from_file('data/important_points.txt')
print('x1='+str(points[0,0])+'y1='+str(points[0,1]))
x1 = points[0,0]
y1 = points[0,1]
x2 = points[1,0]
y2 = points[1,1]
window = helper.read_matrix_from_file('data/window.txt')
wx1 = window[0,0]
wy1 = window[0,1]
wx2 = window[1,0]
wy2 = window[1,1]
def rho(x):
import math
import preparation
import assembly
import solve
import plots
import helper
from helper import Timer
ELEMENTS_PER_LINE = 12
NODES_PER_ELEMENT = 4
STEP_SIZE = 0.2
points = helper.read_matrix_from_file('data/important_points.txt')
print('x1=' + str(points[0, 0]) + 'y1=' + str(points[0, 1]))
x1 = points[0, 0]
y1 = points[0, 1]
x2 = points[1, 0]
y2 = points[1, 1]
window = helper.read_matrix_from_file('data/window.txt')
wx1 = window[0, 0]
wy1 = window[0, 1]
wx2 = window[1, 0]
wy2 = window[1, 1]
#import matplotlib
import matplotlib.pyplot as plt
import numeric
import helper
from helper import Timer
import fem
import assembly
#def rho(x): return 0#math.exp(-0.5/4**2*((x[0]-4)**2+(x[1]-12)**2))*1e-10-math.exp(-0.5/4**2*((x[0]-20)**2+(x[1]-12)**2))*1e-10
#def eps(x): return 1#return if ((x[0] > x1 ) and (x[0] < x2) and (x[1] > y1) and (x[1] < y2))#1
with Timer('importing solution'):
nodes_to_coordinates = helper.read_matrix_from_file('data/nodes_to_coordinates.txt')
local_to_global = helper.read_matrix_from_file('data/local_to_global.txt')
boundary_table = helper.read_matrix_from_file('data/boundary_table.txt')
#print nodes_to_coordinates
u_h = helper.read_matrix_from_file('data/u_h.txt')
with Timer('Calculation 1'):
energy = 0
# phi1 = (1-xx) * (1-yy)
# phi2 = xx * (1-yy)
