[Insight-developers] (no subject)
C. Aaron Cois
cacst11+ at pitt . edu
Thu, 20 Jun 2002 13:28:10 -0400
This is a multi-part message in MIME format.
------=_NextPart_000_0036_01C2185E.521D0060
Content-Type: multipart/alternative;
boundary="----=_NextPart_001_0037_01C2185E.521D0060"
------=_NextPart_001_0037_01C2185E.521D0060
Content-Type: text/plain;
charset="iso-8859-1"
Content-Transfer-Encoding: quoted-printable
We have run across an error that has stumped us for almost a week. It =
seems like a good time to see if any other developers have any thoughts =
or similar experiences. =20
We modified itkLevenbergMarquardtOptimizerTest.cxx to operate on our
particular cost function. The major modifications were:
1. GetValue() modified
2. GetDerivative() modified
3. Added ERF function
4. Changed Max_Iterations to 100
And a few other changes like prints and const. We used data generated =
with
known parameters from excel and the test will fit the parameters almost
perfectly:
----------------
vnl_levenberg_marquardt: 13 iterations, 66 evaluations, 5 residuals. =
RMS
error
start/end 60.749/0.000208907
End condition =3D Failed X Tolerance too small
Number of iters =3D 13
Number of evals =3D 66
Solution =3D (31.9992,255.022,4.0002,2.00008)
Test passed.
----------------
After this solutions is printed out, we get a crash as the smart pointer =
for the cost function goes out of scope...some kind of crash on cleanup. =
=20
I've included the source file (replace
itkLevenbergMarquardtOptimizerTest.cxx and hit 4 when running) and a =
screen
capture of the error.
Any ideas?
Thanks
------=_NextPart_001_0037_01C2185E.521D0060
Content-Type: text/html;
charset="iso-8859-1"
Content-Transfer-Encoding: quoted-printable
<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.0 Transitional//EN">
<HTML><HEAD>
<META content=3D"text/html; charset=3Diso-8859-1" =
http-equiv=3DContent-Type>
<META content=3D"MSHTML 5.00.2920.0" name=3DGENERATOR>
<STYLE></STYLE>
</HEAD>
<BODY bgColor=3D#ffffff>
<DIV><FONT face=3DArial size=3D2>We have run across an error that has =
stumped us for=20
almost a week. It seems like a good time to see if any other =
developers=20
have any thoughts or similar experiences. </FONT></DIV>
<DIV> </DIV>
<DIV><FONT face=3DArial size=3D2>We modified =
itkLevenbergMarquardtOptimizerTest.cxx=20
to operate on our<BR>particular cost function. The major =
modifications=20
were:<BR>1. GetValue() modified<BR>2. GetDerivative() modified<BR>3. =
Added ERF=20
function<BR>4. Changed Max_Iterations to 100<BR><BR>And a few other =
changes like=20
prints and const. We used data generated with<BR>known parameters from =
excel and=20
the test will fit the parameters=20
almost<BR>perfectly:<BR><BR>----------------<BR>vnl_levenberg_marquardt: =
13=20
iterations, 66 evaluations, 5 residuals. RMS<BR>error<BR>start/end =
60.749/0.000208907<BR>End condition =3D Failed X =
Tolerance too=20
small<BR>Number of iters =3D 13<BR>Number of evals =3D=20
66<BR><BR>Solution =3D=20
(31.9992,255.022,4.0002,2.00008)<BR>Test=20
passed.<BR>----------------<BR></FONT></DIV>
<DIV><FONT face=3DArial size=3D2>After this solutions is printed out, we =
get a crash=20
as the smart pointer for the cost function goes out of scope...some kind =
of=20
crash on cleanup. </FONT></DIV>
<DIV><FONT face=3DArial size=3D2><BR>I've included the source file=20
(replace<BR>itkLevenbergMarquardtOptimizerTest.cxx and hit 4 when =
running) and a=20
screen<BR>capture of the error.<BR><BR>Any=20
ideas?<BR><BR>Thanks<BR></DIV></FONT></BODY></HTML>
------=_NextPart_001_0037_01C2185E.521D0060--
------=_NextPart_000_0036_01C2185E.521D0060
Content-Type: application/octet-stream;
name="itkLevenbergMarquardtOptimizer.cxx"
Content-Transfer-Encoding: quoted-printable
Content-Disposition: attachment;
filename="itkLevenbergMarquardtOptimizer.cxx"
/*=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=
=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=
=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D
Program: Insight Segmentation & Registration Toolkit
Module: $RCSfile: itkLevenbergMarquardtOptimizer.cxx,v $
Language: C++
Date: $Date: 2002/03/01 18:31:57 $
Version: $Revision: 1.14 $
Copyright (c) 2002 Insight Consortium. All rights reserved.
See ITKCopyright.txt or http://www.itk.org/HTML/Copyright.htm for =
details.
This software is distributed WITHOUT ANY WARRANTY; without even=20
the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR =
PURPOSE. See the above copyright notices for more information.
=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=
=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=
=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D*/
#ifndef _itkLevenbergMarquardtOptimizer_txx
#define _itkLevenbergMarquardtOptimizer_txx
#include "itkLevenbergMarquardtOptimizer.h"
namespace itk
{
/**
* Constructor
*/
LevenbergMarquardtOptimizer
::LevenbergMarquardtOptimizer()
{
m_OptimizerInitialized =3D false;
m_VnlOptimizer =3D 0;
}
/**
* Destructor
*/
LevenbergMarquardtOptimizer
::~LevenbergMarquardtOptimizer()
{
delete m_VnlOptimizer;
}
/**
* Connect a Cost Function
*/
void
LevenbergMarquardtOptimizer
::SetCostFunction( MultipleValuedCostFunction * costFunction )
{
const unsigned int numberOfParameters =3D=20
costFunction->GetNumberOfParameters();
const unsigned int numberOfValues =3D=20
costFunction->GetNumberOfValues();
CostFunctionAdaptorType * adaptor =3D=20
new CostFunctionAdaptorType( numberOfParameters, numberOfValues =
);
=20
adaptor->SetCostFunction( costFunction );
if( m_OptimizerInitialized )
{=20
delete m_VnlOptimizer;
}
=20
this->SetCostFunctionAdaptor( adaptor );
m_VnlOptimizer =3D new vnl_levenberg_marquardt( *adaptor );
ScalesType scales( numberOfParameters );
scales.Fill( 1.0f );
SetScales( scales );
m_OptimizerInitialized =3D true;
}
/**
* Start the optimization
*/
void
LevenbergMarquardtOptimizer
::StartOptimization( void )
{
=20
ParametersType initialPosition =3D GetInitialPosition();
InternalParametersType parameters( initialPosition.Size() );
CostFunctionAdaptorType::ConvertExternalToInternalParameters(=20
GetInitialPosition(),=20
parameters );
m_VnlOptimizer->minimize( parameters );
// InternalParametersType is different than ParametersType....
ParametersType p(parameters.size());
for(unsigned int i=3D0; i < parameters.size(); ++i)
{
p[i] =3D parameters[i];
}
this->SetCurrentPosition( p );
=20
}
/**
* Get the Optimizer
*/
vnl_levenberg_marquardt *=20
LevenbergMarquardtOptimizer
::GetOptimizer()
{
return m_VnlOptimizer;
}
} // end namespace itk
#endif
------=_NextPart_000_0036_01C2185E.521D0060
Content-Type: image/gif;
name="error.gif"
Content-Transfer-Encoding: base64
Content-Disposition: attachment;
filename="error.gif"
R0lGODdhvAG/APcAAAAAAIAAAACAAICAAAAAgIAAgACAgMDAwMDcwKbK8EAgAGAgAIAgAKAgAMAg
AOAgAABAACBAAEBAAGBAAIBAAKBAAMBAAOBAAABgACBgAEBgAGBgAIBgAKBgAMBgAOBgAACAACCA
AECAAGCAAICAAKCAAMCAAOCAAACgACCgAECgAGCgAICgAKCgAMCgAOCgAADAACDAAEDAAGDAAIDA
AKDAAMDAAODAAADgACDgAEDgAGDgAIDgAKDgAMDgAODgAAAAQCAAQEAAQGAAQIAAQKAAQMAAQOAA
QAAgQCAgQEAgQGAgQIAgQKAgQMAgQOAgQABAQCBAQEBAQGBAQIBAQKBAQMBAQOBAQABgQCBgQEBg
QGBgQIBgQKBgQMBgQOBgQACAQCCAQECAQGCAQICAQKCAQMCAQOCAQACgQCCgQECgQGCgQICgQKCg
QMCgQOCgQADAQCDAQEDAQGDAQIDAQKDAQMDAQODAQADgQCDgQEDgQGDgQIDgQKDgQMDgQODgQAAA
gCAAgEAAgGAAgIAAgKAAgMAAgOAAgAAggCAggEAggGAggIAggKAggMAggOAggABAgCBAgEBAgGBA
gIBAgKBAgMBAgOBAgABggCBggEBggGBggIBggKBggMBggOBggACAgCCAgECAgGCAgICAgKCAgMCA
gOCAgACggCCggECggGCggICggKCggMCggOCggADAgCDAgEDAgGDAgIDAgKDAgMDAgODAgADggCDg
gEDggGDggIDggKDggMDggODggAAAwCAAwEAAwGAAwIAAwKAAwMAAwOAAwAAgwCAgwEAgwGAgwIAg
wKAgwMAgwOAgwABAwCBAwEBAwGBAwIBAwKBAwMBAwOBAwABgwCBgwEBgwGBgwIBgwKBgwMBgwOBg
wACAwCCAwECAwGCAwICAwKCAwMCAwOCAwACgwCCgwECgwGCgwICgwKCgwMCgwOCgwADAwCDAwEDA
wGDAwIDAwKDAwP/78KCgpICAgP8AAAD/AP//AAAA//8A/wD//////yH5BAAAAAAALAAAAAC8Ab8A
AAj/ABEIHEiwoMGDCBMqXMiwocOHECNKnEixosWLGDNq3MixI0IpCP6JHEmypMmTKFOqXMmypcuX
MGPKnEmzps2bOHPq3Mmzp0lSIP95HEq0qNGjSJMqXcq0qUOgIREgmkq1qtVEWLNq3Rqpq9evYMOK
HUu2rFmwktKqXauprdu3cOPKnUu3blttePPq3bt3m99t4wILHky4sOHDiBMrXhx43Tp2kCNLnky5
suXLmDNrZkevs2d69UKLFp2gtOnTqFOrXs26tWvTCKAKtUq76tbbWs/q3s27d6S1wCXZHU68uFy+
yJH/Bcy4ufPn0Ac73ky9uvXrkj97Hs39tffv4MHH/w4qtXZt3Oh9q1/vOzhb4/Djz01OP+/y6Pjz
62/8GLv//wBq1xl3o4Vn4IEvSVHaeFGZdx56t7En4YRjuaeWfBhmWF999+3n4YeJTQfgiCRqJiBo
BIZ24IrfCcUQUAvKVp6DtkHIVW8jUcibSGZZmFaGcYkE5HzaiLRhX3+BqOSS/JXo5JPZCZiiiixW
yZqLC8GYAIOzISKSVV9+aeONu/EYiZlfmckjmuqxeSZJv/loV0l1CQmknW0ZaeSR9iXJ5J/7iQjl
oCOeOGU9ViaaGpakkEJQowJpyWV5I1FVKVVj5rbjP125+WannFLIpptr/QOcW//IZSeeQaZ6p6tu
4f+1J5/adAjordAJSuiu1xk6paLAloZlbI4SK9J4MZI3FUleViomSYmIFC1WcIKqZo5vclotqNxe
u2ao3oIbKrdflXpqnqrCuupIeZLU7rua0IkqrPH+UyRJstpr771GbiOSv+OIFJjAuBZcmK68Jmyi
lL8Gq+iwxDZ6LLJbyrhss2H+02yzWkErbbTajvvpp9eKVe23aYoLFspekaqWqe/Fm+689b67bs01
30zzvPfmy+++ef1LUsD/DFy0wUg3qfDSmfmaosMPGyQxowpWrKyXy2q8ccYd/5PVx9+WTLLKYZ0c
stjhhpwyWielddJbeAppEs5yu1o3vHjX23PPs/7//I+/fws2dNJJI8z04ZE5TSDUiULcKKQDSWox
1mBqzPXXXlObebgjs+xptmTDWfLoIrvs9rkyt0ozq3SnerPOeeu5r+xBz/73v4KPRHjh/SHuO2SK
d8d4lcNODTnFk15c+daZa4655muTTm7Z446q9tjdli5ynKfHTO/OsMN7d85208sq7XwDvafQgee+
O++/xx98gcOzyCgpE0Mu+dVaV5UxIiCT1sc0xzlrkS1b1GvZ2apnthy5qSTckwTMLoSuVrFrXuyS
V70uqMEN1s1vPqvdSADHHKO9D2mGi5/C5kea+q1oWBMTyLH21yAa1ShTidCRDifkI+EM6VV2oRVf
52xFsBMWLIUq5BULqeRCA0EMITTskg2ngkOs7PCK6unhD+XjriAKUS9EPJoRcYXEJBJqiYhqohMf
x8Y2NqpqyZsiFauIxTruRotbzONxvtgnv4xxd2U0I5TQqEYDSeGQiEykIuE4OTlWMYd2jCRZ8KjH
St6Fj3ix1R8NFkhBOomQhWxiHB1JR0maEi1ysqQlMZlJP22Sk73zpBIZ9rRQunBSTsmlLnfJy176
8pfAzIiMgknMYhrzmMhMpjIjR55lOvOZ0IymNKdZkGFS85rYzKY2t7kRa3Lzm+AMpziv6c1xmvOc
6ExnU//Kqc52uvOd8HwIO+NJz3raE53zvKc+98nPaOaznwANqEB5+c+BGvSgCO1IQRPK0IY69EXN
rEgAJkrRiT70ohh950ITQtF8ePSjHqVoRkdK0m9u9CABAKlKVRqApgCgpDCNaUVOSpCJrvSmH7Vo
RgDA05cixKcS6SlPKQJUmRo1nTQVSEpxytSQ7nQgRSVIVB8y1aAe9ar4jChDltrUprb0IkX1aU8F
IlSyQpWsYz2rVBHw0raiVaxsNStW5xrOpHK1qzgNQLEqEta4mhWobv2rWt8q1sLKNa5DPSxdF5tN
mt7Vo/9oamRZuteJ9FWogUXsWQG71s769bOZ/SxjRzvczZM+Nh8jwWlqWWqRyxaEs7AdrGhBO9jQ
VpW0uE2maVVakpX2FqR65ataYyvY4oo2qpw9bGhzy9xn7pa3JAHpb4F7vIiUFaqJvS5b06rd7Y41
uW9V7G2bS15fbvS0kI0uatW70uD6crzlja9JtYoQ9KYXJUx1ry7TKt/+jvO8kj2JV6vr3wJfFcAB
Zm9eCWzgBsMUwUw1SVf1WlkHW3ikEL4pfvPL4At7uKEZhu6Gb0rhD5v4oSH+6G+nS9kKn/jFAn2u
dBXM4pA+DsY4HqiM7ysS3yo4HxR2cY6HXP/PHd9XtT1uMZGXrE/H+lay7b0xk6ccT7vidcBSprKW
1WnlKy+4w1sOMziTigD74jXIYk6zOclcZjOTOMhCVrOcyUnfhpDCzdRl45z3vE02R6qjb4ZznPkJ
X4ggd7ZFKXRCqqpoozRatg0Z76OLXGd53rmiFHWjRroLTUa/1iCc3uynARtqUEP6p6/lL6pNjdZT
W3YhpE7sQSatXNpi19VWVQhmGUJri/j5UW4M9kbAq81Df1qxxxa1Z2Ft3WUzm9XbxXVHwHvbXhv7
uMQ+yqN7PdNKSyTYgxauYAsba7eSO7yytjZ/v0tY7567trQN7HX7it1Yexex6+YussPb7nv/y1Xe
8R5qurOrWX/r+7j8nvdwD95WgY/buA1veL0fbvB0S3XgEef3xeH6b1kf5Nf73SzHC45vU3N85A4h
LsTjDeqDazazKKc2wleeXJTXWuUn9yvMdU5zuNJ7uGvFOcv/TfPi7jznsPW5s21NcqYnW+lOr6a3
j6nwmyu73Tl/9V+FrnJ4x3bnM3e3xK2+cp4/PeAMp3jBfV7Wrns85mIfOsmRfmuiH93ubb96tD27
a4u3uuXrrrfHpR6VTp+96GS3OVVra+6rL1fuja8532d997IrnuusVjpxoe71lhNd75C3NecB/nJl
c57sF4e355F77clDtPDPvDbmE/9515sd1+GVbzrqIw901ePe8aKH9OzNTuqhj17unR258kMv+blD
vPLH77rvif34rPN8+VCcujFbn2+HC574+5Y52jvufZ2HtfuNd7fvK158sGu33Bqf99tBS27Sv7vm
Ap9/d7mL/n6vHebFp35Z53Lnpmrn927RFngJt3DwBXLRxG1JAYFaRxTZBla8JIG6Jm0I5YDLpGq/
5IGONhQGyBEY+FSbNn8ZxYF8toIapX0s+IL7pIIwOIPiJIM0eIN95oI4uINZBXs8+IPuZINAOITK
JIREeITFZIRIuITm/6WDTPiEy6SEUDiFSyGFVHiFRmGFWLiFCuWEXPiFTqGFYDiGFCGG8TR4ISd8
vLZow0aGWeKFH+iBA1d3e/d9uxZpGuh5+zaBs6Vo1TZqy2aAFZhsG7doKBh/g1d/EZddaMhYZugR
1BeAdXiIe/hsiKZrg5hryOaHq+Z44td6swZtntZxnadcigd8ufWI05Z8drd2tXdqU8V/GWdwv2d+
mpd/E4dtBOdv9wZ3gud3gAhtCXeIuSd+zueKoRh1jgiH78WK14eMl6eAoWiMY2d6Afd5h/Zz6dd0
y5WNhLh3+4dqsQh02Md00EdY9JaJc6WKJOiMA4iN6BZp7GaK88iKxeGIjfXIcm43jkH3jWXXh/lY
awJpjsDXd8mojIvFjm0ocnl3fQYJegcJdn2Yec7XfLmnj3rXjb/4c/04kMg3kNRYijM3ighJVwq5
ab1nXCq5fPwokjLXkuXnfohndW4nXh1pexKpjYgmeeX4krJ1ea5Xggx1kjsljbpnW3ZYdRtnW/cn
ik0JjwK4lOQ3cdQHet2XizvJf6mGfVeJiMbWkvFHWkQ5Z0IpjG5IeE90lmWZh2Q4lmf5lszkg3A5
l0fhlnTphnZ5l2OYl3r5hXzZl1v4l4B5hYI5mFNYmIb5hIiZmEv/uJiMeYSO+ZhDGJmS+YOUWZk7
eJmYeYOauZkz2Jme+YKgGZorOJqkuWemeZpylpqqmWas2Zph9pqwqWWyOZtTVpu2uWS4mZtDtpu8
iWO++ZsvtptrSUwkKW6QSIF4yYzNqH4YUWqV2Il4aJaLF50WyInWuXRYiJvliJx/WJ0pd5DNlp3I
aYng6ZfMeYGdh4CyiJOMyHvi9Z6Ad38E6H+92JWF6JwKR58FOIyUGJjpuUsvyZLa2Xxz153GyHwF
+ZMuKY4L+ngzmZN7GaBpqJ+/OImxiH4SeZTrl3T1eI/TF3dSiXnd+KGhR56NSaG5lG2DOIcpKXs3
CZKTB6NCF6LfSFmRooaUEUqQYFibLNqhgfiiUQeTDLqjNSqjVWmk3Lig/1icoqmiK7qJ+dmeQcqL
5beH8rmV/Al/znlr+Jeh+feeUCeLhvWUt/+HnnLZgcXkpOVpgmy5kD0KpVEaTGxqXY1ogdopgsuZ
psIJmMHZpx/2p4B6YYI6qA5WqIZqYIiaqP61qIwqX476qOUVqZLaXJRaqakop5g6mZq6qZbZqZ6a
maAaqpw5qqT6maZ6qk/Kp6pKmKn6TZgVq7I6q7Raq7Z6q7JqED6xqy0BEgXBq8CaEr6qY6/KTQAQ
rMg6ElIwVWlZWsM6EM0qTVBhUJe6adHqXMuqq9w0rQRxrdDErTFWrNt0rNoEFMy6rc8qQ+h6UNW6
U1jCe/BaevE6r/Jar4Flrtqag/larukKUO36nC5Sp0jxUvj6q+tqsPpKreJabN6qiUlRsN0kerAR
m7DEyqqEFrB8pY7UaWjjca6tBYJv6mv92rA/dafr1K/9//SvYEWysKax/vht2Yqw3imdQwGu6kpU
wGSz/rqw2USulzid9ZeL5mayzAax0Iqn9NefGzoROsuyEfldXSlvRPt6Cmux++Sz4kagW4eib+ix
Mwt53cm0I5u1TFpyJWlnKMtPKttaGKt1WnuMOGu0N5uxGBeTMVqGY4uzCvp8IPsUaRuDPItNWGtZ
D4mjXKsQcmu1QKuGNlmWTUu2Snq4iPu3TRa41+SzEOhaOypcieu0nfi2oIi3+5prwyeholu1ngtP
g+uwiyi0qddtXqu3U9lqoSu2o2toYip6u6hxtou6CIW500SwMTuxc4oRjwunTKGzKWu51LS6dnqH
VTi8R8/7sdC7hsKUt0EFtQLrt+zKvMHbtg/YsbdLZ+NLTcqrtt4rTVhLr+xrr+7bvunXuRI7vfza
vYp7T84rrdI7t42Fvf1rv6n7TscKbgRcwAZ8wAgMbvt7v850vPXruwcFAIs0wRRcwRZ8wRhswcya
wBzcwR5swCP7wSI8wh1Mufe0tnSLqyq8wiyswgaRwTAcwzI8wy88wzZ8wzgMwK06oQy8w5CZvj7M
ZCgcxCkIxETcm0Z8xMCZxEo8nEzcxCY2xFA8lE88xYRaxVZ8qFicxYo1ChIk/MVgHMZiPMZkXMZm
fMZonMZqvMZs/MU4/MZwHMdyPMd0XMd2fMd4nMd6vMd83MeKFBAAOw==
------=_NextPart_000_0036_01C2185E.521D0060--