| 1 | // | 
|---|
| 2 | // scextrap.cc | 
|---|
| 3 | // | 
|---|
| 4 | // Copyright (C) 1996 Limit Point Systems, Inc. | 
|---|
| 5 | // | 
|---|
| 6 | // Author: Curtis Janssen <cljanss@limitpt.com> | 
|---|
| 7 | // Maintainer: LPS | 
|---|
| 8 | // | 
|---|
| 9 | // This file is part of the SC Toolkit. | 
|---|
| 10 | // | 
|---|
| 11 | // The SC Toolkit is free software; you can redistribute it and/or modify | 
|---|
| 12 | // it under the terms of the GNU Library General Public License as published by | 
|---|
| 13 | // the Free Software Foundation; either version 2, or (at your option) | 
|---|
| 14 | // any later version. | 
|---|
| 15 | // | 
|---|
| 16 | // The SC Toolkit is distributed in the hope that it will be useful, | 
|---|
| 17 | // but WITHOUT ANY WARRANTY; without even the implied warranty of | 
|---|
| 18 | // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the | 
|---|
| 19 | // GNU Library General Public License for more details. | 
|---|
| 20 | // | 
|---|
| 21 | // You should have received a copy of the GNU Library General Public License | 
|---|
| 22 | // along with the SC Toolkit; see the file COPYING.LIB.  If not, write to | 
|---|
| 23 | // the Free Software Foundation, 675 Mass Ave, Cambridge, MA 02139, USA. | 
|---|
| 24 | // | 
|---|
| 25 | // The U.S. Government is granted a limited license as per AL 91-7. | 
|---|
| 26 | // | 
|---|
| 27 |  | 
|---|
| 28 | #ifdef __GNUC__ | 
|---|
| 29 | #pragma implementation | 
|---|
| 30 | #endif | 
|---|
| 31 |  | 
|---|
| 32 | #include <util/state/stateio.h> | 
|---|
| 33 | #include <math/optimize/scextrap.h> | 
|---|
| 34 |  | 
|---|
| 35 | using namespace sc; | 
|---|
| 36 |  | 
|---|
| 37 | static ClassDesc SCExtrapData_cd( | 
|---|
| 38 | typeid(SCExtrapData),"SCExtrapData",1,"public SavableState", | 
|---|
| 39 | 0, 0, 0); | 
|---|
| 40 |  | 
|---|
| 41 | SCExtrapData::SCExtrapData() | 
|---|
| 42 | { | 
|---|
| 43 | } | 
|---|
| 44 |  | 
|---|
| 45 | SCExtrapData::SCExtrapData(StateIn& s) : | 
|---|
| 46 | SavableState(s) | 
|---|
| 47 | { | 
|---|
| 48 | } | 
|---|
| 49 |  | 
|---|
| 50 | SCExtrapData::~SCExtrapData() | 
|---|
| 51 | { | 
|---|
| 52 | } | 
|---|
| 53 |  | 
|---|
| 54 | void | 
|---|
| 55 | SCExtrapData::save_data_state(StateOut& s) | 
|---|
| 56 | { | 
|---|
| 57 | } | 
|---|
| 58 |  | 
|---|
| 59 | //////////////////////////////////////////////////////////////////////////// | 
|---|
| 60 |  | 
|---|
| 61 | static ClassDesc SCExtrapError_cd( | 
|---|
| 62 | typeid(SCExtrapError),"SCExtrapError",1,"public SavableState", | 
|---|
| 63 | 0, 0, 0); | 
|---|
| 64 |  | 
|---|
| 65 | SCExtrapError::SCExtrapError() | 
|---|
| 66 | { | 
|---|
| 67 | } | 
|---|
| 68 |  | 
|---|
| 69 | SCExtrapError::SCExtrapError(StateIn& s) : | 
|---|
| 70 | SavableState(s) | 
|---|
| 71 | { | 
|---|
| 72 | } | 
|---|
| 73 |  | 
|---|
| 74 | SCExtrapError::~SCExtrapError() | 
|---|
| 75 | { | 
|---|
| 76 | } | 
|---|
| 77 |  | 
|---|
| 78 | void | 
|---|
| 79 | SCExtrapError::save_data_state(StateOut& s) | 
|---|
| 80 | { | 
|---|
| 81 | } | 
|---|
| 82 |  | 
|---|
| 83 | //////////////////////////////////////////////////////////////////////////// | 
|---|
| 84 |  | 
|---|
| 85 | static ClassDesc SelfConsistentExtrapolation_cd( | 
|---|
| 86 | typeid(SelfConsistentExtrapolation),"SelfConsistentExtrapolation",1,"public SavableState", | 
|---|
| 87 | 0, 0, 0); | 
|---|
| 88 |  | 
|---|
| 89 | SelfConsistentExtrapolation::SelfConsistentExtrapolation() | 
|---|
| 90 | { | 
|---|
| 91 | errorset_ = 0; | 
|---|
| 92 | tolerance_ = 1.0e-8; | 
|---|
| 93 | } | 
|---|
| 94 |  | 
|---|
| 95 | SelfConsistentExtrapolation::SelfConsistentExtrapolation(StateIn& s) : | 
|---|
| 96 | SavableState(s) | 
|---|
| 97 | { | 
|---|
| 98 | s.get(error_); | 
|---|
| 99 | s.get(errorset_); | 
|---|
| 100 | s.get(tolerance_); | 
|---|
| 101 | } | 
|---|
| 102 |  | 
|---|
| 103 | SelfConsistentExtrapolation::SelfConsistentExtrapolation( | 
|---|
| 104 | const Ref<KeyVal>&keyval) | 
|---|
| 105 | { | 
|---|
| 106 | errorset_ = 0; | 
|---|
| 107 | tolerance_ = keyval->doublevalue("tolerance"); | 
|---|
| 108 | if (keyval->error() != KeyVal::OK) tolerance_ = 1.0e-8; | 
|---|
| 109 | } | 
|---|
| 110 |  | 
|---|
| 111 | SelfConsistentExtrapolation::~SelfConsistentExtrapolation() | 
|---|
| 112 | { | 
|---|
| 113 | } | 
|---|
| 114 |  | 
|---|
| 115 | void | 
|---|
| 116 | SelfConsistentExtrapolation::save_data_state(StateOut& s) | 
|---|
| 117 | { | 
|---|
| 118 | s.put(error_); | 
|---|
| 119 | s.put(errorset_); | 
|---|
| 120 | s.put(tolerance_); | 
|---|
| 121 | } | 
|---|
| 122 |  | 
|---|
| 123 | void | 
|---|
| 124 | SelfConsistentExtrapolation::start_extrapolation() | 
|---|
| 125 | { | 
|---|
| 126 | } | 
|---|
| 127 |  | 
|---|
| 128 | ///////////////////////////////////////////////////////////////////////////// | 
|---|
| 129 |  | 
|---|
| 130 | // Local Variables: | 
|---|
| 131 | // mode: c++ | 
|---|
| 132 | // c-file-style: "CLJ" | 
|---|
| 133 | // End: | 
|---|