/**************************************************************************/
/**************************************************************************/
/**************************   "steric"   **********************************/
/**************************************************************************/
/*************     Program to calculate ligand cone    ********************/
/*************     angles as a measure of steric size  ********************/
/**************************************************************************/
/**************************************************************************/

/**************************************************************************/
/******************   Additional Aid Functions   **************************/
/**************************************************************************/
/******************        This module is        **************************/
/******************      system independant      **************************/
/**************************************************************************/

#include <stdio.h>
#include <math.h>
#include <string.h>

#include "sterdefn.h" /* declaration of structures, functions and globals */
#include "steraid.h"  /* additional functions needed                      */
#include "stertext.h" /* all text functions for text mode                 */

/**************************************************************************/
/***************************************************************************/

char *Last_Character(char *string, char chr)
{
  int i=0,len;
  char *ptr;
  char newline[LINELEN];
  ptr=string;
  len=strlen(string);
  while((*ptr==chr)&&(i<len))
  {
    i++;
    ptr++;
  }
  strcpy(newline,ptr);
  strcpy(string,newline);
  return(string);
}

/***************************************************************************/

int Alpha_Beta(char ch)
{
  if(((ch>64)&&(ch<91))||((ch>96)&&(ch<123))) return(1);
  else return(0);
}

/**************************************************************************/

int Int_String(char *line)
{
  int i,len=strlen(line);
  char newline[LINELEN];
  strcpy(newline,line);
  for(i=0;i<len;i++) if((newline[i]>57)||(newline[i]<48)) newline[i]=' ';
  i=0;
  if(sscanf(newline,"%d",&i)) return(i);
  return(0);
}

/**************************************************************************/

char *strupr(char *s)
{
  int i;
  for(i=0;s[i]!=0;i++) if((s[i]>96)&&(s[i]<123)) s[i]-=32;
  return(s);
}

/**************************************************************************/

int stricmp(char *s1, char *s2)
{
  char ns1[LINELEN], ns2[LINELEN];
  strcpy(ns1,s1);
  strcpy(ns2,s2);
  strupr(ns1);
  strupr(ns2);
  return(strcmp(ns1,ns2));
}

/**************************************************************************/

double Change_Double(char *title, double val, FILE *input)
{
  char line[LINELEN];

  sprintf(line,"%s (%f):",title,val);
  Get_Input_Line(line,input);
  if (strlen(line)!=0) sscanf(line,"%lf",&val);
  return(val);
}

/**************************************************************************/

int Change_Integer(char *title, int val, FILE *input)
{
  char line[LINELEN];

  sprintf(line,"%s (%d):",title,val);
  Get_Input_Line(line,input);
  if (strlen(line)!=0) sscanf(line,"%d",&val);
  return(val);
}

/**************************************************************************/

long get_next_line(FILE *in, char *string, char *line)
{
  char card[LINELEN];
  int i,len=strlen(string);
  line[0]=0;
  strupr(string);
  while(fgets(line,148,in)!=NULL)
  {
    for(i=0;i<len;i++) card[i]=line[i];
    card[i]=0;
    strupr(card);
    if(strcmp(card,string)==0) return(ftell(in));
  }
  return(0L);
}

/**************************************************************************/

int count_lines(FILE *in, char *string)
{
  int num=0;
  long fpos=0L;
  int i,len=strlen(string);
  char line[LINELEN], card[LINELEN];
  fpos=ftell(in);
  if(fpos<0) fpos=0;
  fseek(in,0L,SEEK_SET);
  strupr(string);
  while(fgets(line,148,in)!=NULL)
  {
    for(i=0;i<len;i++) card[i]=line[i];
    card[i]=0;
    strupr(card);
    if(strcmp(card,string)==0) num++;
  }
  fseek(in,fpos,SEEK_SET);
  return(num);
}

/**************************************************************************/

double Radians(double A)
{
  return(A * PI / 180);
}

/**************************************************************************/

double Degrees(double A)
{
  return(A * 180 / PI);
}

/**************************************************************************/

double fpos(double A)
{
  if(A>0) return(A);
  else return(0.0);
}

/**************************************************************************/

double fmax(double A, double B)
{
  if (A > B) return(A);
  else return(B);
}

/**************************************************************************/

double fmin(double A, double B)
{
  if (A < B) return(A);
  else return(B);
}

/**************************************************************************/

double fmax3(double A, double B, double C)
{
  if(A>B) if(A>C) return(A);
	  else return(C);
  else if(C>B) return (C);
	  else return(B);
}

/**************************************************************************/

double fmin3(double A, double B, double C)
{
  if(A<B) if(A<C) return(A);
	  else return(C);
  else if(C<B) return (C);
	  else return(B);
}

/**************************************************************************/

double StoD(double A)
{
  return(Degrees(2*acos( 1-(A/(2*PI)) )));
}

/**************************************************************************/

double RtoD(double A)
{
  return(360*(A/(2*PI)));
}

/**************************************************************************/

char Small(double D)
{
  return(fabs(D) < 0.05);
}

/**************************************************************************/

char AlmostZero(double D)
{
  return(fabs(D) < 0.00000001);
}

/**************************************************************************/

int Similar(double A, double B, int percent)
{
  if(2*fabs(A-B)/(A+B)<(((double)percent)/100.0)) return(1);
  return(0);
}

/**************************************************************************/

double msqrt(double X, char *name)
{
  if (X<0) {X=0; Error_Message(E_SQRTNG,name);}
  return(sqrt(X));
}

/**************************************************************************/

double arcTan(double Y, double X)
{
  double A;
  if (AlmostZero(X))
  {
    if (Y >= 0) return(PI/2);
    else return(-PI/2);
  }
  else
  {
    A = atan(Y/X);
    if (A < 0) A = A + PI;
    return(A);
  }
}

/**************************************************************************/

double arcCos(double X)
{
  if (X>1)
  {
    printf("\nArcCos domain error x=%f\n",X);
    return(0.0);
  }
  if (X<-1)
  {
    printf("\nArcCos domain error x=%f\n",X);
    return(PI);
  }
  else return(acos(X));
}

/**************************************************************************/

double cosrule_angle(double A, double B, double C)
{
  double cosalpha;

  if ((A==0.0)||(B==0.0)||(C==0.0)) return(0.0);
  cosalpha = ((A*A)+(B*B)-(C*C))/(2*A*B);
  if (cosalpha>1) {cosalpha=1;printf(" coserr ");}
  if (cosalpha<-1) {cosalpha=-1;printf(" coserr ");}
  return (acos(cosalpha));
}

/**************************************************************************/

int power(int X, int Y)
{
  int count;
  int ans;
  ans=1;
  for(count=0;count<Y;count++)
   ans*=X;
  return(ans);
}

/**************************************************************************/
/**************************************************************************/
/*************************  The End ... ***********************************/
/**************************************************************************/