[Insight-users] Re: help on how to add new element to itk

[Luis Ibanez]

This is a multi-part message in MIME format.
--------------060401060303080006030905
Content-Type: text/plain; charset=us-ascii; format=flowed
Content-Transfer-Encoding: 7bit


Hi Luzhong,

Thanks for sending the other files.

Your code works for me, compiles fine in
VisualStudio.Net and when executed it
generates the following output:


Mass matrix:
0.388889 0 0.138889 0 0.111111 0 0.277778 0
0 0.388889 0 0.138889 0 0.111111 0 0.277778
0.138889 0 0.166667 0 0.166667 0 0.111111 0
0 0.138889 0 0.166667 0 0.166667 0 0.111111
0.111111 0 0.166667 0 0.5 0 0.305556 0
0 0.111111 0 0.166667 0 0.5 0 0.305556
0.277778 0 0.111111 0 0.305556 0 0.722222 0
0 0.277778 0 0.111111 0 0.305556 0 0.722222

Material matrix:
30000 0 0 0
0 30000 0 0
0 0 30000 0
0 0 0 30000

#dof per node = 2
Test PASSED!

---------------------------------------------

The only line I changed was on the "pragma"
that was misspelled as "progam" in your file.

I built this against a CVS checkout of ITK
(updated today) so the vnl methods "resize"
were changed to "set_size()".


Please find attached the files.



Regards,



    Luis



-----------------------
luzhong yin wrote:

> Hi, Luis
> 
> Thank you in advance.
> 
> Luzhong
> 
> The rest code is
> 
......
 >


--------------060401060303080006030905
Content-Type: application/x-zip-compressed;
 name="LuzhongYin.zip"
Content-Transfer-Encoding: base64
Content-Disposition: inline;
 filename="LuzhongYin.zip"

UEsDBBQAAAAIAAV+KTC2h8Ow9QAAAHYBAAAZAAAATHV6aG9uZ1lpbi9DTWFrZUxpc3RzLnR4
dHVPQWrDMBA81+A/iJCDezGlP1CkVVAtS64kQ26ipAqExjbUChhC/t6Vc2qhB61WszOzo86a
N2C+EtDGJT2XRVkIqXnoKGvoHirpGwSlyE0Qptccn09SM9VzqLa3DPcOgpAK7jgC5eA3twXn
spOgnqoA1hpLNsgg45TIabqOnzXpLvFjjmSOiWQtl7beZDPN/2zGKI/VmYO5jZXgqrIgZHt7
fCE401u2zu8ZX0WU8wAHYL2nOwXVKQ4ET31cFnJOXyiESxzimF45e1Hv/+HtOkBDT+0efFBS
N1h2lq4psitmRVW+pMmVTcMwjbnT1yF+n4/zmucHUEsDBBQAAAAIAAJdKDBq4TnFTAIAAN4F
AAASAAAATHV6aG9uZ1lpbi9mZW0uY3h4jVRdj5pAFH0n4T9MNdmgcd0B+oRiYoVtTBa1i9s0
2W4IhcvupMxghjGtafrfO3yI4qopDwznnnvP/RhmuoRF6TYG1CHi573ruSlQYMJwZvjhizd8
66hK97zLpzCHkleVu7tTn5pozGOS5YJDSCeFeZsT9opYSCHfhBGgXMSj91apdMaaAB1VOVWF
MIFoSJgmPwYoegs56j+/9FTlj6og+SyyGCxrlUkaOGJ4wPQBMwbMHFV83YllfYVIZHy92wDa
CM3olRkKD4btSmQBv7ReHbYRz/jFxsMD1I8gw7cTH8Qsy3hMWCgg1zbiSFG/qKi3FRvI9OuK
xkVFs61oNIrGdUXzf7tuEjDzqqInLZyE6QNhEHI3DXNBIiJ2h92htRC1L/sel0NvJ58XNj4g
1zaxfIZHpqksFxsHPC8w/thU1f7ZD7UARraqnLL77FUs4FYBUO16MTYND9BNn+HeGUovKf0c
ZZSUcY4yS8qsKVUpOBrQUNjxTp4MEgWRHNL40uD6E+2mT3unbVuWDODkd/nbOwPkQbs5EF6Y
55WP5sFxYSVXJat5Z09H2Vag8Rh1imB5OgvWQp3CBCxOi9WDPXofU4mej3MuhXXjLEGb4ojL
aSG7jKqmJAtdbOkP4MvEgVcOkC+Te7nEGV0BL4fbO1ItbqyExZAgeYcFT74b+N70cR2slvPF
2n30q7wxpCBADmLUwrQN2QnN9BNsnODiVurKQkhS2au3vBstq+l0DblAq6nvu86H76xTK3AQ
W86Q+22+Dvyn2cz1fcn8VZV/UEsDBBQAAAAIAPd9KTD3NlvvWwYAAA4WAAAiAAAATHV6aG9u
Z1lpbi9pdGtGRU1FbGVtZW50MkRDMExRLmN4eMVYbW/aVhT+juT/cNRqkw3mxcbdJiiRMghV
qqTpljabVlXI4AtYNdfUNom9Nv9959xrg20MY8mHJQLse87zPOec+2q/dOcOm8Pk+nY4ubv4
Xam9XAf2YmXDgx1wly9Uxw3tqcd6YP38y08a2hl33LlSwyuXz7yNw+CFG30ZX1xfeGzFeGSO
hp2r31rLF3mPe+618TNZ2dFSmJQat1csXNszBgiHb/mGOVtRQ/Z/77uOUivQK7Ve7w2LLnnE
FoEduT5/77s8OufOH8xdLCN1w0N3wZkD2AqurtTu2Czygw/Jmv0I60iHsefb0Y/woEPB1Q8c
Fmgw83kYKTUMAqDdhvHln9cXPQhsvmAwW7LZFywNhSasWEB740WEz6KRPGR356CKm8Ggo8E3
aRiMJCQX/w219+FRsq6jVsiiSej+zVRT62eNnzqfBwt7E4Yi20+C6/Mn9wd50U+djEqn9tYJ
3ZTag/SR5Soz1atseYLHXe/ky1fdS+82qykLbubl3grVitIrtf+79gGLNkHKUhff/ULChRR7
vd3Aqsj+dmmv2XjDZyQTqnJcQXEsQmm41eHK5cwO4OvGdgLXsyMW2B4wyQ1LO4S5vwkgJG6l
Ns/YAeptItixS5examUjqF2nbzhBwfNntkd9gQVwOTqEKVJtGnrT0HRQd7/0Y3NH2AxNOopQ
0nwW7j3jkhLW1PMQD9RADxH2sGQBg6YBrweALfRjKDUiI4YUL9KALFEweqCqBjQBJyoGRFeh
pkEbLACV+w4WxdDSCGQFcN7AQHqKSbTFidlCdy3zVf+Qnin1Gkf0zJKekeo1nqLXLes19vS6
JT2zUq9xmt5ePff1rJJet7KelXrpdJLA4kw6sLAXJs2IBe49ztJ7dmD64Np+bUeBG8sGoTPa
m1N12M0TOwLHDVh6F4ihG7Ygn+Eot/zquelTh1088NAKWhHi536wV1JZOKDZkKfFWsmx2Dw8
HCo0wv+kYZQ1OqdoVOZhphpmRR5ijDeen8dBDSPVaBYH9lPy6KYa3Yo8zCyPA9Pl1DwOahhl
jafnYaUaVkUe3UKtnp7HQQ0j1WhUj6vtxJ76vld9EriiLWAc+Ks3nj+1veFud6mY20ptIbze
R6AX5rzYSN5H+SlOUYoTHcSGDrGJny5+LB0SvE/wPsH7BO/jGcML+oqnuHQkU2ri1MR1sXOJ
v3Md3iLwLRF10IgfJ8ZPosMUoVMEEMWMp9VNQyod2nbtc9fz1E6rk60ksYF1jJZu2DxbTaja
WMvmGVYoXxING/sYTILOSu3fvQ15AozNMrdxhNvc4670zri7ZW7zCHd3j7vSO+O2ytzdI9zW
Hnelt/E5q/gU2WMauFifBiWCV5YQTsiUkCkhU0KmxMpw+DwigA0CNotAYUvIlpCtjOQ5SYFs
7JA8pymQjR3yHNDYab2iXVVVpaIh9lhMmxC4OzfRFFuCGy1KTZXihqbJkfeWRlgGtiTYKICN
DAwCi3PByrCmxMqoc1irgE1ZTYHtatux3RHRmyL8rHCNLH2pkOR9sgI2sjJkTE6MXtkiQHsa
inYE3knyFoMsSSdFTalTsIYYKBLUsXs1ZFURUsdRIPWnVP7mcR/Rual5xkTi0o430oHnHHje
gWcZpPM/PX+aaCSmNp6Uw69BpKoYax0DFliyYq9hX6IPtUxZYRFJD5WChAsSSEk4kXCRgSAx
BQmXJNtYaFkGN7zEZy6HDaJgw1IDPSblAn2Nx/FWB75/z0d/BrJNqeWi2fc0Mk/xFPBNLqhb
zbnthUykVHzUyhz6hdOheFjBTSuwH3BZYNvHE9yiHHbvYvfg+h+xOIL1aNgi3za9chDvM8YX
15NfP15ejSZ3l7cfz68u/zr/cHnz7uCRk0TU4WhYJy4dbn1vQ3thrzekLUY8rTLcKH2vvOmI
hxljcHwRr4+Gk1sskMydIMlxiFGCbJXMwfElvULpOOSgUndwfIGvUDoOOahkDY4v9xVKxyFF
JYXedsRGY4C9J5yzzr2zvQ1Tq7bVEVsEjN3Mx/jt+Cu1o2mFIHC5exafsccXmyfxGSfH9yy+
ivi6J/GZJ8f3LL6K+KyT+Lonx/csvlJ8NGxxWWmeXfv37IOPe6KO+3R/20wvY6jZ1BOzormr
4+a632zpiVXRnHGLpXT3slb+Pz7SOzNshMLb115vzlZk/wdQSwMEFAAAAAgAoX0pMF1WUDR2
AgAAygUAACAAAABMdXpob25nWWluL2l0a0ZFTUVsZW1lbnQyREMwTFEuaL1U32vbMBB+N/h/
OBooafCaEvqUboPUbkogbbI57WAMgmKdEzFZMrKctpT87z3ZyxKHlL6M2S+29N3dd9/9aIlU
cUxhPhf29/Dm7kZihsr2ovBi/G2+8r0W3QqF7wMIIlQiS45w0oDElp+vThxAsQyLnCUIBIDX
/YMUM3ewfbudju9BB34lkhUFNILVFwsjiO7lJ6U58gAkcWMmgPACEq2sUKUuCyDGwqLvYW0P
QkEvgiriuXPT9b0jAaAPebmQIoFdBp8vg95X3yOK9iVHp9TBHcwwyyWzyKfM0Hno/F75Hukw
H1zHs++DcDYPx4M4bjejBb53zBTOfK9m0d+T5Z89lD10ncZAYt6hXWlegMGKBVgNqszQkAJC
WVwaZoVWNbjreAAgAeAVIkxZKe1oh5oYjga+QA82VzV0LYwtmYS1FhxucR881eR/oPgPFMuV
bZeqEEuF3PeEK1YAj5hYbWak+CnkNoCh1MyewlMAW6gjCNrFJL2o8oU9iLoP9D0Kf19mCzST
9JBG0f7A5/8tgF0hLFFnaM0L6BSYgj9d3KjDNs2dUhCvWI7DUiUutaINVQb7UlJjWTiDY2pV
NWo4iKgN1qTSGreufO+wLHfMGvFc/xfOOHrH/UJr6VpgrBMmh0Znt1IvmAw1qS0UpX48xrJC
TW2zIaRzMrXNSPXicJpGhj1VKm6HX6vqt8gxEakglTmuBa0ety/weSdrS6Ruvt3gXj+MxtH8
cRQ/DMajn4PZaHJfpeNUcv7bYRR2II/CAGItS6dXvx86OlVL0SAUWu41UQsVF6ljuvnbUvS9
gW6XJopDY0H2+7QUq71aWYEDtT5a0773BlBLAwQUAAAACADBfSkwjm5Y1+QBAADVAwAAIwAA
AEx1emhvbmdZaW4vaXRrRkVNRWxlbWVudDJEQzBMUU0uY3h4jVNda9swFH03+D9cUhh28Jo0
LdvwlsJI3RFIShuX7VGo8nVyV1s2stosK/nvu3KSddn8MMsg+Zz76Xt0QnmGOYh5OhFfk4Xv
ndRGLksJa2k06WWQUSMfCozh4v2HdyHzqDPKfY9PpFXxlCH0yD5eJ/OkwBK1HV1NhrO7+emq
54y0LLGppUJgI3j5E8ixdMBhHbv7XhwfI0EIMaRPNRpVyKYJuBb23v6Xq+9B+9xUGU6v7jc1
gj4TURc86obPu+GLV3guLRqSRRxPKt3Y24o0A1CKzroBBgNI0YJdISyxKtEaUrKA2vk1zsCu
qHl7yTYuYQDDyNUM4ccO7ixyhXdzo8hV382dR66FlmtL6rfN9GGqyXIr9BPb+up9L7aCct8l
VA/fUVmQOvM9tUL1yJbS7v3XyK/hxugZ9d8h3Keh5cpC+0NODznzlvmdoCVZhg2wzipjdtkA
fyisLVXa96hhD1Ot9T7EwG2UBxCUgsOMsw2LjZRQsrGflBuK7x2mNCON0iScw5Iiu+lfBm/6
PKsQxmMYQuhCvexG2+bwPafwQ+5vptLLSTtOIa6ns0SISIjZ9MYdesfi+1fHvd0stnv17haH
F5PZ5zQVi+TLNL1PFsGxX/hqu9069fBFhKPrFcd8pRz/C1BLAwQUAAAACAD2XCgwyvRxEJwB
AAC7AwAAIQAAAEx1emhvbmdZaW4vaXRrRkVNRWxlbWVudDJEQzBMUU0uaIVTTUvDQBC9L+x/
GOzFlmq09hRFKE3FQlqV9igsm2RiFzebkGwPQfzvTrL1Ixo1l8Cb92bmvd0dqNQkmIIQyj7f
LFYLjRkaOwnmZ+HDSuw4G1BZGfyDQRxlYr1PEI56OKe7oz8YK8yiUhpsWZwZmWFVyBiBePDy
FUgxawDOvNGIMxjBY6xlVUF3IVeJSkWmpicmTzAB2l9ZBHREcDJlYBLAoXVecqbJpSyJJSur
YmVrKMo8Ik3T0uOsbxr4UOwjrWLOfji66lCvOaPlybmYh7PN5rjbZ/yfesiZG+Q3CQDcztZB
uBCLkKJcb0V4Nws2x8NLV3QBAQURYCr3miznprLlPrZ5CbnRNWVAXiuwOwraWCyN1JxVVJZP
6LRe8+tu2TNg/t4YpPkIOKopV4xVWivzBEXeTqBh+UE0dVhFmgQkZJLKSurT3we3FfrWdJ7L
YFsXCOZcjKEPn4hxH3zRD08/4dVhD99vXd27vaEQ4Hy/XoLnQf+l+3K0Yrlebr+d77DVvzZ6
JM+dW+77aXPH6I1QSaXQkAb/PkvO3gBQSwECFAsUAAAACAAFfikwtofDsPUAAAB2AQAAGQAA
AAAAAAABACAAAAAAAAAATHV6aG9uZ1lpbi9DTWFrZUxpc3RzLnR4dFBLAQIUCxQAAAAIAAJd
KDBq4TnFTAIAAN4FAAASAAAAAAAAAAEAIAAAACwBAABMdXpob25nWWluL2ZlbS5jeHhQSwEC
FAsUAAAACAD3fSkw9zZb71sGAAAOFgAAIgAAAAAAAAABACAAAACoAwAATHV6aG9uZ1lpbi9p
dGtGRU1FbGVtZW50MkRDMExRLmN4eFBLAQIUCxQAAAAIAKF9KTBdVlA0dgIAAMoFAAAgAAAA
AAAAAAEAIAAAAEMKAABMdXpob25nWWluL2l0a0ZFTUVsZW1lbnQyREMwTFEuaFBLAQIUCxQA
AAAIAMF9KTCObljX5AEAANUDAAAjAAAAAAAAAAEAIAAAAPcMAABMdXpob25nWWluL2l0a0ZF
TUVsZW1lbnQyREMwTFFNLmN4eFBLAQIUCxQAAAAIAPZcKDDK9HEQnAEAALsDAAAhAAAAAAAA
AAEAIAAAABwPAABMdXpob25nWWluL2l0a0ZFTUVsZW1lbnQyREMwTFFNLmhQSwUGAAAAAAYA
BgDFAQAA9xAAAAAA
--------------060401060303080006030905--