# 28.05.2009, c
from sfepy import data_dir

filename_mesh = data_dir + '/meshes/3d/elbow.mesh'

fields = {
    'scalar' : ('real', 'scalar', 'Omega', 1),
    'vector' : ('real', 'vector', 'Omega', 1),
}

integrals = {
    'i1' : ('v', 'gauss_o2_d3'),
    'i2' : ('s', 'gauss_o2_d2'),
}

regions = {
    'Omega' : ('all', {}),
    'Gamma' : ('nodes of surface', {'can_cells' : True}),
}

expressions = {
    'volume_p' : 'd_volume.i1.Omega( p )',
    'volume_u' : 'd_volume.i1.Omega( u )',
    'surface_p' : 'd_volume_surface.i2.Gamma( p )',
    'surface_u' : 'd_volume_surface.i2.Gamma( u )',
}

fe = {
    'chunk_size' : 1000
}

import numpy as nm
from sfepy.base.testing import TestCommon
from sfepy.base.base import debug, pause

##
# 10.07.2007, c
class Test( TestCommon ):
    tests = ['test_volume']

    def from_conf( conf, options ):
        from sfepy.fem import ProblemDefinition

        problem = 
ProblemDefinition.from_conf(conf, init_equations=False)