Hi all,
I am working on a 3D model that requires a very large mesh. When I tried a
uniform grid of 80*80*80 cells, I got the following attached memory error. I
wonder if there are ways to overcome this limit. Thank you in advance for
any advice.
Regards,
Zhiwen
python(12494) malloc: *** mmap(size=73728000) failed (error code=12)
*** error: can't allocate region
*** set a breakpoint in malloc_error_break to debug
Traceback (most recent call last):
File "run_pyramid.py", line 55, in <module>
eq_u,eq_v,eq_w=getEquations(u=u,v=v,w=w,c11=c11,c12=c12,c13=c13,c33=c33,c44=
c44,beta11=beta11,beta22=beta22,beta33=beta33,cinn=cinn)
File "/Volumes/RAID2/Users/liangz/DEV/Elastic/Solve/elasticEquations.py",
line
41, in getEquations
graduvw_w=Gradc(c13,Zhat)*Gradc(u,Xhat)+Gradc(c44,Xhat)*Gradc(u,Zhat)+Gradc(
c13,Zhat)*Gradc(v,Yhat)+Gradc(c44,Yhat)*Gradc(v,Zhat)
File "/Volumes/RAID2/Users/liangz/DEV/Elastic/Solve/myGrads.py", line 15,
in G
radc
return _Gradc(var).dot(hat)
File
"/Volumes/RAID2/usr/local/stow/fipy-trunk/lib/python2.5/site-packages/fip
y/variables/meshVariable.py", line 193, in dot
return _MeshVariable.__dot(self, other,
self._OperatorVariableClass(baseClas
s))
File
"/Volumes/RAID2/usr/local/stow/fipy-trunk/lib/python2.5/site-packages/fip
y/variables/meshVariable.py", line 165, in __dot
opShape = numerix._broadcastShape(A[index].shape, B.shape)
File
"/Volumes/RAID2/usr/local/stow/fipy-trunk/lib/python2.5/site-packages/fip
y/variables/variable.py", line 1287, in __getitem__
unit=self.getUnit(),
File
"/Volumes/RAID2/usr/local/stow/fipy-trunk/lib/python2.5/site-packages/fip
y/variables/variable.py", line 249, in getUnit
return self._extractUnit(self.getValue())
File
"/Volumes/RAID2/usr/local/stow/fipy-trunk/lib/python2.5/site-packages/fip
y/variables/variable.py", line 453, in getValue
value = self._calcValue()
File
"/Volumes/RAID2/usr/local/stow/fipy-trunk/lib/python2.5/site-packages/fip
y/variables/gaussCellGradVariable.py", line 93, in _calcValue
return inline._optionalInline(self._calcValueIn, self._calcValuePy, N,
M, id
s, orientations, volumes)
File
"/Volumes/RAID2/usr/local/stow/fipy-trunk/lib/python2.5/site-packages/fip
y/tools/inline/inline.py", line 9, in _optionalInline
return pythonFn(*args)
File
"/Volumes/RAID2/usr/local/stow/fipy-trunk/lib/python2.5/site-packages/fip
y/variables/gaussCellGradVariable.py", line 80, in _calcValuePy
grad = numerix.array(numerix.sum(orientations * contributions, 1))
File "/sw/lib/python2.5/site-packages/numpy/ma/core.py", line 1712, in
__mul__
return multiply(self, other)
File "/sw/lib/python2.5/site-packages/numpy/ma/core.py", line 514, in
__call__
(d1, d2) = (get_data(a), get_data(b))
File "/sw/lib/python2.5/site-packages/numpy/ma/core.py", line 317, in
get_data
data = getattr(a, '_data', np.array(a, subok=subok))
MemoryError
