Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Nonlinear fitting (work in progress) #151

Open
wants to merge 18 commits into
base: develop
Choose a base branch
from

Conversation

PrasadBabarendaGamage
Copy link
Contributor

Addresses #148

Distance vector sign was updated in previous commits by setting:
distanceVector(1:numberOfCoordinates)=interpolatedPoint%values(:,NO_PART_DERIV)-dataPointLocation
and by multiplying gradient and Hessian values by -1 compared with the original implementation. 

However, with the updated code, the gradient values are the same as what the old code produced, however the Hessian isn't. Fix was to revert to the original implementation for the gradient and Hessian and setting the distance vector as:
distanceVector(1:numberOfCoordinates)=dataPointLocation-interpolatedPoint%values(:,NO_PART_DERIV) 
This has been done for all cases of the projection (line, face, element).
…TTING_SUBTYPE

Based this on those used for Finite Elasticity.
…iteElementCalculate

(linear) Fitting_FiniteElementCalculate was originally used for EQUATIONS_SET_DATA_POINT_FITTING_SUBTYPE
Setting FiniteElementJacobianEvaluate to not implemented and instead will be calculated with FD
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

Successfully merging this pull request may close these issues.

2 participants