if f.LorY == 3: print 'Paramenters: E = %f, thickness a = %f' % (membrane_param.mem_E, membrane_param.mem_a) print '\nTotal elements : %i' % param.N print 'Total equations: %i' % param.Neq JtoCal = 4.184 #### Compute interactions print '\nCompute interactions' computeInter(surf_array, field_array, param, membrane_param) #### Generate RHS print '\nGenerate RHS' F, F_sym, X_sym, Nblock = generateRHS(surf_array, field_array, membrane_param, Neq, electricField) print '\nRHS generated...' #### Generate matrix M, M_sym = generateMatrix(surf_array, Neq) print '\nSymbolic system:\n' counter = 0 for i in range(len(M_sym)): for j in range(len(M_sym[i])): counter += 1 buff = '' for k in range(len(M_sym[i][j])): for l in range(len(M_sym[i][j][k])): buff += M_sym[i][j][k][l]
f.E.imag) else: print 'Parameters: kappa: %f, E: %f' % (f.kappa, f.E) print '\nTotal elements : %i' % param.N print 'Total equations: %i' % param.Neq JtoCal = 4.184 #### Compute interactions print '\nCompute interactions' computeInter(surf_array, field_array, param) #### Generate RHS print '\nGenerate RHS' F, F_sym, X_sym, Nblock = generateRHS(surf_array, field_array, Neq, electricField) print '\nRHS generated...' #### Generate matrix M, M_sym = generateMatrix(surf_array, Neq) print '\nSymbolic system:\n' counter = 0 for i in range(len(M_sym)): for j in range(len(M_sym[i])): counter += 1 buff = '' for k in range(len(M_sym[i][j])): for l in range(len(M_sym[i][j][k])): buff += M_sym[i][j][k][l]
else: print 'Parameters: kappa: %f, E: %f'%(f.kappa, f.E) print '\nTotal elements : %i'%param.N print 'Total equations: %i'%param.Neq JtoCal = 4.184 #### Compute interactions print '\nCompute interactions' computeInter(surf_array, field_array, param) #### Generate RHS print '\nGenerate RHS' F, F_sym, X_sym, Nblock = generateRHS(surf_array, field_array, Neq, electricField) print '\nRHS generated...' #### Generate matrix M, M_sym = generateMatrix(surf_array, Neq) print '\nSymbolic system:\n' counter = 0 for i in range(len(M_sym)): for j in range(len(M_sym[i])): counter += 1 buff = '' for k in range(len(M_sym[i][j])): for l in range(len(M_sym[i][j][k])):
print 'Is an inner-most region' print 'Parameters: kappa: %f, E: %f'%(f.kappa, f.E) print '\nTotal elements : %i'%param.N print 'Total equations: %i'%param.Neq JtoCal = 4.184 #### Compute interactions print '\nCompute interactions' computeInter(surf_array, field_array, param) #### Generate RHS print '\nGenerate RHS' F, F_sym, X_sym, Nblock = generateRHS(surf_array, field_array, Neq) print '\nRHS generated...' #### Generate matrix M, M_sym = generateMatrix(surf_array, Neq) print '\nSymbolic system:\n' counter = 0 for i in range(len(M_sym)): for j in range(len(M_sym[i])): counter += 1 buff = '' for k in range(len(M_sym[i][j])): for l in range(len(M_sym[i][j][k])):