Index: sm/trunk-jpl/src/c/analyses/GiaAnalysis.cpp =================================================================== --- /issm/trunk-jpl/src/c/analyses/GiaAnalysis.cpp (revision 21531) +++ (revision ) @@ -1,69 +1,0 @@ -#include "./GiaAnalysis.h" -#include "../toolkits/toolkits.h" -#include "../classes/classes.h" -#include "../shared/shared.h" -#include "../modules/modules.h" - -/*Model processing*/ -void GiaAnalysis::CreateConstraints(Constraints* constraints,IoModel* iomodel){/*{{{*/ - /*No constraints*/ -}/*}}}*/ -void GiaAnalysis::CreateLoads(Loads* loads, IoModel* iomodel){/*{{{*/ - /*No loads*/ -}/*}}}*/ -void GiaAnalysis::CreateNodes(Nodes* nodes,IoModel* iomodel){/*{{{*/ - ::CreateNodes(nodes,iomodel,GiaAnalysisEnum,P1Enum); -}/*}}}*/ -int GiaAnalysis::DofsPerNode(int** doflist,int domaintype,int approximation){/*{{{*/ - return 1; -}/*}}}*/ -void GiaAnalysis::UpdateElements(Elements* elements,IoModel* iomodel,int analysis_counter,int analysis_type){/*{{{*/ - - /*Update elements: */ - int counter=0; - for(int i=0;inumberofelements;i++){ - if(iomodel->my_elements[i]){ - Element* element=(Element*)elements->GetObjectByOffset(counter); - element->Update(i,iomodel,analysis_counter,analysis_type,P1Enum); - counter++; - } - } - - iomodel->FetchDataToInput(elements,"md.geometry.thickness",ThicknessEnum); - iomodel->FetchDataToInput(elements,"md.gia.mantle_viscosity",GiaMantleViscosityEnum); - iomodel->FetchDataToInput(elements,"md.gia.lithosphere_thickness",GiaLithosphereThicknessEnum); - iomodel->FetchDataToInput(elements,"md.mask.ice_levelset",MaskIceLevelsetEnum); -}/*}}}*/ -void GiaAnalysis::UpdateParameters(Parameters* parameters,IoModel* iomodel,int solution_enum,int analysis_enum){/*{{{*/ -}/*}}}*/ - -/*Finite Element Analysis*/ -void GiaAnalysis::Core(FemModel* femmodel){/*{{{*/ - _error_("not implemented"); -}/*}}}*/ -ElementVector* GiaAnalysis::CreateDVector(Element* element){/*{{{*/ - /*Default, return NULL*/ - return NULL; -}/*}}}*/ -ElementMatrix* GiaAnalysis::CreateJacobianMatrix(Element* element){/*{{{*/ -_error_("Not implemented"); -}/*}}}*/ -ElementMatrix* GiaAnalysis::CreateKMatrix(Element* element){/*{{{*/ - _error_("not implemented yet"); -}/*}}}*/ -ElementVector* GiaAnalysis::CreatePVector(Element* element){/*{{{*/ -_error_("not implemented yet"); -}/*}}}*/ -void GiaAnalysis::GetSolutionFromInputs(Vector* solution,Element* element){/*{{{*/ - _error_("not implemented yet"); -}/*}}}*/ -void GiaAnalysis::GradientJ(Vector* gradient,Element* element,int control_type,int control_index){/*{{{*/ - _error_("Not implemented yet"); -}/*}}}*/ -void GiaAnalysis::InputUpdateFromSolution(IssmDouble* solution,Element* element){/*{{{*/ - _error_("not implemented yet"); -}/*}}}*/ -void GiaAnalysis::UpdateConstraints(FemModel* femmodel){/*{{{*/ - /*Default, do nothing*/ - return; -}/*}}}*/ Index: sm/trunk-jpl/src/c/analyses/GiaAnalysis.h =================================================================== --- /issm/trunk-jpl/src/c/analyses/GiaAnalysis.h (revision 21531) +++ (revision ) @@ -1,33 +1,0 @@ -/*! \file GiaAnalysis.h - * \brief: header file for generic external result object - */ - -#ifndef _GiaAnalysis_ -#define _GiaAnalysis_ - -/*Headers*/ -#include "./Analysis.h" - -class GiaAnalysis: public Analysis{ - - public: - /*Model processing*/ - void CreateConstraints(Constraints* constraints,IoModel* iomodel); - void CreateLoads(Loads* loads, IoModel* iomodel); - void CreateNodes(Nodes* nodes,IoModel* iomodel); - int DofsPerNode(int** doflist,int domaintype,int approximation); - void UpdateElements(Elements* elements,IoModel* iomodel,int analysis_counter,int analysis_type); - void UpdateParameters(Parameters* parameters,IoModel* iomodel,int solution_enum,int analysis_enum); - - /*Finite element Analysis*/ - void Core(FemModel* femmodel); - ElementVector* CreateDVector(Element* element); - ElementMatrix* CreateJacobianMatrix(Element* element); - ElementMatrix* CreateKMatrix(Element* element); - ElementVector* CreatePVector(Element* element); - void GetSolutionFromInputs(Vector* solution,Element* element); - void GradientJ(Vector* gradient,Element* element,int control_type,int control_index); - void InputUpdateFromSolution(IssmDouble* solution,Element* element); - void UpdateConstraints(FemModel* femmodel); -}; -#endif Index: /issm/trunk-jpl/src/c/analyses/GiaIvinsAnalysis.cpp =================================================================== --- /issm/trunk-jpl/src/c/analyses/GiaIvinsAnalysis.cpp (revision 21532) +++ /issm/trunk-jpl/src/c/analyses/GiaIvinsAnalysis.cpp (revision 21532) @@ -0,0 +1,69 @@ +#include "./GiaIvinsAnalysis.h" +#include "../toolkits/toolkits.h" +#include "../classes/classes.h" +#include "../shared/shared.h" +#include "../modules/modules.h" + +/*Model processing*/ +void GiaIvinsAnalysis::CreateConstraints(Constraints* constraints,IoModel* iomodel){/*{{{*/ + /*No constraints*/ +}/*}}}*/ +void GiaIvinsAnalysis::CreateLoads(Loads* loads, IoModel* iomodel){/*{{{*/ + /*No loads*/ +}/*}}}*/ +void GiaIvinsAnalysis::CreateNodes(Nodes* nodes,IoModel* iomodel){/*{{{*/ + ::CreateNodes(nodes,iomodel,GiaIvinsAnalysisEnum,P1Enum); +}/*}}}*/ +int GiaIvinsAnalysis::DofsPerNode(int** doflist,int domaintype,int approximation){/*{{{*/ + return 1; +}/*}}}*/ +void GiaIvinsAnalysis::UpdateElements(Elements* elements,IoModel* iomodel,int analysis_counter,int analysis_type){/*{{{*/ + + /*Update elements: */ + int counter=0; + for(int i=0;inumberofelements;i++){ + if(iomodel->my_elements[i]){ + Element* element=(Element*)elements->GetObjectByOffset(counter); + element->Update(i,iomodel,analysis_counter,analysis_type,P1Enum); + counter++; + } + } + + iomodel->FetchDataToInput(elements,"md.geometry.thickness",ThicknessEnum); + iomodel->FetchDataToInput(elements,"md.giaivins.mantle_viscosity",GiaIvinsMantleViscosityEnum); + iomodel->FetchDataToInput(elements,"md.giaivins.lithosphere_thickness",GiaIvinsLithosphereThicknessEnum); + iomodel->FetchDataToInput(elements,"md.mask.ice_levelset",MaskIceLevelsetEnum); +}/*}}}*/ +void GiaIvinsAnalysis::UpdateParameters(Parameters* parameters,IoModel* iomodel,int solution_enum,int analysis_enum){/*{{{*/ +}/*}}}*/ + +/*Finite Element Analysis*/ +void GiaIvinsAnalysis::Core(FemModel* femmodel){/*{{{*/ + _error_("not implemented"); +}/*}}}*/ +ElementVector* GiaIvinsAnalysis::CreateDVector(Element* element){/*{{{*/ + /*Default, return NULL*/ + return NULL; +}/*}}}*/ +ElementMatrix* GiaIvinsAnalysis::CreateJacobianMatrix(Element* element){/*{{{*/ +_error_("Not implemented"); +}/*}}}*/ +ElementMatrix* GiaIvinsAnalysis::CreateKMatrix(Element* element){/*{{{*/ + _error_("not implemented yet"); +}/*}}}*/ +ElementVector* GiaIvinsAnalysis::CreatePVector(Element* element){/*{{{*/ +_error_("not implemented yet"); +}/*}}}*/ +void GiaIvinsAnalysis::GetSolutionFromInputs(Vector* solution,Element* element){/*{{{*/ + _error_("not implemented yet"); +}/*}}}*/ +void GiaIvinsAnalysis::GradientJ(Vector* gradient,Element* element,int control_type,int control_index){/*{{{*/ + _error_("Not implemented yet"); +}/*}}}*/ +void GiaIvinsAnalysis::InputUpdateFromSolution(IssmDouble* solution,Element* element){/*{{{*/ + _error_("not implemented yet"); +}/*}}}*/ +void GiaIvinsAnalysis::UpdateConstraints(FemModel* femmodel){/*{{{*/ + /*Default, do nothing*/ + return; +}/*}}}*/ Index: /issm/trunk-jpl/src/c/analyses/GiaIvinsAnalysis.h =================================================================== --- /issm/trunk-jpl/src/c/analyses/GiaIvinsAnalysis.h (revision 21532) +++ /issm/trunk-jpl/src/c/analyses/GiaIvinsAnalysis.h (revision 21532) @@ -0,0 +1,33 @@ +/*! \file GiaIvinsAnalysis.h + * \brief: header file for generic external result object + */ + +#ifndef _GiaIvinsAnalysis_ +#define _GiaIvinsAnalysis_ + +/*Headers*/ +#include "./Analysis.h" + +class GiaIvinsAnalysis: public Analysis{ + + public: + /*Model processing*/ + void CreateConstraints(Constraints* constraints,IoModel* iomodel); + void CreateLoads(Loads* loads, IoModel* iomodel); + void CreateNodes(Nodes* nodes,IoModel* iomodel); + int DofsPerNode(int** doflist,int domaintype,int approximation); + void UpdateElements(Elements* elements,IoModel* iomodel,int analysis_counter,int analysis_type); + void UpdateParameters(Parameters* parameters,IoModel* iomodel,int solution_enum,int analysis_enum); + + /*Finite element Analysis*/ + void Core(FemModel* femmodel); + ElementVector* CreateDVector(Element* element); + ElementMatrix* CreateJacobianMatrix(Element* element); + ElementMatrix* CreateKMatrix(Element* element); + ElementVector* CreatePVector(Element* element); + void GetSolutionFromInputs(Vector* solution,Element* element); + void GradientJ(Vector* gradient,Element* element,int control_type,int control_index); + void InputUpdateFromSolution(IssmDouble* solution,Element* element); + void UpdateConstraints(FemModel* femmodel); +}; +#endif Index: sm/trunk-jpl/src/c/classes/GiaDeflectionCoreArgs.h =================================================================== --- /issm/trunk-jpl/src/c/classes/GiaDeflectionCoreArgs.h (revision 21531) +++ (revision ) @@ -1,40 +1,0 @@ -/*!\file: GiaDeflectionCoreArgs.h - * \brief place holder for arguments to the GiaDeflectionCoreArgs routine - */ - -#ifndef _GIADEFLECTIONCOREARGS_H_ -#define _GIADEFLECTIONCOREARGS_H_ - -struct GiaDeflectionCoreArgs{ - - /*inputs: */ - IssmDouble currenttime; - IssmDouble* hes; //loading history (in ice thickness) - int numtimes; //loading history length - IssmDouble ri; //radial distance from center of disk to vertex i - IssmDouble re; //radius of disk - IssmDouble* times; //loading history times - - /*gia material parameters: */ - IssmDouble lithosphere_density; - IssmDouble lithosphere_shear_modulus; - IssmDouble lithosphere_thickness; - IssmDouble mantle_density; - IssmDouble mantle_shear_modulus; - IssmDouble mantle_viscosity; - - /*gia solution parameters: */ - int iedge; - - /*ice properties: */ - IssmDouble rho_ice; - - /*constants: */ - IssmDouble yts; - - /*debug info: */ - int idisk; //id of the element we are running the gia code in. - -}; - -#endif Index: /issm/trunk-jpl/src/c/classes/GiaIvinsDeflectionCoreArgs.h =================================================================== --- /issm/trunk-jpl/src/c/classes/GiaIvinsDeflectionCoreArgs.h (revision 21532) +++ /issm/trunk-jpl/src/c/classes/GiaIvinsDeflectionCoreArgs.h (revision 21532) @@ -0,0 +1,40 @@ +/*!\file: GiaIvinsDeflectionCoreArgs.h + * \brief place holder for arguments to the GiaIvinsDeflectionCoreArgs routine + */ + +#ifndef _GIAIVINSDEFLECTIONCOREARGS_H_ +#define _GIAIVINSDEFLECTIONCOREARGS_H_ + +struct GiaIvinsDeflectionCoreArgs{ + + /*inputs: */ + IssmDouble currenttime; + IssmDouble* hes; //loading history (in ice thickness) + int numtimes; //loading history length + IssmDouble ri; //radial distance from center of disk to vertex i + IssmDouble re; //radius of disk + IssmDouble* times; //loading history times + + /*gia material parameters: */ + IssmDouble lithosphere_density; + IssmDouble lithosphere_shear_modulus; + IssmDouble lithosphere_thickness; + IssmDouble mantle_density; + IssmDouble mantle_shear_modulus; + IssmDouble mantle_viscosity; + + /*gia solution parameters: */ + int iedge; + + /*ice properties: */ + IssmDouble rho_ice; + + /*constants: */ + IssmDouble yts; + + /*debug info: */ + int idisk; //id of the element we are running the gia code in. + +}; + +#endif Index: sm/trunk-jpl/src/c/cores/gia_core.cpp =================================================================== --- /issm/trunk-jpl/src/c/cores/gia_core.cpp (revision 21531) +++ (revision ) @@ -1,60 +1,0 @@ -/*!\file: gia_core.cpp - * \brief: core of the GIA solution - */ - -#include "./cores.h" -#include "../toolkits/toolkits.h" -#include "../classes/classes.h" -#include "../shared/shared.h" -#include "../modules/modules.h" -#include "../solutionsequences/solutionsequences.h" -void gia_core(FemModel* femmodel){ - - Vector *wg = NULL; - Vector *dwdtg = NULL; - IssmDouble *x = NULL; - IssmDouble *y = NULL; - - /*parameters: */ - bool save_results; - int gsize; - int configuration_type; - - /*Recover some parameters: */ - femmodel->parameters->FindParam(&save_results,SaveResultsEnum); - femmodel->parameters->FindParam(&configuration_type,ConfigurationTypeEnum); - - if(VerboseSolution()) _printf0_(" computing GIA\n"); - - /*Call on core computations: */ - femmodel->SetCurrentConfiguration(GiaAnalysisEnum); - - /*Figure out size of g-set deflection vector and allocate solution vector: */ - gsize = femmodel->nodes->NumberOfDofs(configuration_type,GsetEnum); - wg = new Vector(gsize); - dwdtg = new Vector(gsize); - - /*first, recover x and y vectors from vertices: */ - VertexCoordinatesx(&x,&y,NULL,femmodel->vertices); //no need for z coordinate - - /*call the main module: */ - femmodel->Deflection(wg,dwdtg,x,y); - - /*assemble vector: */ - wg->Assemble(); - dwdtg->Assemble(); - - InputUpdateFromVectorx(femmodel,wg,GiaWEnum,VertexSIdEnum); - InputUpdateFromVectorx(femmodel,dwdtg,GiadWdtEnum,VertexSIdEnum); - - if(save_results){ - if(VerboseSolution()) _printf0_(" saving results\n"); - int outputs[2] = {GiaWEnum,GiadWdtEnum}; - femmodel->RequestedOutputsx(&femmodel->results,&outputs[0],2); - } - - xDelete(x); - xDelete(y); - - -} Index: /issm/trunk-jpl/src/c/cores/giaivins_core.cpp =================================================================== --- /issm/trunk-jpl/src/c/cores/giaivins_core.cpp (revision 21532) +++ /issm/trunk-jpl/src/c/cores/giaivins_core.cpp (revision 21532) @@ -0,0 +1,60 @@ +/*!\file: giaivins_core.cpp + * \brief: core of the GIAIVINS solution + */ + +#include "./cores.h" +#include "../toolkits/toolkits.h" +#include "../classes/classes.h" +#include "../shared/shared.h" +#include "../modules/modules.h" +#include "../solutionsequences/solutionsequences.h" +void giaivins_core(FemModel* femmodel){ + + Vector *wg = NULL; + Vector *dwdtg = NULL; + IssmDouble *x = NULL; + IssmDouble *y = NULL; + + /*parameters: */ + bool save_results; + int gsize; + int configuration_type; + + /*Recover some parameters: */ + femmodel->parameters->FindParam(&save_results,SaveResultsEnum); + femmodel->parameters->FindParam(&configuration_type,ConfigurationTypeEnum); + + if(VerboseSolution()) _printf0_(" computing GIAIVINS\n"); + + /*Call on core computations: */ + femmodel->SetCurrentConfiguration(GiaIvinsAnalysisEnum); + + /*Figure out size of g-set deflection vector and allocate solution vector: */ + gsize = femmodel->nodes->NumberOfDofs(configuration_type,GsetEnum); + wg = new Vector(gsize); + dwdtg = new Vector(gsize); + + /*first, recover x and y vectors from vertices: */ + VertexCoordinatesx(&x,&y,NULL,femmodel->vertices); //no need for z coordinate + + /*call the main module: */ + femmodel->Deflection(wg,dwdtg,x,y); + + /*assemble vector: */ + wg->Assemble(); + dwdtg->Assemble(); + + InputUpdateFromVectorx(femmodel,wg,GiaIvinsWEnum,VertexSIdEnum); + InputUpdateFromVectorx(femmodel,dwdtg,GiaIvinsdWdtEnum,VertexSIdEnum); + + if(save_results){ + if(VerboseSolution()) _printf0_(" saving results\n"); + int outputs[2] = {GiaIvinsWEnum,GiaIvinsdWdtEnum}; + femmodel->RequestedOutputsx(&femmodel->results,&outputs[0],2); + } + + xDelete(x); + xDelete(y); + + +} Index: /issm/trunk-jpl/src/c/modules/GiaIvinsDeflectionCorex/.deps/libISSMCore_la-GiaIvinsDeflectionCorex.Plo =================================================================== --- /issm/trunk-jpl/src/c/modules/GiaIvinsDeflectionCorex/.deps/libISSMCore_la-GiaIvinsDeflectionCorex.Plo (revision 21532) +++ /issm/trunk-jpl/src/c/modules/GiaIvinsDeflectionCorex/.deps/libISSMCore_la-GiaIvinsDeflectionCorex.Plo (revision 21532) @@ -0,0 +1,1243 @@ +./modules/GiaIvinsDeflectionCorex/libISSMCore_la-GiaIvinsDeflectionCorex.lo: \ + modules/GiaIvinsDeflectionCorex/GiaIvinsDeflectionCorex.cpp \ + modules/GiaIvinsDeflectionCorex/./GiaIvinsDeflectionCorex.h \ + modules/GiaIvinsDeflectionCorex/../../classes/classes.h \ + modules/GiaIvinsDeflectionCorex/../../classes/./Contour.h \ + modules/GiaIvinsDeflectionCorex/../../classes/../shared/shared.h \ + modules/GiaIvinsDeflectionCorex/../../classes/../shared/./FSanalyticals/fsanalyticals.h \ + modules/GiaIvinsDeflectionCorex/../../classes/../shared/./FSanalyticals/../Numerics/types.h \ + ../.././config.h /usr/include/stdio.h /usr/include/sys/cdefs.h \ + /usr/include/sys/_symbol_aliasing.h \ + /usr/include/sys/_posix_availability.h /usr/include/Availability.h \ + /usr/include/AvailabilityInternal.h /usr/include/_types.h \ + /usr/include/sys/_types.h /usr/include/machine/_types.h \ + /usr/include/i386/_types.h /usr/include/sys/_pthread/_pthread_types.h \ + /usr/include/sys/_types/_va_list.h /usr/include/sys/_types/_size_t.h \ + /usr/include/sys/_types/_null.h /usr/include/sys/stdio.h \ + /usr/include/sys/_types/_off_t.h /usr/include/sys/_types/_ssize_t.h \ + modules/GiaIvinsDeflectionCorex/../../classes/../shared/./MemOps/MemOps.h \ + /Applications/Xcode.app/Contents/Developer/Toolchains/XcodeDefault.xctoolchain/usr/bin/../include/c++/v1/cassert \ + /Applications/Xcode.app/Contents/Developer/Toolchains/XcodeDefault.xctoolchain/usr/bin/../include/c++/v1/__config \ + /usr/include/unistd.h /usr/include/sys/unistd.h \ + /usr/include/sys/_types/_posix_vdisable.h \ + /usr/include/sys/_types/_seek_set.h /usr/include/_types/_uint64_t.h \ + /usr/include/sys/_types/_uid_t.h /usr/include/sys/_types/_gid_t.h \ + /usr/include/sys/_types/_intptr_t.h /usr/include/sys/_types/_pid_t.h \ + /usr/include/sys/_types/_useconds_t.h /usr/include/sys/select.h \ + /usr/include/sys/appleapiopts.h /usr/include/sys/_types/_fd_def.h \ + /usr/include/sys/_types/_timespec.h /usr/include/sys/_types/_timeval.h \ + /usr/include/sys/_types/_time_t.h \ + /usr/include/sys/_types/_suseconds_t.h \ + /usr/include/sys/_types/_sigset_t.h \ + /usr/include/sys/_types/_fd_setsize.h \ + /usr/include/sys/_types/_fd_set.h /usr/include/sys/_types/_fd_clr.h \ + /usr/include/sys/_types/_fd_isset.h /usr/include/sys/_types/_fd_zero.h \ + /usr/include/sys/_types/_fd_copy.h /usr/include/sys/_select.h \ + /usr/include/sys/_types/_dev_t.h /usr/include/sys/_types/_mode_t.h \ + /usr/include/sys/_types/_uuid_t.h /usr/include/gethostuuid.h \ + /usr/include/assert.h /usr/include/stdlib.h /usr/include/sys/wait.h \ + /usr/include/sys/_types/_id_t.h /usr/include/sys/signal.h \ + /usr/include/machine/signal.h /usr/include/i386/signal.h \ + /usr/include/machine/_mcontext.h /usr/include/i386/_mcontext.h \ + /usr/include/mach/i386/_structs.h \ + /usr/include/sys/_pthread/_pthread_attr_t.h \ + /usr/include/sys/_types/_sigaltstack.h \ + /usr/include/sys/_types/_ucontext.h /usr/include/sys/resource.h \ + /Applications/Xcode.app/Contents/Developer/Toolchains/XcodeDefault.xctoolchain/usr/bin/../lib/clang/7.0.2/include/stdint.h \ + /usr/include/stdint.h /usr/include/sys/_types/_int8_t.h \ + /usr/include/sys/_types/_int16_t.h /usr/include/sys/_types/_int32_t.h \ + /usr/include/sys/_types/_int64_t.h /usr/include/_types/_uint8_t.h \ + /usr/include/_types/_uint16_t.h /usr/include/_types/_uint32_t.h \ + /usr/include/sys/_types/_uintptr_t.h /usr/include/_types/_intmax_t.h \ + /usr/include/_types/_uintmax_t.h /usr/include/machine/endian.h \ + /usr/include/i386/endian.h /usr/include/sys/_endian.h \ + /usr/include/libkern/_OSByteOrder.h \ + /usr/include/libkern/i386/_OSByteOrder.h /usr/include/alloca.h \ + /usr/include/sys/_types/_ct_rune_t.h /usr/include/sys/_types/_rune_t.h \ + /usr/include/sys/_types/_wchar_t.h /usr/include/machine/types.h \ + /usr/include/i386/types.h \ + modules/GiaIvinsDeflectionCorex/../../classes/../shared/./Bamg/shared.h \ + modules/GiaIvinsDeflectionCorex/../../classes/../shared/./Bamg/Abs.h \ + modules/GiaIvinsDeflectionCorex/../../classes/../shared/./Bamg/BigPrimeNumber.h \ + modules/GiaIvinsDeflectionCorex/../../classes/../shared/./Bamg/BinaryRand.h \ + /Applications/Xcode.app/Contents/Developer/Toolchains/XcodeDefault.xctoolchain/usr/bin/../include/c++/v1/cstdlib \ + modules/GiaIvinsDeflectionCorex/../../classes/../shared/./Bamg/Exchange.h \ + modules/GiaIvinsDeflectionCorex/../../classes/../shared/./Bamg/extrema.h \ + modules/GiaIvinsDeflectionCorex/../../classes/../shared/./Bamg/HeapSort.h \ + modules/GiaIvinsDeflectionCorex/../../classes/../shared/./Bamg/OppositeAngle.h \ + modules/GiaIvinsDeflectionCorex/../../classes/../shared/./Bamg/../Numerics/constants.h \ + modules/GiaIvinsDeflectionCorex/../../classes/../shared/./Elements/elements.h \ + modules/GiaIvinsDeflectionCorex/../../classes/../shared/./Enum/Enum.h \ + modules/GiaIvinsDeflectionCorex/../../classes/../shared/./Enum/./EnumDefinitions.h \ + modules/GiaIvinsDeflectionCorex/../../classes/../shared/./Exceptions/exceptions.h \ + /Applications/Xcode.app/Contents/Developer/Toolchains/XcodeDefault.xctoolchain/usr/bin/../include/c++/v1/exception \ + /Applications/Xcode.app/Contents/Developer/Toolchains/XcodeDefault.xctoolchain/usr/bin/../include/c++/v1/cstddef \ + /Applications/Xcode.app/Contents/Developer/Toolchains/XcodeDefault.xctoolchain/usr/bin/../lib/clang/7.0.2/include/stddef.h \ + /Applications/Xcode.app/Contents/Developer/Toolchains/XcodeDefault.xctoolchain/usr/bin/../include/c++/v1/type_traits \ + /Applications/Xcode.app/Contents/Developer/Toolchains/XcodeDefault.xctoolchain/usr/bin/../include/c++/v1/string \ + /Applications/Xcode.app/Contents/Developer/Toolchains/XcodeDefault.xctoolchain/usr/bin/../include/c++/v1/iosfwd \ + /usr/include/wchar.h /usr/include/sys/_types/_mbstate_t.h \ + /Applications/Xcode.app/Contents/Developer/Toolchains/XcodeDefault.xctoolchain/usr/bin/../lib/clang/7.0.2/include/stdarg.h \ + /usr/include/time.h /usr/include/sys/_types/_clock_t.h \ + /usr/include/_wctype.h /usr/include/sys/_types/_wint_t.h \ + /usr/include/_types/_wctype_t.h /usr/include/ctype.h \ + /usr/include/runetype.h \ + /Applications/Xcode.app/Contents/Developer/Toolchains/XcodeDefault.xctoolchain/usr/bin/../include/c++/v1/cstring \ + /usr/include/string.h /usr/include/sys/_types/_rsize_t.h \ + /usr/include/sys/_types/_errno_t.h /usr/include/strings.h \ + /Applications/Xcode.app/Contents/Developer/Toolchains/XcodeDefault.xctoolchain/usr/bin/../include/c++/v1/cstdio \ + /Applications/Xcode.app/Contents/Developer/Toolchains/XcodeDefault.xctoolchain/usr/bin/../include/c++/v1/cwchar \ + /Applications/Xcode.app/Contents/Developer/Toolchains/XcodeDefault.xctoolchain/usr/bin/../include/c++/v1/cwctype \ + /Applications/Xcode.app/Contents/Developer/Toolchains/XcodeDefault.xctoolchain/usr/bin/../include/c++/v1/cctype \ + /usr/include/wctype.h /usr/include/_types/_wctrans_t.h \ + /Applications/Xcode.app/Contents/Developer/Toolchains/XcodeDefault.xctoolchain/usr/bin/../include/c++/v1/algorithm \ + /Applications/Xcode.app/Contents/Developer/Toolchains/XcodeDefault.xctoolchain/usr/bin/../include/c++/v1/initializer_list \ + /Applications/Xcode.app/Contents/Developer/Toolchains/XcodeDefault.xctoolchain/usr/bin/../include/c++/v1/utility \ + /Applications/Xcode.app/Contents/Developer/Toolchains/XcodeDefault.xctoolchain/usr/bin/../include/c++/v1/__tuple \ + /Applications/Xcode.app/Contents/Developer/Toolchains/XcodeDefault.xctoolchain/usr/bin/../include/c++/v1/__tuple_03 \ + /Applications/Xcode.app/Contents/Developer/Toolchains/XcodeDefault.xctoolchain/usr/bin/../include/c++/v1/memory \ + /Applications/Xcode.app/Contents/Developer/Toolchains/XcodeDefault.xctoolchain/usr/bin/../include/c++/v1/typeinfo \ + /Applications/Xcode.app/Contents/Developer/Toolchains/XcodeDefault.xctoolchain/usr/bin/../include/c++/v1/cstdint \ + /Applications/Xcode.app/Contents/Developer/Toolchains/XcodeDefault.xctoolchain/usr/bin/../include/c++/v1/new \ + /Applications/Xcode.app/Contents/Developer/Toolchains/XcodeDefault.xctoolchain/usr/bin/../include/c++/v1/limits \ + /Applications/Xcode.app/Contents/Developer/Toolchains/XcodeDefault.xctoolchain/usr/bin/../include/c++/v1/__undef_min_max \ + /Applications/Xcode.app/Contents/Developer/Toolchains/XcodeDefault.xctoolchain/usr/bin/../include/c++/v1/iterator \ + /Applications/Xcode.app/Contents/Developer/Toolchains/XcodeDefault.xctoolchain/usr/bin/../include/c++/v1/__functional_base \ + /Applications/Xcode.app/Contents/Developer/Toolchains/XcodeDefault.xctoolchain/usr/bin/../include/c++/v1/__functional_base_03 \ + /Applications/Xcode.app/Contents/Developer/Toolchains/XcodeDefault.xctoolchain/usr/bin/../include/c++/v1/__debug \ + /Applications/Xcode.app/Contents/Developer/Toolchains/XcodeDefault.xctoolchain/usr/bin/../include/c++/v1/tuple \ + /Applications/Xcode.app/Contents/Developer/Toolchains/XcodeDefault.xctoolchain/usr/bin/../include/c++/v1/stdexcept \ + /Applications/Xcode.app/Contents/Developer/Toolchains/XcodeDefault.xctoolchain/usr/bin/../include/c++/v1/iostream \ + /Applications/Xcode.app/Contents/Developer/Toolchains/XcodeDefault.xctoolchain/usr/bin/../include/c++/v1/ios \ + /Applications/Xcode.app/Contents/Developer/Toolchains/XcodeDefault.xctoolchain/usr/bin/../include/c++/v1/__locale \ + /Applications/Xcode.app/Contents/Developer/Toolchains/XcodeDefault.xctoolchain/usr/bin/../include/c++/v1/mutex \ + /Applications/Xcode.app/Contents/Developer/Toolchains/XcodeDefault.xctoolchain/usr/bin/../include/c++/v1/__mutex_base \ + /Applications/Xcode.app/Contents/Developer/Toolchains/XcodeDefault.xctoolchain/usr/bin/../include/c++/v1/chrono \ + /Applications/Xcode.app/Contents/Developer/Toolchains/XcodeDefault.xctoolchain/usr/bin/../include/c++/v1/ctime \ + /Applications/Xcode.app/Contents/Developer/Toolchains/XcodeDefault.xctoolchain/usr/bin/../include/c++/v1/ratio \ + /Applications/Xcode.app/Contents/Developer/Toolchains/XcodeDefault.xctoolchain/usr/bin/../include/c++/v1/climits \ + /Applications/Xcode.app/Contents/Developer/Toolchains/XcodeDefault.xctoolchain/usr/bin/../lib/clang/7.0.2/include/limits.h \ + /usr/include/limits.h /usr/include/machine/limits.h \ + /usr/include/i386/limits.h /usr/include/i386/_limits.h \ + /usr/include/sys/syslimits.h \ + /Applications/Xcode.app/Contents/Developer/Toolchains/XcodeDefault.xctoolchain/usr/bin/../include/c++/v1/system_error \ + /Applications/Xcode.app/Contents/Developer/Toolchains/XcodeDefault.xctoolchain/usr/bin/../include/c++/v1/cerrno \ + /usr/include/errno.h /usr/include/sys/errno.h /usr/include/pthread.h \ + /usr/include/pthread/pthread_impl.h /usr/include/pthread/sched.h \ + /usr/include/sys/_pthread/_pthread_cond_t.h \ + /usr/include/sys/_pthread/_pthread_condattr_t.h \ + /usr/include/sys/_pthread/_pthread_key_t.h \ + /usr/include/sys/_pthread/_pthread_mutex_t.h \ + /usr/include/sys/_pthread/_pthread_mutexattr_t.h \ + /usr/include/sys/_pthread/_pthread_once_t.h \ + /usr/include/sys/_pthread/_pthread_rwlock_t.h \ + /usr/include/sys/_pthread/_pthread_rwlockattr_t.h \ + /usr/include/sys/_pthread/_pthread_t.h /usr/include/pthread/qos.h \ + /usr/include/sys/qos.h /usr/include/sys/_types/_mach_port_t.h \ + /Applications/Xcode.app/Contents/Developer/Toolchains/XcodeDefault.xctoolchain/usr/bin/../include/c++/v1/functional \ + /Applications/Xcode.app/Contents/Developer/Toolchains/XcodeDefault.xctoolchain/usr/bin/../include/c++/v1/__functional_03 \ + /usr/include/locale.h /usr/include/_locale.h /usr/include/xlocale.h \ + /usr/include/_xlocale.h /usr/include/xlocale/_ctype.h \ + /usr/include/xlocale/__wctype.h /usr/include/xlocale/_stdio.h \ + /usr/include/xlocale/_stdlib.h /usr/include/xlocale/_string.h \ + /usr/include/xlocale/_time.h /usr/include/xlocale/_wchar.h \ + /usr/include/xlocale/_wctype.h \ + /Applications/Xcode.app/Contents/Developer/Toolchains/XcodeDefault.xctoolchain/usr/bin/../include/c++/v1/streambuf \ + /Applications/Xcode.app/Contents/Developer/Toolchains/XcodeDefault.xctoolchain/usr/bin/../include/c++/v1/istream \ + /Applications/Xcode.app/Contents/Developer/Toolchains/XcodeDefault.xctoolchain/usr/bin/../include/c++/v1/ostream \ + /Applications/Xcode.app/Contents/Developer/Toolchains/XcodeDefault.xctoolchain/usr/bin/../include/c++/v1/locale \ + /usr/include/nl_types.h /usr/include/sys/types.h \ + /usr/include/sys/_types/_blkcnt_t.h \ + /usr/include/sys/_types/_blksize_t.h \ + /usr/include/sys/_types/_in_addr_t.h \ + /usr/include/sys/_types/_in_port_t.h /usr/include/sys/_types/_ino_t.h \ + /usr/include/sys/_types/_ino64_t.h /usr/include/sys/_types/_key_t.h \ + /usr/include/sys/_types/_nlink_t.h \ + /usr/include/sys/_types/_fsblkcnt_t.h \ + /usr/include/sys/_types/_fsfilcnt_t.h /usr/include/_types/_nl_item.h \ + /Applications/Xcode.app/Contents/Developer/Toolchains/XcodeDefault.xctoolchain/usr/bin/../include/c++/v1/bitset \ + /Applications/Xcode.app/Contents/Developer/Toolchains/XcodeDefault.xctoolchain/usr/bin/../include/c++/v1/__bit_reference \ + /Applications/Xcode.app/Contents/Developer/Toolchains/XcodeDefault.xctoolchain/usr/bin/../include/c++/v1/sstream \ + /Applications/Xcode.app/Contents/Developer/Toolchains/XcodeDefault.xctoolchain/usr/bin/../include/c++/v1/iomanip \ + modules/GiaIvinsDeflectionCorex/../../classes/../shared/./io/io.h \ + modules/GiaIvinsDeflectionCorex/../../classes/../shared/./io/./Comm/IssmComm.h \ + modules/GiaIvinsDeflectionCorex/../../classes/../shared/./io/./Comm/../../../toolkits/mpi/issmmpi.h \ + /Users/adhikari/issm/trunk-jpl-sesaw/externalpackages/mpich/install/include/mpi.h \ + /Users/adhikari/issm/trunk-jpl-sesaw/externalpackages/mpich/install/include/mpio.h \ + /Users/adhikari/issm/trunk-jpl-sesaw/externalpackages/mpich/install/include/mpicxx.h \ + modules/GiaIvinsDeflectionCorex/../../classes/../shared/./io/./Comm/../../../toolkits/mpi/./commops/commops.h \ + modules/GiaIvinsDeflectionCorex/../../classes/../shared/./io/./Comm/../../../toolkits/mpi/./commops/../../../shared/io/Comm/IssmComm.h \ + modules/GiaIvinsDeflectionCorex/../../classes/../shared/./io/./Comm/../../../toolkits/mpi/./commops/../../mpi/issmmpi.h \ + modules/GiaIvinsDeflectionCorex/../../classes/../shared/./io/./Disk/diskio.h \ + modules/GiaIvinsDeflectionCorex/../../classes/../shared/./io/./Print/Print.h \ + modules/GiaIvinsDeflectionCorex/../../classes/../shared/./io/./Marshalling/Marshalling.h \ + modules/GiaIvinsDeflectionCorex/../../classes/../shared/./io/./Marshalling/IoCodeConversions.h \ + modules/GiaIvinsDeflectionCorex/../../classes/../shared/./Exp/exp.h \ + modules/GiaIvinsDeflectionCorex/../../classes/../shared/./Exp/../Numerics/recast.h \ + modules/GiaIvinsDeflectionCorex/../../classes/../shared/./Matrix/matrix.h \ + modules/GiaIvinsDeflectionCorex/../../classes/../shared/./Numerics/numerics.h \ + modules/GiaIvinsDeflectionCorex/../../classes/../shared/./FSanalyticals/../Numerics/./Verbosity.h \ + modules/GiaIvinsDeflectionCorex/../../classes/../shared/./FSanalyticals/../Numerics/./GaussPoints.h \ + modules/GiaIvinsDeflectionCorex/../../classes/../shared/./FSanalyticals/../Numerics/./isnan.h \ + /Applications/Xcode.app/Contents/Developer/Toolchains/XcodeDefault.xctoolchain/usr/bin/../include/c++/v1/cmath \ + /usr/include/math.h \ + modules/GiaIvinsDeflectionCorex/../../classes/../shared/./FSanalyticals/../Numerics/./OptPars.h \ + modules/GiaIvinsDeflectionCorex/../../classes/../shared/./Sorting/sorting.h \ + modules/GiaIvinsDeflectionCorex/../../classes/../shared/./String/sharedstring.h \ + modules/GiaIvinsDeflectionCorex/../../classes/../shared/./Threads/issm_threads.h \ + modules/GiaIvinsDeflectionCorex/../../classes/../shared/./TriMesh/trimesh.h \ + modules/GiaIvinsDeflectionCorex/../../classes/../shared/./LatLong/latlong.h \ + modules/GiaIvinsDeflectionCorex/../../classes/../datastructures/datastructures.h \ + modules/GiaIvinsDeflectionCorex/../../classes/../datastructures/./DataSet.h \ + /Applications/Xcode.app/Contents/Developer/Toolchains/XcodeDefault.xctoolchain/usr/bin/../include/c++/v1/vector \ + /Applications/Xcode.app/Contents/Developer/Toolchains/XcodeDefault.xctoolchain/usr/bin/../include/c++/v1/__split_buffer \ + modules/GiaIvinsDeflectionCorex/../../classes/../datastructures/./Object.h \ + modules/GiaIvinsDeflectionCorex/../../classes/./Vertices.h \ + modules/GiaIvinsDeflectionCorex/../../classes/./Vertex.h \ + modules/GiaIvinsDeflectionCorex/../../classes/./classes.h \ + modules/GiaIvinsDeflectionCorex/../../classes/../toolkits/toolkits.h \ + modules/GiaIvinsDeflectionCorex/../../classes/../toolkits/./petsc/petscincludes.h \ + /Users/adhikari/issm/trunk-jpl-sesaw/externalpackages/petsc/install/include/petscmat.h \ + /Users/adhikari/issm/trunk-jpl-sesaw/externalpackages/petsc/install/include/petscvec.h \ + /Users/adhikari/issm/trunk-jpl-sesaw/externalpackages/petsc/install/include/petscis.h \ + /Users/adhikari/issm/trunk-jpl-sesaw/externalpackages/petsc/install/include/petscsys.h \ + /Users/adhikari/issm/trunk-jpl-sesaw/externalpackages/petsc/install/include/petscconf.h \ + /Users/adhikari/issm/trunk-jpl-sesaw/externalpackages/petsc/install/include/petscfix.h \ + /Users/adhikari/issm/trunk-jpl-sesaw/externalpackages/petsc/install/include/petscversion.h \ + /Users/adhikari/issm/trunk-jpl-sesaw/externalpackages/petsc/install/include/petscmath.h \ + /Applications/Xcode.app/Contents/Developer/Toolchains/XcodeDefault.xctoolchain/usr/bin/../include/c++/v1/complex \ + /Users/adhikari/issm/trunk-jpl-sesaw/externalpackages/petsc/install/include/petscerror.h \ + /Users/adhikari/issm/trunk-jpl-sesaw/externalpackages/petsc/install/include/petscviewertypes.h \ + /Users/adhikari/issm/trunk-jpl-sesaw/externalpackages/petsc/install/include/petscoptions.h \ + /Users/adhikari/issm/trunk-jpl-sesaw/externalpackages/petsc/install/include/petsclog.h \ + /Applications/Xcode.app/Contents/Developer/Toolchains/XcodeDefault.xctoolchain/usr/bin/../lib/clang/7.0.2/include/xmmintrin.h \ + /Applications/Xcode.app/Contents/Developer/Toolchains/XcodeDefault.xctoolchain/usr/bin/../lib/clang/7.0.2/include/mmintrin.h \ + /Applications/Xcode.app/Contents/Developer/Toolchains/XcodeDefault.xctoolchain/usr/bin/../lib/clang/7.0.2/include/mm_malloc.h \ + /Applications/Xcode.app/Contents/Developer/Toolchains/XcodeDefault.xctoolchain/usr/bin/../lib/clang/7.0.2/include/emmintrin.h \ + /usr/include/sys/param.h /usr/include/machine/param.h \ + /usr/include/i386/param.h /usr/include/i386/_param.h \ + /Users/adhikari/issm/trunk-jpl-sesaw/externalpackages/petsc/install/include/petscsftypes.h \ + /Users/adhikari/issm/trunk-jpl-sesaw/externalpackages/petsc/install/include/petscistypes.h \ + /Users/adhikari/issm/trunk-jpl-sesaw/externalpackages/petsc/install/include/petscviewer.h \ + /Users/adhikari/issm/trunk-jpl-sesaw/externalpackages/petsc/install/include/petscdrawtypes.h \ + /Users/adhikari/issm/trunk-jpl-sesaw/externalpackages/petsc/install/include/petscksp.h \ + /Users/adhikari/issm/trunk-jpl-sesaw/externalpackages/petsc/install/include/petscpc.h \ + /Users/adhikari/issm/trunk-jpl-sesaw/externalpackages/petsc/install/include/petscpctypes.h \ + /Users/adhikari/issm/trunk-jpl-sesaw/externalpackages/petsc/install/include/petscdmtypes.h \ + modules/GiaIvinsDeflectionCorex/../../classes/../toolkits/./petsc/patches/petscpatches.h \ + modules/GiaIvinsDeflectionCorex/../../classes/../toolkits/./petsc/patches/./SolverEnum.h \ + modules/GiaIvinsDeflectionCorex/../../classes/../toolkits/./petsc/patches/../../toolkitsenums.h \ + modules/GiaIvinsDeflectionCorex/../../classes/../toolkits/./petsc/objects/petscobjects.h \ + modules/GiaIvinsDeflectionCorex/../../classes/../toolkits/./petsc/objects/./PetscMat.h \ + modules/GiaIvinsDeflectionCorex/../../classes/../toolkits/./petsc/objects/../petscincludes.h \ + modules/GiaIvinsDeflectionCorex/../../classes/../toolkits/./petsc/objects/./PetscVec.h \ + modules/GiaIvinsDeflectionCorex/../../classes/../toolkits/./petsc/objects/./PetscSolver.h \ + modules/GiaIvinsDeflectionCorex/../../classes/../toolkits/./metis/metisincludes.h \ + /Users/adhikari/issm/trunk-jpl-sesaw/externalpackages/petsc/install/include/metis.h \ + /usr/include/inttypes.h /usr/include/xlocale/_inttypes.h \ + modules/GiaIvinsDeflectionCorex/../../classes/../toolkits/./metis/patches/metispatches.h \ + modules/GiaIvinsDeflectionCorex/../../classes/../toolkits/./triangle/triangleincludes.h \ + modules/GiaIvinsDeflectionCorex/../../classes/../toolkits/./objects/toolkitobjects.h \ + modules/GiaIvinsDeflectionCorex/../../classes/../toolkits/./objects/./Vector.h \ + modules/GiaIvinsDeflectionCorex/../../classes/../toolkits/./objects/../issm/issmtoolkit.h \ + modules/GiaIvinsDeflectionCorex/../../classes/../toolkits/./objects/../issm/./IssmAbsMat.h \ + modules/GiaIvinsDeflectionCorex/../../classes/../toolkits/./objects/../issm/./IssmAbsVec.h \ + modules/GiaIvinsDeflectionCorex/../../classes/../toolkits/./objects/../issm/./IssmDenseMat.h \ + modules/GiaIvinsDeflectionCorex/../../classes/../toolkits/./objects/../issm/./IssmSeqVec.h \ + modules/GiaIvinsDeflectionCorex/../../classes/../toolkits/./objects/../issm/../gsl/gslincludes.h \ + modules/GiaIvinsDeflectionCorex/../../classes/../toolkits/./objects/../issm/./IssmMat.h \ + modules/GiaIvinsDeflectionCorex/../../classes/../toolkits/./objects/../issm/../ToolkitOptions.h \ + modules/GiaIvinsDeflectionCorex/../../classes/../toolkits/./objects/../issm/./IssmToolkitUtils.h \ + modules/GiaIvinsDeflectionCorex/../../classes/../toolkits/./objects/../issm/./IssmVec.h \ + modules/GiaIvinsDeflectionCorex/../../classes/../toolkits/./objects/../issm/./IssmSolver.h \ + modules/GiaIvinsDeflectionCorex/../../classes/../toolkits/./objects/../issm/./IssmMpiDenseMat.h \ + modules/GiaIvinsDeflectionCorex/../../classes/../toolkits/./objects/../issm/../mumps/mumpsincludes.h \ + modules/GiaIvinsDeflectionCorex/../../classes/../toolkits/./objects/../issm/./Bucket.h \ + modules/GiaIvinsDeflectionCorex/../../classes/../toolkits/./objects/../issm/./IssmMpiVec.h \ + modules/GiaIvinsDeflectionCorex/../../classes/../toolkits/./objects/../issm/./IssmMpiSparseMat.h \ + modules/GiaIvinsDeflectionCorex/../../classes/../toolkits/./objects/../issm/./SparseRow.h \ + modules/GiaIvinsDeflectionCorex/../../classes/../toolkits/./objects/./Matrix.h \ + modules/GiaIvinsDeflectionCorex/../../classes/../toolkits/./objects/./Solver.h \ + modules/GiaIvinsDeflectionCorex/../../classes/./Nodes.h \ + modules/GiaIvinsDeflectionCorex/../../classes/./Contours.h \ + modules/GiaIvinsDeflectionCorex/../../classes/./Node.h \ + modules/GiaIvinsDeflectionCorex/../../classes/./DofIndexing.h \ + modules/GiaIvinsDeflectionCorex/../../classes/./Update.h \ + modules/GiaIvinsDeflectionCorex/../../classes/./Profiler.h \ + modules/GiaIvinsDeflectionCorex/../../classes/./DependentObject.h \ + modules/GiaIvinsDeflectionCorex/../../classes/./Segment.h \ + modules/GiaIvinsDeflectionCorex/../../classes/./Massfluxatgate.h \ + modules/GiaIvinsDeflectionCorex/../../classes/./Definition.h \ + modules/GiaIvinsDeflectionCorex/../../classes/./Elements/Element.h \ + modules/GiaIvinsDeflectionCorex/../../classes/./Elements/Elements.h \ + modules/GiaIvinsDeflectionCorex/../../classes/./FemModel.h \ + modules/GiaIvinsDeflectionCorex/../../classes/./Misfit.h \ + modules/GiaIvinsDeflectionCorex/../../classes/../modules/SurfaceAreax/SurfaceAreax.h \ + modules/GiaIvinsDeflectionCorex/../../classes/../classes/Params/Parameters.h \ + modules/GiaIvinsDeflectionCorex/../../classes/../classes/Inputs/Input.h \ + modules/GiaIvinsDeflectionCorex/../../classes/../classes/gauss/Gauss.h \ + modules/GiaIvinsDeflectionCorex/../../classes/./Nodalvalue.h \ + modules/GiaIvinsDeflectionCorex/../../classes/./Masscon.h \ + modules/GiaIvinsDeflectionCorex/../../classes/./Massconaxpby.h \ + modules/GiaIvinsDeflectionCorex/../../classes/./Constraints/Constraints.h \ + modules/GiaIvinsDeflectionCorex/../../classes/./Constraints/Constraint.h \ + modules/GiaIvinsDeflectionCorex/../../classes/./Constraints/SpcStatic.h \ + modules/GiaIvinsDeflectionCorex/../../classes/./Constraints/SpcTransient.h \ + modules/GiaIvinsDeflectionCorex/../../classes/./Constraints/SpcDynamic.h \ + modules/GiaIvinsDeflectionCorex/../../classes/./Loads/Loads.h \ + modules/GiaIvinsDeflectionCorex/../../classes/./Loads/Load.h \ + modules/GiaIvinsDeflectionCorex/../../classes/./Loads/Friction.h \ + modules/GiaIvinsDeflectionCorex/../../classes/./Loads/Numericalflux.h \ + modules/GiaIvinsDeflectionCorex/../../classes/./Loads/Neumannflux.h \ + modules/GiaIvinsDeflectionCorex/../../classes/./Loads/Riftfront.h \ + modules/GiaIvinsDeflectionCorex/../../classes/./Loads/Penpair.h \ + modules/GiaIvinsDeflectionCorex/../../classes/./Loads/Pengrid.h \ + modules/GiaIvinsDeflectionCorex/../../classes/./Loads/Moulin.h \ + modules/GiaIvinsDeflectionCorex/../../classes/./Elements/Penta.h \ + modules/GiaIvinsDeflectionCorex/../../classes/./Elements/./ElementHook.h \ + modules/GiaIvinsDeflectionCorex/../../classes/./Elements/./PentaRef.h \ + modules/GiaIvinsDeflectionCorex/../../classes/./Elements/Seg.h \ + modules/GiaIvinsDeflectionCorex/../../classes/./Elements/./SegRef.h \ + modules/GiaIvinsDeflectionCorex/../../classes/./Elements/Tria.h \ + modules/GiaIvinsDeflectionCorex/../../classes/./Elements/./TriaRef.h \ + modules/GiaIvinsDeflectionCorex/../../classes/./Elements/Tetra.h \ + modules/GiaIvinsDeflectionCorex/../../classes/./Elements/./TetraRef.h \ + modules/GiaIvinsDeflectionCorex/../../classes/./Options/Option.h \ + modules/GiaIvinsDeflectionCorex/../../classes/./Options/Options.h \ + modules/GiaIvinsDeflectionCorex/../../classes/./Options/./GenericOption.h \ + modules/GiaIvinsDeflectionCorex/../../classes/./Options/./OptionUtilities.h \ + modules/GiaIvinsDeflectionCorex/../../classes/./Inputs/Inputs.h \ + modules/GiaIvinsDeflectionCorex/../../classes/./Inputs/BoolInput.h \ + modules/GiaIvinsDeflectionCorex/../../classes/./Inputs/DoubleInput.h \ + modules/GiaIvinsDeflectionCorex/../../classes/./Inputs/DoubleArrayInput.h \ + modules/GiaIvinsDeflectionCorex/../../classes/./Inputs/IntInput.h \ + modules/GiaIvinsDeflectionCorex/../../classes/./Inputs/TetraInput.h \ + modules/GiaIvinsDeflectionCorex/../../classes/./Inputs/PentaInput.h \ + modules/GiaIvinsDeflectionCorex/../../classes/./Inputs/TriaInput.h \ + modules/GiaIvinsDeflectionCorex/../../classes/./Inputs/SegInput.h \ + modules/GiaIvinsDeflectionCorex/../../classes/./Inputs/ControlInput.h \ + modules/GiaIvinsDeflectionCorex/../../classes/./Inputs/DatasetInput.h \ + modules/GiaIvinsDeflectionCorex/../../classes/./Inputs/TransientInput.h \ + modules/GiaIvinsDeflectionCorex/../../classes/./ExternalResults/Results.h \ + modules/GiaIvinsDeflectionCorex/../../classes/./ExternalResults/ExternalResult.h \ + modules/GiaIvinsDeflectionCorex/../../classes/./ExternalResults/GenericExternalResult.h \ + modules/GiaIvinsDeflectionCorex/../../classes/./Materials/Materials.h \ + modules/GiaIvinsDeflectionCorex/../../classes/./Materials/Material.h \ + modules/GiaIvinsDeflectionCorex/../../classes/./Materials/Matice.h \ + modules/GiaIvinsDeflectionCorex/../../classes/./Materials/../Hook.h \ + modules/GiaIvinsDeflectionCorex/../../classes/./Materials/Matestar.h \ + modules/GiaIvinsDeflectionCorex/../../classes/./Materials/Matpar.h \ + modules/GiaIvinsDeflectionCorex/../../classes/./Params/GenericParam.h \ + modules/GiaIvinsDeflectionCorex/../../classes/../classes/Params/./Param.h \ + modules/GiaIvinsDeflectionCorex/../../classes/./Params/BoolParam.h \ + modules/GiaIvinsDeflectionCorex/../../classes/./Params/DoubleMatParam.h \ + modules/GiaIvinsDeflectionCorex/../../classes/./Params/DoubleTransientMatParam.h \ + modules/GiaIvinsDeflectionCorex/../../classes/./Params/DoubleMatArrayParam.h \ + modules/GiaIvinsDeflectionCorex/../../classes/./Params/DoubleParam.h \ + modules/GiaIvinsDeflectionCorex/../../classes/./Params/DoubleVecParam.h \ + modules/GiaIvinsDeflectionCorex/../../classes/./Params/IntParam.h \ + modules/GiaIvinsDeflectionCorex/../../classes/./Params/IntVecParam.h \ + modules/GiaIvinsDeflectionCorex/../../classes/./Params/IntMatParam.h \ + modules/GiaIvinsDeflectionCorex/../../classes/./Params/FileParam.h \ + modules/GiaIvinsDeflectionCorex/../../classes/./Params/MatrixParam.h \ + modules/GiaIvinsDeflectionCorex/../../classes/./Params/VectorParam.h \ + modules/GiaIvinsDeflectionCorex/../../classes/./Params/StringArrayParam.h \ + modules/GiaIvinsDeflectionCorex/../../classes/./Params/StringParam.h \ + modules/GiaIvinsDeflectionCorex/../../classes/./Params/TransientParam.h \ + modules/GiaIvinsDeflectionCorex/../../classes/./Params/DataSetParam.h \ + modules/GiaIvinsDeflectionCorex/../../classes/./matrix/matrixobjects.h \ + modules/GiaIvinsDeflectionCorex/../../classes/./matrix/./ElementMatrix.h \ + modules/GiaIvinsDeflectionCorex/../../classes/./matrix/./ElementVector.h \ + modules/GiaIvinsDeflectionCorex/../../classes/./gauss/gaussobjects.h \ + modules/GiaIvinsDeflectionCorex/../../classes/../classes/gauss/./GaussSeg.h \ + modules/GiaIvinsDeflectionCorex/../../classes/../classes/gauss/./GaussTria.h \ + modules/GiaIvinsDeflectionCorex/../../classes/../classes/gauss/./GaussTetra.h \ + modules/GiaIvinsDeflectionCorex/../../classes/../classes/gauss/./GaussPenta.h \ + modules/GiaIvinsDeflectionCorex/../../classes/./kriging/krigingobjects.h \ + modules/GiaIvinsDeflectionCorex/../../classes/./kriging/./Variogram.h \ + modules/GiaIvinsDeflectionCorex/../../classes/./kriging/./GaussianVariogram.h \ + modules/GiaIvinsDeflectionCorex/../../classes/./kriging/./ExponentialVariogram.h \ + modules/GiaIvinsDeflectionCorex/../../classes/./kriging/./SphericalVariogram.h \ + modules/GiaIvinsDeflectionCorex/../../classes/./kriging/./PowerVariogram.h \ + modules/GiaIvinsDeflectionCorex/../../classes/./kriging/./Quadtree.h \ + modules/GiaIvinsDeflectionCorex/../../classes/./kriging/./Covertree.h \ + /Applications/Xcode.app/Contents/Developer/Toolchains/XcodeDefault.xctoolchain/usr/bin/../include/c++/v1/map \ + /Applications/Xcode.app/Contents/Developer/Toolchains/XcodeDefault.xctoolchain/usr/bin/../include/c++/v1/__tree \ + modules/GiaIvinsDeflectionCorex/../../classes/./kriging/./Observation.h \ + modules/GiaIvinsDeflectionCorex/../../classes/./kriging/./Observations.h \ + modules/GiaIvinsDeflectionCorex/../../classes/./Dakota/IssmDirectApplicInterface.h \ + modules/GiaIvinsDeflectionCorex/../../classes/./Dakota/IssmParallelDirectApplicInterface.h \ + modules/GiaIvinsDeflectionCorex/../../classes/./IoModel.h \ + modules/GiaIvinsDeflectionCorex/../../classes/./GiaIvinsDeflectionCoreArgs.h \ + modules/GiaIvinsDeflectionCorex/../../classes/./RiftStruct.h \ + modules/GiaIvinsDeflectionCorex/../InputUpdateFromConstantx/InputUpdateFromConstantx.h + +modules/GiaIvinsDeflectionCorex/./GiaIvinsDeflectionCorex.h: + +modules/GiaIvinsDeflectionCorex/../../classes/classes.h: + +modules/GiaIvinsDeflectionCorex/../../classes/./Contour.h: + +modules/GiaIvinsDeflectionCorex/../../classes/../shared/shared.h: + +modules/GiaIvinsDeflectionCorex/../../classes/../shared/./FSanalyticals/fsanalyticals.h: + +modules/GiaIvinsDeflectionCorex/../../classes/../shared/./FSanalyticals/../Numerics/types.h: + +../.././config.h: + +/usr/include/stdio.h: + +/usr/include/sys/cdefs.h: + +/usr/include/sys/_symbol_aliasing.h: + +/usr/include/sys/_posix_availability.h: + +/usr/include/Availability.h: + +/usr/include/AvailabilityInternal.h: + +/usr/include/_types.h: + +/usr/include/sys/_types.h: + +/usr/include/machine/_types.h: + +/usr/include/i386/_types.h: + +/usr/include/sys/_pthread/_pthread_types.h: + +/usr/include/sys/_types/_va_list.h: + +/usr/include/sys/_types/_size_t.h: + +/usr/include/sys/_types/_null.h: + +/usr/include/sys/stdio.h: + +/usr/include/sys/_types/_off_t.h: + +/usr/include/sys/_types/_ssize_t.h: + +modules/GiaIvinsDeflectionCorex/../../classes/../shared/./MemOps/MemOps.h: + +/Applications/Xcode.app/Contents/Developer/Toolchains/XcodeDefault.xctoolchain/usr/bin/../include/c++/v1/cassert: + +/Applications/Xcode.app/Contents/Developer/Toolchains/XcodeDefault.xctoolchain/usr/bin/../include/c++/v1/__config: + +/usr/include/unistd.h: + +/usr/include/sys/unistd.h: + +/usr/include/sys/_types/_posix_vdisable.h: + +/usr/include/sys/_types/_seek_set.h: + +/usr/include/_types/_uint64_t.h: + +/usr/include/sys/_types/_uid_t.h: + +/usr/include/sys/_types/_gid_t.h: + +/usr/include/sys/_types/_intptr_t.h: + +/usr/include/sys/_types/_pid_t.h: + +/usr/include/sys/_types/_useconds_t.h: + +/usr/include/sys/select.h: + +/usr/include/sys/appleapiopts.h: + +/usr/include/sys/_types/_fd_def.h: + +/usr/include/sys/_types/_timespec.h: + +/usr/include/sys/_types/_timeval.h: + +/usr/include/sys/_types/_time_t.h: + +/usr/include/sys/_types/_suseconds_t.h: + +/usr/include/sys/_types/_sigset_t.h: + +/usr/include/sys/_types/_fd_setsize.h: + +/usr/include/sys/_types/_fd_set.h: + +/usr/include/sys/_types/_fd_clr.h: + +/usr/include/sys/_types/_fd_isset.h: + +/usr/include/sys/_types/_fd_zero.h: + +/usr/include/sys/_types/_fd_copy.h: + +/usr/include/sys/_select.h: + +/usr/include/sys/_types/_dev_t.h: + +/usr/include/sys/_types/_mode_t.h: + +/usr/include/sys/_types/_uuid_t.h: + +/usr/include/gethostuuid.h: + +/usr/include/assert.h: + +/usr/include/stdlib.h: + +/usr/include/sys/wait.h: + +/usr/include/sys/_types/_id_t.h: + +/usr/include/sys/signal.h: + +/usr/include/machine/signal.h: + +/usr/include/i386/signal.h: + +/usr/include/machine/_mcontext.h: + +/usr/include/i386/_mcontext.h: + +/usr/include/mach/i386/_structs.h: + +/usr/include/sys/_pthread/_pthread_attr_t.h: + +/usr/include/sys/_types/_sigaltstack.h: + +/usr/include/sys/_types/_ucontext.h: + +/usr/include/sys/resource.h: + +/Applications/Xcode.app/Contents/Developer/Toolchains/XcodeDefault.xctoolchain/usr/bin/../lib/clang/7.0.2/include/stdint.h: + +/usr/include/stdint.h: + +/usr/include/sys/_types/_int8_t.h: + +/usr/include/sys/_types/_int16_t.h: + +/usr/include/sys/_types/_int32_t.h: + +/usr/include/sys/_types/_int64_t.h: + +/usr/include/_types/_uint8_t.h: + +/usr/include/_types/_uint16_t.h: + +/usr/include/_types/_uint32_t.h: + +/usr/include/sys/_types/_uintptr_t.h: + +/usr/include/_types/_intmax_t.h: + +/usr/include/_types/_uintmax_t.h: + +/usr/include/machine/endian.h: + +/usr/include/i386/endian.h: + +/usr/include/sys/_endian.h: + +/usr/include/libkern/_OSByteOrder.h: + +/usr/include/libkern/i386/_OSByteOrder.h: + +/usr/include/alloca.h: + +/usr/include/sys/_types/_ct_rune_t.h: + +/usr/include/sys/_types/_rune_t.h: + +/usr/include/sys/_types/_wchar_t.h: + +/usr/include/machine/types.h: + +/usr/include/i386/types.h: + +modules/GiaIvinsDeflectionCorex/../../classes/../shared/./Bamg/shared.h: + +modules/GiaIvinsDeflectionCorex/../../classes/../shared/./Bamg/Abs.h: + +modules/GiaIvinsDeflectionCorex/../../classes/../shared/./Bamg/BigPrimeNumber.h: + +modules/GiaIvinsDeflectionCorex/../../classes/../shared/./Bamg/BinaryRand.h: + +/Applications/Xcode.app/Contents/Developer/Toolchains/XcodeDefault.xctoolchain/usr/bin/../include/c++/v1/cstdlib: + +modules/GiaIvinsDeflectionCorex/../../classes/../shared/./Bamg/Exchange.h: + +modules/GiaIvinsDeflectionCorex/../../classes/../shared/./Bamg/extrema.h: + +modules/GiaIvinsDeflectionCorex/../../classes/../shared/./Bamg/HeapSort.h: + +modules/GiaIvinsDeflectionCorex/../../classes/../shared/./Bamg/OppositeAngle.h: + +modules/GiaIvinsDeflectionCorex/../../classes/../shared/./Bamg/../Numerics/constants.h: + +modules/GiaIvinsDeflectionCorex/../../classes/../shared/./Elements/elements.h: + +modules/GiaIvinsDeflectionCorex/../../classes/../shared/./Enum/Enum.h: + +modules/GiaIvinsDeflectionCorex/../../classes/../shared/./Enum/./EnumDefinitions.h: + +modules/GiaIvinsDeflectionCorex/../../classes/../shared/./Exceptions/exceptions.h: + +/Applications/Xcode.app/Contents/Developer/Toolchains/XcodeDefault.xctoolchain/usr/bin/../include/c++/v1/exception: + +/Applications/Xcode.app/Contents/Developer/Toolchains/XcodeDefault.xctoolchain/usr/bin/../include/c++/v1/cstddef: + +/Applications/Xcode.app/Contents/Developer/Toolchains/XcodeDefault.xctoolchain/usr/bin/../lib/clang/7.0.2/include/stddef.h: + +/Applications/Xcode.app/Contents/Developer/Toolchains/XcodeDefault.xctoolchain/usr/bin/../include/c++/v1/type_traits: + +/Applications/Xcode.app/Contents/Developer/Toolchains/XcodeDefault.xctoolchain/usr/bin/../include/c++/v1/string: + +/Applications/Xcode.app/Contents/Developer/Toolchains/XcodeDefault.xctoolchain/usr/bin/../include/c++/v1/iosfwd: + +/usr/include/wchar.h: + +/usr/include/sys/_types/_mbstate_t.h: + +/Applications/Xcode.app/Contents/Developer/Toolchains/XcodeDefault.xctoolchain/usr/bin/../lib/clang/7.0.2/include/stdarg.h: + +/usr/include/time.h: + +/usr/include/sys/_types/_clock_t.h: + +/usr/include/_wctype.h: + +/usr/include/sys/_types/_wint_t.h: + +/usr/include/_types/_wctype_t.h: + +/usr/include/ctype.h: + +/usr/include/runetype.h: + +/Applications/Xcode.app/Contents/Developer/Toolchains/XcodeDefault.xctoolchain/usr/bin/../include/c++/v1/cstring: + +/usr/include/string.h: + +/usr/include/sys/_types/_rsize_t.h: + +/usr/include/sys/_types/_errno_t.h: + +/usr/include/strings.h: + +/Applications/Xcode.app/Contents/Developer/Toolchains/XcodeDefault.xctoolchain/usr/bin/../include/c++/v1/cstdio: + +/Applications/Xcode.app/Contents/Developer/Toolchains/XcodeDefault.xctoolchain/usr/bin/../include/c++/v1/cwchar: + +/Applications/Xcode.app/Contents/Developer/Toolchains/XcodeDefault.xctoolchain/usr/bin/../include/c++/v1/cwctype: + +/Applications/Xcode.app/Contents/Developer/Toolchains/XcodeDefault.xctoolchain/usr/bin/../include/c++/v1/cctype: + +/usr/include/wctype.h: + +/usr/include/_types/_wctrans_t.h: + +/Applications/Xcode.app/Contents/Developer/Toolchains/XcodeDefault.xctoolchain/usr/bin/../include/c++/v1/algorithm: + +/Applications/Xcode.app/Contents/Developer/Toolchains/XcodeDefault.xctoolchain/usr/bin/../include/c++/v1/initializer_list: + +/Applications/Xcode.app/Contents/Developer/Toolchains/XcodeDefault.xctoolchain/usr/bin/../include/c++/v1/utility: + +/Applications/Xcode.app/Contents/Developer/Toolchains/XcodeDefault.xctoolchain/usr/bin/../include/c++/v1/__tuple: + +/Applications/Xcode.app/Contents/Developer/Toolchains/XcodeDefault.xctoolchain/usr/bin/../include/c++/v1/__tuple_03: + +/Applications/Xcode.app/Contents/Developer/Toolchains/XcodeDefault.xctoolchain/usr/bin/../include/c++/v1/memory: + +/Applications/Xcode.app/Contents/Developer/Toolchains/XcodeDefault.xctoolchain/usr/bin/../include/c++/v1/typeinfo: + +/Applications/Xcode.app/Contents/Developer/Toolchains/XcodeDefault.xctoolchain/usr/bin/../include/c++/v1/cstdint: + +/Applications/Xcode.app/Contents/Developer/Toolchains/XcodeDefault.xctoolchain/usr/bin/../include/c++/v1/new: + +/Applications/Xcode.app/Contents/Developer/Toolchains/XcodeDefault.xctoolchain/usr/bin/../include/c++/v1/limits: + +/Applications/Xcode.app/Contents/Developer/Toolchains/XcodeDefault.xctoolchain/usr/bin/../include/c++/v1/__undef_min_max: + +/Applications/Xcode.app/Contents/Developer/Toolchains/XcodeDefault.xctoolchain/usr/bin/../include/c++/v1/iterator: + +/Applications/Xcode.app/Contents/Developer/Toolchains/XcodeDefault.xctoolchain/usr/bin/../include/c++/v1/__functional_base: + +/Applications/Xcode.app/Contents/Developer/Toolchains/XcodeDefault.xctoolchain/usr/bin/../include/c++/v1/__functional_base_03: + +/Applications/Xcode.app/Contents/Developer/Toolchains/XcodeDefault.xctoolchain/usr/bin/../include/c++/v1/__debug: + +/Applications/Xcode.app/Contents/Developer/Toolchains/XcodeDefault.xctoolchain/usr/bin/../include/c++/v1/tuple: + +/Applications/Xcode.app/Contents/Developer/Toolchains/XcodeDefault.xctoolchain/usr/bin/../include/c++/v1/stdexcept: + +/Applications/Xcode.app/Contents/Developer/Toolchains/XcodeDefault.xctoolchain/usr/bin/../include/c++/v1/iostream: + +/Applications/Xcode.app/Contents/Developer/Toolchains/XcodeDefault.xctoolchain/usr/bin/../include/c++/v1/ios: + +/Applications/Xcode.app/Contents/Developer/Toolchains/XcodeDefault.xctoolchain/usr/bin/../include/c++/v1/__locale: + +/Applications/Xcode.app/Contents/Developer/Toolchains/XcodeDefault.xctoolchain/usr/bin/../include/c++/v1/mutex: + +/Applications/Xcode.app/Contents/Developer/Toolchains/XcodeDefault.xctoolchain/usr/bin/../include/c++/v1/__mutex_base: + +/Applications/Xcode.app/Contents/Developer/Toolchains/XcodeDefault.xctoolchain/usr/bin/../include/c++/v1/chrono: + +/Applications/Xcode.app/Contents/Developer/Toolchains/XcodeDefault.xctoolchain/usr/bin/../include/c++/v1/ctime: + +/Applications/Xcode.app/Contents/Developer/Toolchains/XcodeDefault.xctoolchain/usr/bin/../include/c++/v1/ratio: + +/Applications/Xcode.app/Contents/Developer/Toolchains/XcodeDefault.xctoolchain/usr/bin/../include/c++/v1/climits: + +/Applications/Xcode.app/Contents/Developer/Toolchains/XcodeDefault.xctoolchain/usr/bin/../lib/clang/7.0.2/include/limits.h: + +/usr/include/limits.h: + +/usr/include/machine/limits.h: + +/usr/include/i386/limits.h: + +/usr/include/i386/_limits.h: + +/usr/include/sys/syslimits.h: + +/Applications/Xcode.app/Contents/Developer/Toolchains/XcodeDefault.xctoolchain/usr/bin/../include/c++/v1/system_error: + +/Applications/Xcode.app/Contents/Developer/Toolchains/XcodeDefault.xctoolchain/usr/bin/../include/c++/v1/cerrno: + +/usr/include/errno.h: + +/usr/include/sys/errno.h: + +/usr/include/pthread.h: + +/usr/include/pthread/pthread_impl.h: + +/usr/include/pthread/sched.h: + +/usr/include/sys/_pthread/_pthread_cond_t.h: + +/usr/include/sys/_pthread/_pthread_condattr_t.h: + +/usr/include/sys/_pthread/_pthread_key_t.h: + +/usr/include/sys/_pthread/_pthread_mutex_t.h: + +/usr/include/sys/_pthread/_pthread_mutexattr_t.h: + +/usr/include/sys/_pthread/_pthread_once_t.h: + +/usr/include/sys/_pthread/_pthread_rwlock_t.h: + +/usr/include/sys/_pthread/_pthread_rwlockattr_t.h: + +/usr/include/sys/_pthread/_pthread_t.h: + +/usr/include/pthread/qos.h: + +/usr/include/sys/qos.h: + +/usr/include/sys/_types/_mach_port_t.h: + +/Applications/Xcode.app/Contents/Developer/Toolchains/XcodeDefault.xctoolchain/usr/bin/../include/c++/v1/functional: + +/Applications/Xcode.app/Contents/Developer/Toolchains/XcodeDefault.xctoolchain/usr/bin/../include/c++/v1/__functional_03: + +/usr/include/locale.h: + +/usr/include/_locale.h: + +/usr/include/xlocale.h: + +/usr/include/_xlocale.h: + +/usr/include/xlocale/_ctype.h: + +/usr/include/xlocale/__wctype.h: + +/usr/include/xlocale/_stdio.h: + +/usr/include/xlocale/_stdlib.h: + +/usr/include/xlocale/_string.h: + +/usr/include/xlocale/_time.h: + +/usr/include/xlocale/_wchar.h: + +/usr/include/xlocale/_wctype.h: + +/Applications/Xcode.app/Contents/Developer/Toolchains/XcodeDefault.xctoolchain/usr/bin/../include/c++/v1/streambuf: + +/Applications/Xcode.app/Contents/Developer/Toolchains/XcodeDefault.xctoolchain/usr/bin/../include/c++/v1/istream: + +/Applications/Xcode.app/Contents/Developer/Toolchains/XcodeDefault.xctoolchain/usr/bin/../include/c++/v1/ostream: + +/Applications/Xcode.app/Contents/Developer/Toolchains/XcodeDefault.xctoolchain/usr/bin/../include/c++/v1/locale: + +/usr/include/nl_types.h: + +/usr/include/sys/types.h: + +/usr/include/sys/_types/_blkcnt_t.h: + +/usr/include/sys/_types/_blksize_t.h: + +/usr/include/sys/_types/_in_addr_t.h: + +/usr/include/sys/_types/_in_port_t.h: + +/usr/include/sys/_types/_ino_t.h: + +/usr/include/sys/_types/_ino64_t.h: + +/usr/include/sys/_types/_key_t.h: + +/usr/include/sys/_types/_nlink_t.h: + +/usr/include/sys/_types/_fsblkcnt_t.h: + +/usr/include/sys/_types/_fsfilcnt_t.h: + +/usr/include/_types/_nl_item.h: + +/Applications/Xcode.app/Contents/Developer/Toolchains/XcodeDefault.xctoolchain/usr/bin/../include/c++/v1/bitset: + +/Applications/Xcode.app/Contents/Developer/Toolchains/XcodeDefault.xctoolchain/usr/bin/../include/c++/v1/__bit_reference: + +/Applications/Xcode.app/Contents/Developer/Toolchains/XcodeDefault.xctoolchain/usr/bin/../include/c++/v1/sstream: + +/Applications/Xcode.app/Contents/Developer/Toolchains/XcodeDefault.xctoolchain/usr/bin/../include/c++/v1/iomanip: + +modules/GiaIvinsDeflectionCorex/../../classes/../shared/./io/io.h: + +modules/GiaIvinsDeflectionCorex/../../classes/../shared/./io/./Comm/IssmComm.h: + +modules/GiaIvinsDeflectionCorex/../../classes/../shared/./io/./Comm/../../../toolkits/mpi/issmmpi.h: + +/Users/adhikari/issm/trunk-jpl-sesaw/externalpackages/mpich/install/include/mpi.h: + +/Users/adhikari/issm/trunk-jpl-sesaw/externalpackages/mpich/install/include/mpio.h: + +/Users/adhikari/issm/trunk-jpl-sesaw/externalpackages/mpich/install/include/mpicxx.h: + +modules/GiaIvinsDeflectionCorex/../../classes/../shared/./io/./Comm/../../../toolkits/mpi/./commops/commops.h: + +modules/GiaIvinsDeflectionCorex/../../classes/../shared/./io/./Comm/../../../toolkits/mpi/./commops/../../../shared/io/Comm/IssmComm.h: + +modules/GiaIvinsDeflectionCorex/../../classes/../shared/./io/./Comm/../../../toolkits/mpi/./commops/../../mpi/issmmpi.h: + +modules/GiaIvinsDeflectionCorex/../../classes/../shared/./io/./Disk/diskio.h: + +modules/GiaIvinsDeflectionCorex/../../classes/../shared/./io/./Print/Print.h: + +modules/GiaIvinsDeflectionCorex/../../classes/../shared/./io/./Marshalling/Marshalling.h: + +modules/GiaIvinsDeflectionCorex/../../classes/../shared/./io/./Marshalling/IoCodeConversions.h: + +modules/GiaIvinsDeflectionCorex/../../classes/../shared/./Exp/exp.h: + +modules/GiaIvinsDeflectionCorex/../../classes/../shared/./Exp/../Numerics/recast.h: + +modules/GiaIvinsDeflectionCorex/../../classes/../shared/./Matrix/matrix.h: + +modules/GiaIvinsDeflectionCorex/../../classes/../shared/./Numerics/numerics.h: + +modules/GiaIvinsDeflectionCorex/../../classes/../shared/./FSanalyticals/../Numerics/./Verbosity.h: + +modules/GiaIvinsDeflectionCorex/../../classes/../shared/./FSanalyticals/../Numerics/./GaussPoints.h: + +modules/GiaIvinsDeflectionCorex/../../classes/../shared/./FSanalyticals/../Numerics/./isnan.h: + +/Applications/Xcode.app/Contents/Developer/Toolchains/XcodeDefault.xctoolchain/usr/bin/../include/c++/v1/cmath: + +/usr/include/math.h: + +modules/GiaIvinsDeflectionCorex/../../classes/../shared/./FSanalyticals/../Numerics/./OptPars.h: + +modules/GiaIvinsDeflectionCorex/../../classes/../shared/./Sorting/sorting.h: + +modules/GiaIvinsDeflectionCorex/../../classes/../shared/./String/sharedstring.h: + +modules/GiaIvinsDeflectionCorex/../../classes/../shared/./Threads/issm_threads.h: + +modules/GiaIvinsDeflectionCorex/../../classes/../shared/./TriMesh/trimesh.h: + +modules/GiaIvinsDeflectionCorex/../../classes/../shared/./LatLong/latlong.h: + +modules/GiaIvinsDeflectionCorex/../../classes/../datastructures/datastructures.h: + +modules/GiaIvinsDeflectionCorex/../../classes/../datastructures/./DataSet.h: + +/Applications/Xcode.app/Contents/Developer/Toolchains/XcodeDefault.xctoolchain/usr/bin/../include/c++/v1/vector: + +/Applications/Xcode.app/Contents/Developer/Toolchains/XcodeDefault.xctoolchain/usr/bin/../include/c++/v1/__split_buffer: + +modules/GiaIvinsDeflectionCorex/../../classes/../datastructures/./Object.h: + +modules/GiaIvinsDeflectionCorex/../../classes/./Vertices.h: + +modules/GiaIvinsDeflectionCorex/../../classes/./Vertex.h: + +modules/GiaIvinsDeflectionCorex/../../classes/./classes.h: + +modules/GiaIvinsDeflectionCorex/../../classes/../toolkits/toolkits.h: + +modules/GiaIvinsDeflectionCorex/../../classes/../toolkits/./petsc/petscincludes.h: + +/Users/adhikari/issm/trunk-jpl-sesaw/externalpackages/petsc/install/include/petscmat.h: + +/Users/adhikari/issm/trunk-jpl-sesaw/externalpackages/petsc/install/include/petscvec.h: + +/Users/adhikari/issm/trunk-jpl-sesaw/externalpackages/petsc/install/include/petscis.h: + +/Users/adhikari/issm/trunk-jpl-sesaw/externalpackages/petsc/install/include/petscsys.h: + +/Users/adhikari/issm/trunk-jpl-sesaw/externalpackages/petsc/install/include/petscconf.h: + +/Users/adhikari/issm/trunk-jpl-sesaw/externalpackages/petsc/install/include/petscfix.h: + +/Users/adhikari/issm/trunk-jpl-sesaw/externalpackages/petsc/install/include/petscversion.h: + +/Users/adhikari/issm/trunk-jpl-sesaw/externalpackages/petsc/install/include/petscmath.h: + +/Applications/Xcode.app/Contents/Developer/Toolchains/XcodeDefault.xctoolchain/usr/bin/../include/c++/v1/complex: + +/Users/adhikari/issm/trunk-jpl-sesaw/externalpackages/petsc/install/include/petscerror.h: + +/Users/adhikari/issm/trunk-jpl-sesaw/externalpackages/petsc/install/include/petscviewertypes.h: + +/Users/adhikari/issm/trunk-jpl-sesaw/externalpackages/petsc/install/include/petscoptions.h: + +/Users/adhikari/issm/trunk-jpl-sesaw/externalpackages/petsc/install/include/petsclog.h: + +/Applications/Xcode.app/Contents/Developer/Toolchains/XcodeDefault.xctoolchain/usr/bin/../lib/clang/7.0.2/include/xmmintrin.h: + +/Applications/Xcode.app/Contents/Developer/Toolchains/XcodeDefault.xctoolchain/usr/bin/../lib/clang/7.0.2/include/mmintrin.h: + +/Applications/Xcode.app/Contents/Developer/Toolchains/XcodeDefault.xctoolchain/usr/bin/../lib/clang/7.0.2/include/mm_malloc.h: + +/Applications/Xcode.app/Contents/Developer/Toolchains/XcodeDefault.xctoolchain/usr/bin/../lib/clang/7.0.2/include/emmintrin.h: + +/usr/include/sys/param.h: + +/usr/include/machine/param.h: + +/usr/include/i386/param.h: + +/usr/include/i386/_param.h: + +/Users/adhikari/issm/trunk-jpl-sesaw/externalpackages/petsc/install/include/petscsftypes.h: + +/Users/adhikari/issm/trunk-jpl-sesaw/externalpackages/petsc/install/include/petscistypes.h: + +/Users/adhikari/issm/trunk-jpl-sesaw/externalpackages/petsc/install/include/petscviewer.h: + +/Users/adhikari/issm/trunk-jpl-sesaw/externalpackages/petsc/install/include/petscdrawtypes.h: + +/Users/adhikari/issm/trunk-jpl-sesaw/externalpackages/petsc/install/include/petscksp.h: + +/Users/adhikari/issm/trunk-jpl-sesaw/externalpackages/petsc/install/include/petscpc.h: + +/Users/adhikari/issm/trunk-jpl-sesaw/externalpackages/petsc/install/include/petscpctypes.h: + +/Users/adhikari/issm/trunk-jpl-sesaw/externalpackages/petsc/install/include/petscdmtypes.h: + +modules/GiaIvinsDeflectionCorex/../../classes/../toolkits/./petsc/patches/petscpatches.h: + +modules/GiaIvinsDeflectionCorex/../../classes/../toolkits/./petsc/patches/./SolverEnum.h: + +modules/GiaIvinsDeflectionCorex/../../classes/../toolkits/./petsc/patches/../../toolkitsenums.h: + +modules/GiaIvinsDeflectionCorex/../../classes/../toolkits/./petsc/objects/petscobjects.h: + +modules/GiaIvinsDeflectionCorex/../../classes/../toolkits/./petsc/objects/./PetscMat.h: + +modules/GiaIvinsDeflectionCorex/../../classes/../toolkits/./petsc/objects/../petscincludes.h: + +modules/GiaIvinsDeflectionCorex/../../classes/../toolkits/./petsc/objects/./PetscVec.h: + +modules/GiaIvinsDeflectionCorex/../../classes/../toolkits/./petsc/objects/./PetscSolver.h: + +modules/GiaIvinsDeflectionCorex/../../classes/../toolkits/./metis/metisincludes.h: + +/Users/adhikari/issm/trunk-jpl-sesaw/externalpackages/petsc/install/include/metis.h: + +/usr/include/inttypes.h: + +/usr/include/xlocale/_inttypes.h: + +modules/GiaIvinsDeflectionCorex/../../classes/../toolkits/./metis/patches/metispatches.h: + +modules/GiaIvinsDeflectionCorex/../../classes/../toolkits/./triangle/triangleincludes.h: + +modules/GiaIvinsDeflectionCorex/../../classes/../toolkits/./objects/toolkitobjects.h: + +modules/GiaIvinsDeflectionCorex/../../classes/../toolkits/./objects/./Vector.h: + +modules/GiaIvinsDeflectionCorex/../../classes/../toolkits/./objects/../issm/issmtoolkit.h: + +modules/GiaIvinsDeflectionCorex/../../classes/../toolkits/./objects/../issm/./IssmAbsMat.h: + +modules/GiaIvinsDeflectionCorex/../../classes/../toolkits/./objects/../issm/./IssmAbsVec.h: + +modules/GiaIvinsDeflectionCorex/../../classes/../toolkits/./objects/../issm/./IssmDenseMat.h: + +modules/GiaIvinsDeflectionCorex/../../classes/../toolkits/./objects/../issm/./IssmSeqVec.h: + +modules/GiaIvinsDeflectionCorex/../../classes/../toolkits/./objects/../issm/../gsl/gslincludes.h: + +modules/GiaIvinsDeflectionCorex/../../classes/../toolkits/./objects/../issm/./IssmMat.h: + +modules/GiaIvinsDeflectionCorex/../../classes/../toolkits/./objects/../issm/../ToolkitOptions.h: + +modules/GiaIvinsDeflectionCorex/../../classes/../toolkits/./objects/../issm/./IssmToolkitUtils.h: + +modules/GiaIvinsDeflectionCorex/../../classes/../toolkits/./objects/../issm/./IssmVec.h: + +modules/GiaIvinsDeflectionCorex/../../classes/../toolkits/./objects/../issm/./IssmSolver.h: + +modules/GiaIvinsDeflectionCorex/../../classes/../toolkits/./objects/../issm/./IssmMpiDenseMat.h: + +modules/GiaIvinsDeflectionCorex/../../classes/../toolkits/./objects/../issm/../mumps/mumpsincludes.h: + +modules/GiaIvinsDeflectionCorex/../../classes/../toolkits/./objects/../issm/./Bucket.h: + +modules/GiaIvinsDeflectionCorex/../../classes/../toolkits/./objects/../issm/./IssmMpiVec.h: + +modules/GiaIvinsDeflectionCorex/../../classes/../toolkits/./objects/../issm/./IssmMpiSparseMat.h: + +modules/GiaIvinsDeflectionCorex/../../classes/../toolkits/./objects/../issm/./SparseRow.h: + +modules/GiaIvinsDeflectionCorex/../../classes/../toolkits/./objects/./Matrix.h: + +modules/GiaIvinsDeflectionCorex/../../classes/../toolkits/./objects/./Solver.h: + +modules/GiaIvinsDeflectionCorex/../../classes/./Nodes.h: + +modules/GiaIvinsDeflectionCorex/../../classes/./Contours.h: + +modules/GiaIvinsDeflectionCorex/../../classes/./Node.h: + +modules/GiaIvinsDeflectionCorex/../../classes/./DofIndexing.h: + +modules/GiaIvinsDeflectionCorex/../../classes/./Update.h: + +modules/GiaIvinsDeflectionCorex/../../classes/./Profiler.h: + +modules/GiaIvinsDeflectionCorex/../../classes/./DependentObject.h: + +modules/GiaIvinsDeflectionCorex/../../classes/./Segment.h: + +modules/GiaIvinsDeflectionCorex/../../classes/./Massfluxatgate.h: + +modules/GiaIvinsDeflectionCorex/../../classes/./Definition.h: + +modules/GiaIvinsDeflectionCorex/../../classes/./Elements/Element.h: + +modules/GiaIvinsDeflectionCorex/../../classes/./Elements/Elements.h: + +modules/GiaIvinsDeflectionCorex/../../classes/./FemModel.h: + +modules/GiaIvinsDeflectionCorex/../../classes/./Misfit.h: + +modules/GiaIvinsDeflectionCorex/../../classes/../modules/SurfaceAreax/SurfaceAreax.h: + +modules/GiaIvinsDeflectionCorex/../../classes/../classes/Params/Parameters.h: + +modules/GiaIvinsDeflectionCorex/../../classes/../classes/Inputs/Input.h: + +modules/GiaIvinsDeflectionCorex/../../classes/../classes/gauss/Gauss.h: + +modules/GiaIvinsDeflectionCorex/../../classes/./Nodalvalue.h: + +modules/GiaIvinsDeflectionCorex/../../classes/./Masscon.h: + +modules/GiaIvinsDeflectionCorex/../../classes/./Massconaxpby.h: + +modules/GiaIvinsDeflectionCorex/../../classes/./Constraints/Constraints.h: + +modules/GiaIvinsDeflectionCorex/../../classes/./Constraints/Constraint.h: + +modules/GiaIvinsDeflectionCorex/../../classes/./Constraints/SpcStatic.h: + +modules/GiaIvinsDeflectionCorex/../../classes/./Constraints/SpcTransient.h: + +modules/GiaIvinsDeflectionCorex/../../classes/./Constraints/SpcDynamic.h: + +modules/GiaIvinsDeflectionCorex/../../classes/./Loads/Loads.h: + +modules/GiaIvinsDeflectionCorex/../../classes/./Loads/Load.h: + +modules/GiaIvinsDeflectionCorex/../../classes/./Loads/Friction.h: + +modules/GiaIvinsDeflectionCorex/../../classes/./Loads/Numericalflux.h: + +modules/GiaIvinsDeflectionCorex/../../classes/./Loads/Neumannflux.h: + +modules/GiaIvinsDeflectionCorex/../../classes/./Loads/Riftfront.h: + +modules/GiaIvinsDeflectionCorex/../../classes/./Loads/Penpair.h: + +modules/GiaIvinsDeflectionCorex/../../classes/./Loads/Pengrid.h: + +modules/GiaIvinsDeflectionCorex/../../classes/./Loads/Moulin.h: + +modules/GiaIvinsDeflectionCorex/../../classes/./Elements/Penta.h: + +modules/GiaIvinsDeflectionCorex/../../classes/./Elements/./ElementHook.h: + +modules/GiaIvinsDeflectionCorex/../../classes/./Elements/./PentaRef.h: + +modules/GiaIvinsDeflectionCorex/../../classes/./Elements/Seg.h: + +modules/GiaIvinsDeflectionCorex/../../classes/./Elements/./SegRef.h: + +modules/GiaIvinsDeflectionCorex/../../classes/./Elements/Tria.h: + +modules/GiaIvinsDeflectionCorex/../../classes/./Elements/./TriaRef.h: + +modules/GiaIvinsDeflectionCorex/../../classes/./Elements/Tetra.h: + +modules/GiaIvinsDeflectionCorex/../../classes/./Elements/./TetraRef.h: + +modules/GiaIvinsDeflectionCorex/../../classes/./Options/Option.h: + +modules/GiaIvinsDeflectionCorex/../../classes/./Options/Options.h: + +modules/GiaIvinsDeflectionCorex/../../classes/./Options/./GenericOption.h: + +modules/GiaIvinsDeflectionCorex/../../classes/./Options/./OptionUtilities.h: + +modules/GiaIvinsDeflectionCorex/../../classes/./Inputs/Inputs.h: + +modules/GiaIvinsDeflectionCorex/../../classes/./Inputs/BoolInput.h: + +modules/GiaIvinsDeflectionCorex/../../classes/./Inputs/DoubleInput.h: + +modules/GiaIvinsDeflectionCorex/../../classes/./Inputs/DoubleArrayInput.h: + +modules/GiaIvinsDeflectionCorex/../../classes/./Inputs/IntInput.h: + +modules/GiaIvinsDeflectionCorex/../../classes/./Inputs/TetraInput.h: + +modules/GiaIvinsDeflectionCorex/../../classes/./Inputs/PentaInput.h: + +modules/GiaIvinsDeflectionCorex/../../classes/./Inputs/TriaInput.h: + +modules/GiaIvinsDeflectionCorex/../../classes/./Inputs/SegInput.h: + +modules/GiaIvinsDeflectionCorex/../../classes/./Inputs/ControlInput.h: + +modules/GiaIvinsDeflectionCorex/../../classes/./Inputs/DatasetInput.h: + +modules/GiaIvinsDeflectionCorex/../../classes/./Inputs/TransientInput.h: + +modules/GiaIvinsDeflectionCorex/../../classes/./ExternalResults/Results.h: + +modules/GiaIvinsDeflectionCorex/../../classes/./ExternalResults/ExternalResult.h: + +modules/GiaIvinsDeflectionCorex/../../classes/./ExternalResults/GenericExternalResult.h: + +modules/GiaIvinsDeflectionCorex/../../classes/./Materials/Materials.h: + +modules/GiaIvinsDeflectionCorex/../../classes/./Materials/Material.h: + +modules/GiaIvinsDeflectionCorex/../../classes/./Materials/Matice.h: + +modules/GiaIvinsDeflectionCorex/../../classes/./Materials/../Hook.h: + +modules/GiaIvinsDeflectionCorex/../../classes/./Materials/Matestar.h: + +modules/GiaIvinsDeflectionCorex/../../classes/./Materials/Matpar.h: + +modules/GiaIvinsDeflectionCorex/../../classes/./Params/GenericParam.h: + +modules/GiaIvinsDeflectionCorex/../../classes/../classes/Params/./Param.h: + +modules/GiaIvinsDeflectionCorex/../../classes/./Params/BoolParam.h: + +modules/GiaIvinsDeflectionCorex/../../classes/./Params/DoubleMatParam.h: + +modules/GiaIvinsDeflectionCorex/../../classes/./Params/DoubleTransientMatParam.h: + +modules/GiaIvinsDeflectionCorex/../../classes/./Params/DoubleMatArrayParam.h: + +modules/GiaIvinsDeflectionCorex/../../classes/./Params/DoubleParam.h: + +modules/GiaIvinsDeflectionCorex/../../classes/./Params/DoubleVecParam.h: + +modules/GiaIvinsDeflectionCorex/../../classes/./Params/IntParam.h: + +modules/GiaIvinsDeflectionCorex/../../classes/./Params/IntVecParam.h: + +modules/GiaIvinsDeflectionCorex/../../classes/./Params/IntMatParam.h: + +modules/GiaIvinsDeflectionCorex/../../classes/./Params/FileParam.h: + +modules/GiaIvinsDeflectionCorex/../../classes/./Params/MatrixParam.h: + +modules/GiaIvinsDeflectionCorex/../../classes/./Params/VectorParam.h: + +modules/GiaIvinsDeflectionCorex/../../classes/./Params/StringArrayParam.h: + +modules/GiaIvinsDeflectionCorex/../../classes/./Params/StringParam.h: + +modules/GiaIvinsDeflectionCorex/../../classes/./Params/TransientParam.h: + +modules/GiaIvinsDeflectionCorex/../../classes/./Params/DataSetParam.h: + +modules/GiaIvinsDeflectionCorex/../../classes/./matrix/matrixobjects.h: + +modules/GiaIvinsDeflectionCorex/../../classes/./matrix/./ElementMatrix.h: + +modules/GiaIvinsDeflectionCorex/../../classes/./matrix/./ElementVector.h: + +modules/GiaIvinsDeflectionCorex/../../classes/./gauss/gaussobjects.h: + +modules/GiaIvinsDeflectionCorex/../../classes/../classes/gauss/./GaussSeg.h: + +modules/GiaIvinsDeflectionCorex/../../classes/../classes/gauss/./GaussTria.h: + +modules/GiaIvinsDeflectionCorex/../../classes/../classes/gauss/./GaussTetra.h: + +modules/GiaIvinsDeflectionCorex/../../classes/../classes/gauss/./GaussPenta.h: + +modules/GiaIvinsDeflectionCorex/../../classes/./kriging/krigingobjects.h: + +modules/GiaIvinsDeflectionCorex/../../classes/./kriging/./Variogram.h: + +modules/GiaIvinsDeflectionCorex/../../classes/./kriging/./GaussianVariogram.h: + +modules/GiaIvinsDeflectionCorex/../../classes/./kriging/./ExponentialVariogram.h: + +modules/GiaIvinsDeflectionCorex/../../classes/./kriging/./SphericalVariogram.h: + +modules/GiaIvinsDeflectionCorex/../../classes/./kriging/./PowerVariogram.h: + +modules/GiaIvinsDeflectionCorex/../../classes/./kriging/./Quadtree.h: + +modules/GiaIvinsDeflectionCorex/../../classes/./kriging/./Covertree.h: + +/Applications/Xcode.app/Contents/Developer/Toolchains/XcodeDefault.xctoolchain/usr/bin/../include/c++/v1/map: + +/Applications/Xcode.app/Contents/Developer/Toolchains/XcodeDefault.xctoolchain/usr/bin/../include/c++/v1/__tree: + +modules/GiaIvinsDeflectionCorex/../../classes/./kriging/./Observation.h: + +modules/GiaIvinsDeflectionCorex/../../classes/./kriging/./Observations.h: + +modules/GiaIvinsDeflectionCorex/../../classes/./Dakota/IssmDirectApplicInterface.h: + +modules/GiaIvinsDeflectionCorex/../../classes/./Dakota/IssmParallelDirectApplicInterface.h: + +modules/GiaIvinsDeflectionCorex/../../classes/./IoModel.h: + +modules/GiaIvinsDeflectionCorex/../../classes/./GiaIvinsDeflectionCoreArgs.h: + +modules/GiaIvinsDeflectionCorex/../../classes/./RiftStruct.h: + +modules/GiaIvinsDeflectionCorex/../InputUpdateFromConstantx/InputUpdateFromConstantx.h: Index: /issm/trunk-jpl/src/c/modules/GiaIvinsDeflectionCorex/GiaIvinsDeflectionCorex.cpp =================================================================== --- /issm/trunk-jpl/src/c/modules/GiaIvinsDeflectionCorex/GiaIvinsDeflectionCorex.cpp (revision 21532) +++ /issm/trunk-jpl/src/c/modules/GiaIvinsDeflectionCorex/GiaIvinsDeflectionCorex.cpp (revision 21532) @@ -0,0 +1,171 @@ +/*!\file GiaIvinsDeflectionCorex + * \brief: GIA solution from Erik Ivins. + * Compute deflection wi from a single disk of radius re, load history hes for + * numtimes time steps. + */ + +#include "./GiaIvinsDeflectionCorex.h" + +#include "../../classes/classes.h" +#include "../../shared/shared.h" +#include "../../toolkits/toolkits.h" +#include "../InputUpdateFromConstantx/InputUpdateFromConstantx.h" + +/*External blocks: {{{*/ +struct blockp{ + double pset[7]; +}; + +struct blocko{ + double rhoi; +}; + +struct blockrad{ + double distrad; +}; + +struct blocks{ + double aswokm_w; + double aswokm_dwdt; +}; + +extern "C" { + int distme_(int* pNtime,int* pNtimp,int* pNtimm,double* time,double* bi,double* dmi,double* zhload); + + int what0_(int* piedge,int* pNtimp,int* pNtimm,double* time,double* bi,double* dmi); + extern struct blockp blockp_; + extern struct blocko blocko_; + extern struct blockrad blockrad_; + extern struct blocks blocks_; +} + +/*}}}*/ + +void GiaIvinsDeflectionCorex( IssmDouble* pwi, IssmDouble* pdwidt, GiaIvinsDeflectionCoreArgs* arguments){ + + /*intermediary: */ + int i; + + /*output: */ + IssmDouble wi=0; + IssmDouble dwidt=0; + + /*inputs: {{{*/ + /*constant: */ + int idisk=1; // disk # + IssmDouble yts; + + /*coming from the model structure, runtime configurable:*/ + /*gia solution parameters: */ + int iedge=0; + + /*gia inputs: */ + IssmDouble ri; //radial distance from center of disk to vertex i + IssmDouble re; //radius of disk + IssmDouble* hes; //loading history (in ice thickness) + IssmDouble* times; //loading history times + int numtimes; //loading history length + IssmDouble currenttime; + int Ntime; // number of times with load history + int Ntimm; // Ntime-1 : for slope/y-cept of load segments + int Ntimp; // Ntime+1 : for evaluation time + IssmDouble* blockt_time=NULL; + IssmDouble* blockt_bi=NULL; + IssmDouble* blockt_dmi=NULL; + IssmDouble* blocky_zhload=NULL; + + /*gia material parameters: */ + IssmDouble lithosphere_shear_modulus; + IssmDouble lithosphere_density; + IssmDouble mantle_shear_modulus; + IssmDouble mantle_viscosity; + IssmDouble mantle_density; + IssmDouble lithosphere_thickness; + + /*ice properties: */ + IssmDouble rho_ice; + + /*some debug info: */ + int disk_id; + +/*}}}*/ + + /*Recover material parameters and loading history: see GiaIvinsDeflectionCoreArgs for more details {{{*/ + ri = arguments->ri; + re = arguments->re; + hes = arguments->hes; + times = arguments->times; + numtimes = arguments->numtimes; + currenttime = arguments->currenttime; + lithosphere_shear_modulus = arguments->lithosphere_shear_modulus; + lithosphere_density = arguments->lithosphere_density; + mantle_shear_modulus = arguments->mantle_shear_modulus; + mantle_viscosity = arguments->mantle_viscosity; + mantle_density = arguments->mantle_density; + lithosphere_thickness = arguments->lithosphere_thickness; + rho_ice = arguments->rho_ice; + disk_id = arguments->idisk; + iedge = arguments->iedge; + yts = arguments->yts; + + /*}}}*/ + + /*Modify inputs to match naruse code: */ + Ntime=numtimes; + Ntimm=Ntime-1; + Ntimp=Ntime+1; + + /*Prepare block inputs for fortran distme and what0 routines of the naruse code: {{{*/ + /*Now, let's set pset from the data that we got in input to GiaIvinsDeflectionCorex: */ + blockp_.pset[0]=lithosphere_thickness; + blockp_.pset[1]=mantle_viscosity; + blockp_.pset[2]=lithosphere_shear_modulus; + blockp_.pset[3]=mantle_shear_modulus; + blockp_.pset[4]=lithosphere_density; + blockp_.pset[5]=mantle_density; + blockp_.pset[6]=re; + blocko_.rhoi=rho_ice; + + /*loading history: */ + blocky_zhload=xNew(Ntime); + for(i=0;i(Ntimp); + for (i=0;i(Ntimm); + + /*dmi: */ + blockt_dmi=xNew(Ntimm); + + /*radial distance of i-th element: */ + blockrad_.distrad=ri/1000.0; // in km + /*}}}*/ + + /*Call distme driver: */ + distme_(&Ntime,&Ntimp,&Ntimm,blockt_time,blockt_bi,blockt_dmi,blocky_zhload); + + /*Call what0 driver: */ + what0_(&iedge,&Ntimp,&Ntimm,blockt_time,blockt_bi,blockt_dmi); + + /*output solution: */ + wi = blocks_.aswokm_w; + dwidt = blocks_.aswokm_dwdt; + *pwi=wi; + *pdwidt=dwidt; + + /*Free ressources: */ + xDelete(blockt_time); + xDelete(blockt_bi); + xDelete(blockt_dmi); + xDelete(blocky_zhload); + +} Index: /issm/trunk-jpl/src/c/modules/GiaIvinsDeflectionCorex/GiaIvinsDeflectionCorex.h =================================================================== --- /issm/trunk-jpl/src/c/modules/GiaIvinsDeflectionCorex/GiaIvinsDeflectionCorex.h (revision 21532) +++ /issm/trunk-jpl/src/c/modules/GiaIvinsDeflectionCorex/GiaIvinsDeflectionCorex.h (revision 21532) @@ -0,0 +1,13 @@ +/*!\file: GiaIvinsDeflectionCorex.h + * \brief header file for ... + */ + +#ifndef _GIAIVINSDEFLECTIONCOREX_H +#define _GIAIVINSDEFLECTIONCOREX_H + +#include "../../classes/classes.h" + +/* local prototypes: */ +void GiaIvinsDeflectionCorex( IssmDouble* pwi, IssmDouble* pdwidt, GiaIvinsDeflectionCoreArgs* arguments); + +#endif /* _GIAIVINSDEFLECTIONCOREX_H */ Index: /issm/trunk-jpl/src/c/modules/GiaIvinsDeflectionCorex/distme.f =================================================================== --- /issm/trunk-jpl/src/c/modules/GiaIvinsDeflectionCorex/distme.f (revision 21532) +++ /issm/trunk-jpl/src/c/modules/GiaIvinsDeflectionCorex/distme.f (revision 21532) @@ -0,0 +1,73 @@ + subroutine distme(Ntime,Ntimp,Ntimm,time,bi,dmi,zhload) + implicit double precision (a-h,o-y) + integer Ntime,Ntimp,Ntimm + parameter (Nafter=1) + double precision pset(7) + double precision time(Ntimp),dmi(Ntimm),bi(Ntimm),dumbt(Ntimp) + double precision hload(Ntime),qpat(Ntime),qt(Ntime) + double precision zhload(Ntime),rhoi,distrad +c :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: + common /blockp/ pset + common /blockrad/ distrad + common /blocko/ rhoi + data g /9.832186d0/, yearco /3.15576d7/, eradm/6.371d6/ + data dpi /3.1415926535897932d0/, dzero/0.0d0/ +c :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: +c The units of time(Ntimp) are ka and the height of the load in meters. +c The slope, then for example, is in units of meters per ka. +c Note that "dumbt( )" is designed to perserve the initial "time( )" variable. +c :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: + do 776 k = 1, Ntimp + dumbt(k) = time(k) + 776 continue +c :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: + do 39 itime = 1, Ntime + hload(itime) = dble( zhload(itime) ) + 39 continue +c :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: +c now set up a piece-wise history: bi() = y-intercept +c dmi() = slope +c :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: + do 70 i = 2, Ntime + dmi(i-1) = ( hload(i) - hload(i-1) )/( dumbt(i) - dumbt(i-1) ) + bi(i-1) = hload(i-1) - ( dmi(i)*dumbt(i-1) ) + 70 continue +c write(6,*) zhload(1,1), zhload(1,2) +c :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: +c With pset(6) in mks units, lets convert the piecewise linear formulas +c for the time-dependent ice load heights to dimensionless values w.r.t. time. +c (tfact is in seconds) +c :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: + tfact = pset(2)/pset(4) +c :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: +c get all times as dimensionless +c :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: + do 20 jt = 1, Nafter + time(Ntime + jt) = ( dumbt(Ntime + jt) * yearco * 1.0d3 ) / tfact + 20 continue + do 75 ind = 1, Ntimm + dmi(ind) = dmi(ind) / (( yearco * 1.0d3 ) / tfact ) + 75 continue + do 77 j = 1, Ntime + time(j) = ( dumbt(j) * yearco * 1.0d3 ) / tfact + 77 continue +c :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: +c create an incremental load in Pa and non-dimensionalized: +c :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: + do 80 iq = 1, Ntime + qpat(iq) = hload(iq)*rhoi*g + qt(iq) = qpat(iq) / pset(4) + 80 continue +c :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: +c As the final step in this routine, create a dimensionless stress load from +c qp. Here we'll use bi( ) and dmi( ) vectors with dimensionless time. Then +c qp (and it's piece-wise decomposition) is ready for the direct dimensionless +c integrals for the inverse Laplace transform and inverse Hankel transform +c without further mutiplicative factors. +c :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: + do 85 i = 2, Ntime + dmi(i-1) = ( qt(i) - qt(i-1) )/( time(i) - time(i-1) ) + bi(i-1) = qt(i-1) - ( dmi(i-1)*time(i-1) ) + 85 continue + 999 return + end Index: /issm/trunk-jpl/src/c/modules/GiaIvinsDeflectionCorex/freed.f =================================================================== --- /issm/trunk-jpl/src/c/modules/GiaIvinsDeflectionCorex/freed.f (revision 21532) +++ /issm/trunk-jpl/src/c/modules/GiaIvinsDeflectionCorex/freed.f (revision 21532) @@ -0,0 +1,222 @@ + subroutine freed(r2,u2,r1,u1,h,zk,e1,e2,e4,b0,b1,a2,a1,a0,decay + 1,amps) + implicit double precision (a-h,o-z) + double precision decay(2),amps(5) + double precision ac0,ac1,ac2,ac3,ac4,ac5,ac6,ac7,ac8,ac9,ac10, + 1ac11 + common /blockz/ zkp + data zero /0.0d0/, g /9.832186d0/ +c :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: +c Given the inputs to this subroutine(r2 through zk in the call +c statement above), the outputs are coefficients of "s" that +c are crucial to the Laplace transform inversion. From b0 and b1 +c we can compute the decay poles (or eigenvalues). +c +c This is NOT true in our case, though. - SA +c NOTE IN THE CODE THAT A CALL TO THIS SUBROUTINE NEED NOT +c BE MADE AT EACH TIME STEP --- BUT WILL HAVE TO BE CALLED +c IN THE NUMERICAL INTEGRATION FOR COMPUTING THE INVERSE HANKEL +c TRANSFORM **** +c +c Each term should be returned as dimensionless +c h => length u2 => stress taumx2 => time +c :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: + zkp2 = zkp*zkp + ur = u1/u2 + ghu2 = (g*h) / u2 +c :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: +c ac0 dimensional units are stress times l^-2 +c :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: + ac0 = 4.0d0*ur*zkp2*( 1.0d0 + e4 + + 1 2.0d0*e2*(1.0d0 + 2.0d0*zkp2) ) +c DIMESIONLESS +c :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: +c ac1 dimensional units are stress times l^-2 +c :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: + ac1 = 2.0d0*r1*ghu2*zkp*(1.0d0 - e4 + 4.0d0*zkp*e2) +c DIMESIONLESS +c :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: +c ac2 dimensional units are stress^2 times l^-2 +c :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: + ac2 = 8.0d0*ur*ur*(-1.0d0 + e1)* + 1 (1.0d0 + e1)*(1.0d0 + e2)*zkp2 +c DIMESIONLESS +c :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: +c ac3 dimensional units are stress^2 times l^-2 +c :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: + ac3 = + 1 2.0d0*zkp*ghu2*ur*((r1 + r2)*(1.0d0 + e4) + + 2 2.0d0*(r2 - r1)*e2*( 1.0d0 + 2.0d0*zkp2 )) +c DIMESIONLESS +c :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: +c ac4 dimensional units are stress^2 times l^-2 +c :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: + ac4 = ghu2*ghu2*r1*(r2 - r1)* + 1 (1.0d0 - e4 + 4.0d0*zkp*e2) +c DIMESIONLESS +c :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: +c ac5 dimensional units are stress^3 times l^-2 +c :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: + ac5 = + 1 4.0d0*zkp2*ur*ur*ur*(1.0d0 - e2 - 2.0d0*e1*zkp)* + 2 (1.0d0 - e2 + 2.0d0*e1*zkp) +c DIMESIONLESS +c :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: +c ac6 dimensional units are stress^3 times l^-2 +c :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: + ac6 = + 1 2.0d0*zkp*ur*ur*(1.0d0 - e4 + 2 - 4.0d0*e2*zkp)*ghu2*r2 +c DIMESIONLESS +c :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: +c ac7 dimensional units are stress^3 times l^-2 +c :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: + ac7 = + 1 ur * ( ( (1.0d0 - e1)*(1.0d0 + e1) )**2)*r1*(r2 - r1) + 2 * ( ghu2*ghu2 ) +c DIMESIONLESS +c :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: +c ac8 dimensional units are stress^0 times l^-1 +c :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: + ac8 = -2.0d0*zkp*(1.0d0 + e2*(1.0d0 + 2.0d0*zkp*(1.0d0 + zkp))) +c DIMESIONLESS +c :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: +c ac9 dimensional units are stress^1 times l^-1 +c :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: + ac9 = + 1 ( 4.0d0*zkp*u1 - + 2 g*h*(r2 - r1)*(1.0d0 + e2*(1.0d0 + 2.0d0*zkp*(1.0d0 + zkp))) + 3 ) / u2 +c DIMESIONLESS +c :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: +c ac10 dimensional units are stress^2 times l^-1 +c :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: + ac10 = + 1 -2.0d0*zkp*ur*ur*( 1.0d0 - e2 + 2 - 2.0d0*zkp*e2*(1.0d0 + zkp) ) +c DIMESIONLESS +c :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: +c ac11 dimensional units are stress^2 times l^-1 +c :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: + ac11 = + 1 ghu2*ur*(r2 - r1)*(1.0d0 - e2*(1.0d0 + 2.0d0*zkp)) +c DIMESIONLESS +c :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: +c Other functions may be found in file "apcw0.record" +c (Nov. 9 1996) +c The following is a Mathematica version of the isolation of the +c coefficeints of the L transform variable s in the denominator. +c Here is where the set-up is performed to obtain the "free decay" +c times (with the Hankel transform variable "zk" embedded. Note that +c some greater efficency could be achieved by further simplifying the +c combinations of "acn" functions which are now a series of function +c subroutines in the fortran code. The corresponding Mathematica +c session is "twolayer.Linversion" dated Nov. 23, 1996. +c +c In[59]:= +c Together[%] +c Out[59]= +c ac2 + ac3 - ac4 + 2 ac5 - 2 ac6 + 2 ac7 +c --------------------------------------------- +c ac0 - ac1 + ac2 + ac3 - ac4 + ac5 - ac6 + ac7 +c In[61]:= +c Simplify[Coefficient[els,s^2]] +c Out[61]= +c 1 +c In[65]:= +c eslnos = +c ac5/(ac0 - ac1 + ac2 + ac3 - ac4 + ac5 - ac6 + ac7) - +c ac6/(ac0 - ac1 + ac2 + ac3 - ac4 + ac5 - ac6 + ac7) + +c +c ac7/(ac0 - ac1 + ac2 + ac3 - ac4 + ac5 - ac6 + ac7) +c +c Out[65]= +c ac5 +c --------------------------------------------- - +c ac0 - ac1 + ac2 + ac3 - ac4 + ac5 - ac6 + ac7 +c +c ac6 +c --------------------------------------------- + +c ac0 - ac1 + ac2 + ac3 - ac4 + ac5 - ac6 + ac7 +c +c ac7 +c --------------------------------------------- +c ac0 - ac1 + ac2 + ac3 - ac4 + ac5 - ac6 + ac7 +c In[66]:= +c Simplify[%] +c Out[66]= +c ac5 - ac6 + ac7 +c --------------------------------------------- +c ac0 - ac1 + ac2 + ac3 - ac4 + ac5 - ac6 + ac7 +c +c :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: +c a common denominator factor is: bc +c :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: + bc = ( ac0 - ac1 + + 1 ac2 + ac3 + 2 - ac4 + ac5 - + 3 ac6 + ac7 ) +c :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: +c b1: Denominator coefficent of s: +c :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: + b1 = + 1 ( ac2 + ac3 + 2 - ac4 + ( 2.0d0 * ac5 ) + 3 - ( 2.0d0 * ac6 ) + 4 + ( 2.0d0 * ac7 ) ) / bc +c :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: +c b0: Denominator coefficent of s^0: +c :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: + b0 = + 1 ( ac5 - ac6 + + 2 ac7 ) / bc +c :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: +c the eigenvaules are just the solution of the quadratic in s: +c so return as "decay" +c *** Note that the decay times are defined as positive *** +c if a negative inverse decay time is returned there is an error! +c :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: + sb1 = b1*b1 + fb0 = 4.0d0*b0 + diff = sb1 - fb0 + if(diff.le.zero) go to 25 + rs = dsqrt( diff ) + decay(1) = -( - b1 - rs ) / 2.0d0 + decay(2) = -( - b1 + rs ) / 2.0d0 + go to 26 + 25 idgen = 100 + go to 9990 +c :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: +c and for the numerator part of the quadratic s dependence +c the Mathematica session is: +c +c Out[14]= +c 2 +c ac10 + ac11 + (2 ac10 + 2 ac11 + ac9) s + (ac10 + ac11 + ac8 + ac9) s +c---------------------------------------------------------------------- +c ac0 - ac1 + ac2 + ac3 - ac4 + ac5 - ac6 + ac7 +c +c :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: + 26 a0 = ( ac10 + ac11 ) / bc +c :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: + a1 = ( 2.0d0*( ac10 + ac11 ) + 1 + ac9 ) / bc +c :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: + a2 = (ac10 + ac11 + 1 + ac8 + ac9) / bc +c :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: +c and the following terms are the amplitudes of the inverse Laplace +c transform solution for the non-q part. (See the boxed equation on +c page 4 of the Nov. 23 1996 notes.) +c :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: + decdif = 1.0d0/(decay(2) - decay(1)) + amps(1) = -decdif*( decay(1) * ( a1 - a2*decay(1) ) - a0 ) + amps(2) = decdif*( decay(2) * ( a1 - a2*decay(2) ) - a0 ) + amps(3) = a2 + amps(4) = - decay(1) * amps(1) + amps(5) = - decay(2) * amps(2) + go to 999 + 9990 write(6,998) idgen + 998 format(' idgen val ** fatal error ** degenerate e.v.'/1h ,1p,1i12) + 999 return + end Index: /issm/trunk-jpl/src/c/modules/GiaIvinsDeflectionCorex/ojrule.f =================================================================== --- /issm/trunk-jpl/src/c/modules/GiaIvinsDeflectionCorex/ojrule.f (revision 21532) +++ /issm/trunk-jpl/src/c/modules/GiaIvinsDeflectionCorex/ojrule.f (revision 21532) @@ -0,0 +1,63 @@ + subroutine ojrule(dk,bcin_w,bcin_dwdt) + implicit double precision(a-h,o-z) + parameter (nhank = 1024) + double precision yvalue_w(nhank),yvalue_dwdt(nhank) + double precision bcin_w(nhank),bcin_dwdt(nhank) + double precision wok_w,wok_dwdt,rpos + double precision swok_w,swok_dwdt + double precision pset(7) + double precision aswokm_w,aswokm_dwdt,distrad +c :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: + common /blockrad/ distrad + common /blockp/ pset + common /blocks/ aswokm_w,aswokm_dwdt + data zero /0.0d0/, one /1.0d0/, two /2.0d0/, three /3.0d0/, + 1rescal/ 1.0d0/ + data yearco /3.15576d7/ +c :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: + bath = dk / three +c :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: +c rpos should be normalized wrt lithosphere thickness +c give r is normalized dist_rad :: r == dist_rad / h +c :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: + r = distrad / (pset(1) / 1.0d3) + rpos = r + ak = zero +c :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: +c form the yvalue's for the Simpson's rule formulas +c :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: + do 425 ik = 1, nhank + ak = ak + dk + rak = ak * r + rarg = dbesj0( rak ) + yvalue_w(ik) = bcin_w(ik) * rarg + yvalue_dwdt(ik) = bcin_dwdt(ik) * rarg + 425 continue +c :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: +c correct to end point val. in Simp. Rule +c yvalue(nhank) = bcin(nhank) * rarg / two +c find the area under the curve using the Simpson's rule formulas +c :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: + sumde_w = zero + sumde_dwdt = zero + do 300 int = 1, nhank + intp1 = int + 1 + ide = 2 + ( (-1)**intp1 + 1 ) + fide = dfloat(ide) + sumde_w = ( fide * yvalue_w(int) ) + sumde_w + sumde_dwdt = ( fide * yvalue_dwdt(int) ) + sumde_dwdt + 300 continue + wok_w = bath * sumde_w + wok_dwdt = bath * sumde_dwdt + +c :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: + hscale = sngl(pset(1)) + hsckm = hscale / 1.0e3 + swok_w = hscale * sngl(wok_w) + aswokm_w = swok_w + swok_dwdt = (hscale * yearco * 1.0e3 * sngl(wok_dwdt)) + 1 * ( sngl(pset(4))/ sngl(pset(2)) ) + aswokm_dwdt = swok_dwdt + return + + end Index: /issm/trunk-jpl/src/c/modules/GiaIvinsDeflectionCorex/pwise.f =================================================================== --- /issm/trunk-jpl/src/c/modules/GiaIvinsDeflectionCorex/pwise.f (revision 21532) +++ /issm/trunk-jpl/src/c/modules/GiaIvinsDeflectionCorex/pwise.f (revision 21532) @@ -0,0 +1,42 @@ + subroutine pwise(t,ta,tb,xi1,xi2,xi3,xi4,slope,ycept,decay, + 1bhaq_w,bhaq_dwdt) + implicit double precision (a-h,o-z) + double precision decay(2) +c :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: +c This subroutine retrieves the convolution for the ith linear piece-wise +c q hat function (the load shape or Bessel function part having +c been removed) with the free-decay solution. (see notes of +c 12-31-96 "Convolution in time"). The convolution is returned as "bhaq". +c :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: + tbt = tb - t + tat = ta - t + gat1 = tat * decay(1) + gat2 = tat * decay(2) + gbt1 = tbt * decay(1) + gbt2 = tbt * decay(2) + ea1 = dexp(gat1) + ea2 = dexp(gat2) + eb1 = dexp(gbt1) + eb2 = dexp(gbt2) +c :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: +c define xit1 term: +c :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: + xit1 =(ycept/decay(1)) * (eb1 - ea1) - + 1(slope/(decay(1)*decay(1))) * + 2 ( (1.0d0 - tb*decay(1))*eb1 + 3 - (1.0d0 - ta*decay(1))*ea1 ) +c :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: +c define xit2 term: +c :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: + xit2 =(ycept/decay(2)) * (eb2 - ea2) - + 1(slope/(decay(2)*decay(2))) * + 2 ( (1.0d0 - tb*decay(2))*eb2 + 3 - (1.0d0 - ta*decay(2))*ea2 ) +c :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: +c add terms for the i-th interval contribution. +c ABOVE IS THE NON-DEGENERATE CASE +c :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: + bhaq_w = (xi1 * xit1) + (xi2 * xit2) + bhaq_dwdt = (xi3 * xit1) + (xi4 * xit2) + return + end Index: /issm/trunk-jpl/src/c/modules/GiaIvinsDeflectionCorex/qwise.f =================================================================== --- /issm/trunk-jpl/src/c/modules/GiaIvinsDeflectionCorex/qwise.f (revision 21532) +++ /issm/trunk-jpl/src/c/modules/GiaIvinsDeflectionCorex/qwise.f (revision 21532) @@ -0,0 +1,65 @@ + subroutine qwise(t,ta,qjadon,xi0,xi1,xi2,xi3,xi4,slope,ycept, + 1decay,bhaq_w,bhaq_dwdt) + implicit double precision (a-h,o-z) + double precision decay(2) +c :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: +c This subroutine retrieves the convolution for the J-th linear piece-wise +c q hat function (the load shape or Bessel function part having been removed) +c with the free-decay solution. (see notes of 3-27-97 "convo.ice" Mathematica +c session). The convolution is returned as "bhaq". +c +c THIS ROUTINE REPLACES pwise.f ONLY FOR t < time(Ntime) * +c (such that the load is still in place at time t). * +c +c Note irate = 1 case has to be applied to the linear term only (freed.f applies +c this correction to exponential terms) +c :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: + xg1 = xi1/(decay(1)*decay(1)) + xg2 = xi2/(decay(2)*decay(2)) + xg3 = xi3/(decay(1)*decay(1)) + xg4 = xi4/(decay(2)*decay(2)) + gb1 = decay(1)*ycept + gb2 = decay(2)*ycept +c :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: +c define xit0 term: +c :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: + xit0_w = (xi0 + qjadon) * ( ( slope * t ) + ycept ) +c :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: +c define xit1 term: +c :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: + xit1_w = xg1 * ( + 1 gb1 + slope * ( ( t * decay(1) ) - 1.0d0 ) + 2 - ( gb1 + slope * ( ( ta * decay(1) ) - 1.0d0 )) + 3 * dexp( decay(1) * (ta - t) ) + 4 ) +c :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: +c define xit2 term: +c :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: + xit2_w = xg2 * ( + 1 gb2 + slope * ( ( t * decay(2) ) - 1.0d0 ) + 2 - ( gb2 + slope * ( ( ta * decay(2) ) - 1.0d0 ) ) + 3 * dexp( decay(2) * (ta - t) ) + 4 ) +c :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: +c And the rate equivalents: +c (sign switch due to freed.f already +c having corrected in x1t, x2t pass). +c :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: + xit0_dwdt = (xi0 + qjadon) * slope + xit1_dwdt =-xg3 * ( + 1 slope + 2 + ( gb1 + slope * ( ( ta * decay(1) ) - 1.0d0 )) + 3 * dexp( decay(1) * (ta - t) ) + 4 ) + xit2_dwdt =-xg4 * ( + 1 slope + 2 + ( gb2 + slope * ( ( ta * decay(2) ) - 1.0d0 )) + 3 * dexp( decay(2) * (ta - t) ) + 4 ) +c :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: +c add terms for the J-th (and final) interval contribution. +c :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: + bhaq_w = xit0_w + xit1_w + xit2_w + bhaq_dwdt = xit0_dwdt + xit1_dwdt + xit2_dwdt + return + end Index: /issm/trunk-jpl/src/c/modules/GiaIvinsDeflectionCorex/stot.f =================================================================== --- /issm/trunk-jpl/src/c/modules/GiaIvinsDeflectionCorex/stot.f (revision 21532) +++ /issm/trunk-jpl/src/c/modules/GiaIvinsDeflectionCorex/stot.f (revision 21532) @@ -0,0 +1,79 @@ + subroutine stot(ikval,qjadon,fltng_w,fltng_dwdt,Ntimp,Ntimm, + 1time,bi,dmi) + implicit double precision (a-h,o-z) + integer Ntimp,Ntimm + parameter (Nafter = 1) + parameter (nhank = 1024) + double precision decay(2) + double precision pset(7) + double precision time(Ntimp),bi(Ntimm),dmi(Ntimm) + double precision dekay1(nhank),dekay2(nhank),amp0(nhank), + 1amp1(nhank),amp2(nhank),amp3(nhank),amp4(nhank) +c :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: + common /blockm/ dekay1,dekay2,amp0,amp1,amp2,amp3,amp4 +c :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: +c This subroutine returns the inverse Laplace transform to the +c time-domain for the vertical displacement at time t for Hankel wavenumber +c ikval. (In general this routine needs to be called nhank times). +c The main derivation uses the Faltung theorem of Laplace transforms. +c (1-1-97) NEW CASE OF 3-27-97 IS FOR t(Ntime + i) < t(Ntime) OR IN OTHER +c WORDS, THE LOAD STILL IN PLACE AT t. OPTION CALL to qwise.f +c PERFORMS THIS. +c +c A theory for the degenerate case was worked out but has been removed +c as an option from this code. +c +c Definition of tspan: nondimensional time span backwards form present +c when this routine is first called the dimensional +c equivalent might be say tspan = 12 ka, then 11 and +c then finally tspan = 0. +c Additional note for r.s.l calculations: the routines qwise and pwise +c are identical to the previous case for computations of present-day only +c vertical deformation field. +c :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: + decay(1) = dekay1(ikval) + decay(2) = dekay2(ikval) + xi0 = amp0(ikval) + xi1 = amp1(ikval) + xi2 = amp2(ikval) + xi3 = amp3(ikval) + xi4 = amp4(ikval) + +c :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: +c note that tspan must be updated in the calling routine "what0.f" +c :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: + t = time(Ntimp) + sumb_w = 0.0d0 + sumb_dwdt = 0.0d0 + ta = time(1) + tb = time(2) +c :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: + do 97 i = 1,Ntimm +c :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: +c note that this "if" prevents adding load +c segments of "future" times when computing +c an r.s.l. history (10-06-98). +c :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: + if(t.lt.ta) go to 97 + slope=dmi(i) + ycept=bi(i) +c :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: + if( t . gt . ta . and . t . le . tb) go to 38 + call pwise(t,ta,tb,xi1,xi2,xi3,xi4,slope,ycept,decay, + 1bhaq_w,bhaq_dwdt) + go to 39 +c :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: +c Note that qwise is employed only for the J-th Q hat term when t for +c evaluation still has to consider the load itself +c :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: + 38 call qwise(t,ta,qjadon,xi0,xi1,xi2,xi3,xi4,slope,ycept,decay, + 1bhaq_w,bhaq_dwdt) + 39 sumb_w = bhaq_w + sumb_w + sumb_dwdt = bhaq_dwdt + sumb_dwdt + ta = time(i + 1) + tb = time(i + 2) + 97 continue + fltng_w = sumb_w + fltng_dwdt = sumb_dwdt + return + end Index: /issm/trunk-jpl/src/c/modules/GiaIvinsDeflectionCorex/what0.f =================================================================== --- /issm/trunk-jpl/src/c/modules/GiaIvinsDeflectionCorex/what0.f (revision 21532) +++ /issm/trunk-jpl/src/c/modules/GiaIvinsDeflectionCorex/what0.f (revision 21532) @@ -0,0 +1,127 @@ + subroutine what0(iedge,Ntimp,Ntimm,time,bi,dmi) + implicit double precision (a-h,o-z) + integer Ntimp,Ntimm + parameter (nhank = 1024) + parameter (N = nhank/2) + double precision dekay1(nhank),dekay2(nhank),amp0(nhank), + 1amp1(nhank),amp2(nhank),amp3(nhank),amp4(nhank), + 1zksam(nhank),zksamp(nhank) + double precision decay(2),pset(7),amps(5), + 1decta(2),dyri1(nhank),dyri2(nhank),sna(nhank) + double precision cinner_w(nhank),cinner_dwdt(nhank) + double precision bcin_w(nhank),bcin_dwdt(nhank) + double precision time(Ntimp),bi(Ntimm),dmi(Ntimm) + integer maxk +c :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: + common /blockp/ pset + common /blockz/ zkp + common /blockm/ dekay1,dekay2,amp0,amp1,amp2,amp3,amp4 + data yearco /3.15576d7/, pi /3.1415926535897932384d0/ + data g /9.832186d0/, four /4.d0/, two /2.0d0/, + 1 one /1.0d0/, zero/0.0d0/ , maxk/64/ +c :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: + twopi = two * pi + r2 = pset(6) + u2 = pset(4) + r1 = pset(5) + u1 = pset(3) + h = pset(1) + urat = u1/u2 +c :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: +c alphap is dimensionless disk radius +c :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: + alphap = pset(7)/pset(1) + twoap = two * alphap + rghm = ( r1 * g * h * alphap ) / (two * u2) + taumx = pset(2)/pset(4) + tmxyr = taumx / yearco +c + dfac = dfloat(nhank)/dfloat(maxk) + endk = dfloat(nhank)/dfac + dk = endk/dfloat(nhank) +c + ak = zero + do 7000 ik = 1,nhank + ak = ak + dk + zkp = ak + pikn = (6.371d6 * zkp) / h + zkd = pikn / 6.371d6 +c + zkp2 = 2.0d0 * zkp + zkp4 = 4.0d0 * zkp + e1 = dexp(zkp) + e2 = dexp(zkp2) + e4 = dexp(zkp4) +c :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: + call freed(r2,u2,r1,u1,h,zkd,e1,e2,e4,b0,b1,a2,a1,a0,decay,amps) +c :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: + decta(1) = decay(1)/tmxyr + decta(2) = decay(2)/tmxyr + dyri1(ik) = decta(1) + dyri2(ik) = decta(2) + sna(ik) = pikn +c :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: +c Form vectors for full construction in pwise.f and stot.f +c Note that freed will produce decay spectra defined as positive, +c ie. negative decay must reinsert a minus sign. +c :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: + dekay1(ik) = decay(1) + dekay2(ik) = decay(2) + amp0(ik) = amps(3) + amp1(ik) = amps(1) + amp2(ik) = amps(2) + amp3(ik) = amps(4) + amp4(ik) = amps(5) + zksam(ik) = zkd + zksamp(ik) = zkp + 7000 continue +c :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: +c The following looped call sets up the free solution convolved with the +c load function q hat. Note that the returned vector set "cinner" is the +c inner-most part of the arguement of the inverse Hankel trans. integral. +c It is time-dependent and the loop is for the k-dependancy. The time for +c calculation is given in the vector "time(Ntimp)" in the routine stot.f that is +c called below. Note that the sign on cinner(ik) below is for a load directed +c downward. ** For iedge = 1 assume sq. edge load and for iedge = 2 assume an +c elliptical cross section. Note loops 8500,8000 and 9500,9000 for the two +c cases, respectively. +c :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: + 49 go to (8499,9499), iedge + 8499 do 8000 ik = 1, nhank + xakap = zksamp(ik)*alphap + diku = xakap * urat + pref = diku / ( diku + rghm ) + qjadon = one / ( four * zksamp(ik) * urat ) + call stot(ik,qjadon,fltng_w,fltng_dwdt,Ntimp,Ntimm,time,bi,dmi) + cinner_w(ik) = - fltng_w * pref * twoap + cinner_dwdt(ik) = - fltng_dwdt * pref * twoap + bcin_w(ik) = cinner_w(ik) * dbesj1(xakap) + bcin_dwdt(ik) = cinner_dwdt(ik) * dbesj1(xakap) + 8000 continue +c :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: +c "ojrule.f" computes the inverse Hankel trasform with a simple +c Simpson's rule. The routine "ojrule" is buliding a set of solutions stored +c in common "blocks" in r or "asrpos(nrv) ", and computed rate or displacement +c for each of N3G disks in "aswokm(nrv,N3G)" . +c :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: + call ojrule(dk,bcin_w,bcin_dwdt) + go to 999 +c :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: + 9499 do 9000 ik = 1, nhank + xakap = zksamp(ik)*alphap + oxakap = one/xakap + diku = xakap * urat + pref = diku / ( diku + rghm ) + qjadon = one / ( four * zksamp(ik) * urat ) + call stot(ik,qjadon,fltng_w,fltng_dwdt,Ntimp,Ntimm,time,bi,dmi) + cinner_w(ik) = - fltng_w * pref * twoap + cinner_dwdt(ik) = - fltng_dwdt * pref * twoap + bcin_w(ik) = cinner_w(ik) * oxakap * ( dsin(xakap) * oxakap + 1 - dcos(xakap) ) + bcin_dwdt(ik) = cinner_dwdt(ik) * oxakap * ( dsin(xakap) * oxakap + 1 - dcos(xakap) ) + 9000 continue + call ojrule(dk,bcin_w,bcin_dwdt) + 999 return + + end