C言語 float型の桁落ち回避について

Question

C言語の二次方程式の桁落ち回避の課題です．
エラーは出ないしおよそ値は合っているのですが，課題の指示でfloat型では通常の解の公式では精度が著しく落ちることを確認せよ，とあるのに値が変わりません．これはどこかミスしているのでしょうか？
以下がコードです．a=1,b=200,c=1で確認せよ，との指示です．

#include<stdio.h>

int func(int a,int b, int c, float *p01, float *p02, float *p11,float *p12);

int main (){
	float ans01, ans02, ans11, ans12;
	printf("float________________________________");
	func(1,200,1, &ans01, &ans02, &ans11, &ans12); //ans[0通常,1桁落ち回避][解の番号]
	printf("\n\n ans01=%f, ans02=%f \n",ans01,ans02);	
	printf("\n\n ans11=%f, ans12=%f \n\n",ans11,ans12);
	return 0;
	}
	
int func(int a,int b, int c, float *p01, float *p02, float *p11,float *p12){
	float ans01,ans02,ans11,ans12;
	//通常
	ans01=(-b+sqrt(b*b-4*a*c))/(2*a);
	ans02=(-b-sqrt(b*b-4*a*c))/(2*a);
	*p01=ans01;
	*p02=ans02;
	//桁落ち回避
	ans11=(-2*c)/(b+sqrt(b*b-4*a*c));
	ans12=(-b-sqrt(b*b-4*a*c))/(2*a);
	*p11=ans11;
	*p12=ans12;
	return 0;	
	}
	
コード

Answer 1

こんにちは。

とあるのに値が変わりません．これはどこかミスしているのでしょうか？

その通りです。floatで計算せず、doubleで計算しているところがあります。
ちょうどそこがfloatでは大きく桁落ちする場所なので、doubleでやった結果桁落ちが見えなくなってます。

具体的に書いてしまうと答えになってしまうので、ヒントだけ。

Comments

[north_redwings]

sqrtの部分がdoubleで計算されている，ということでしょうか？
一応(float)でキャスティングしたところ，値が変わったので合っていると思いたいのですが...．

[Chironian]

sqrtの説明を確認するとよいですよ。http://www.c-tipsref.com/reference/math/sqrt.html

標準ライブラリを全て暗記するのは現実的ではありませんから、疑問を感じた時は調べてみると意外な発見がありますよ。

[north_redwings]

迅速な回答ありがとうございました．
C言語は他の言語より難しいですね...
もっと勉強したいと思います．