Python generateRHS Exemples

Exemple #1

    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]

Exemple #2

Fichier : main_matrix.py Projet : lizhen-dlut/pygbe

                                                    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]

Exemple #3

Fichier : main_matrix.py Projet : LEONOB2014/pygbe

    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])):

Exemple #4

        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])):