二分法(2変数)

Question

###前提・実現したいこと
xを0から2πまで動かしたとき、cos(x)+cos(y)=0 をyについて解きたいです。
###発生している問題・エラーメッセージ

xを動かしてもyが動きません。

###該当のソースコード

#include <iostream>
#include <iomanip>
#include <cmath>
#include <complex>
using namespace std;

double f(double x,double y){
    return cos(x)+cos(y);
}

class bisection{
    int i;
public:
    bisection();
    double m;
    double x;
    double operator()(double(*func)(double ,double),double a,double b,double eps=1.e-12,int Nmax=50){
        for ( i= 1;i<=Nmax;i++) {
            m = (a + b) /2.;
            if (func(x,m) > 0)
                b = m;
            else
                a = m;
            if (func(x,m) == 0 || fabs(a - b) < eps) {
                return m;
                break;
            }
        }
        
        if (i > Nmax)
            cout << "bisection cannot slove it" << endl;
        return m;
    }
};

int main(){
    bisection ob;
    int N=100;
    double dx=6.28/N;
    for(int i=0;i<N;i++){
        double x=i*dx;
    ob(f,0.,6.28);
    cout << x << " " <<ob.m<< endl;
    }
    return 0;
}

###試したこと
課題に対してアプローチしたことを記載してください

###補足情報(言語/FW/ツール等のバージョンなど)
コードが汚く見にくいかと思いますが、回答宜しくお願いします。
より洗練されたコーディングの仕方があればそちらも教えていただけるとありがたいです。

Answer 1

より洗練されたコーディングの仕方があればそちらも教えていただけるとありがたいです。

少なくともメンバ変数を晒してない版:

#include <iostream>
#include <functional>
#include <cmath>

class bisection {
private:
    double eps_;
    int nmax_;
public:
    bisection(double eps = 1.0e-12, int nmax =50) : eps_(eps), nmax_(nmax) {};
    double operator()(std::function<double(double,double)> func, double x, double a,double b) {
      int i;
      double m;
      for ( i = 0; i < nmax_; ++i ) {
        m = (a + b) / 2.;
        if ( func(x,m) > 0 ) b = m; else a = m;
        if (func(x,m) == 0 || fabs(a - b) < eps_) {
          break;
        }
      }

      if ( i == nmax_ )
        std::cout << "bisection cannot slove it" << std::endl;
      return m;
    }
};

int main(){
  bisection ob;
  const int     N = 100;
  const double dx = 6.28/N;
  auto f = [](double x,double y) -> double { return cos(x) + cos(y); };

  for ( int i = 0; i < N; ++i) {
    double x = i * dx;
    double y = ob(f, x, 0.,6.28);
    std::cout << x << " " << y << std::endl;
  }
}

Answer 2

for(int i=0;i<N;i++){
  double x=i*dx;
  ob(f,0.,6.28);
  cout << x << " " <<ob.m<< endl;
}

これじゃダメでしょ。ob.x を書き換えないと。
正しくは:

for(int i=0;i<N;i++){
    ob.x=i*dx;
    ob(f,0.,6.28);
    cout << ob.x << " " <<ob.m<< endl;
}

Comments

[sukiyaki]

解決しました。ありがとうございます。