| 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:
 | 
|---|